data-raw/rds-for-each-project-in-meta.R

meta <- readRDS("inst/META/AllExperimentMetaData.RDS")
meta <- janitor::clean_names(meta, case = "big_camel")
saveRDS(meta, file = "inst/META/AllExperimentMetaData.RDS")

length(meta)

real <- readRDS("inst/RAW/AMP1046-001.RDS")
length(real)

demo <- readRDS("inst/DATA/DemoInput.RDS")
length(demo)

demo2 <- demo[-5]
# X <- lapply(demo2, colnames)
# Y <- lapply(real, colnames)
#
# purrr::map2(X, Y, \(x,y) x %in% y)
#

demo_1 <- demo_2 <- demo2
demo_1[[5]] <- demo_1[[5]][1:11,]

demo_2[[5]] <- demo2[[5]][12:17,]

demo_1.1 <- lapply(demo_1, janitor::clean_names, case = "big_camel")
demo_2.1 <- lapply(demo_2, janitor::clean_names, case = "big_camel")
real_1.1 <- lapply(real  , janitor::clean_names, case = "big_camel")

demo_2.2 <- lapply(demo_2.1, function(x) {
  if ("ExperimentId" %in% colnames(x)) {
    x$ExperimentId <- "ANALYSIS02"
  }

  if ("AnalysisType"  %in% colnames(x)) {
    x$AnalysisType <- "ANALYSIS02TypeB"
  }
  return(x)
})

table_names <- c(
  "QCTable",
  "EditingFreq",
  "IndelFreq",
  "CRISPRSummary",
  "SampleMeta",
  "AlleleAlignFreq"
)

names(demo_1.1) <- table_names
names(demo_2.2) <- table_names
names(real_1.1) <- table_names

real_1.1$EditingFreq$SampleName <- rep(real_1.1$SampleMeta$SampleName, 2)[-1]

saveRDS(object = real_1.1, file = "inst/DATA/AMP1046-001.RDS")
saveRDS(object = demo_1.1, file = "inst/DATA/ANALYSIS01.RDS")
saveRDS(object = demo_2.2, file = "inst/DATA/ANALYSIS02.RDS")
teofiln/gene.editing.dash documentation built on Feb. 21, 2022, 12:59 a.m.