[BioC] Agilent duplicate spots and multiple probes per gene
Nathan S. Watson-Haigh
nathan.watson-haigh at csiro.au
Wed May 13 02:51:05 CEST 2009
-----BEGIN PGP SIGNED MESSAGE-----
Hash: SHA1
I'm a bit confused/worried about duplicated spots on the Agilent bovine array,
in particular how best to handle them in limma as I have technical rep dye-swaps.
<code>
> #Get all non-control spots:
> genes.idx <- which(RG$genes$ControlType == 0)
> MA.none <- normalizeWithinArrays(RG[genes.idx,], method="none")
> table(table(MA.none$genes$ProbeUID))
2 6
21465 10
</code>
So 21465 probes are duplicated twice and 10 probes are duplicated 6 times. Does
anyone have suggestions on how best to proceed given that I have technical rep
dye-swaps?
Also, how do I best handle those genes with multiple different probes? Simply
average?
Cheers,
Nathan
- --
- --------------------------------------------------------
Dr. Nathan S. Watson-Haigh
OCE Post Doctoral Fellow
CSIRO Livestock Industries
Queensland Bioscience Precinct
St Lucia, QLD 4067
Australia
Tel: +61 (0)7 3214 2922
Fax: +61 (0)7 3214 2900
Web: http://www.csiro.au/people/Nathan.Watson-Haigh.html
- --------------------------------------------------------
-----BEGIN PGP SIGNATURE-----
Version: GnuPG v1.4.9 (MingW32)
Comment: Using GnuPG with Mozilla - http://enigmail.mozdev.org
iEYEARECAAYFAkoKGXkACgkQ9gTv6QYzVL4RyACeNQUXJABLApt0qF8fU6Yn90On
QY4AoJ+rYzYa8MGgmMt3mzzsBZ68+k8m
=1+Hx
-----END PGP SIGNATURE-----
More information about the Bioconductor
mailing list