R14

raw_PDUI = read.table("header.txt", header = TRUE, stringsAsFactors = F)
source("./Processing_raw_PDUI_matrix.R")
raw_PDUI_process.res = Processing_raw_PDUI_matrix(raw_PDUI,filter_gene_thre = 0.2,filter_cell_thre = 0.1)
raw_PDUI_matrix_sc = raw_PDUI_process.res$raw_PDUI_matrix_sc
raw_PDUI_matrix_sNA = raw_PDUI_process.res$raw_PDUI_matrix_sNA
source("./Find_Neighboring_Cells.R")
Neighboring_Cell.res = Find_Neighboring_Cells(raw_PDUI_sNA = raw_PDUI_matrix_sNA)
cluster_info = Neighboring_Cell.res$cluster_info
num_clusters = Neighboring_Cell.res$num_clusters
outliers = Neighboring_Cell.res$outliers
source("./imputation_model.R")
Imputation.res = Impuation_model(raw_PDUI_sNA = raw_PDUI_matrix_sNA,raw_PDUI_sc = raw_PDUI_matrix_sc,num_clusters = num_clusters,cluster_info = cluster_info)
PDUI_sc_imputed = Imputation.res
write.table(PDUI_sc_imputed, paste0("./scDaPars", "/scDaPars_imputed_results.txt"), quote = FALSE)