Displaying 2 results from an estimated 2 matches for "yijk".
Did you mean:
dijk
2009 Sep 23
1
Best way to arrange data
...rray(0,c(G,T,length(nj)))
for (b in 1:length(gene)) {
for (i in 1:max(gene)) {
for (j in 1:max(groups)) {
for (k in 1:max(reps)) {
if ((gene[b]==i) && (groups[b]==j) && (replicate[b]==k))
{
yijk_array[i,j,k]=yijk[b]
}
}
}
}
}
But this gives an array that has a lot of zeros in places where I don't have
any data....
2009 Aug 21
0
data layout for crossed factors w/interaction in linear mix models
Dear All,
I am trying to fit a simple linear mixed model (see below this paragraph) arising from a crossed factorial design with 2 factors and ubalanced number of replicates (from two to five) in each cell, but I keep getting an error message (see bottom of message). The model is:
yijk = intercept + ai + bj + abij + ejik, where:
"intercept" is fixed, and the crosss factors, ai, i = 1,..,10, and bj, j= 1,..,10, are random. I am interested in estimating the variance components of these factors AND their interaction. I have tried:
fm1 <- lmer(formula = V1~1 + (1|V2)...