Walmes,
Agora a dúvida é com o uso da função contrast(), apesar de
serem 2 fatores e 2 níveis gostaria de aprender a usar a função para ocasiões
futuras, então após a escolha do modelo reduzido fiz:
> ##Fazendo os contrastes
> require(contrast)
> contr<-contrast(m4,
list(Altitude=levels(as.factor(H.termes$Altitude)),
+
group="pature"),list(Altitude=levels(as.factor(H.termes$Altitude)),group="foret"))
Erro em gendata.default(fit = list(coefficients
= c(-2.17718382645277, :
factor(s) not in design: group
Esta dando esse que não estou compreendendo o significado,
poderia me ajudar,
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: r-br-bounces@listas.c3sl.ufpr.br
[mailto:r-br-bounces@listas.c3sl.ufpr.br] On Behalf Of Walmes Zeviani
Sent: segunda-feira, 25 de abril de 2011 17:58
To: r-br@listas.c3sl.ufpr.br
Subject: Re: [R-br] Dúvida contrastes glm binomial
Alexandre,
Você precisa saber as restrições usadas para obter uma solução do sistema de
equações indeterminado. Mas ao mesmo tempo você não precisa saber porque as funções
estimáveis (contrastes) independem dessas restrições. Nesse caso
Estimate Std. Error z value
Pr(>|z|)
(Intercept)
-2.20627 0.23768 -9.283 < 2e-16 *** (média de
a1 e v1)
Altitudebaixa
1.91926 0.21846 8.785 < 2e-16 ***
(diferença p a2 e v1)
Vegetationpature
-0.83836 0.36768 -2.280 0.02260 *
(diferença p a1 e v2)
Parcelle
0.01484 0.06585 0.225
0.82173 (ef linear de parcela??)
Altitudebaixa:Vegetationpature
-1.10337 0.41007 -2.691 0.00713 **
(diferença pa a2 e v2)
Use a
contrast::contrast() para obter o valor desses contrastes, embora com 2 fatores
de 2 níveis isso seja fácil na mão. Porque "parcelle" entrou como
regressora/contínua?
À 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: walmes@ufpr.br
twitter: @walmeszeviani
homepage: http://www.leg.ufpr.br/~walmes
linux user number: 531218
==========================================================================