[R-sig-ME] mcmcsamp with Poisson distribution
bbolker at gmail.com
Sun Sep 9 23:03:50 CEST 2012
Nate Fronk <nrf5017 at ...> writes:
> I am using a Before-After Control-Impact study
> design to look at the effect of natural gas activity on bird species.The
> interaction in my equation tests for the effect of well pads being placed in
> the forest.
> I'm aware I
> could take the Bayesian route but I just wanted to see if there may be some
> code out there somewhere to get around this. Below is an example of what my
> data looks like and my code. Any input would be greatly appreciated.
> Block Pads Time Abundance
> 34B25 2 1 8
> 35C74 4 1 15
> Revi <- read.table("revibaci2.txt", header = TRUE)
The following should both be unnecessary, although they don't hurt ...
mcmcsamp has never (as far as I know) worked for GLMMs (i.e.
family not equal to the default "gaussian")
> mms1<- mcmcsamp(mm1,1000)
An alternative to try:
mm2 <- glmmadmb(Abundance~Pads*Time+(1|Block),
More information about the R-sig-mixed-models