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: [email protected] [mailto:[email protected]] On Behalf Of Walmes Zeviani Sent: segunda-feira, 25 de abril de 2011 17:58 To: [email protected] 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: [email protected] twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
_______________________________________________ R-br mailing list [email protected] https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
