Table 2

Simulations with K=3 longitudinal markers with a mixture of distributions (continuous, count, and binary)

Approach:R-INLA 1
R-INLA 2
JMbayes2
rstanarm 1
rstanarm 2
(Empirical Bayes)
(Full Bayesian)
(1 chain / 1000 iter.)
(4 chains / 2000 iter.)
True valueBias(SD)CP (%)Bias(SD)CPBias(SD)CPBias(SD)CPBias(SD)CP
β10=0.20.001(0.049)950.001(0.049)960.002(0.049)950.006(0.051)93–0.001(0.048)96
β11=0.10.003(0.019)960.003(0.019)950.002(0.019)950.003(0.021)93–0.001(0.019)94
β12=0.10(0.04)950(0.04)94–0.001(0.04)94–0.004(0.041)930.004(0.039)97
β13=0.2–0.001(0.034)970(0.034)980(0.034)98–0.003(0.035)96–0.007(0.033)98
σε1=0.40(0.003)940(0.003)940(0.003)940.002(0.01)920.001(0.004)95
β20=3–0.002(0.053)96–0.002(0.053)950.002(0.053)960.001(0.049)98–0.008(0.047)98
β21=0.1–0.002(0.02)96–0.002(0.02)96–0.001(0.02)96–0.002(0.019)97–0.003(0.019)98
β22=0.10.002(0.043)960.002(0.043)950(0.043)960.001(0.041)950.005(0.041)95
β23=0.20.004(0.045)940.004(0.045)940.003(0.045)950.002(0.045)940.009(0.04)96
β30=1–0.014(0.092)95–0.014(0.091)95–0.005(0.093)940.001(0.096)940.026(0.102)97
β31=10.008(0.031)930.009(0.031)93–0.001(0.031)95–0.001(0.032)94–0.007(0.031)96
β32=1–0.006(0.074)96–0.007(0.073)950.004(0.075)97–0.002(0.076)94–0.008(0.071)98
β33=10.01(0.07)950.011(0.071)95–0.001(0.071)95–0.003(0.068)96–0.018(0.077)85
σb102=0.160.003(0.016)950.002(0.014)94–0.001(0.014)99–0.007(0.029)94–0.008(0.021)96
σb112=0.090.003(0.01)940.004(0.009)950(0.009)930.008(0.038)880.011(0.025)94
σb202=0.250.002(0.019)950.001(0.017)960.002(0.017)950.004(0.017)960.002(0.016)95
σb212=0.160.002(0.014)950.002(0.013)950.003(0.014)940.002(0.012)970.007(0.016)82
σb302=0.25–0.01(0.036)96–0.009(0.037)960(0.041)930.004(0.041)960.002(0.037)97
covb10,b11=0.03–0.002(0.009)96–0.004(0.009)94–0.001(0.007)97–0.001(0.01)96–0.003(0.009)96
covb10,b20=0.020(0.011)950(0.011)950(0.01)96–0.001(0.01)94–0.003(0.011)97
covb10,b21=0.04–0.001(0.01)96–0.001(0.01)960(0.009)96–0.002(0.011)940(0.009)97
covb10,b30=00(0.019)96–0.003(0.017)970.001(0.016)95–0.002(0.016)94–0.001(0.017)95
covb11,b20=0.03–0.001(0.01)96–0.001(0.01)96–0.001(0.009)960(0.011)94–0.002(0.01)97
covb11,b21=00(0.008)950(0.008)960(0.008)950.002(0.01)920(0.008)96
covb11,b30=0.060.005(0.014)920.006(0.013)940.002(0.014)930(0.015)95–0.001(0.014)97
covb20,b21=0.080(0.012)950(0.012)95–0.001(0.011)940(0.011)960.002(0.011)96
covb20,b30=0.05–0.001(0.019)95–0.002(0.019)95–0.003(0.018)94–0.001(0.019)960.001(0.018)96
covb21,b30=0.04–0.001(0.017)95–0.001(0.017)96–0.003(0.016)940.009(0.019)920.011(0.018)95
φ1=0.5–0.003(0.101)95–0.003(0.102)940(0.103)980.009(0.113)940.021(0.114)96
φ2=0.20.006(0.074)930.005(0.073)940.005(0.074)970.002(0.077)930.006(0.072)94
φ3=0.3–0.001(0.094)940.001(0.092)94–0.004(0.092)94–0.006(0.093)94–0.003(0.087)95
Conv. rate1110.510.64
Comp. time (s)43.55 (5.58)69.8 (44.59)295.73 (20.22)3766.41 (1429.06)10,861.22 (3735.02)
Approach:R-INLA 1
R-INLA 2
JMbayes2
rstanarm 1
rstanarm 2
(Empirical Bayes)
(Full Bayesian)
(1 chain / 1000 iter.)
(4 chains / 2000 iter.)
True valueBias(SD)CP (%)Bias(SD)CPBias(SD)CPBias(SD)CPBias(SD)CP
β10=0.20.001(0.049)950.001(0.049)960.002(0.049)950.006(0.051)93–0.001(0.048)96
β11=0.10.003(0.019)960.003(0.019)950.002(0.019)950.003(0.021)93–0.001(0.019)94
β12=0.10(0.04)950(0.04)94–0.001(0.04)94–0.004(0.041)930.004(0.039)97
β13=0.2–0.001(0.034)970(0.034)980(0.034)98–0.003(0.035)96–0.007(0.033)98
σε1=0.40(0.003)940(0.003)940(0.003)940.002(0.01)920.001(0.004)95
β20=3–0.002(0.053)96–0.002(0.053)950.002(0.053)960.001(0.049)98–0.008(0.047)98
β21=0.1–0.002(0.02)96–0.002(0.02)96–0.001(0.02)96–0.002(0.019)97–0.003(0.019)98
β22=0.10.002(0.043)960.002(0.043)950(0.043)960.001(0.041)950.005(0.041)95
β23=0.20.004(0.045)940.004(0.045)940.003(0.045)950.002(0.045)940.009(0.04)96
β30=1–0.014(0.092)95–0.014(0.091)95–0.005(0.093)940.001(0.096)940.026(0.102)97
β31=10.008(0.031)930.009(0.031)93–0.001(0.031)95–0.001(0.032)94–0.007(0.031)96
β32=1–0.006(0.074)96–0.007(0.073)950.004(0.075)97–0.002(0.076)94–0.008(0.071)98
β33=10.01(0.07)950.011(0.071)95–0.001(0.071)95–0.003(0.068)96–0.018(0.077)85
σb102=0.160.003(0.016)950.002(0.014)94–0.001(0.014)99–0.007(0.029)94–0.008(0.021)96
σb112=0.090.003(0.01)940.004(0.009)950(0.009)930.008(0.038)880.011(0.025)94
σb202=0.250.002(0.019)950.001(0.017)960.002(0.017)950.004(0.017)960.002(0.016)95
σb212=0.160.002(0.014)950.002(0.013)950.003(0.014)940.002(0.012)970.007(0.016)82
σb302=0.25–0.01(0.036)96–0.009(0.037)960(0.041)930.004(0.041)960.002(0.037)97
covb10,b11=0.03–0.002(0.009)96–0.004(0.009)94–0.001(0.007)97–0.001(0.01)96–0.003(0.009)96
covb10,b20=0.020(0.011)950(0.011)950(0.01)96–0.001(0.01)94–0.003(0.011)97
covb10,b21=0.04–0.001(0.01)96–0.001(0.01)960(0.009)96–0.002(0.011)940(0.009)97
covb10,b30=00(0.019)96–0.003(0.017)970.001(0.016)95–0.002(0.016)94–0.001(0.017)95
covb11,b20=0.03–0.001(0.01)96–0.001(0.01)96–0.001(0.009)960(0.011)94–0.002(0.01)97
covb11,b21=00(0.008)950(0.008)960(0.008)950.002(0.01)920(0.008)96
covb11,b30=0.060.005(0.014)920.006(0.013)940.002(0.014)930(0.015)95–0.001(0.014)97
covb20,b21=0.080(0.012)950(0.012)95–0.001(0.011)940(0.011)960.002(0.011)96
covb20,b30=0.05–0.001(0.019)95–0.002(0.019)95–0.003(0.018)94–0.001(0.019)960.001(0.018)96
covb21,b30=0.04–0.001(0.017)95–0.001(0.017)96–0.003(0.016)940.009(0.019)920.011(0.018)95
φ1=0.5–0.003(0.101)95–0.003(0.102)940(0.103)980.009(0.113)940.021(0.114)96
φ2=0.20.006(0.074)930.005(0.073)940.005(0.074)970.002(0.077)930.006(0.072)94
φ3=0.3–0.001(0.094)940.001(0.092)94–0.004(0.092)94–0.006(0.093)94–0.003(0.087)95
Conv. rate1110.510.64
Comp. time (s)43.55 (5.58)69.8 (44.59)295.73 (20.22)3766.41 (1429.06)10,861.22 (3735.02)
Table 2

