Question: multifactor microarray data sample analysis using limma package
0
5.8 years ago by
India/Delhi/DIPAS-DRDO
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.

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?

9.5 Interaction Models: 2  2 Factorial Designs

can you help me to understand above part ?

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

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.

0
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).

Thank you Sir.

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

I know nothing about the maSigPro package.

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

Hi Devon,

Kindly shed some light on following post, thanks.

Microarray analyis of different platform (GPL96 vs. GPL3921)