library(ggplot2)
library(scales)
library(reshape2)
# This function prints the mean and variance of Median Absolute Deviation (MAD), Mean Absolute Error
# (MAE) and Root Mean Square Error (RMSE) of each of the three main simualtion parameters,
# including library size, bin mean, and non-zero cell proportion, over all cell types in a
# dataset.
# Inputs:
# name: The name of the dataset ("Buenrostro2018", "Cusanovich2018", "PBMCs").
# Output: -
analyze_MAD_MAE_RMSE_stats <- function(name){
print(paste("Results/", name, "/", name, "_stats_MAD_MAE_RMSE_revision.txt", sep = ""))
data <- read.table(paste("Results/", name, "/", name, "_stats_MAD_MAE_RMSE_revision.txt", sep = ""), header = FALSE, dec = ".")
print(data)
data <- data[grep("_0_0_1$", data$V2),]
data.mean <- data[which(data$V3 == "bin_mean"),]
data.lib <- data[which(data$V3 == "lib_size"),]
data.nzp <- data[which(data$V3 == "non_zero_proportion"),]
mat <- matrix(0, nrow = 3, ncol = 3)
mat[1,] <- c(mean(data.nzp[,4]), mean(data.nzp[,5]), mean(data.nzp[,6]))
mat[2,] <- c(mean(data.mean[,4]), mean(data.mean[,5]), mean(data.mean[,6]))
mat[3,] <- c(mean(data.lib[,4]), mean(data.lib[,5]), mean(data.lib[,6]))
sd <- matrix(0, nrow = 3, ncol = 3)
sd[1,] <- c(sd(data.nzp[,4]), sd(data.nzp[,5]), sd(data.nzp[,6]))
sd[2,] <- c(sd(data.mean[,4]), sd(data.mean[,5]), sd(data.mean[,6]))
sd[3,] <- c(sd(data.lib[,4]), sd(data.lib[,5]), sd(data.lib[,6]))
colnames(mat) <- c("MAD", "MAE", "RMSE")
rownames(mat) <- c("Non-zero proportion", "Bin mean", "Library size")
colnames(sd) <- c("MAD", "MAE", "RMSE")
rownames(sd) <- c("Non-zero proportion", "Bin mean", "Library size")
print('Mean:')
print(mat)
print('standard deviation:')
print(sd)
print(scientific_10x(sd))
}
analyze_MAD_MAE_RMSE_stats("Buenrostro2018")
analyze_MAD_MAE_RMSE_stats("Cusanovich2018")
analyze_MAD_MAE_RMSE_stats("PBMCs")
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.