inst/appliedPAClassificationScripts/runMethods/run_fast_methods.R

#options(warn = 2, error = recover)

for(data_set in c("Mannini", "SasakiFreeLiving", "SasakiLab")) {
#for(data_set in c("SasakiFreeLiving")) {
  for(location in c("ankle", "wrist")) {
    if(data_set %in% c("SasakiFreeLiving", "SasakiLab")) {
      class_var_levels <- c("y_category3", "y_intensity")
    } else {
      class_var_levels <- c("y_category4", "y_intensity")
    }
    for(class_var in class_var_levels) {
       for(fit_method in c("RF", "normalHMM", "normalFMM")) {
#      for(fit_method in c("normalFMM")) {
        setwd(file.path("C:/Stat/HMM/PACwithDDM/inst/results", data_set, location, class_var, fit_method))
        
        for(reduced_trans_mat_parameterization in FALSE) {
          if(identical(data_set, "SasakiFreeLiving")) {
            N <- 15
          } else if(identical(data_set, "SasakiLab")) {
            if(identical(location, "ankle")) {
              N <- 34
            } else {
              N <- 35
            }
          } else if(identical(data_set, "Mannini")) {
            N <- 33
          }
          
          for(subject in seq_len(N)) {
#          for(subject in 1L) {
            system_cmd <- paste0("R --vanilla --args ", subject, " ", data_set, " ", location, " ", class_var, " ", fit_method, " ", reduced_trans_mat_parameterization,
              " < C:/Stat/HMM/PACwithDDM/inst/appliedPAClassificationScripts/runMethods/LOSO_crossval.R ")
            
            system(system_cmd, intern = TRUE)
          }
        }
      }
    }
  }
}
elray1/PACwithDDM documentation built on May 13, 2019, 8:37 a.m.