x <- read.table("../../true.dat",header=TRUE)
pars <- x$value
names(pars) <- as.character(x$parameter)
load("../../../DATA/Owls.rda") ## check whether this works in sim context ...
owlsim <- Owls
Owls$ArrivalTime <- scale(Owls$ArrivalTime,center=TRUE,scale=FALSE)

X <- model.matrix(~(FoodTreatment+ArrivalTime)*SexParent,data=Owls)
eta0 <- X %*% pars[1:6]        ## fixed effects
eta1 <- eta0 + Owls$logBroodSize  ## add offset
b <- rnorm(length(levels(Owls$Nest)),sd=sqrt(pars["nestvar"]))
eta2 <- eta1 + b[Owls$Nest]       ## add RE
pz <- plogis(pars["logitpz"])
owlsim$SiblingNegotiation <- ifelse(runif(nrow(Owls))<pz,0,
                                 rpois(nrow(Owls),lambda=exp(eta2)))
Owls <- owlsim  ## now overwrite original
save("Owls",file="Owls.rda")

## testing
## Owls <- rename(Owls,c(SiblingNegotiation="NCalls"))
## Owls$ArrivalTime <- scale(Owls$ArrivalTime,center=TRUE,scale=FALSE)
## ofit_glmmadmb <-  glmmadmb(NCalls~(FoodTreatment+ArrivalTime)*SexParent+
##                            offset(logBroodSize),
##                            random=~1|Nest,
##                            data=Owls,
##                            zeroInflation=TRUE,
##                            family="poisson")
## coef(ofit_glmmadmb)
## x