[R-br] fatorial_duplo_regressao

[Wagner Pinheiro]

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>:

> 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 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
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140312/904d4ea3/attachment.html>