Model Bayesian latent profile plus proportional odds model for Beat the Blues Data (M=3) ++++++++++++++++++++++++ model { for(i in 1:N){ T[i] ~ dcat( P[i, 1:3] ) # for(j in 1:nr[i]){ cfb[i, j] <- bdi[ne[i]+j-1] - base[ne[i]+j-1] cfb[i, j] ~ dnorm(mu[T[i], visit[ne[i]+j-1]], tau) } # P[i,1] <- 1/(1+exp( - (theta11 + gam1*treatment[ne[i]] + gam2*(base[ne[i]] - base.bar) + gam3*drug[ne[i]] + gam4*length[ne[i]] ))) P[i,2] <- 1/(1+exp( - (theta22 + gam1*treatment[ne[i]] + gam2*(base[ne[i]] - base.bar) + gam3*drug[ne[i]] + gam4*length[ne[i]] ))) - P[i,1] P[i,3] <- 1 - P[i,1] - P[i,2] } # for (k in 1:4){ mu[1, time[k]]< -beta11* time[k] + beta12 * time[k]*time[k] mu[2, time[k]]<- (beta11+alp21)* time[k] + (beta12+alp22) * time[k]*time[k] mu[3, time[k]]<- 0 } base.bar <- mean(base[]) # priors beta11 ~ dnorm(0.0,1.0E-4) beta12 ~ dnorm(0.0,1.0E-4) theta11 ~ dnorm(0.0,1.0E-4) theta22 ~ dnorm(0.0,1.0E-4) alp21 ~ dnorm(0.0,1.0E-4) alp22 ~ dnorm(0.0,1.0E-4) gam1 ~ dnorm(0.0,1.0E-4) gam2 ~ dnorm(0.0,1.0E-4) gam3 ~ dnorm(0.0,1.0E-4) gam4 ~ dnorm(0.0,1.0E-4) tau ~ dgamma(1.0E-3,1.0E-3) sigma <- 1.0/ sqrt(tau) # beta21<-beta11+alp21 beta22<-beta12+alp22 theta1<-theta11 - gam2*base.bar theta2<-theta22 - gam2*base.bar } Initial values list(tau =0.5, beta11=-8, beta12=0, alp21=5, alp22=0.5, theta11=1, theta22=2, gam1=1, gam2=1, gam3=0, gam4=0) list(tau =0.2, beta11=-6, beta12=1, alp21=3, alp22=0.8, theta11=1, theta22=2, gam1=0.5, gam2=0.5, gam3=1, gam4=0.5) Data list(N=97, time=c(2,3,5,8), nr=c(2, 1, 4, 4, 4, 2, 4, 4, 4, 4, 2, 1, 1, 1, 3, 4, 1, 4, 4, 4, 4, 4, 4, 1, 4, 2, 2, 2, 4, 4, 3, 3, 3, 4, 1, 2, 1, 4, 1, 4, 4, 4, 4, 4, 4, 4, 4, 1, 4, 4, 4, 1, 4, 2, 4, 4, 1, 2, 4, 4, 4, 4, 1, 4, 2, 1, 1, 4, 1, 4, 3, 1, 1, 2, 4, 1, 1, 1, 2, 1, 4, 2, 4, 4, 4, 1, 2, 4, 1, 4, 4, 2, 3, 4, 4, 4, 4)) no[ ] subject[ ] drug[ ] length[ ] treatment[ ] visit[ ] bdi[ ] base[ ] id[ ] 1 1 0 1 0 2 2 29 1 2 1 0 1 0 3 2 29 1 3 3 1 0 0 2 20 25 2 4 7 1 0 0 2 7 17 3 5 7 1 0 0 3 7 17 3 6 7 1 0 0 5 3 17 3 7 7 1 0 0 8 7 17 3 8 8 0 1 0 2 20 20 4 9 8 0 1 0 3 21 20 4 10 8 0 1 0 5 19 20 4 11 8 0 1 0 8 13 20 4 12 11 0 1 0 2 32 30 5 13 11 0 1 0 3 24 30 5 14 11 0 1 0 5 12 30 5 15 11 0 1 0 8 2 30 5 16 13 0 1 0 2 27 26 6 17 13 0 1 0 3 23 26 6 18 14 1 1 0 2 26 30 7 19 14 1 1 0 3 36 30 7 20 14 1 1 0 5 27 30 7 21 14 1 1 0 8 22 30 7 22 16 0 0 0 2 13 16 8 23 16 0 0 0 3 3 16 8 24 16 0 0 0 5 2 16 8 25 16 0 0 0 8 0 16 8 26 19 0 1 0 2 30 37 9 27 19 0 1 0 3 33 37 9 28 19 0 1 0 5 31 37 9 29 19 0 1 0 8 22 37 9 30 22 0 0 0 2 17 26 10 31 22 0 0 0 3 17 26 10 32 22 0 0 0 5 20 26 10 33 22 0 0 0 8 12 26 10 34 23 0 1 0 2 22 29 11 35 23 0 1 0 3 10 29 11 36 24 0 1 0 2 21 20 12 37 25 0 1 0 2 23 33 13 38 27 1 0 0 2 15 12 14 39 28 1 1 0 2 36 47 15 40 28 1 1 0 3 49 47 15 41 28 1 1 0 5 34 47 15 42 31 0 0 0 2 7 27 16 43 31 0 0 0 3 15 27 16 44 31 0 0 0 5 16 27 16 45 31 0 0 0 8 0 27 16 46 36 0 0 0 2 38 38 17 47 37 0 0 0 2 14 21 18 48 37 0 0 0 3 20 21 18 49 37 0 0 0 5 1 21 18 50 37 0 0 0 8 8 21 18 51 38 1 1 0 2 17 34 19 52 38 1 1 0 3 8 34 19 53 38 1 1 0 5 9 34 19 54 38 1 1 0 8 13 34 19 55 40 1 1 0 2 27 38 20 56 40 1 1 0 3 19 38 20 57 40 1 1 0 5 20 38 20 58 40 1 1 0 8 30 38 20 59 42 0 0 0 2 19 20 21 60 42 0 0 0 3 18 20 21 61 42 0 0 0 5 19 20 21 62 42 0 0 0 8 18 20 21 63 43 1 1 0 2 29 17 22 64 43 1 1 0 3 2 17 22 65 43 1 1 0 5 0 17 22 66 43 1 1 0 8 0 17 22 67 50 0 1 0 2 42 36 23 68 50 0 1 0 3 49 36 23 69 50 0 1 0 5 47 36 23 70 50 0 1 0 8 40 36 23 71 51 0 0 0 2 30 35 24 72 53 0 1 0 2 48 31 25 73 53 0 1 0 3 38 31 25 74 53 0 1 0 5 38 31 25 75 53 0 1 0 8 37 31 25 76 55 1 0 0 2 21 23 26 77 55 1 0 0 3 26 23 26 78 57 0 0 0 2 13 14 27 79 57 0 0 0 3 14 14 27 80 58 0 0 0 2 36 40 28 81 58 0 0 0 3 33 40 28 82 61 0 1 0 2 20 22 29 83 61 0 1 0 3 16 22 29 84 61 0 1 0 5 24 22 29 85 61 0 1 0 8 16 22 29 86 62 0 1 0 2 23 23 30 87 62 0 1 0 3 15 23 30 88 62 0 1 0 5 25 23 30 89 62 0 1 0 8 17 23 30 90 63 0 0 0 2 7 15 31 91 63 0 0 0 3 13 15 31 92 63 0 0 0 5 13 15 31 93 64 0 1 0 2 12 8 32 94 64 0 1 0 3 11 8 32 95 64 0 1 0 5 26 8 32 96 66 0 1 0 2 6 7 33 97 66 0 1 0 3 2 7 33 98 66 0 1 0 5 1 7 33 99 67 1 0 0 2 9 17 34 100 67 1 0 0 3 3 17 34 101 67 1 0 0 5 1 17 34 102 67 1 0 0 8 0 17 34 103 69 0 0 0 2 20 27 35 104 73 0 1 0 2 8 11 36 105 73 0 1 0 3 7 11 36 106 74 0 0 0 2 8 9 37 107 75 0 1 0 2 22 14 38 108 75 0 1 0 3 21 14 38 109 75 0 1 0 5 24 14 38 110 75 0 1 0 8 19 14 38 111 79 0 0 0 2 9 23 39 112 80 0 1 0 2 22 21 40 113 80 0 1 0 3 24 21 40 114 80 0 1 0 5 23 21 40 115 80 0 1 0 8 22 21 40 116 81 0 1 0 2 31 27 41 117 81 0 1 0 3 28 27 41 118 81 0 1 0 5 22 27 41 119 81 0 1 0 8 14 27 41 120 83 0 1 0 2 13 10 42 121 83 0 1 0 3 12 10 42 122 83 0 1 0 5 8 10 42 123 83 0 1 0 8 20 10 42 124 84 1 0 0 2 9 21 43 125 84 1 0 0 3 6 21 43 126 84 1 0 0 5 7 21 43 127 84 1 0 0 8 1 21 43 128 88 1 1 0 2 0 35 44 129 88 1 1 0 3 6 35 44 130 88 1 1 0 5 0 35 44 131 88 1 1 0 8 1 35 44 132 99 0 0 0 2 5 13 45 133 99 0 0 0 3 5 13 45 134 99 0 0 0 5 0 13 45 135 99 0 0 0 8 6 13 45 136 2 1 1 1 2 16 32 46 137 2 1 1 1 3 24 32 46 138 2 1 1 1 5 17 32 46 139 2 1 1 1 8 20 32 46 140 4 0 1 1 2 17 21 47 141 4 0 1 1 3 16 21 47 142 4 0 1 1 5 10 21 47 143 4 0 1 1 8 9 21 47 144 5 1 1 1 2 23 26 48 145 6 1 0 1 2 0 7 49 146 6 1 0 1 3 0 7 49 147 6 1 0 1 5 0 7 49 148 6 1 0 1 8 0 7 49 149 9 1 0 1 2 13 18 50 150 9 1 0 1 3 14 18 50 151 9 1 0 1 5 20 18 50 152 9 1 0 1 8 11 18 50 153 10 1 1 1 2 5 20 51 154 10 1 1 1 3 5 20 51 155 10 1 1 1 5 8 20 51 156 10 1 1 1 8 12 20 51 157 12 1 0 1 2 35 49 52 158 15 1 1 1 2 13 23 53 159 15 1 1 1 3 13 23 53 160 15 1 1 1 5 12 23 53 161 15 1 1 1 8 23 23 53 162 17 0 1 1 2 30 30 54 163 17 0 1 1 3 29 30 54 164 18 0 0 1 2 8 13 55 165 18 0 0 1 3 8 13 55 166 18 0 0 1 5 7 13 55 167 18 0 0 1 8 6 13 55 168 20 1 0 1 2 12 35 56 169 20 1 0 1 3 10 35 56 170 20 1 0 1 5 8 35 56 171 20 1 0 1 8 10 35 56 172 21 0 1 1 2 6 21 57 173 26 0 1 1 2 12 19 58 174 26 0 1 1 3 13 19 58 175 29 1 1 1 2 6 36 59 176 29 1 1 1 3 0 36 59 177 29 1 1 1 5 0 36 59 178 29 1 1 1 8 2 36 59 179 30 0 0 1 2 8 10 60 180 30 0 0 1 3 6 10 60 181 30 0 0 1 5 3 10 60 182 30 0 0 1 8 3 10 60 183 32 0 0 1 2 10 18 61 184 32 0 0 1 3 10 18 61 185 32 0 0 1 5 6 18 61 186 32 0 0 1 8 8 18 61 187 33 1 0 1 2 8 11 62 188 33 1 0 1 3 3 11 62 189 33 1 0 1 5 2 11 62 190 33 1 0 1 8 15 11 62 191 34 1 0 1 2 7 6 63 192 35 1 1 1 2 24 44 64 193 35 1 1 1 3 20 44 64 194 35 1 1 1 5 29 44 64 195 35 1 1 1 8 14 44 64 196 39 1 0 1 2 7 9 65 197 39 1 0 1 3 1 9 65 198 41 1 0 1 2 40 46 66 199 44 0 1 1 2 20 18 67 200 45 1 1 1 2 1 42 68 201 45 1 1 1 3 8 42 68 202 45 1 1 1 5 10 42 68 203 45 1 1 1 8 6 42 68 204 46 0 0 1 2 30 30 69 205 47 1 0 1 2 27 33 70 206 47 1 0 1 3 16 33 70 207 47 1 0 1 5 30 33 70 208 47 1 0 1 8 15 33 70 209 48 0 0 1 2 1 12 71 210 48 0 0 1 3 0 12 71 211 48 0 0 1 5 0 12 71 212 49 1 0 1 2 5 2 72 213 52 0 0 1 2 20 23 73 214 54 1 0 1 2 5 8 74 215 54 1 0 1 3 7 8 74 216 56 1 0 1 2 7 7 75 217 56 1 0 1 3 5 7 75 218 56 1 0 1 5 4 7 75 219 56 1 0 1 8 0 7 75 220 59 1 0 1 2 30 23 76 221 60 0 1 1 2 3 14 77 222 65 0 1 1 2 18 12 78 223 68 1 0 1 2 18 33 79 224 68 1 0 1 3 16 33 79 225 70 0 0 1 2 30 27 80 226 71 0 0 1 2 6 9 81 227 71 0 0 1 3 10 9 81 228 71 0 0 1 5 1 9 81 229 71 0 0 1 8 0 9 81 230 72 0 1 1 2 30 40 82 231 72 0 1 1 3 12 40 82 232 76 1 1 1 2 9 28 83 233 76 1 1 1 3 20 28 83 234 76 1 1 1 5 18 28 83 235 76 1 1 1 8 13 28 83 236 77 0 1 1 2 9 15 84 237 77 0 1 1 3 13 15 84 238 77 0 1 1 5 14 15 84 239 77 0 1 1 8 10 15 84 240 78 1 1 1 2 10 22 85 241 78 1 1 1 3 5 22 85 242 78 1 1 1 5 5 22 85 243 78 1 1 1 8 12 22 85 244 82 1 1 1 2 15 14 86 245 85 1 1 1 2 36 46 87 246 85 1 1 1 3 53 46 87 247 86 0 1 1 2 14 36 88 248 86 0 1 1 3 7 36 88 249 86 0 1 1 5 15 36 88 250 86 0 1 1 8 15 36 88 251 87 1 1 1 2 17 23 89 252 89 1 0 1 2 13 33 90 253 89 1 0 1 3 13 33 90 254 89 1 0 1 5 10 33 90 255 89 1 0 1 8 8 33 90 256 90 0 0 1 2 4 19 91 257 90 0 0 1 3 27 19 91 258 90 0 0 1 5 1 19 91 259 90 0 0 1 8 2 19 91 260 92 1 0 1 2 26 30 92 261 92 1 0 1 3 28 30 92 262 93 1 0 1 2 8 17 93 263 93 1 0 1 3 7 17 93 264 93 1 0 1 5 12 17 93 265 94 0 1 1 2 4 19 94 266 94 0 1 1 3 3 19 94 267 94 0 1 1 5 3 19 94 268 94 0 1 1 8 3 19 94 269 95 0 1 1 2 11 16 95 270 95 0 1 1 3 4 16 95 271 95 0 1 1 5 2 16 95 272 95 0 1 1 8 3 16 95 273 96 1 1 1 2 16 16 96 274 96 1 1 1 3 10 16 96 275 96 1 1 1 5 10 16 96 276 96 1 1 1 8 8 16 96 277 98 0 1 1 2 22 11 97 278 98 0 1 1 3 9 11 97 279 98 0 1 1 5 11 11 97 280 98 0 1 1 8 11 11 97 END ne[ ] iid[ ] 1 1 3 2 4 3 8 4 12 5 16 6 18 7 22 8 26 9 30 10 34 11 36 12 37 13 38 14 39 15 42 16 46 17 47 18 51 19 55 20 59 21 63 22 67 23 71 24 72 25 76 26 78 27 80 28 82 29 86 30 90 31 93 32 96 33 99 34 103 35 104 36 106 37 107 38 111 39 112 40 116 41 120 42 124 43 128 44 132 45 136 46 140 47 144 48 145 49 149 50 153 51 157 52 158 53 162 54 164 55 168 56 172 57 173 58 175 59 179 60 183 61 187 62 191 63 192 64 196 65 198 66 199 67 200 68 204 69 205 70 209 71 212 72 213 73 214 74 216 75 220 76 221 77 222 78 223 79 225 80 226 81 230 82 232 83 236 84 240 85 244 86 245 87 247 88 251 89 252 90 256 91 260 92 262 93 265 94 269 95 273 96 277 97 END T[] 1 NA 2 3 2 NA NA 2 2 2 2 NA NA NA 2 2 NA 2 1 2 3 2 3 NA 3 NA NA NA NA NA NA 3 NA 2 NA NA NA 3 NA 3 NA 3 2 1 2 2 2 NA 2 NA 2 NA 2 NA NA 1 NA 2 1 NA 2 NA NA 1 NA NA NA 1 NA 2 2 NA NA NA NA NA NA NA NA NA NA 1 2 NA 2 NA NA 1 NA 1 2 NA 2 2 2 NA 3 END