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
Odirley Campos, bom dia. Recomendo que use o pacote ExpDes lá você encontrará o que precisa e tudo isso em apenas uma linha de comando. Acesse o link e confira. http://cran.r-project.org/web/packages/ExpDes/ExpDes.pdf 2014-03-11 22:52 GMT-03:00 Odirley Campos <camposagro@yahoo.com.br>:
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.
-- Wagner Rogério Ferreira Pinheiro Bacharel em Estatística (UFPA) Mestre em Estatística Aplicada e Biometria (UFV) --------------------------------------------------------------------- *Contatos:* Fone: (31) 9258-2495 (31) 7143-5968 Skype: wagnerpinheiro01
neste link o Walmes da algumas dicas sobre desdobramento de interações: http://ridiculas.wordpress.com/2011/05/21/fatorial-com-desdobramento-da-inte.... Quanto ao erro que você obteve, esta acontecendo o seguinte. Fatores servem para identificar tratamentos. Você não pode fazer operações matemáticas com fatores. Eles aceitam somente operações relacionais (<,>;<=,>=,!=). A expreção I(fcr^2), não é permitido na função aov. Quanto ao Erro em `contrasts<-`(`*tmp*`, value = contr.funs[1 + isOF[nn]]) : #contrastes podem ser aplicados apenas a fatores com 2 ou mais níveis; Deve estar relacionado com o valor -1 na expressão m1<-aov(ph~-1+ dos/(fcr+I(fcr^2)), data=b), Você está confundindo o uso dessa função com a função lm{stats} Coloque um CMR (código mínimo REPRODUZÍVEL), para que possa lhe ajudar melhor; Utilize a função dput dos dados necessários para o exemplo para facilitar a execução de seus códigos por outros. att Em Qua 12 Mar 2014 08:36:35 BRT, Wagner Pinheiro escreveu:
Odirley Campos, bom dia.
Recomendo que use o pacote ExpDes lá você encontrará o que precisa e tudo isso em apenas uma linha de comando. Acesse o link e confira.
http://cran.r-project.org/web/packages/ExpDes/ExpDes.pdf
2014-03-11 22:52 GMT-03:00 Odirley Campos <camposagro@yahoo.com.br <mailto:camposagro@yahoo.com.br>>:
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 <tel:%5B16%5D%202000%202000> 2000 2500 2500 2500 3000 3000 3000 1500 1500 1500 2000 2000 2000 [31] 2500 2500 <tel:%5B31%5D%202500%202500> 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 valuePr(>F) #dos245.2522.625 137.133 7.37e-14 *** #fcr32.440.8124.9210.00836 ** #dos:fcr614.082.34714.223 7.24e-07 *** #Residuals243.960.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)))) #DfSum SqMeam SqF valuePr(>F) #dos245.1022.55043.910 9.56e-09 *** #dos:fcr98.690.9661.8810.10465 #dos:fcr: dos030.420.1400.2720.84507 #dos:fcr: dos137.502.4994.8660.00877 ** #dos:fcr: dos230.780.2590.503 0.68359 #Residuals2412.330.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 <mailto: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.
-- Wagner Rogério Ferreira Pinheiro
Bacharel em Estatística (UFPA) Mestre em Estatística Aplicada e Biometria (UFV) --------------------------------------------------------------------- *Contatos:* Fone: (31) 9258-2495 (31) 7143-5968 Skype: wagnerpinheiro01
_______________________________________________ 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.
participantes (3)
-
Fernando Souza -
Odirley Campos -
Wagner Pinheiro