Guest User

Untitled

a guest
Jan 18th, 2018
80
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.36 KB | None | 0 0
  1. library(lme4)
  2. library(runjags)
  3.  
  4. model <- template.jags(total.fruits ~ reg + amd*gen + (1 | reg), data=Arabidopsis, family = "poisson")
  5.  
  6. # can ignore warning for this demo
  7.  
  8. ml.jags <- run.jags(model)
  9. summary(ml.jags)
  10.  
  11. > results[row.names(results) == "amd_effect[2]", "Mean"] * results[row.names(results) == "gen_coefficient_amd_level[2]", "Mean"]
  12. [1] 0.0001745244
Add Comment
Please, Sign In to add comment