From 38fe56c726c6d76d07be82159313138f28f77cbf Mon Sep 17 00:00:00 2001 From: SiWen314 Date: Fri, 16 Aug 2024 09:22:59 -0400 Subject: [PATCH] Modify the modality-case interaction term in sim.NormalIG.Hierarchical.R for the within-modality situation --- Rpackage/iMRMC/R/sim.NormalIG.Hierarchical.R | 19 +++++++++++++++---- 1 file changed, 15 insertions(+), 4 deletions(-) diff --git a/Rpackage/iMRMC/R/sim.NormalIG.Hierarchical.R b/Rpackage/iMRMC/R/sim.NormalIG.Hierarchical.R index 92c85977..475ad332 100644 --- a/Rpackage/iMRMC/R/sim.NormalIG.Hierarchical.R +++ b/Rpackage/iMRMC/R/sim.NormalIG.Hierarchical.R @@ -148,15 +148,26 @@ sim.NormalIG.Hierarchical = function(config,R = NULL,AR = NULL,BR = NULL,is.with if(C_dist == 'normal'){ C = matrix(rep(rnorm(nC,0,config$sigma_C),nR),nR,nC,byrow =TRUE) + ## Mod A/ first replicate tauC.A = matrix(rep(rnorm(nC,0,config$sigma_C.A),nR),nR,nC,byrow =TRUE) - tauC.B = matrix(rep(rnorm(nC,0,config$sigma_C.B),nR),nR,nC,byrow =TRUE) - + ## Mod B/ second replicate + if(is.within){ + tauC.B = tauC.A + }else{ + tauC.B = matrix(rep(rnorm(nC,0,config$sigma_C.B),nR),nR,nC,byrow =TRUE) + } + }else if(C_dist == 'beta'){ C = matrix(rep(rbeta(nC,config$a_C,config$b_C),nR),nR,nC,byrow =TRUE) + ## Mod A/ first replicate tauC.A = matrix(rep(rbeta(nC,config$a_C.A,config$b_C.A),nR),nR,nC,byrow =TRUE) - tauC.B = matrix(rep(rbeta(nC,config$a_C.B,config$b_C.B),nR),nR,nC,byrow =TRUE) - + ## Mod B/ second replicate + if(is.within){ + tauC.B = tauC.A + }else{ + tauC.B = matrix(rep(rbeta(nC,config$a_C.B,config$b_C.B),nR),nR,nC,byrow =TRUE) + } } # RC