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@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 ==========================================================================