Displaying 1 result from an estimated 1 matches for "fijos2".
Did you mean:
fijos
2008 Apr 30
1
error with lme within a loop
...=B%*%U%*%Delta
X=NULL
for(i in 0:(pord-1)){X=cbind(X,x^i)}
### This X and Z matrix are the matrices in the mixed model
representation of a spline
n<-nrow(X)
Id<- factor(rep(1,n))
Z.block<-list(list(Id=pdIdent(~Z-1)),list(barrio=pdIdent(~1)))
Z.block<-unlist(Z.block,recursive=FALSE)
fijos2<-c(-326.9178203,0.1645375 )
ef.spline<-c(-0.2214524, -0.1649163, -0.1649163, 0.8878776, -0.2214524)
set.seed(48)
ef.error<-rnorm(135,0,sqrt(0.03073974))
where
ef.spline are standardized random effects, and errores are the residuals.
My loop looks like this
B<-1000
y.boot<...