In dshenker/RFSLAM: What the Package Does (One Line, Title Case)

knitr::opts_chunk$set(echo = TRUE)

df <- read_csv("~/workspace/SAFE/code/SLAM_ML/dshenker/RFSLAM_Package/RFSLAM/vignettes/sim_data.csv")
head(df)

df[df$subjid == "a101083",]

df[df$visit_type == "M24","subjid"]

df[df$visit_type == "Exit",]

unique(df[df$visit_type == "Exit", "subjid"])

unique(df[df$visit_type == "M24", "subjid"])

response_cols <- c("final_cohaoutb___1", "final_cohaoutb___2", "final_cohaoutb___3", "final_cohaoutb___4", "final_cohaoutb___5", "final_cohaoutb___6", "final_cohaoutb___7", "final_cohaoutb___8", "final_cohaoutb___9", "final_cohaoutb___10", "final_cohaoutb___11")
ids <- unique(df$subjid)
for (i in 1:length(ids)) {
  id = ids[i]
  print(sum(rowSums(is.na(df[df$subjid == id, response_cols])) != 11))
}

response_cols <- c("final_cohaoutb___1", "final_cohaoutb___2", "final_cohaoutb___3", "final_cohaoutb___4", "final_cohaoutb___5", "final_cohaoutb___6", "final_cohaoutb___7", "final_cohaoutb___8", "final_cohaoutb___9", "final_cohaoutb___10", "final_cohaoutb___11")
ids <- unique(df$subjid)
for (i in 1:length(ids)) {
  id = ids[i]
  print(id)
  print(rowSums(df[df$subjid == id, response_cols], na.rm = TRUE))
}

df[df$subjid == "a106084", c("visit_type", "bl_visdat", response_cols)]

dshenker/RFSLAM documentation built on Sept. 18, 2022, 3:26 a.m.