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
"parc" "ntree" "time" "ndeath" "nhole"
1 20 4 0 9
2 25 4 5 14
3 24 4 5 36
4 21 4 4 21
5 14 4 0 1
6 19 4 0 5
7 21 4 1 5
8 11 4 1 3
9 23 4 3 35
10 24 4 2 26
11 18 4 3 27
12 18 4 0 11
13 21 4 4 4
14 10 4 6 3
15 24 4 5 7
16 26 4 8 16
17 19 4 7 39
18 10 4 4 7
19 18 4 2 2
20 19 4 4 1
21 23 4 2 3
22 17 4 1 5
23 21 4 5 21
24 17 4 4 8
25 11 4 3 5
26 13 4 6 2
27 11 4 5 1
28 16 4 1 8
29 13 4 4 7
30 24 4 3 5
31 19 4 3 2
32 24 4 4 16
33 11 4 3 3
34 13 4 1 2
35 25 4 5 5
36 24 4 3 2
37 19 4 2 3
38 23 4 5 2
39 12 4 0 0
40 21 4 4 2
41 22 4 1 3
42 24 4 5 26
43 20 4 5 16
44 25 4 1 5
45 24 4 3 2
46 20 4 3 2
47 21 4 0 1
48 20 4 3 3
49 25 4 2 1
50 27 4 0 1
51 19 4 3 3
52 17 4 0 1
53 24 4 1 2
54 21 4 0 2
55 22 4 3 2
56 13 4 6 3
57 27 4 5 5
58 17 4 2 1
59 25 4 4 3
60 31 4 0 2
61 31 4 1 1
62 33 4 1 1
63 35 4 3 1
64 20 4 2 1
65 26 4 1 2
66 35 4 6 14
67 31 4 6 5
68 32 4 4 4
69 32 4 4 2
70 34 4 1 4
71 35 4 2 3
72 24 4 2 2
73 36 4 4 3
74 33 4 7 1
75 33 4 5 4
76 32 4 0 2
77 33 4 2 1
78 31 4 1 1
79 28 4 1 1
80 30 4 0 0
81 25 4 0 0
82 29 4 6 4
83 28 4 3 2
84 29 4 4 3
85 29 4 2 4
86 31 4 5 6
87 33 4 4 2
88 18 4 5 3
89 20 4 0 8
90 25 4 2 6
91 30 4 0 0
92 28 4 0 0
93 31 4 1 1
94 33 4 1 1
95 28 4 5 2
96 25 4 3 2
97 28 4 1 1
98 30 4 2 1
99 20 4 0 0
100 22 4 1 1
101 29 4 2 1
102 34 4 0 0
103 33 4 0 0
104 39 4 3 2
105 28 4 3 4
106 31 4 3 2
107 30 4 0 1
108 30 4 6 3
_______________________________________________
R-br mailing list
[email protected]
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.
