how to add contrast for lmFit in limma?
0
0
Entering edit mode
Jianping Jin ▴ 890
@jianping-jin-1212
Last seen 10.2 years ago
Dear list: I am confused in how to add contrasts in the design table. I added "batch", RNA amplification batches, as a factor to see and or to try avoid batch effect on the results. The three factors are normal/disease, smoke/non-smoke, and batch 1/batch 2. I am interested to know if these 3 variables have any effect on gene expression. It appears to me that there is no way to use factor effects model (with intercept), which could make design easier. For the cell means model design as shown below there are so many combinations that needs to build contrasts and I am not sure of an efficient way. I will greatly appreciate it if anyone could help me out! The design is shown with the column notations as the following: n: normal p: patient N: non-smoke Y: smoke 1 and 2: batch > design nN1 nN2 nY1 nY2 pN1 pN2 pY1 PY2 1 1 0 0 0 0 0 0 0 2 0 0 1 0 0 0 0 0 3 0 0 1 0 0 0 0 0 4 1 0 0 0 0 0 0 0 5 0 1 0 0 0 0 0 0 6 0 1 0 0 0 0 0 0 7 0 0 0 1 0 0 0 0 8 0 0 0 1 0 0 0 0 9 0 0 0 0 0 0 1 0 10 0 0 0 0 0 0 1 0 11 0 0 0 0 1 0 0 0 12 0 0 0 0 0 0 0 1 13 0 0 0 0 0 1 0 0 14 0 0 0 0 0 0 1 0 15 0 0 0 0 0 1 0 0 16 0 0 0 0 0 1 0 0 17 0 0 0 0 1 0 0 0 18 0 0 0 0 0 0 0 1 attr(,"assign") [1] 1 1 1 1 1 1 1 1 attr(,"contrasts") attr(,"contrasts")$treatment [1] "contr.treatment" For your convenience I listed the target below as well: >> target > Sample Expt Smoke Batch > 3 LCM_06CR.CEL n N 1 > 5 LCM_11A.CEL n Y 1 > 7 LCM_14A.CEL n Y 1 > 8 LCM_15BR.CEL n N 1 > 14 LCM_30A.CEL n N 2 > 15 LCM_32A.CEL n N 2 > 16 LCM_33B.CEL n Y 2 > 17 LCM_34A.CEL n Y 2 > 1 LCM_04A.CEL p Y 1 > 2 LCM_05A.CEL p Y 1 > 4 LCM_08BR.CEL p N 1 > 6 LCM_13B.CEL p Y 2 > 9 LCM_16B1.CEL p N 2 > 10 LCM_19A.CEL p Y 1 > 11 LCM_22B.CEL p N 2 > 12 LCM_24B.CEL p N 2 > 13 LCM_27A.CEL p N 1 > 18 LCM_37B.CEL p Y 2 Jianping xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx x Jianping Jin Ph.D. x x Bioinformatics scientist x x Center for bioinformatics x x 3133 Bioinformatics Building x x CB# 7104 x x University of North Carolina x x Chapel Hill, NC 27599 x x Tel: (919)843-6105 x x Fax: (919)843-3103 x x E-mail: jjin at email.unc.edu x xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
• 863 views
ADD COMMENT

Login before adding your answer.

Traffic: 471 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6