Description Usage Arguments Value Author(s) Examples
Call the function to print the output by BUSgibbs.
1 2 |
x |
The BUSfits object output by the function |
... |
not used. |
print the results from the output by BUSgibbs.
Xiangyu Luo
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 | rm(list = ls(all = TRUE))
set.seed(123)
#a toy example, there are 6 samples and 20 genes in each batch
example_Data <- list()
#batch 1
example_Data[[1]] <- rbind(matrix(c(1,1,5,5,10,10,
3,3,7,7,12,12), ncol=6, byrow=TRUE), matrix(c(1,2),nrow=18, ncol=6))
#batch 2
batch2_effect <- c(2,2,2,1,1)
example_Data[[2]] <- rbind(matrix(c(1,1,5,5,10,10,
3,3,7,7,12,12), ncol=6, byrow=TRUE), matrix(c(1,2),nrow=18, ncol=6)) + batch2_effect
#batch 3
batch3_effect <- c(3,2,1,1,2)
example_Data[[3]] <- rbind(matrix(c(1,1,5,5,10,10,
3,3,7,7,12,12), ncol=6, byrow=TRUE), matrix(c(1,2),nrow=18, ncol=6)) + batch3_effect
set.seed(123)
BUSfits <- BUSgibbs(example_Data, n.subtypes = 3, n.iterations = 100, showIteration = FALSE)
print(BUSfits)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.