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

_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.