Obrigado Prof. Paulo e Walmes deu certo,
Respondendo a pergunta Walmes, já tinha retirado a Parcelle do
modelo esta errado mesmo, e agora declarei as.factor() para as variáveis
antes da glm, obrigado facilitou bem as coisas, ficando:
> m4
Call: glm(formula = Infestacao ~ Altitude + Vegetation +
Altitude:Vegetation, family = "binomial", data = H.termes)
Coefficients:
(Intercept)
Altitudebaixa Vegetationpature
-2.1772
1.9240 -0.8366
Altitudebaixa:Vegetationpature
-1.1144
Degrees of Freedom: 1532 Total (i.e. Null); 1529 Residual
Null Deviance: 1579
Residual Deviance: 1323 AIC: 1331
> ##Fazendo os contrastes
> require(contrast)
> contr<-contrast(m4,
list(Altitude=levels(H.termes$Altitude),Vegetation="pature"),
+ list(Altitude=levels(H.termes$Altitude),
Vegetation="foret"),type="average")
> print(contr,X=TRUE)
glm model parameter contrast
Contrast S.E. Lower Upper t df Pr(>|t|)
1 -1.393758 0.2035796 -1.793083 -0.9944332 -6.85 1529 0
Contrast coefficients:
(Intercept) Altitudebaixa Vegetationpature Altitudebaixa:Vegetationpature
1 0 0 1 0.5
>
Obrigado,
Alexandre dos Santos
Ingenieur forestier, Msc.
INRA- Biostatistique et Processus Spatiaux (BioSP)
Domaine Saint-Paul
Site Agroparc
84914 - Avignon - France
Tél. : +33 (0)6 87 95 16 29
From: [email protected]
[mailto:[email protected]] On Behalf Of Walmes Zeviani
Sent: terça-feira, 26 de abril de 2011 09:11
To: [email protected]
Subject: Re: [R-br] Dúvida contrast() glm binomial
Alexandre,
Recomendo que passe as variáveis para fator antes da glm(), isso para evitar
que elas sejam convertidas para fator dentro da função e você tenha que
ficar usando as.factor() dentro da contrast(). O fator "group" não existe
nas análises que você postou o summary aqui, então por isso que a função
retorna "group não faz parte do delineamento", porque ela não é termo
declarado na formula do modelo, apenas foi declarado "Altitute",
"Vegetation" e "Parcelle". Parcelle não é uma variável categórica? Você usou
ela como variável regressora (contínua) no modelo. Isso está correto?
À disposição.
Walmes.
==========================================================================
Walmes Marques Zeviani
LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)
Departamento de Estatística - Universidade Federal do Paraná
fone: (+55) 41 3361 3573
VoIP: (3361 3600) 1053 1173
e-mail: [email protected]
twitter: @walmeszeviani
homepage: http://www.leg.ufpr.br/~walmes <http://www.leg.ufpr.br/%7Ewalmes>
linux user number: 531218
==========================================================================
_______________________________________________
R-br mailing list
[email protected]
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
