Simplify your model. Does your TandemRepeat have a lot of levels? Or is
your sample size very small?
Alain
Dear all,
I am
fitting a hurdle model in the following way:
HNB <-
hurdle(chro ~ as.factor(TandemRepeat)| as.factor(TandemRepeat), data
=data_negbin_fin,
dist = "negbin")
But the std.
error for log(theta) = NA
Count model
coefficients (truncated negbin with log link):
Estimate Std. Error z value
Pr(>|z|)
Log(theta) 13.5062 NA NA NA
And it
gives the following error:
In
sqrt(diag(vc_count)[kx + 1]) : NaNs
Can
somebody help me please.
Thanks so
much
Ana
--
Dr. Alain F. Zuur
First author of:
1. Analysing Ecological Data (2007).
Zuur, AF, Ieno, EN and Smith, GM. Springer. 680 p.
URL: www.springer.com/0-387-45967-7
2. Mixed effects models and extensions in ecology with R. (2009).
Zuur, AF, Ieno, EN, Walker, N, Saveliev, AA, and Smith, GM. Springer.
http://www.springer.com/life+sci/ecology/book/978-0-387-87457-9
3. A Beginner's Guide to R (2009).
Zuur, AF, Ieno, EN, Meesters, EHWG. Springer
http://www.springer.com/statistics/computational/book/978-0-387-93836-3
4. Zero Inflated Models and Generalized Linear Mixed Models with R. (2012) Zuur,
Saveliev, Ieno.
http://www.highstat.com/book4.htm
Other books: http://www.highstat.com/books.htm
Statistical consultancy, courses, data analysis and software
Highland Statistics Ltd.
6 Laverock road
UK - AB41 6FN Newburgh
Tel: 0044 1358 788177
Email: highstat@highstat.com
URL: www.highstat.com
URL: www.brodgar.com
[[alternative HTML version deleted]]
