<div dir="ltr"><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">Isso é decorrente do que mencionei na mensagem anterior. Você só tem plantio direto em um lugar. O efeito está confundido. Por isso o rank da matriz é deficiente. Isso ocorre com toda classe de modelo pois de fato não existe informação nos dados que permita identificar os parâmetros em questão. Na lm() um NA aparece nas estimativas dos coeficientes.<br>
<br><span style="font-family:courier new,monospace">> da <- data.frame(local=gl(4,5), sistema="PC", stringsAsFactors=FALSE)<br>> da$sistema[da$local=="4"] <- "PD"<br>> da <- as.data.frame(da, stringsAsFactors=TRUE)<br>
> da$y <- rnorm(nrow(da))<br>> str(da)<br>'data.frame':    20 obs. of  3 variables:<br> $ local  : Factor w/ 4 levels "1","2","3","4": 1 1 1 1 1 2 2 2 2 2 ...<br> $ sistema: chr  "PC" "PC" "PC" "PC" ...<br>
 $ y      : num  0.375 -0.646 -1.727 -1.764 -0.285 ...<br>> xtabs(~local+sistema, da)<br>     sistema<br>local PC PD<br>    1  5  0<br>    2  5  0<br>    3  5  0<br>    4  0  <b>5</b><br>> m0 <- lm(y~local+sistema, data=da)<br>
> coef(m0)<br>(Intercept)      local2      local3      local4   <b>sistemaPD</b> <br> -0.8091214   0.6408831   0.5092444   0.2868672          <b>NA </b><br>> <br></span><br>À disposição.<br>Walmes.<br><br></div></div>