<div><font face="trebuchet ms,sans-serif">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.</font></div>
<div><font face="trebuchet ms,sans-serif"><br></font></div><div><div><font face="'courier new', monospace">sexo <- as.factor(sample(c("M","F","I"), 10, repl=TRUE))</font></div><div>
<font face="'courier new', monospace">recod <- c(-1,0,1)</font></div><div><font face="'courier new', monospace">sexo.rec <- recod[match(sexo, c("F","I","M"))]</font></div>
<div><font face="'courier new', monospace">data.frame(sexo, sexo.rec)</font></div><div><font face="'courier new', monospace"><br></font></div><div><font face="'courier new', monospace">recod <- c(0,0.5,1)</font></div>
<div><font face="'courier new', monospace">sexo.rec <- recod[match(sexo, c("F","I","M"))]</font></div><div><font face="'courier new', monospace">data.frame(sexo, sexo.rec)</font></div>
</div><div><font face="trebuchet ms,sans-serif"><br></font></div><div><font face="trebuchet ms,sans-serif">A interpretação dos coeficientes, bem, teu preditor linear seria</font></div><div><font face="trebuchet ms,sans-serif"><br>
</font></div><div><font face="'courier new', monospace">b0+b1*sexo,</font></div><div><font face="trebuchet ms,sans-serif"><br></font></div><div><font face="trebuchet ms,sans-serif">sexo só assume valores (-1,0,1), então você teria as configurações</font></div>
<div><font face="trebuchet ms,sans-serif"><br></font></div><div><font face="'courier new', monospace">b0-b1  # femea</font></div><div><font face="'courier new', monospace">b0     # indeterminado</font></div>
<div><font face="'courier new', monospace">b0+b1  # macho</font></div><div><font face="trebuchet ms,sans-serif"><br></font></div><div><font face="trebuchet ms,sans-serif">Mas você pode usar (0, 0.5, 1), b0 é a fêmea, b1 é a diferença para macho.</font></div>
<div><font face="trebuchet ms,sans-serif"><br></font></div><div><font face="trebuchet ms,sans-serif">À disposição.</font></div><div><font face="trebuchet ms,sans-serif">Walmes.</font></div><div><font face="'trebuchet ms', sans-serif"><br>
</font></div><span style="font-family:trebuchet ms,sans-serif">==========================================================================</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">Walmes Marques Zeviani</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">Departamento de Estatística - Universidade Federal do Paraná</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">fone: (+55) 41 3361 3573</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">VoIP: (3361 3600) 1053 1173</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">e-mail: <a href="mailto:walmes@ufpr.br" target="_blank">walmes@ufpr.br</a></span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">twitter: @walmeszeviani</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">homepage: <a href="http://www.leg.ufpr.br/%7Ewalmes" target="_blank">http://www.leg.ufpr.br/~walmes</a></span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">linux user number: 531218</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">==========================================================================</span><br>