[BioC] CellMix package use with own data matrix

Renaud Gaujoux renaud at mancala.cbio.uct.ac.za
Fri May 9 12:41:55 CEST 2014


Hi Maria,

note that you do not need to use an ExpressionSet object. Plain matrix
objects should work as well, as long as the rownames match:

x <- as.matrix(read.table('mixed_
data.tsv', sep = "\t"))
y <- as.matrix(read.table('celltype_data.tsv', sep = "\t"))
res <- gedProportions(x, y, rescale = FALSE, log = FALSE, normalize = FALSE)

# unscaled proportions
coef(res)
# rescaled proportions (sum up to one)
scoef(res)

If you later want to create an ExpressionSet then it just about
wrapping the matrix object:

eset <- ExpressionSet(x)

but this does not make much sense if you do not have extra phenotypic
data or gene annotation package name to associate your gene expression
data with.

Hope this helps.
Please let me know if this works fine for you.

Bests,
Renaud



More information about the Bioconductor mailing list