search

claycantrell / PartInvShinyUI

MMM Lab Partial Invariance Shiny UI
GNU General Public License v3.0
0 stars 2 forks source link

Test equivalence of `PartInvMulti_we()` and `PartInv()` for single dimensions #11

Closed marklhc closed 2 years ago

YichiZhang2024 commented 3 years ago

Tested the two functions with the NSI example. The two functions gave the same result.

PartInv(cut_z = 9, kappa_r = 0, kappa_f = 0, phi_r = 1, phi_f = 1, lambda_r = c(.6788576, .8644396,  .8817195), lambda_f = c( .7467853,  .8586934,  .8284677), tau_r = c(1.492308, 1.492308, 1.612088), tau_f = c(1.32531, 1.323569,  1.332582), Theta_r = c(0.4770053, 0.3356521, 0.4490181), Theta_f = c(.3906694, .3017572, .4293952))
PartInvMulti_we(cut_z = 9,
                weights_item = c(1, 1, 1),
                weights_latent = 1,
                alpha_r = 0,
                alpha_f = 0,
                psi_r = 1,
                psi_f = 1,
                lambda_r = c(.6788576, .8644396,  .8817195),
                lambda_f = c( .7467853,  .8586934,  .8284677),
                nu_r = c(1.492308, 1.492308, 1.612088),
                nu_f = c(1.32531, 1.323569,  1.332582),
                Theta_r = c(0.4770053, 0.3356521, 0.4490181),
                Theta_f = c(.3906694, .3017572, .4293952),
                plot_contour = FALSE)