R/temp/replicbase.R

In ORIONZ.G: EAP Scoring in Exploratory FA Solutions with Correlated Residuals

replicbase<-function(LAM, Ree, PSI, N,display = TRUE){


  ## 1 replic Continuous Omega


  lam <- LAM

  siz <- size(lam)
  n <- siz[1]
  m <- siz[2]


  th <- transpose(runif(N))
  L <- chol(Ree)

  TMP <- matrix(runif(N*n),ncol=n)

  E <- TMP %*% L
  EI <- TMP

  ##Note: Z1 is the data matrix with local dependencies according to Ree
  #       Z2 is the data matrix if the item were locally independent

  Z1 <- transpose((lam %*% transpose(th)) + (PSI %*% transpose(E)))
  Z2 <- transpose((lam %*% transpose(th)) + (PSI %*% transpose(EI)))

  # Note: the omegas in the simulation are obtained as the squared correlations between
  # the sum scores and the true theta levels. This is possible here because the true theta levels are known.


  sumsco1 <- Z1 %*% matrix(1,n,1)
  tmp <- cor(sumsco1,th)
  omrep <- tmp * tmp
  sumsco2 <- Z2 %*% matrix(1,n,1)
  tmp2 <- cor(sumsco2,th)
  omrepli <- tmp2 * tmp2
  relefrep <- (omrep * (1-omrepli)) / (omrepli * (1-omrep))



  OUT<-list('omrep'=omrep,'omrepli'=omrepli, 'relefrep'=relefrep)

  invisible(OUT)


}