Nothing
data_dir = system.file("extdata/MAE.rds", package = "animalcules")
MAE = readRDS(data_dir)
alpha.methods <- c("inverse_simpson", "gini_simpson", "shannon")
beta.methods <- c("wUniFrac", "bray")
tax.default <- "genus"
# These are options for rendering datatables
dtopts <- list(scrollX=TRUE, paging=TRUE)
tax.name <- colnames(rowData(MAE[['MicrobeGenetics']]))
sam.name <- rownames(colData(MAE[['MicrobeGenetics']]))
org.name <- rownames(as.data.frame(assays(MAE[['MicrobeGenetics']])))
measure.type <- c('Final Guess', 'Final Best Hit', 'Final High Confidence Hit')
minbatch <- function(batch1){
batch2 <- as.factor(batch1)
batch3 <- split(batch1,batch2)
return(min(unlist(lapply(seq_len(length(batch3)),
function(x) length(batch3[[x]])))))
}
covariates = colnames(colData(MAE))
sam_temp <- as.data.frame(colData(MAE[['MicrobeGenetics']]))
num_select <- lapply(covariates, function(x) is_categorical(unlist(sam_temp[,x])))
num_covariates <- covariates[!unlist(num_select)]
# choose the covariates that has less than 8 levels
covariates.colorbar <- c()
for (i in seq_len(length(covariates))){
num.levels <- length(unique(colData(MAE)[[covariates[i]]]))
if (num.levels < 8){
covariates.colorbar <- c(covariates.colorbar, covariates[i])
}
}
# choose the covariates that has 2 levels
covariates.two.levels <- c()
for (i in seq_len(length(covariates))){
num.levels <- length(unique(colData(MAE)[[covariates[i]]]))
if (num.levels == 2){
covariates.two.levels <- c(covariates.two.levels, covariates[i])
}
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.