compare_RML_HPS.Rd
Simulates data from a simple linear mixed effects model, then calculates REML and HPS effect size estimators as described in Pustejovsky, Hedges, & Shadish (2014).
compare_RML_HPS(iterations, beta, rho, phi, design, m, n, MB = TRUE)
number of independent iterations of the simulation
vector of fixed effect parameters
intra-class correlation parameter
autocorrelation parameter
design matrix. If not specified, it will be calculated based on m
, n
, and MB
.
number of cases. Not used if design
is specified.
number of measurement occasions. Not used if design
is specified.
If true, a multiple baseline design will be used; otherwise, an AB design will be used. Not used if design
is specified.
A matrix reporting the mean and variance of the effect size estimates and various associated statistics.
Pustejovsky, J. E., Hedges, L. V., & Shadish, W. R. (2014). Design-comparable effect sizes in multiple baseline designs: A general modeling framework. Journal of Educational and Behavioral Statistics, 39(4), 211-227. doi:10.3102/1076998614547577
compare_RML_HPS(iterations=10, beta = c(0,1,0,0), rho = 0.3,
phi = 0.5, design=design_matrix(m=3,n=8))
#> Warning: 'g_REML()' is deprecated and may be removed in a later version of the package. Please use 'g_mlm()' instead.
#> Warning: 'g_REML()' is deprecated and may be removed in a later version of the package. Please use 'g_mlm()' instead.
#> Warning: 'g_REML()' is deprecated and may be removed in a later version of the package. Please use 'g_mlm()' instead.
#> Warning: 'g_REML()' is deprecated and may be removed in a later version of the package. Please use 'g_mlm()' instead.
#> Warning: 'g_REML()' is deprecated and may be removed in a later version of the package. Please use 'g_mlm()' instead.
#> Warning: 'g_REML()' is deprecated and may be removed in a later version of the package. Please use 'g_mlm()' instead.
#> Warning: 'g_REML()' is deprecated and may be removed in a later version of the package. Please use 'g_mlm()' instead.
#> Warning: 'g_REML()' is deprecated and may be removed in a later version of the package. Please use 'g_mlm()' instead.
#> Warning: 'g_REML()' is deprecated and may be removed in a later version of the package. Please use 'g_mlm()' instead.
#> Warning: 'g_REML()' is deprecated and may be removed in a later version of the package. Please use 'g_mlm()' instead.
#> mean var
#> p_beta 0.9847397 0.073613086
#> r_theta 1.2720722 0.460753406
#> delta_AB 0.9289647 0.079624297
#> nu 6.9054883 10.142207845
#> kappa 0.3545964 0.006256692
#> g_AB 0.8093198 0.079686515
#> V_g_AB 0.2467520 0.002796612
#> cnvg_warn 0.0000000 0.000000000
#> sigma_sq 0.8119100 0.103050307
#> phi 0.4551261 0.077315335
#> Tau.id.var(constant) 0.4601622 0.537825084
#> rho 0.2786148 0.081782119
#> RML_coverage1 1.0000000 0.000000000
#> RML_coverage2 1.0000000 0.000000000
#> D_bar 1.1252940 0.087464761
#> S_sq 1.2659887 0.438082999
#> delta_hat_unadj 1.0529873 0.058389798
#> phi_hat 0.1862142 0.038234799
#> rho_hat 0.4724283 0.071940262
#> kappa_hat 0.3213402 0.008799495
#> df 6.5016596 15.608659515
#> delta_hat 0.8851809 0.045195879
#> V_delta_hat 0.3008532 0.009866625
#> HPS_coverage1 1.0000000 0.000000000
#> HPS_coverage2 1.0000000 0.000000000
#> cov NA 0.036446512