Simulations with K=3 longitudinal markers with a mixture of distributions (continuous, count, and binary)

Approach:R-INLA 1
R-INLA 2
JMbayes2
rstanarm 1
rstanarm 2
(Empirical Bayes)
(Full Bayesian)
(1 chain / 1000 iter.)
(4 chains / 2000 iter.)
True valueBias(SD)CP (%)Bias(SD)CPBias(SD)CPBias(SD)CPBias(SD)CP
β10=0.20.001(0.049)950.001(0.049)960.002(0.049)950.006(0.051)93–0.001(0.048)96
β11=0.10.003(0.019)960.003(0.019)950.002(0.019)950.003(0.021)93–0.001(0.019)94
β12=0.10(0.04)950(0.04)94–0.001(0.04)94–0.004(0.041)930.004(0.039)97
β13=0.2–0.001(0.034)970(0.034)980(0.034)98–0.003(0.035)96–0.007(0.033)98
σε1=0.40(0.003)940(0.003)940(0.003)940.002(0.01)920.001(0.004)95
β20=3–0.002(0.053)96–0.002(0.053)950.002(0.053)960.001(0.049)98–0.008(0.047)98
β21=0.1–0.002(0.02)96–0.002(0.02)96–0.001(0.02)96–0.002(0.019)97–0.003(0.019)98
β22=0.10.002(0.043)960.002(0.043)950(0.043)960.001(0.041)950.005(0.041)95
β23=0.20.004(0.045)940.004(0.045)940.003(0.045)950.002(0.045)940.009(0.04)96
β30=1–0.014(0.092)95–0.014(0.091)95–0.005(0.093)940.001(0.096)940.026(0.102)97
β31=10.008(0.031)930.009(0.031)93–0.001(0.031)95–0.001(0.032)94–0.007(0.031)96
β32=1–0.006(0.074)96–0.007(0.073)950.004(0.075)97–0.002(0.076)94–0.008(0.071)98
β33=10.01(0.07)950.011(0.071)95–0.001(0.071)95–0.003(0.068)96–0.018(0.077)85
σb102=0.160.003(0.016)950.002(0.014)94–0.001(0.014)99–0.007(0.029)94–0.008(0.021)96
σb112=0.090.003(0.01)940.004(0.009)950(0.009)930.008(0.038)880.011(0.025)94
σb202=0.250.002(0.019)950.001(0.017)960.002(0.017)950.004(0.017)960.002(0.016)95
σb212=0.160.002(0.014)950.002(0.013)950.003(0.014)940.002(0.012)970.007(0.016)82
σb302=0.25–0.01(0.036)96–0.009(0.037)960(0.041)930.004(0.041)960.002(0.037)97
covb10,b11=0.03–0.002(0.009)96–0.004(0.009)94–0.001(0.007)97–0.001(0.01)96–0.003(0.009)96
covb10,b20=0.020(0.011)950(0.011)950(0.01)96–0.001(0.01)94–0.003(0.011)97
covb10,b21=0.04–0.001(0.01)96–0.001(0.01)960(0.009)96–0.002(0.011)940(0.009)97
covb10,b30=00(0.019)96–0.003(0.017)970.001(0.016)95–0.002(0.016)94–0.001(0.017)95
covb11,b20=0.03–0.001(0.01)96–0.001(0.01)96–0.001(0.009)960(0.011)94–0.002(0.01)97
covb11,b21=00(0.008)950(0.008)960(0.008)950.002(0.01)920(0.008)96
covb11,b30=0.060.005(0.014)920.006(0.013)940.002(0.014)930(0.015)95–0.001(0.014)97
covb20,b21=0.080(0.012)950(0.012)95–0.001(0.011)940(0.011)960.002(0.011)96
covb20,b30=0.05–0.001(0.019)95–0.002(0.019)95–0.003(0.018)94–0.001(0.019)960.001(0.018)96
covb21,b30=0.04–0.001(0.017)95–0.001(0.017)96–0.003(0.016)940.009(0.019)920.011(0.018)95
φ1=0.5–0.003(0.101)95–0.003(0.102)940(0.103)980.009(0.113)940.021(0.114)96
φ2=0.20.006(0.074)930.005(0.073)940.005(0.074)970.002(0.077)930.006(0.072)94
φ3=0.3–0.001(0.094)940.001(0.092)94–0.004(0.092)94–0.006(0.093)94–0.003(0.087)95
Conv. rate1110.510.64
Comp. time (s)43.55 (5.58)69.8 (44.59)295.73 (20.22)3766.41 (1429.06)10,861.22 (3735.02)
Approach:R-INLA 1
R-INLA 2
JMbayes2
rstanarm 1
rstanarm 2
(Empirical Bayes)
(Full Bayesian)
(1 chain / 1000 iter.)
(4 chains / 2000 iter.)
True valueBias(SD)CP (%)Bias(SD)CPBias(SD)CPBias(SD)CPBias(SD)CP
β10=0.20.001(0.049)950.001(0.049)960.002(0.049)950.006(0.051)93–0.001(0.048)96
β11=0.10.003(0.019)960.003(0.019)950.002(0.019)950.003(0.021)93–0.001(0.019)94
β12=0.10(0.04)950(0.04)94–0.001(0.04)94–0.004(0.041)930.004(0.039)97
β13=0.2–0.001(0.034)970(0.034)980(0.034)98–0.003(0.035)96–0.007(0.033)98
σε1=0.40(0.003)940(0.003)940(0.003)940.002(0.01)920.001(0.004)95
β20=3–0.002(0.053)96–0.002(0.053)950.002(0.053)960.001(0.049)98–0.008(0.047)98
β21=0.1–0.002(0.02)96–0.002(0.02)96–0.001(0.02)96–0.002(0.019)97–0.003(0.019)98
β22=0.10.002(0.043)960.002(0.043)950(0.043)960.001(0.041)950.005(0.041)95
β23=0.20.004(0.045)940.004(0.045)940.003(0.045)950.002(0.045)940.009(0.04)96
β30=1–0.014(0.092)95–0.014(0.091)95–0.005(0.093)940.001(0.096)940.026(0.102)97
β31=10.008(0.031)930.009(0.031)93–0.001(0.031)95–0.001(0.032)94–0.007(0.031)96
β32=1–0.006(0.074)96–0.007(0.073)950.004(0.075)97–0.002(0.076)94–0.008(0.071)98
β33=10.01(0.07)950.011(0.071)95–0.001(0.071)95–0.003(0.068)96–0.018(0.077)85
σb102=0.160.003(0.016)950.002(0.014)94–0.001(0.014)99–0.007(0.029)94–0.008(0.021)96
σb112=0.090.003(0.01)940.004(0.009)950(0.009)930.008(0.038)880.011(0.025)94
σb202=0.250.002(0.019)950.001(0.017)960.002(0.017)950.004(0.017)960.002(0.016)95
σb212=0.160.002(0.014)950.002(0.013)950.003(0.014)940.002(0.012)970.007(0.016)82
σb302=0.25–0.01(0.036)96–0.009(0.037)960(0.041)930.004(0.041)960.002(0.037)97
covb10,b11=0.03–0.002(0.009)96–0.004(0.009)94–0.001(0.007)97–0.001(0.01)96–0.003(0.009)96
covb10,b20=0.020(0.011)950(0.011)950(0.01)96–0.001(0.01)94–0.003(0.011)97
covb10,b21=0.04–0.001(0.01)96–0.001(0.01)960(0.009)96–0.002(0.011)940(0.009)97
covb10,b30=00(0.019)96–0.003(0.017)970.001(0.016)95–0.002(0.016)94–0.001(0.017)95
covb11,b20=0.03–0.001(0.01)96–0.001(0.01)96–0.001(0.009)960(0.011)94–0.002(0.01)97
covb11,b21=00(0.008)950(0.008)960(0.008)950.002(0.01)920(0.008)96
covb11,b30=0.060.005(0.014)920.006(0.013)940.002(0.014)930(0.015)95–0.001(0.014)97
covb20,b21=0.080(0.012)950(0.012)95–0.001(0.011)940(0.011)960.002(0.011)96
covb20,b30=0.05–0.001(0.019)95–0.002(0.019)95–0.003(0.018)94–0.001(0.019)960.001(0.018)96
covb21,b30=0.04–0.001(0.017)95–0.001(0.017)96–0.003(0.016)940.009(0.019)920.011(0.018)95
φ1=0.5–0.003(0.101)95–0.003(0.102)940(0.103)980.009(0.113)940.021(0.114)96
φ2=0.20.006(0.074)930.005(0.073)940.005(0.074)970.002(0.077)930.006(0.072)94
φ3=0.3–0.001(0.094)940.001(0.092)94–0.004(0.092)94–0.006(0.093)94–0.003(0.087)95
Conv. rate1110.510.64
Comp. time (s)43.55 (5.58)69.8 (44.59)295.73 (20.22)3766.41 (1429.06)10,861.22 (3735.02)
Close
This Feature Is Available To Subscribers Only

Sign In or Create an Account

Close

This PDF is available to Subscribers Only

View Article Abstract & Purchase Options

For full access to this pdf, sign in to an existing account, or purchase an annual subscription.

Close