[R-sig-ME] confusing error message
bbolker at gmail.com
Mon Aug 6 06:19:22 CEST 2012
Victoria Buller <Victoria.Buller at ...> writes:
> Can someone explain to me what this error message means? -
> () : Downdated VtV is not positive definite
> I keep getting it after using this code for my data-
> > m1<-glmer(a.martiensseni~distance.tree+herbangular.trans+
> I get that it is to do with my random effects... but I don't
> understand why there is a problem. 'Site' is just an area known as
> A,B,C,D or E and location is the location within the site the frog
> was found, so either transect 1,2, or 3. I am trying to analyse frog
> presence or absence data with habitat characteristics... can anyone
Site is just on the edge of having enough levels to estimate
a random effect reliably, and transect has not enough (the rule
of thumb is a *minimum* of 5-6: see e.g. http://glmm.wikidot.com/faq ).
If you have three transects within each site, I suspect you want
(1|site/location) or (equivalently) (1|site)+(1|site/location).
I hope you have a pretty big data set: one rule of thumb is
that you need the effective size of the data set to about
10 times the number of parameters; for binary data 'effective size'
is (approximately) the minimum of the number of presences and
the number of absences, so the minimum of # presences/# absences
should be at least 200 or so ...
Frank Harrell's book is a good reference for what to do
if you have more predictors than your data can support.
More information about the R-sig-mixed-models