Hi!

First I would like to thank you for the help with my first question. However,I am afraid I am still into trouble. I think I have properly used the getGropus commad but seems to be a fatal error when calculating the fixed effects (see bellow) and I cannot understand why. Any inputs? 

Thanks in advance,

/Jose

hongos1=read.table("c:/temp/hongos1.txt", header=T)
> attach(hongos1)
> Tra=factor(Tra)
> Ind=factor(Ind)
> Rep=factor(Rep)
> names(hongos1)
[1] "Esp" "Tra" "Ind" "Rep" "Are" "Cre"
> hongos1$TraInd= getGroups(hongos1, form = ~1|Tra/Ind, level = 2)
> model=lme(Cre~Tra,random=~1|Tra/Ind/Rep)
> summary(model)
Linear mixed-effects model fit by REML
 Data: NULL 
       AIC     BIC    logLik
  30.88467 41.5423 -7.442334

Random effects:
 Formula: ~1 | Tra
        (Intercept)
StdDev:  0.07849531

 Formula: ~1 | Ind %in% Tra
        (Intercept)
StdDev:   0.1536001

 Formula: ~1 | Rep %in% Ind %in% Tra
        (Intercept)   Residual
StdDev:   0.2242143 0.08325685

Fixed effects: Cre ~ Tra 
Error in if (any(wchLv <- (as.double(levels(xtTab[, wchPval])) == 0))) { : 
        missing value where logical needed
In addition: Warning message: 
NaNs produced in: pt(q, df, lower.tail, log.p) 


	[[alternate HTML version deleted]]


