[R-br] fatorial_duplo_regressao

[Odirley Campos]

Srs(as), boa noite.
Sou iniciante no R. Como faço para ajustar um modelo de 2º grau para os dados abaixo. Antecipademante agradeco

a<-read.table("dose_fcr_ph2.txt",
h=T)
names(a)
tra<-as.factor(a$trat) #usado apenas
para verificar pressupostos da anova
fcr<-as.factor(a$fcr) #4 niveis
(1500, 2000, 2500 e 3000)
dos<-as.factor(a$dose) #3 niveis (0,
1, 2), tratado como qualitativo
ph<-as.numeric(a$ph)
b<-data.frame(tra,fcr,dos,ph)
b
>
fcr
 [1] 1500 1500 1500 2000 2000 2000 2500 2500
2500 3000 3000 3000 1500 1500 1500
[16] 2000 2000
2000 2500 2500 2500 3000 3000 3000 1500 1500 1500 2000 2000 2000
[31] 2500 2500
2500 3000 3000 3000
>
dos
 [1] 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1
1 1 1 2 2 2 2 2 2 2 2 2 2 2 2
>
ph
 [1] 5.86 6.19 6.87 5.87 6.30 6.30 6.05 6.37
7.03 5.69 5.72 6.52 6.05 6.02 6.09
[16] 6.19 6.73
6.80 7.08 7.08 7.08 9.25 8.81 9.06 9.25 9.24 8.81 9.13 8.97 9.36
[31] 9.05 8.99
8.93 9.21 8.90 7.42
ane2<-aov(ph~dos*fcr,
data=b)
summary(ane2)#tenho interesse no
desdobramento
#            Df
Sum Sq Mean Sq F value   Pr(>F)    
#dos          2  45.25  22.625 137.133 7.37e-14 ***
#fcr          3   2.44   0.812   4.921  0.00836 ** 
#dos:fcr      6  14.08   2.347  14.223 7.24e-07 ***
#Residuals   24   3.96   0.165                     
ane3<-aov(ph~dos/fcr,
data=b)#guarda efeitos aninhados
summary(ane3,
split=list("dos:fcr"=list("dos0"=c(1,4,7), "dos1"=c(2,5,8),
"dos2"=c(3,6,9))))
#              Df  Sum Sq     Meam Sq     F value   Pr(>F)    
#dos            2    45.10      22.550    43.910    9.56e-09 ***
#dos:fcr        9    8.69        0.966    1.881     0.10465    
#dos:fcr:
dos0  3    0.42        0.140    0.272     0.84507    
#dos:fcr:
dos1  3    7.50        2.499    4.866     0.00877 ** 
#dos:fcr:
dos2  3    0.78        0.259    0.503     0.68359    
#Residuals       24   12.33      0.514                      
#Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
#---------------------------------------------------------

#Ajustando modelo de 2 grau
> m1<-aov(ph~-1+ dos/(fcr+I(fcr^2)), data=b)# aparece a mesagem de erro abaixo 
#Erro em `contrasts<-`(`*tmp*`, value = contr.funs[1 + isOF[nn]]) : 
#contrastes podem ser aplicados apenas a fatores com 2 ou mais níveis
#Além disso: Mensagens de aviso perdidas:
#In Ops.factor(fcr, 2) : ^ not meaningful for factors

 # está faltando alguma informação no "comando"?



Odirley R. Campos
Engenheiro Agrônomo UFV/MG
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140311/b0f76dd9/attachment-0001.html>