Bom dia Pessoal,

          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:

require(gamlss)

mort_syn<-read.table(file="m_syn.txt", header=T)

dd=NULL

dd$y=cbind(ndeath,ntree-ndeath)

dd$y

dd$nhole=nhole

mm=gamlss(formula=y~1+nhole,nu.formula=~1+nhole,sigma.formula=~1,

data=dd,

      family=ZIBB(mu.link="logit",nu.link="logit",sigma.link="log"))

Mais o objeto mm que deveria estimar o parâmetro de Nu esta dando NA, sem erros, sendo:

> mm

Family:  c("ZIBB", "Zero Inflated Beta Binomial")

Fitting method: RS()

Call:  gamlss(formula = y ~ 1 + nhole, sigma.formula = ~1, nu.formula = ~1 +      nhole, family = ZIBB(mu.link = "logit", nu.link = "logit", 

    sigma.link = "log"), data = dd)

Mu Coefficients:

(Intercept)        nhole 

   -2.10530      0.02949 

Sigma Coefficients:

(Intercept) 

     -3.422 

Nu Coefficients:

(Intercept)        nhole 

     -2.367           NA

 Degrees of Freedom for the fit: 4 Residual Deg. of Freedom   104

Global Deviance:     420.808

            AIC:     428.808

            SBC:     439.537

  Alguém que utiliza ou já utilizou o pacote gamlss, saberia me dizer o que esta acontecendo, os dados que utilizei estão anexo,

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