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

[Paulo Justiniano]

Veja a mensagem de erro...
Voce está passando uma string para o argumento group
e nao um objeto com um fator


On Tue, 26 Apr 2011, Alexandre dos Santos wrote:

> 
> 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 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
> linux user number: 531218
> ==========================================================================
> 
> 
>