Re: [R] negative binomial GAMM with variance structures

2011-09-22 Thread Simon Wood
gamm can't estimate the theta parameter for a negative binomial automatically. It can only work with fixed user supplied values for theta (i.e. negbin(2.3) should work, but negbin(c(1,10)) won't). Is negative binomial the only thing you can use here (it doesn't seem like the most natural choice

[R] negative binomial GAMM with variance structures

2011-09-22 Thread Meredith Jantzen
Hello, I am having some difficulty converting my gam code to a correct gamm code, and I'm really hoping someone will be able to help me. I was previously using this script for my overdispersed gam data: M30 <-gam(efuscus~s(mic, k=7) +temp +s(date)+s(For3k, k=7) + pressure+ humidity, fam