search for: dati1

...e of R in order to estimate heritability. I want to estimate the confidence interval of the variance coefficient and so I should use a bootstrap simulation. The pedigree file has 1386 subjects so I create a kinship matrix [1386*1386].This is the code of R I use: kfit2 <- lmekin(IT~1+AGE +(1|ID), dati1, varlist=list(kmat))kfit2kfit2$vcoef library(boot) var <- function(formula,data,indices,varlist) { d <- data[indices,] # allows boot to select sample kfit2 <- lmekin(formula, data=d, varlist=list(kmat)) return(kfit2$vcoef)} # bootstrapping with 1000 replications results <- boot(da...

I'm fitting a Poisson gam model, say model<-gam(a65tm~as.factor(day.week )+as.factor(week)+offset(log(pop65))+s(time,k=10,bs="cr",fx=FALSE,by=NA,m=1),sp=c( 0.001),data=dati1,family=poisson) Currently I've difficulties in obtaining right predictions by using gam.predict function with MGCV package in R version 2.2.1 (see below my syntax). previ<-predict.gam(model,dati2,type="response") I expect to get the predicted values in the same scale of the re...