Thank you for your quick reply! 
 
When I run the code you provide I get this output:
 
LCOVER
LOCS      1      2      3      4      5      6      7      9
   0 214507  79939  69803 778359  22932  32391  99630   8082
   1     15      7      1     32      0      0      0      0
   2      2      1      0      0      0      0      0      0
   3      0      0      0      1      0      0      0      0
 
Nate
 
 
 

________________________________

From: David Winsemius [mailto:dwinsemius@comcast.net]
Sent: Thu 5/24/2012 1:54 PM
To: Nathan Svoboda
Cc: r-help@r-project.org
Subject: Re: [R] R Error: System is computationally singular




On May 24, 2012, at 1:57 PM, Nathan Svoboda wrote:

> Greetings,
>
> I am trying to fit a zero-inflated Poisson model using zeroinfl() 
> from the
> pscl library. I have 5 covariates (4 continuous, 1 categorical); the
> categorical variable has 7 levels.  I have had success fitting 
> models that
> contain only the continuous covariates; however, when I add the 
> categorical
> variable to any of the models (or if I run it by itself) I get the 
> following
> error:
>
> Error in solve.default(as.matrix(fit$hessian)) :
>
>  system is computationally singular: reciprocal condition number =
> 3.46934e-20
>
> The code I am using is:
>
> library(pscl)
> f1 <- formula(LOCS ~ as.factor(LCOVER) + D_ROADS + D_WATER + D_EDGE +
> D_GRASS)
> ZIP1 <- zeroinfl(f1, dist="poisson", link = "logit", data = FAWNS)
>
> There is no correlation between my covariates. Also, I tried 
> reducing my
> categorical covariate to 3 levels and still receive the same error. 
> Can
> anyone suggest why I may be getting this error when I add the 
> categorical
> covariate?
>

What does this show:

with( FAWNS, table(LOCS, LCOVER) )

--
David Winsemius, MD
West Hartford, CT




	[[alternative HTML version deleted]]

