[BioC] limma matrix design for biological replicates.

Gordon K Smyth smyth at wehi.EDU.AU
Sat Jun 26 08:53:03 CEST 2010


Dear Neeraj Rana,

You targets file is ok, although your samples are really normal and tumour 
rather than wt and mutant.  You can use

   design <- model.matrix(~factor(targets$Cy5))
   colnames(design) <- c("NodNeg","PosvsNeg")
   ...
   topTable(fit, coef="PosvsNeg")

Best wishes
Gordon


> Date: Fri, 25 Jun 2010 13:55:46 +0530
> From: neeraj <kushrn at gmail.com>
> To: Bioconductor at stat.math.ethz.ch
> Subject: [BioC] limma matrix design for biological replicates.
>
> hi,
>
> i am doing analysis for 10 breast cancer arrays (two color Agilent4x44
> arrays) with limma,where Cy3 is normal and Cy5 is tumor.And all are
> biological replicates ,means all 10 arrays have been prepared from 10
> different patients.
> there are two categories as nod negative and nod positive.Five arrays from
> nod negative(5 patients) and five are from nod positive(another 5
> patients).I want to see the differentially regulated genes between two
> categories(NOD NEGATIVE vs NOD POSITIVE).
> i designed the target file as given below..I want to make it sure whether it
> is correct or not.
>
> SampleNumber                 FileName         Cy3              Cy5
> 1                                      1135_NN.txt
> wt1                 mu1
> 2                                      2157_NN.txt              wt1
>        mu1
> 3                                      3159_NN.txt
> wt1                 mu1
> 4                                      4171_NN.txt
> wt1                 mu1
> 5                                      5179_NN.txt              wt1
>                mu1
> 6                                      628_NP.txt
> wt2                 mu2
> 7                                      758_NP.txt
> wt2                 mu2
> 8                                      880_NP.txt
> wt2                 mu2
> 9                                      993_NP.txt
> wt2                 mu2
> 10                                   1096_NP.txt
> wt2                 mu2
>
> 1 to 5 are nod negative and ,6 to 10 are nod positive.
>
> thanx.
> NEERAJ RANA
> JRF
> IISC BANGLORE(INDIA)

______________________________________________________________________
The information in this email is confidential and intend...{{dropped:4}}



More information about the Bioconductor mailing list