Displaying 1 result from an estimated 1 matches for "densitylow".
2012 Jun 21
1
lme random effects in additive models with interaction
...del. But, I have two
data sets (Replication 1 and 2) and I want that Replication is random
effect. Here is my code:
knots <- default.knots(Time)
z <- outer(Time, knots, "-")
z <- z * (z > 0)
z<-z^2
i.size50 <- I(Size==50)
i.chargepos <- I(Charge=="+")
i.densitylow <- I(Density==0.001)
X <- cbind( I(Time^2),Time*i.size50,Time*i.chargepos,Time*i.densitylow)
Z <- cbind( z, z*i.size50, z*i.chargepos,z*i.densitylow)
K <- length(knots)
block.ind <- list(1:K, (K+1):(2*K),(2*K+1):(3*K),(3*K+1):(4*K))
Z.block <- list()
for (i in 1:length(block.in...