<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv=Content-Type content="text/html; charset=iso-8859-1">
<meta name=Generator content="Microsoft Word 12 (filtered medium)">
<style>
<!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri","sans-serif";
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;}
@page Section1
{size:612.0pt 792.0pt;
margin:70.85pt 3.0cm 70.85pt 3.0cm;}
div.Section1
{page:Section1;}
-->
</style>
<!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang=PT-BR link=blue vlink=purple>
<div class=Section1>
<p class=MsoNormal> Bom dia Pessoal,<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal> Estou tentando modelizar a variável resposta ndeath
(plantas mortas)em função da variável nhole (numeros de insetos), com o uso de
um modelo hierárquico de Bernoulli e Beta binomial, para isso usei o pacote gamlss,
para estimar os parâmetros sendo:<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal><span lang=EN-US>require(gamlss)<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>mort_syn<-read.table(file="m_syn.txt",
header=T)<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>dd=NULL<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>dd$y=cbind(ndeath,ntree-ndeath)<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>dd$y <o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>dd$nhole=nhole<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US>mm=gamlss(formula=y~1+nhole,nu.formula=~1+nhole,sigma.formula=~1,<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> </span>data=dd,<o:p></o:p></p>
<p class=MsoNormal><span lang=EN-US>
family=ZIBB(mu.link="logit",nu.link="logit",sigma.link="log"))<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US><o:p> </o:p></span></p>
<p class=MsoNormal>Mais o objeto mm que deveria estimar o parâmetro de Nu esta
dando NA, sem erros, sendo:<o:p></o:p></p>
<p class=MsoNormal><span lang=EN-US>> mm<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US>Family: c("ZIBB", "Zero
Inflated Beta Binomial") <o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>Fitting method: RS() <o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US>Call: gamlss(formula = y ~ 1 + nhole,
sigma.formula = ~1, nu.formula = ~1 + nhole, family = ZIBB(mu.link =
"logit", nu.link = "logit", <o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> sigma.link = "log"), data =
dd) <o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US>Mu Coefficients:<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>(Intercept) nhole <o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> -2.10530 0.02949 <o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>Sigma Coefficients:<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>(Intercept) <o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> -3.422 <o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>Nu Coefficients:<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>(Intercept) nhole <o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> -2.367 <b><span
style='color:red'> NA</span></b> <o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US> Degrees of Freedom for the fit: 4 Residual
Deg. of Freedom 104 <o:p></o:p></span></p>
<p class=MsoNormal>Global Deviance: 420.808 <o:p></o:p></p>
<p class=MsoNormal> AIC: 428.808 <o:p></o:p></p>
<p class=MsoNormal> SBC: 439.537<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal> Alguém que utiliza ou já utilizou o pacote gamlss, saberia
me dizer o que esta acontecendo, os dados que utilizei estão anexo,<o:p></o:p></p>
<p class=MsoNormal>Obrigado,<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal><b><span style='font-size:10.0pt;font-family:"Arial","sans-serif";
color:black'>Alexandre dos Santos</span></b><span style='font-size:10.0pt;
color:black'><o:p></o:p></span></p>
<p class=MsoNormal><b><span style='font-size:10.0pt;font-family:"Arial","sans-serif";
color:black'>Ingenieur forestier, Msc.</span></b><span style='font-size:10.0pt;
color:black'><o:p></o:p></span></p>
<p class=MsoNormal><b><span style='font-size:10.0pt;font-family:"Arial","sans-serif";
color:black'>INRA- Biostatistique et Processus Spatiaux (BioSP)</span></b><span
style='font-size:10.0pt;color:black'><o:p></o:p></span></p>
<p class=MsoNormal><b><span lang=EN-US style='font-size:10.0pt;font-family:
"Arial","sans-serif";color:black'>Domaine Saint-Paul<br>
Site Agroparc <br>
84914 - Avignon - France<br>
Tél. : +33 (0)6 87 95 16 29</span></b><span lang=EN-US style='font-size:10.0pt;
color:black'><o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US><o:p> </o:p></span></p>
</div>
</body>
</html>