[BioC] Limma: technical replicates in multi-group design

Wei Shi shi at wehi.EDU.AU
Thu Mar 25 23:37:26 CET 2010


Dear Adam:

    You can use the following commands to create a design matrix for 
your experiment:

s <- factor(c("A","A","B","B","C","C"))
design <- model.matrix(~0+s)
colnames(design) <- levels(s)
fit <- lmFit(x,design)

    To do a pairwise comparison, you can use the makeConstrasts function:

contr <- makeContrasts(A-B,A-C,B-C,levels=design)
fit.contr <- eBayes(contrasts.fit(fit,contr))

Cheers,
Wei

Adam Kiezun wrote:
> Hi,
> How do I write the design matrix to specify technical replicates in
> multi-group experiment on a single-channel array (data coming from
> lumi)?
>
> I have 3 groups of samples (A, B, C), each group has 2 technical
> replicates (so I have 6 expression vectors: A1, A2, B1, B2, B3). How
> do I write the design matrix to tell limma that A1 and A2 etc are
> technical replicates and that I want to do all pairwise comparisons
> between groups (ie. AvsB, AvsC, BvsC)?
>
> The limma manual covers technical replicates for two-channel arrays,
> and the multi-group experiments with biological replicates. I'd like
> to know how to deal with a one-channel arrays, multi-group experiment
> with technical replicates.
>
> regards
> Adam
>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor
>   

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



More information about the Bioconductor mailing list