[BioC] Help with

Guillaume Meurice guillaume.meurice at igr.fr
Fri May 20 15:44:45 CEST 2011


Sorry, I made some mistake in my previous mail, into the contrast matrix (bad copy paste from another project), so here I have corrected them.

> I have a question regarding the way to properly design biological replicate and technical replicates.
> 
> In the projet, we have two groups of sample : R (respond to the treatment), NR (no response).
> For each groups, there is several replicates (3) that come from different patients (so not strictely biological).
> Each patient have two sample : one at day one(d1), the second after one week (d7)
> The hybridization are planned to be performed in dual color, with dye-sap.
> 
> Here is the target file :
> 
> Cy3	Cy5	Patient	Issu	TechnicalReplicat
> d1	d7	A	R	1
> d7	d1	A	R	1
> d1	d7	B	R	2
> d7	d1	B	R	2
> d1	d7	C	R	3
> d7	d1	C	R	3
> d1	d7	D	NR	4
> d7	d1	D	NR	4
> d1	d7	E	NR	5
> d7	d1	E	NR	5
> d1	d7	F	NR	6
> d7	d1	F	NR	6
> 
> 
> my question is how to properly design the matrix design and the contrast matrix to answer the difference between the two groups ?
> Which column should I use as biological replicat (so as I can use the duplicateCorrelation function) ?
> 
> here is how I've started, but I can't figure out how to use duplicateCorrelation with this design
> 
> design <- cbind(
> 		R_D7vsD1  = c(-1,1,-1,1,-1,1,  0,0, 0,0, 0,0),
> 		NR_D7vsD1 = c( 0,0,0,0,0,0,-1,1,-1,1,-1,1)
> )
> fit <- lmFit(MAn, design, weights = NULL)
> cont.mat = makeContrasts(NRvsR = NR_D7sD1 - R_D7vsD1, levels = design)
> fit2 <- contrasts.fit(fit, cont.mat)
> fit2 <- eBayes(fit2)
> 
> 
> 
> 
> Many thanks by advance for any help.
> 
> --
> Guillaume
> 	[[alternative HTML version deleted]]
> 
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at r-project.org
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor



More information about the Bioconductor mailing list