[BioC] a question on creating model.matrix using limma
Srinivas Iyyer
srini_iyyer_bio at yahoo.com
Wed Jul 12 00:43:24 CEST 2006
Dear Dr. Smith and group,
apologies for simple question.
I have data with simple plan. There are 8 samples (4
normal and 4 tumor). These tissue samples are from
different patients. No replicate experiments were
done.
I wish to see genes differentially expressed between
Tumors and normals.
In such cases how can I create model.matrix
Following is taken from Limma manual:
design <- model.matrix(~ -1 +
factor(c(1,1,1,2,2,3,3,3)))
my question:
1. what is '~ -1' is for?
2. Will this work for my study structure. That means
I have only two types of RNA 1. tumor and 2. Normal.
design <- model.matrix(~ -1 +
factor(c(1,1,1,1,2,2,2,2))).
Thanks
sri
More information about the Bioconductor
mailing list