<div dir="ltr"><div class="gmail_default" style="font-family:trebuchet ms,sans-serif"><span style="font-family:courier new,monospace"><span style="font-family:trebuchet ms,sans-serif">Baseado nas discussões, segue um exemplo de análise.</span><br>
</span></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif"><span style="font-family:courier new,monospace"><br>##-----------------------------------------------------------------------------<br><br>
pkg <- c("lattice", "nlme", "doBy")<br>sapply(pkg, require, character.only=TRUE)<br><br>data <- structure(list(N = c(0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 100L, 100L,<br>100L, 100L, 100L, 100L, 100L, 100L), K = c(0L, 0L, 0L, 0L, 50L,<br>
50L, 50L, 50L, 0L, 0L, 0L, 0L, 50L, 50L, 50L, 50L), Tempo = c(10L,<br>10L, 30L, 30L, 10L, 10L, 30L, 30L, 10L, 10L, 30L, 30L, 10L, 10L,<br>30L, 30L), resp = c(24.30206, 32.96457, 36.58875, 42.1957, 43.33543,<br>44.00179, 45.53739, 51.26205, 54.23789, 54.94291, 55.0155, 57.36688,<br>
58.63317, 66.21197, 68.81468, 69.32965)), .Names = c("N", "K",<br>"Tempo", "resp"), class = "data.frame", row.names = c(NA, -16L<br>))<br><br>str(data)<br><br>## Gráfico.<br>
xyplot(resp~Tempo|N, groups=K, data=data, type=c("p","a"))<br><br>## Tabela de frequência.<br>ftable(formula=Tempo~N+K, data=data)<br><br>## Os fatores tem 2 níveis, não requer conveter para fator.<br>
<br>data<br><br>## Não informação de unidades experimentais, então assumir que está na<br>## ordem correta e colocar o indicador de repetição.<br><br>data$rept <- 1:2<br>data$parcela <- with(data, interaction(N, K, rept))<br>
str(data)<br><br>## Especificação e ajuste.<br>m0 <- lme(resp~N*K*Tempo, random=~1|parcela, data=data)<br><br>## Dignóstico simples.<br>r <- residuals(m0)<br>f <- fitted(m0)<br>plot(r~f)<br>plot(sqrt(abs(r))~f)<br>
qqnorm(r)<br><br>## Teste de Wald para os termos de efeito fixo.<br>anova(m0)<br><br>## Quadro de estimativas e medidas de ajuste.<br>summary(m0)<br><br>## Predição.<br>pred <- data.frame(N=c(0,100,50,50,50,50),<br>                   K=c(25,25,0,50,25,25),<br>
                   Tempo=c(20,20,20,20,10,30))<br><br>pred$y <- predict(m0, newdata=pred, level=0)<br>pred<br><br>##-----------------------------------------------------------------------------</span><br><br></div>​<div class="gmail_default" style="font-family:trebuchet ms,sans-serif;display:inline">
​À disposição.<br>Walmes.​</div></div>