diff options
author | Matth?us G. Chajdas <dev@anteru.net> | 2019-11-10 13:56:53 +0100 |
---|---|---|
committer | Matth?us G. Chajdas <dev@anteru.net> | 2019-11-10 13:56:53 +0100 |
commit | 1dd3124a9770e11b6684e5dd1e6bc15a0aa3bc67 (patch) | |
tree | 87a171383266dd1f64196589af081bc2f8e497c3 /tests/examplefiles/example.jag | |
parent | f1c080e184dc1bbc36eaa7cd729ff3a499de568a (diff) | |
download | pygments-master.tar.gz |
Diffstat (limited to 'tests/examplefiles/example.jag')
-rw-r--r-- | tests/examplefiles/example.jag | 48 |
1 files changed, 0 insertions, 48 deletions
diff --git a/tests/examplefiles/example.jag b/tests/examplefiles/example.jag deleted file mode 100644 index cae34beb..00000000 --- a/tests/examplefiles/example.jag +++ /dev/null @@ -1,48 +0,0 @@ -# lsat.jags example from classic-bugs examples in JAGS -# See http://sourceforge.net/projects/mcmc-jags/files/Examples/2.x/ -var - response[R,T], m[R], culm[R], alpha[T], a[T], theta[N], r[N,T], - p[N,T], beta, theta.new, p.theta[T], p.item[R,T], P.theta[R]; -data { - for (j in 1:culm[1]) { - r[j, ] <- response[1, ]; - } - for (i in 2:R) { - for (j in (culm[i - 1] + 1):culm[i]) { - r[j, ] <- response[i, ]; - } - } -} -model { - # 2-parameter Rasch model - for (j in 1:N) { - for (k in 1:T) { - probit(p[j,k]) <- delta[k]*theta[j] - eta[k]; - r[j,k] ~ dbern(p[j,k]); - } - theta[j] ~ dnorm(0,1); - } - - # Priors - for (k in 1:T) { - eta[k] ~ dnorm(0,0.0001); - e[k] <- eta[k] - mean(eta[]); # sum-to-zero constraint - - delta[k] ~ dnorm(0,1) T(0,); # constrain variance to 1, slope +ve - d[k] <- delta[k]/pow(prod(delta), 1/T); # PRODUCT_k (d_k) = 1 - - g[k] <- e[k]/d[k]; # equivalent to B&A's threshold parameters - } - - # Compute probability of response pattern i, for later use in computing G^2 - theta.new ~ dnorm(0,1); # ability parameter for random student - for(k in 1:T) { - probit(p.theta[k]) <- delta[k]*theta.new - eta[k]; - for(i in 1:R) { - p.item[i,k] <- p.theta[k]^response[i,k] * (1-p.theta[k])^(1-response[i,k]); - } - } - for(i in 1:R) { - P.theta[i] <- prod(p.item[i,]) - } -} |