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]]

