[R-br] fatorial_duplo_regressao

Quarta Março 12 17:34:27 BRT 2014

neste link o Walmes da algumas dicas sobre desdobramento de interações: 
http://ridiculas.wordpress.com/2011/05/21/fatorial-com-desdobramento-da-interacao-grafico-e-tabela-de-medias/.

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 em yahoo.com.br
> <mailto:camposagro em 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 em listas.c3sl.ufpr.br <mailto:R-br em 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 em 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.