Skip to content

Instantly share code, notes, and snippets.

Embed
What would you like to do?
rlogis, R and JAGS
x = rlogis(10000,0,10)
mod = "model{
y ~ dlogis(0, 1/10)
}
"
library(rjags)
m = rjags::jags.model(file = textConnection(mod))
s = rjags::coda.samples(m,n.iter = 10000,variable.names = "y")
x1 = s[[1]][,"y"]
qqplot(x1,x)
abline(0,1)
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment