Hello,
I'm a bit new to the world of R so forgive my ignorance. I'm trying to do a
zero-inflated negative binomial regression and have received an error message
and i'm not sure what it means. I'm running R 2.5.1 on XP. I have just tried
a really simple version of the model to see if it would run before I put all
the variables in. I have attached all the variables to the object alan.
Here is the table:
Date Location Deer Code Sex Age Pel Per Lac Scars Emb Lar Nym Bee Mi Mass
1 12/06/2007 Ballysallagh 0 B1 FALSE 1 A 1 1 4 0 0 0
0 0 21.59
2 15/06/2007 Ballysallagh 0 B10 FALSE 1 A 1 1 5 0 0 0
0 2 26.19
3 15/06/2007 Ballysallagh 0 B12 FALSE 1 A 1 1 0 5 0 0
0 1 28.55
4 15/06/2007 Ballysallagh 0 B13 FALSE 1 A 0 1 5 0 0 0
4 45 23.93
5 15/06/2007 Ballysallagh 0 B16 FALSE 1 A 1 1 0 4 0 0
0 1 34.19
6 15/06/2007 Ballysallagh 0 B17 FALSE 1 A 1 1 5 0 0 0
0 0 25.02
7 15/06/2007 Ballysallagh 0 B18 FALSE 1 A 1 1 5 0 0 0
0 7 33.06
8 12/06/2007 Ballysallagh 0 B5 FALSE 1 A 1 1 4 0 0 0
0 2 23.55
9 12/06/2007 Ballysallagh 0 B6 FALSE 1 A 1 1 4 0 0 0
0 2 22.67
10 12/06/2007 Ballysallagh 0 B7 FALSE 1 A 1 1 5 0 0 0
0 2 22.57
11 12/06/2007 Ballysallagh 0 B8 FALSE 1 A 1 1 5 0 0 0
0 10 24.01
12 24/06/2007 Caledon 1 C12 FALSE 1 A 1 0 0 0 1 0
0 3 21.68
13 25/06/2007 Caledon 1 C17 FALSE 1 A 1 0 4 0 1 0
0 23 20.56
14 25/06/2007 Caledon 1 C18 FALSE 1 A 1 0 0 0 0 0
2 2 17.47
15 25/06/2007 Caledon 1 C20 FALSE 1 A 1 0 8 0 0 0
0 3 19.97
16 25/06/2007 Caledon 1 C21 FALSE 1 A 1 0 3 0 2 0
0 20 20.58
17 25/06/2007 Caledon 1 C24 FALSE 1 A 1 0 0 4 2 0
0 7 17.37
18 25/06/2007 Caledon 1 C27 FALSE 1 A 1 0 0 0 3 0
0 6 24.14
19 25/06/2007 Caledon 1 C28 FALSE 1 A 1 0 0 5 1 0
0 6 21.58
20 25/06/2007 Caledon 1 C33 FALSE 1 A 0 0 0 0 1 0
0 0 19.06
21 25/06/2007 Caledon 1 C35 FALSE 1 A 1 0 3 0 2 0
0 20 24.55
22 24/06/2007 Caledon 1 C4 FALSE 1 A 1 1 4 3 0 0
0 5 22.65
23 24/06/2007 Caledon 1 C8 FALSE 1 A 0 0 0 0 3 0
0 3 23.08
24 01/06/2007 Hillsborough 0 H11 FALSE 1 A 1 1 4 0 0 0
0 0 23.01
25 01/06/2007 Hillsborough 0 H16 FALSE 1 A 1 1 5 0 0 0
0 0 18.18
26 01/06/2007 Hillsborough 0 H17 FALSE 1 A 1 1 4 0 0 0
0 0 24.45
27 01/06/2007 Hillsborough 0 H19 FALSE 1 A 0 1 0 0 0 0
0 3 21.01
28 01/06/2007 Hillsborough 0 H2 FALSE 1 A 1 1 5 0 0 0
0 1 24.86
29 01/06/2007 Hillsborough 0 H21 FALSE 1 A 1 1 5 0 0 0
0 0 20.48
30 01/06/2007 Hillsborough 0 H25 FALSE 1 A 1 1 5 0 0 0
0 0 22.09
31 01/06/2007 Hillsborough 0 H26 FALSE 1 A 1 1 5 0 0 0
0 0 20.18
32 01/06/2007 Hillsborough 0 H3 FALSE 1 A 0 1 5 0 0 0
0 0 25.82
33 01/06/2007 Hillsborough 0 H34 FALSE 1 A 1 0 0 4 0 0
0 1 18.78
34 01/06/2007 Hillsborough 0 H35 FALSE 1 A 0 1 4 0 0 0
0 0 20.02
35 01/06/2007 Hillsborough 0 H5 FALSE 1 A 0 1 5 0 0 0
0 0 19.58
36 01/06/2007 Hillsborough 0 H6 FALSE 1 A 1 1 5 0 0 0
0 6 28.84
37 16/06/2007 LoughNavar 1 LN1 FALSE 1 A 1 1 5 0 3 0
0 1 22.71
38 16/06/2007 LoughNavar 1 LN10 FALSE 1 A 1 1 5 0 0 0
0 0 23.53
39 16/06/2007 LoughNavar 1 LN11 FALSE 1 A 1 1 4 0 0 0
0 2 33.24
40 16/06/2007 LoughNavar 1 LN13 FALSE 1 A 1 1 5 0 0 0
0 0 20.20
41 16/06/2007 LoughNavar 1 LN14 FALSE 1 A 1 1 0 7 0 0
1 1 31.35
42 16/06/2007 LoughNavar 1 LN16 FALSE 1 A 1 1 0 5 0 0
0 2 19.43
43 16/06/2007 LoughNavar 1 LN18 FALSE 1 A 1 0 6 0 0 0
0 12 24.30
44 16/06/2007 LoughNavar 1 LN23 FALSE 1 A 1 1 5 0 0 0
0 1 24.34
45 16/06/2007 LoughNavar 1 LN7 FALSE 1 A 1 1 5 0 1 0
0 0 28.01
46 16/06/2007 LoughNavar 1 LN9 FALSE 1 A 1 1 4 0 8 0
0 5 25.43
47 11/07/2007 Mt.Stewart 0 M1 FALSE 1 A 1 0 5 0 0 0
0 1 25.31
48 11/07/2007 Mt.Stewart 0 M18 FALSE 1 A 1 1 0 0 0 0
0 2 26.66
49 11/07/2007 Mt.Stewart 0 M21 FALSE 1 A 0 0 3 0 0 0
0 7 23.70
50 11/07/2007 Mt.Stewart 0 M4 FALSE 1 A 1 0 4 0 0 0
0 1 27.05
51 05/06/2007 Portglenone 0 P10 FALSE 1 A 1 0 5 0 0 0
1 1 21.08
52 05/06/2007 Portglenone 0 P13 FALSE 1 A 1 1 0 5 0 0
0 23 34.07
53 05/06/2007 Portglenone 0 P14 FALSE 1 A 1 1 4 0 0 0
0 5 29.00
54 05/06/2007 Portglenone 0 P17 FALSE 1 A 1 1 6 0 0 0
0 1 23.61
55 11/06/2007 Portglenone 0 P18 FALSE 1 A 1 0 4 0 0 0
1 25 19.76
56 11/06/2007 Portglenone 0 P19 FALSE 1 A 1 0 5 0 0 0
0 25 19.13
57 05/06/2007 Portglenone 0 P5 FALSE 1 A 1 0 4 0 0 0
0 11 19.10
58 05/06/2007 Portglenone 0 P7 FALSE 1 A 1 1 2 0 0 0
0 1 26.02
59 29/05/2007 Randalstown 1 R5 FALSE 1 A 1 1 0 0 1 0
3 0 26.70
60 30/06/2007 Somerset 0 S14 FALSE 1 A 1 0 0 4 0 0
0 3 22.20
61 30/06/2007 Somerset 0 S15 FALSE 1 A 1 1 4 0 0 0
1 11 33.36
62 30/06/2007 Somerset 0 S17 FALSE 1 A 1 0 3 0 0 0
0 1 27.10
63 30/06/2007 Somerset 0 S2 FALSE 1 A 1 0 0 0 0 0
2 4 27.87
64 03/07/2007 Somerset 0 S25 FALSE 1 A 1 0 6 0 0 0
0 2 27.34
65 03/07/2007 Somerset 0 S28 FALSE 1 A 1 0 5 0 0 0
0 4 23.83
66 03/07/2007 Somerset 0 S29 FALSE 1 A 1 0 4 4 0 0
0 4 21.51
67 03/07/2007 Somerset 0 S31 FALSE 1 A 1 0 4 0 0 0
0 17 27.01
68 03/07/2007 Somerset 0 S35 FALSE 1 A 0 0 0 0 0 0
0 2 15.60
69 30/06/2007 Somerset 0 S7 FALSE 1 A 1 1 0 0 0 0
0 5 25.58
70 30/06/2007 Somerset 0 S9 FALSE 1 A 1 1 4 0 0 0
0 3 31.87
71 29/05/2007 Tollymore 1 T13 FALSE 1 A 1 0 0 0 1 0
0 0 20.24
72 29/05/2007 Tollymore 1 T16 FALSE 1 A 0 0 0 4 9 0
0 1 17.43
73 29/05/2007 Tollymore 1 T19 FALSE 1 A 1 0 5 0 3 0
0 0 19.50
74 29/05/2007 Tollymore 1 T23 FALSE 1 A 1 0 0 5 6 0
0 1 18.67
75 29/05/2007 Tollymore 1 T24 FALSE 1 A 1 0 0 0 0 0
0 2 16.55
76 29/05/2007 Tollymore 1 T30 FALSE 1 A 1 0 0 0 2 0
0 2 20.88
77 29/05/2007 Tollymore 1 T4 FALSE 1 A 1 0 0 0 3 0
2 0 16.43
78 29/05/2007 Tollymore 1 T6 FALSE 1 A 1 0 0 5 28 0
4 0 17.50
79 29/05/2007 Tollymore 1 T9 FALSE 1 A 1 0 4 0 3 0
0 0 20.58
80 15/07/2007 Tullychurry 1 TC1 FALSE 1 A 1 0 0 0 9 0
0 0 25.27
81 16/07/2007 Tullychurry 1 TC12 FALSE 1 A 1 0 5 0 44 1
0 4 24.83
82 16/07/2007 Tullychurry 1 TC15 FALSE 1 A 1 0 5 0 5 0
0 14 22.28
83 16/07/2007 Tullychurry 1 TC18 FALSE 1 A/J 0 0 0 0 15 0
9 0 13.87
84 16/07/2007 Tullychurry 1 TC23 FALSE 1 A 1 0 0 4 23 0
0 1 27.00
85 15/07/2007 Tullychurry 1 TC3 FALSE 1 A 0 1 6 0 75 0
0 4 22.43
86 15/07/2007 Tullychurry 1 TC4 FALSE 1 A 1 0 0 5 16 0
0 1 26.73
87 15/07/2007 Tullychurry 1 TC6 FALSE 1 A 1 0 2 0 5 0
0 8 21.28
88 15/07/2007 Tullychurry 1 TC7 FALSE 1 A 1 1 6 0 42 0
1 1 28.55
89 15/07/2007 Tullychurry 1 TC8 FALSE 1 A 1 0 3 0 3 0
0 0 20.13
90 15/06/2007 Ballysallagh 0 B11 FALSE 0 A/J 0 0 0 0 0 0
0 0 13.82
91 24/06/2007 Caledon 1 C16 FALSE 0 J 0 0 0 0 0 0
0 3 8.40
92 20/06/2007 Caledon 1 C2 FALSE 0 J 0 0 0 0 0 0
0 0 8.08
93 16/06/2007 LoughNavar 1 LN2 FALSE 0 A/J 0 0 0 0 0 0
0 2 8.71
94 18/06/2007 LoughNavar 1 LN27 FALSE 0 A 0 0 0 0 0 0
0 1 16.90
95 18/06/2007 LoughNavar 1 LN34 FALSE 0 J 0 0 0 0 0 0
1 2 9.38
96 11/07/2007 Mt.Stewart 0 M10 FALSE 0 J 0 0 0 0 0 0
0 2 12.52
97 11/07/2007 Mt.Stewart 0 M12 FALSE 0 A 0 0 0 0 0 0
0 6 16.18
98 11/07/2007 Mt.Stewart 0 M3 FALSE 0 J 0 0 0 0 0 0
0 5 11.42
99 15/05/2007 Randalstown 1 R1 FALSE 0 A 0 0 0 4 44 0
4 0 15.72
100 11/06/2007 Randalstown 1 R19 FALSE 0 A 0 0 0 0 2 0
0 7 15.20
101 30/06/2007 Somerset 0 S12 FALSE 0 A 0 0 0 0 0 0
1 7 14.04
102 03/07/2007 Somerset 0 S19 FALSE 0 A 0 0 0 0 0 0
0 3 15.00
103 29/05/2007 Tollymore 1 T11 FALSE 0 A 1 0 0 0 2 0
3 0 16.30
104 29/05/2007 Tollymore 1 T28 FALSE 0 A 0 0 0 0 0 0
0 0 14.49
And here is the syntax is used to run the ZINB with the error message, any help
greatly appreciated.
> zinb.zc<- zicounts(resp=Scars~.,x=~Age,z=~Age, data=alan)
Error in solve.default(z0$hessian) : system is computationally singular:
reciprocal condition number = 3.44599e-18
Cheers
Alan
Alan Harrison
Quercus
Queen's University Belfast
MBC, 97 Lisburn Road
Belfast
BT9 7BL
T: 02890 972219
M: 07798615682
[[alternative HTML version deleted]]
______________________________________________
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
