swh:1:snp:fafdc4d55567af0f04b069bd7bca576c0cea9371
Raw File
Tip revision: c20430824c5a603eec695ee167e32e4cd1b6c024 authored by Aria on 20 December 2020, 07:48:44 UTC
Update README.md
Tip revision: c204308
run_saver.R
library(SAVER)

transform_by_saver <- function(mtx){
  mtx <- t(mtx)
  mtx.saver <- saver(mtx, ncores = 20,estimates.only = TRUE,size.factor = 1)
  return(mtx.saver)
}

write_data <- function(mtx,out){
  write.csv(t(mtx),out,row.names = FALSE)
}

run_main <- function(prefix){
  sample_hela_select <- read.csv(paste0('../result/scRNA_seq_2020_11_1-7/',prefix,'_hela_select.csv'))
  sample_hela_select_saver <- transform_by_saver(sample_hela_select)
  write_data(sample_hela_select_saver,paste0('../result/scRNA_seq_2020_11_1-7/',prefix,'_hela_select_saver.csv'))

  sample_f9_select <- read.csv(paste0('../result/scRNA_seq_2020_11_1-7/',prefix,'_f9_select.csv'))
  sample_f9_select_saver <- transform_by_saver(sample_f9_select)
  write_data(sample_f9_select_saver,paste0('../result/scRNA_seq_2020_11_1-7/',prefix,'_f9_select_saver.csv'))
}

# there is a side effect for saver in R, which is the '-' in gene symbol will be replaced with '.'
names = c('A1A2','B1A3','C1C2')
lapply(names,run_main)
back to top