[R-br] Dúvida contrast() glm binomial

[Alexandre dos Santos]

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="fo
ret"))

 

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 em listas.c3sl.ufpr.br
[mailto:r-br-bounces em listas.c3sl.ufpr.br] On Behalf Of Walmes Zeviani
Sent: segunda-feira, 25 de abril de 2011 17:58
To: r-br em 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 em ufpr.br
twitter: @walmeszeviani
homepage: http://www.leg.ufpr.br/~walmes <http://www.leg.ufpr.br/%7Ewalmes> 
linux user number: 531218
==========================================================================

-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110426/775f4011/attachment.html>