From b752a21c8c5a98cf532e51bb6e8862c298ecc0cc Mon Sep 17 00:00:00 2001
From: corentinchoisy <corentin.choisy@proton.me>
Date: Thu, 1 May 2025 01:10:31 +0200
Subject: [PATCH] Modified JAGS file architecture

---
 src/bpcm      | 23 -----------------------
 src/bpcm_beta | 25 -------------------------
 src/bpcm_dif  | 46 ----------------------------------------------
 3 files changed, 94 deletions(-)
 delete mode 100644 src/bpcm
 delete mode 100644 src/bpcm_beta
 delete mode 100644 src/bpcm_dif

diff --git a/src/bpcm b/src/bpcm
deleted file mode 100644
index 8706406..0000000
--- a/src/bpcm
+++ /dev/null
@@ -1,23 +0,0 @@
-model {
-    for (i in 1:n) {
-        for (j in 1:p) {
-                Y[i,j] ~ dcat(prob[i,j,1:K[j]])
-                }
-                theta[i] ~ dnorm(0,1)
-    }
-    for (i in 1:n) {
-      for (j in 1:p) {
-            for (k in 1:K[j] ) {
-                eta[i,j,k] <- theta[i] - delta[j,k]
-                psum[i,j,k] <- sum(eta[i,j,1:k])
-                exp.psum[i,j,k] <- exp(psum[i,j,k])
-                prob[i,j,k] <- exp.psum[i,j,k] / sum(exp.psum[i,j,1:K[j]])
-    } } }
-    for (j in 1:p) {
-        delta[j,1] <- 0.0
-        for (k in 2:K[j]) {
-                delta[j,k] ~ dnorm(m.delta, pr.delta)
-        }
-    }
-    pr.delta <- pow(s.delta, -2)
-}
diff --git a/src/bpcm_beta b/src/bpcm_beta
deleted file mode 100644
index ffc441e..0000000
--- a/src/bpcm_beta
+++ /dev/null
@@ -1,25 +0,0 @@
-model {
-    for (i in 1:n) {
-        for (j in 1:p) {
-                Y[i,j] ~ dcat(prob[i,j,1:K[j]])
-                }
-                theta[i] ~ dnorm(0,1)
-    }
-    for (i in 1:n) {
-      for (j in 1:p) {
-            for (k in 1:K[j] ) {
-                eta[i,j,k] <- theta[i] + (beta * Z[i,1]) - delta[j,k]
-                psum[i,j,k] <- sum(eta[i,j,1:k])
-                exp.psum[i,j,k] <- exp(psum[i,j,k])
-                prob[i,j,k] <- exp.psum[i,j,k] / sum(exp.psum[i,j,1:K[j]])
-    } } }
-    for (j in 1:p) {
-        delta[j,1] <- 0.0
-        for (k in 2:K[j]) {
-                delta[j,k] ~ dnorm(m.delta, pr.delta)
-        }
-    }
-    beta ~ dnorm(m.beta,pr.beta)
-    pr.delta <- pow(s.delta, -2)
-    pr.beta <- pow(s.beta,-2)
-}
diff --git a/src/bpcm_dif b/src/bpcm_dif
deleted file mode 100644
index 14a13c2..0000000
--- a/src/bpcm_dif
+++ /dev/null
@@ -1,46 +0,0 @@
-model {
-    for (i in 1:n) {
-        for (j in 1:p) {
-                Y[i,j] ~ dcat(prob[i,j,1:K[j]])
-                }
-                theta[i] ~ dnorm(0, 1)
-    }
-    for (i in 1:n) {
-      for (j in 1:p) {
-            for (k in 1:K[j] ) {
-                eta[i,j,k] <- theta[i] + (beta * Z[i,1]) - ( delta[j,k] + Z[i,1]*ifelse(k==1,0,gamma[j,k])  )
-                psum[i,j,k] <- sum(eta[i,j,1:k])
-                exp.psum[i,j,k] <- exp(psum[i,j,k])
-                prob[i,j,k] <- exp.psum[i,j,k] / sum(exp.psum[i,j,1:K[j]])
-    } } }
-    for (j in 1:pnodif) {
-        delta[j,1] <- 0.0
-        gamma[j,1] <- 0.0
-        for (k in 2:K[j]) {
-                delta[j,k] ~ dnorm(m.delta, pr.delta)
-                gamma[j,k] <- 0.0
-        }
-    }
-    for (j in pnodif1:pnounif) {
-        delta[j,1] <- 0.0
-        gamma[j,1] <- 0.0
-        for (k in 2:K[j]) {
-                delta[j,k] ~ dnorm(m.delta, pr.delta)
-                gamma[j,k] ~ dnorm(m.gamma,s.gamma)
-        }
-    }
-    for (j in pnounif1:p) {
-        delta[j,1] <- 0.0
-        gamma[j,1] <- 0.0
-        gamma[j,2] ~ dnorm(m.gamma,pr.gamma)
-        delta[j,2] ~ dnorm(m.delta, pr.delta)
-        for (k in 3:K[j]) {
-                delta[j,k] ~ dnorm(m.delta, pr.delta)
-                gamma[j,k] <- gamma[j,2]
-        }
-    }
-    beta ~ dnorm(m.beta,pr.beta)
-    pr.delta <- pow(s.delta, -2)
-    pr.beta <- pow(s.beta,-2)
-    pr.gamma <- pow(s.gamma,-2)
-}