11.07.2015 Views

statisticalrethinkin..

statisticalrethinkin..

statisticalrethinkin..

SHOW MORE
SHOW LESS
  • No tags were found...

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

11.2. POISSON 293logit(p) ~ a + bdB*deptB + bdC*deptC + bdD*deptD +bdE*deptE + bdF*deptF +bm*male ,a ~ dnorm(0,10) ,bm ~ dnorm(0,10)) ,data=d ,start=list(a=0,bm=0,bdB=0,bdC=0,bdD=0,bdE=0,bdF=0) ,iter=1e4 , warmup=1000 )precis(m10.7stan)Mean StdDev lower 0.95 upper 0.95a 0.68 0.10 0.49 0.88bm -0.10 0.08 -0.27 0.05bdB -0.04 0.11 -0.26 0.16bdC -1.26 0.11 -1.47 -1.05bdD -1.30 0.11 -1.49 -1.08bdE -1.74 0.13 -2.00 -1.49bdF -3.32 0.17 -3.64 -2.99Go ahead and extract the samples and take a look at the pairs plot. You’ll see that again thequadratic approximation has done a very good job.11.1.5. Fitting binomial models with glm. Gives same results as the map approach in previoussection:m0g

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!