glm_lpdf

Author: mcdonohue

src/stan_files/model/ltjmm_lp.stan

@@ -327,43 +327,43 @@
 
   // Log-likelihoods
   if (prior_PD == 0) {
-    glm_lp(yReal1, yInt1, yEta1, yAux1, family[1], link[1], sum_log_y1, sqrt_y1, log_y1, weights_y1);
+    target += glm_lpdf(yReal1 | yInt1, yEta1, yAux1, family[1], link[1], sum_log_y1, sqrt_y1, log_y1);
     if (M > 1)
-      glm_lp(yReal2, yInt2, yEta2, yAux2, family[2], link[2], sum_log_y2, sqrt_y2, log_y2, weights_y2);
+      target += glm_lpdf(yReal2 | yInt2, yEta2, yAux2, family[2], link[2], sum_log_y2, sqrt_y2, log_y2);
     if (M > 2)
-      glm_lp(yReal3, yInt3, yEta3, yAux3, family[3], link[3], sum_log_y3, sqrt_y3, log_y3, weights_y3);
+      target += glm_lpdf(yReal3 | yInt3, yEta3, yAux3, family[3], link[3], sum_log_y3, sqrt_y3, log_y3);
     if (M > 3)
-      glm_lp(yReal4, yInt4, yEta4, yAux4, family[4], link[4], sum_log_y4, sqrt_y4, log_y4, weights_y4);
+      target += glm_lpdf(yReal4 | yInt4, yEta4, yAux4, family[4], link[4], sum_log_y4, sqrt_y4, log_y4);
     if (M > 4)
-      glm_lp(yReal5, yInt5, yEta5, yAux5, family[5], link[5], sum_log_y5, sqrt_y5, log_y5, weights_y5);
+      target += glm_lpdf(yReal5 | yInt5, yEta5, yAux5, family[5], link[5], sum_log_y5, sqrt_y5, log_y5);
     if (M > 5)
-      glm_lp(yReal6, yInt6, yEta6, yAux6, family[6], link[6], sum_log_y6, sqrt_y6, log_y6, weights_y6);
+      target += glm_lpdf(yReal6 | yInt6, yEta6, yAux6, family[6], link[6], sum_log_y6, sqrt_y6, log_y6);
     if (M > 6)
-      glm_lp(yReal7, yInt7, yEta7, yAux7, family[7], link[7], sum_log_y7, sqrt_y7, log_y7, weights_y7);
+      target += glm_lpdf(yReal7 | yInt7, yEta7, yAux7, family[7], link[7], sum_log_y7, sqrt_y7, log_y7);
     if (M > 7)
-      glm_lp(yReal8, yInt8, yEta8, yAux8, family[8], link[8], sum_log_y8, sqrt_y8, log_y8, weights_y8);
+      target += glm_lpdf(yReal8 | yInt8, yEta8, yAux8, family[8], link[8], sum_log_y8, sqrt_y8, log_y8);
     if (M > 8)
-      glm_lp(yReal9, yInt9, yEta9, yAux9, family[9], link[9], sum_log_y9, sqrt_y9, log_y9, weights_y9);
+      target += glm_lpdf(yReal9 | yInt9, yEta9, yAux9, family[9], link[9], sum_log_y9, sqrt_y9, log_y9);
     if (M > 9)
-      glm_lp(yReal10, yInt10, yEta10, yAux10, family[10], link[10], sum_log_y10, sqrt_y10, log_y10, weights_y10);
+      target += glm_lpdf(yReal10 | yInt10, yEta10, yAux10, family[10], link[10], sum_log_y10, sqrt_y10, log_y10);
     if (M > 10)
-      glm_lp(yReal11, yInt11, yEta11, yAux11, family[11], link[11], sum_log_y11, sqrt_y11, log_y11, weights_y11);
+      target += glm_lpdf(yReal11 | yInt11, yEta11, yAux11, family[11], link[11], sum_log_y11, sqrt_y11, log_y11);
     if (M > 11)
-      glm_lp(yReal12, yInt12, yEta12, yAux12, family[12], link[12], sum_log_y12, sqrt_y12, log_y12, weights_y12);
+      target += glm_lpdf(yReal12 | yInt12, yEta12, yAux12, family[12], link[12], sum_log_y12, sqrt_y12, log_y12);
     if (M > 12)
-      glm_lp(yReal13, yInt13, yEta13, yAux13, family[13], link[13], sum_log_y13, sqrt_y13, log_y13, weights_y13);
+      target += glm_lpdf(yReal13 | yInt13, yEta13, yAux13, family[13], link[13], sum_log_y13, sqrt_y13, log_y13);
     if (M > 13)
-      glm_lp(yReal14, yInt14, yEta14, yAux14, family[14], link[14], sum_log_y14, sqrt_y14, log_y14, weights_y14);
+      target += glm_lpdf(yReal14 | yInt14, yEta14, yAux14, family[14], link[14], sum_log_y14, sqrt_y14, log_y14);
     if (M > 14)
-      glm_lp(yReal15, yInt15, yEta15, yAux15, family[15], link[15], sum_log_y15, sqrt_y15, log_y15, weights_y15);
+      target += glm_lpdf(yReal15 | yInt15, yEta15, yAux15, family[15], link[15], sum_log_y15, sqrt_y15, log_y15);
     if (M > 15)
-      glm_lp(yReal16, yInt16, yEta16, yAux16, family[16], link[16], sum_log_y16, sqrt_y16, log_y16, weights_y16);
+      target += glm_lpdf(yReal16 | yInt16, yEta16, yAux16, family[16], link[16], sum_log_y16, sqrt_y16, log_y16);
     if (M > 16)
-      glm_lp(yReal17, yInt17, yEta17, yAux17, family[17], link[17], sum_log_y17, sqrt_y17, log_y17, weights_y17);
+      target += glm_lpdf(yReal17 | yInt17, yEta17, yAux17, family[17], link[17], sum_log_y17, sqrt_y17, log_y17);
     if (M > 17)
-      glm_lp(yReal18, yInt18, yEta18, yAux18, family[18], link[18], sum_log_y18, sqrt_y18, log_y18, weights_y18);
+      target += glm_lpdf(yReal18 | yInt18, yEta18, yAux18, family[18], link[18], sum_log_y18, sqrt_y18, log_y18);
     if (M > 18)
-      glm_lp(yReal19, yInt19, yEta19, yAux19, family[19], link[19], sum_log_y19, sqrt_y19, log_y19, weights_y19);
+      target += glm_lpdf(yReal19 | yInt19, yEta19, yAux19, family[19], link[19], sum_log_y19, sqrt_y19, log_y19);
     if (M > 19)
-      glm_lp(yReal20, yInt20, yEta20, yAux20, family[20], link[20], sum_log_y20, sqrt_y20, log_y20, weights_y20);  
+      target += glm_lpdf(yReal20 | yInt20, yEta20, yAux20, family[20], link[20], sum_log_y20, sqrt_y20, log_y20);
   }