# Genererer tabeller for htlm-docs

#library(R2HTML)

if(F)
{
f2 = par_read("epil2.par")

admb2 = c(f2$b,f2$log_sigma1,f2$log_sigma2,f2$log_alpha)

tmp = read.table("sasoutput1.txt")
sas2 = c(tmp[c(1:6,8:9),2],-log(tmp[7,2]))

tab1 = matrix(
	c(admb2,sas2),
	ncol=2,
	dimnames=list(c(rep("b",6),"log(sigma_1)","log(sigma_2)","log(alpha)"),
	c("ADMB-RE","SAS"))
	)

HTML2clip(tab1,
	align = "left",
	digits=4,
	caption="Comparison of parameter estimates for SAS NLMIXED and ADMB-RE.")
}

if(T)
{
f3 = par_read("epil3b.par")

admb3 = c(f3$b,f3$log_sigma,f3$w,f3$log_alpha)

tmp = read.table("sasoutput3.txt")
sas3 = c(tmp[1:9,2],-log(tmp[10,2]))

tab3 = matrix(
	c(admb3,sas3),
	ncol=2,
	dimnames=list(c(rep("b",3),"log(sigma_1)","log(sigma_2)","log(sigma_3)",
		"w1","w2","w3","log(alpha)"),
	c("ADMB-RE","SAS"))
	)

HTML2clip(tab3,
	align = "left",
	digits=4,
	caption="Comparison of parameter estimates for SAS NLMIXED and ADMB-RE.")

}