R/DIF_designAnalysisFunctions.R

In lucasjfriesen/jamoviPsychoPDA: Psychometrics & Post-Data Analysis

designAnalysis.nagR2 <-
  function(designList,
           Data,
           group,
           match,
           bootSims,
           type,
           hypTrueEff,
           alpha,
           difFlagScale,
           sigOnly) {
    if (hypTrueEff == "") {
      if (difFlagScale == "zt") {
        hypTrueEff <- c(0, 0.13, 0.26)
      } else {
        hypTrueEff <- c(0, 0.035, 0.07)
      }
      labels <- c("A", "B", "C")
    } else {
      hypTrueEff <- as.numeric(hypTrueEff)
      labels <- c("Custom Hyp.")
    }
    
    GC <-
      data.frame(
        "label" = rep(as.character(), times = length(designList) * length(hypTrueEff)),
        "item" = rep(as.character(), times = length(designList) * length(hypTrueEff)),
        "obsEff" = rep(as.numeric(), times = length(designList) *
                         length(hypTrueEff)),
        "hypTrueEff" = rep(as.numeric(), times = length(designList) *
                             length(hypTrueEff)),
        "typeM" = rep(as.numeric(), times = length(designList) *
                        length(hypTrueEff)),
        "power" = rep(as.numeric(), times = length(designList) *
                        length(hypTrueEff)),
        "bootSE" = rep(as.numeric(), times = length(designList) *
                         length(hypTrueEff)),
        stringsAsFactors = FALSE
      )
    
    
    for (item in 1:length(designList)) {
      empDATA <-
        cbind(Item = jmvcore::toNumeric(Data[, designList[item]]),
              as.factor(group),
              match)
      colnames(empDATA) <-
        c(colnames(Data)[item], "GROUP", "SCORES")
      
      myBoot <-
        boot.empDist(empDATA,
                     R = bootSims,
                     type = type,
                     coefficients = FALSE)
      
      tick <- length(hypTrueEff) - 1
      
      for (hypInd in 1:length(hypTrueEff)) {
        retroDesignRes <-
          retroDesign.nagR2(
            hypTrueEff = hypTrueEff[hypInd],
            myBoot,
            alpha = alpha,
            sigOnly = sigOnly
          )
        
        GC[item * length(hypTrueEff) - tick, 1] <-
          labels[hypInd]
        GC[item * length(hypTrueEff) - tick, 2] <-
          designList[item]
        GC[item * length(hypTrueEff) - tick, 3] <-
          as.numeric(retroDesignRes[5])
        GC[item * length(hypTrueEff) - tick, 4:7] <-
          as.numeric(retroDesignRes[1:4])
        tick <- tick - 1
      }
    }
    
    return(GC)
  }

designAnalysis.coefficients <- function(designList,
                                        hypTrueEff,
                                        coefficient,
                                        coefficientsSE,
                                        alpha,
                                        df,
                                        difFlagScale,
                                        sigOnly) {
  if (hypTrueEff == "") {
    hypTrueEff <- 0
    labels <- "Null"
  } else {
    hypTrueEff <- as.numeric(hypTrueEff)
    labels <- c("Custom Hyp.")
  }
  retroDesignRes <- list()
  
  for (item in 1:length(designList)) {
    retroDesignRes[[item]] =
      retroDesign.coefficients(
        hypTrueEff = hypTrueEff,
        coefficient = coefficient[item,],
        coefficientsSE = coefficientsSE[item,],
        alpha = alpha,
        df = df,
        sigOnly = sigOnly
      )
  }
      return(retroDesignRes)
}