Ficou muito bom assim :) Exatamente o que eu precisava. Obrigado Walmes e Heloise pela atenção. Em 31 de agosto de 2012 22:12, Walmes Zeviani <walmeszeviani@gmail.com>escreveu:
Aqui vai uma dica para eliminar os for() que você usou, pode-se usar o match() para fazer essas recodificações. Deve ter uma forma ainda melhor, mas atualmente eu uso essa.
sexo <- as.factor(sample(c("M","F","I"), 10, repl=TRUE)) recod <- c(-1,0,1) sexo.rec <- recod[match(sexo, c("F","I","M"))] data.frame(sexo, sexo.rec)
recod <- c(0,0.5,1) sexo.rec <- recod[match(sexo, c("F","I","M"))] data.frame(sexo, sexo.rec)
A interpretação dos coeficientes, bem, teu preditor linear seria
b0+b1*sexo,
sexo só assume valores (-1,0,1), então você teria as configurações
b0-b1 # femea b0 # indeterminado b0+b1 # macho
Mas você pode usar (0, 0.5, 1), b0 é a fêmea, b1 é a diferença para macho.
À 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 ==========================================================================
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
-- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org Lattes: http://lattes.cnpq.br/7355685961127056