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 NANA 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: highs...@highstat.com
URL: www.highstat.com
URL: www.brodgar.com
[[alternative HTML version deleted]]
__
R-help@r-project.org 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.