[BioC] Duplicated spots in limmaGUI
João Fadista
Joao.Fadista at agrsci.dk
Thu Apr 19 08:31:22 CEST 2007
Hi,
You can combine duplicate spots by ID with a script similar like this (assuming that your duplicated clones have a smaller ID than the ones that are not duplicated):
index <- order(MA$genes$ID)
MAsorted <- MA[index,]
MAsorted$printer$spacing <- 1
And then average the clones in the subset of the MAList where you have duplicated clones. Let say that your duplicated clones have an ID below 200000:
Subset <- MAsorted$genes$ID < 200000
MAsubset <- MA[Subset,]
If your clones that are duplicated have not a smaller ID than the ones that are not duplicated, you can make an extra column in your MAList where you can put a dummy ID where this is true.
Best regards
João Fadista
Ph.d. student
UNIVERSITY OF AARHUS
Faculty of Agricultural Sciences
Dept. of Genetics and Biotechnology
Blichers Allé 20, P.O. BOX 50
DK-8830 Tjele
Phone: +45 8999 1900
Direct: +45 8999 1900
E-mail: Joao.Fadista at agrsci.dk
Web: http://www.agrsci.org
This email may contain information that is confidential.
Any use or publication of this email without written permission from Faculty of Agricultural Sciences is not allowed.
If you are not the intended recipient, please notify Faculty of Agricultural Sciences immediately and delete this email.
-----Original Message-----
From: bioconductor-bounces at stat.math.ethz.ch [mailto:bioconductor-bounces at stat.math.ethz.ch] On Behalf Of Wonjong Moon
Sent: Thursday, April 19, 2007 2:46 AM
To: bioconductor at stat.math.ethz.ch
Subject: [BioC] Duplicated spots in limmaGUI
Our slide design had some spots were duplicated and some were not.
I am using limmaGUI, but duplicated spots were recognized by space between spots.
How can I handle the spots that were not duplicated?
Should add flags in spottype file and exclude those?
Is there any way to combine duplicated spots by ID or name in Limma?
Thank you.
Wonjong.
Bioinformatician
SBRI.
_______________________________________________
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
More information about the Bioconductor
mailing list