Question: multifactor microarray data sample analysis using limma package
gravatar for jaichandpatel2
5.8 years ago by
jaichandpatel210 wrote:

i had microarray data sample files (.CEL files). i had used limma package & found out topTable of differentially expressed genes in R . But now i have some samples in which more than one factors have been applied i.e. strain, time, treatment etc. So i am not able to construct design matrix. Although in limma user guide has given an example of estrogen datasets but i didn't get it.

please help me if anyone have any idea regarding object of experimental design. How can we add factor with raw data i.e. .CEL files. please give the answer with any example if you can.

ADD COMMENTlink modified 5.8 years ago by Devon Ryan97k • written 5.8 years ago by jaichandpatel210

What part of that user guide didn't you get exactly?

ADD REPLYlink written 5.8 years ago by Devon Ryan97k

9.5 Interaction Models: 2  2 Factorial Designs

can you help me to understand above part ?

ADD REPLYlink written 5.8 years ago by jaichandpatel210

Yes, what part of that section don't you understand?

ADD REPLYlink written 5.8 years ago by Devon Ryan97k

i have a data GSE35819. i want to analyse this data.

In this data i have 18 samples (18 .CEL files)

3 cell line, 3 time points. can yoy tell me how to design matrix, contrast matrix.

i have some more data of this kind.

plz help as much you can.


ADD REPLYlink modified 5.8 years ago • written 5.8 years ago by jaichandpatel210
gravatar for Devon Ryan
5.8 years ago by
Devon Ryan97k
Freiburg, Germany
Devon Ryan97k wrote:

The general idea with GSE35819 is:

design <- data.frame(cellLine=factor(c(rep("H9",3), rep("HS401",6), rep("HS360",6), rep("H9",3))),
    timePoint = factor(c(24,7,7,2,2,24,24,7,7,2,2,24,24,7,7,2,2,24)),
    treatment=factor(c(rep(c("Normoxia","Hypoxia"), 12)))

The model matrix is then dependent on what you want to test. Perhaps you just care about "treatment" while controlling for everything else (excluding interactions):

mm <- model.matrix(~cellLine+timePoint+treatment, design)

You don't always need contrast matrices, it depends on the model and the question. In general, you should either take a class on linear models or just talk to a local statistician about individual designs. Providing a general tutorial on linear model is well beyond the scope of this site (that's why there are classes and books dedicated to the subject).

ADD COMMENTlink written 5.8 years ago by Devon Ryan97k

Thank you Sir.

ADD REPLYlink written 5.8 years ago by jaichandpatel210

Can i use maSigPro package for this. As i used maSigPro package and got this error. What is meaning of this error, should i change samething in this R-script.

If you have any idea then suggest me.

> edesign
        Time Replicate Control hypoxic
Array1     3         1       1       2
Array2     2         1       2       1
Array3     2         1       1       2
Array4     1         2       2       1
Array5     1         2       1       2
Array6     3         2       2       1
Array7     3         2       1       2
Array8     2         2       2       1
Array9     2         2       1       2
Array10    1         3       2       1
Array11    1         3       1       2
Array12    3         3       2       1
Array13    3         3       1       2
Array14    2         3       2       1
Array15    2         3       1       2
Array16    1         1       2       1
Array17    1         1       1       2
Array18    3         1       2       1


fit <- p.vector(eset, design, Q = 0.05, MT.adjust = "BH", min.obs = 20)
Error in dat[, as.character(rownames(dis))] : subscript out of bounds

ADD REPLYlink written 5.7 years ago by jaichandpatel210

I know nothing about the maSigPro package.

ADD REPLYlink written 5.7 years ago by Devon Ryan97k

Hello Sir

I have got p-adjusted value i.e. q-value using limma package, but how can i get the gene list of different coefficients in an ascending or descending as well in .xls format. what should be the cut off in case of q-value to select differentially expressed genes (upregulated & downregulated ), Will it be same as in p-value (0.05 or 0.01) or something else? Please Sir suggest me, what should i do ?

Thankyou very much

ADD REPLYlink modified 5.7 years ago • written 5.7 years ago by jaichandpatel210

Hi Devon,

Kindly shed some light on following post, thanks.

Microarray analyis of different platform (GPL96 vs. GPL3921)

ADD REPLYlink written 5.3 years ago by Bioinformatist Newbie250
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1222 users visited in the last hour