[R] Newbie
Stefan Grosse
singularitaet at gmx.net
Thu Aug 16 17:12:33 CEST 2007
> I'm a bit new to the world of R so forgive my ignorance.
That has nothing to do with the knowledge of R but with the model. Age
has only 2 different values: 0 and 1 and if it is 0 there is no scars,
so what exactly have you expected from the model?
I would say if you just want to prove that older deer have more scars
try the Mann Whitney non parametric test...
Stefan
> 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 at 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.
>
>
More information about the R-help
mailing list