Entering edit mode
Laurent David
▴
10
@laurent-david-4759
Last seen 9.6 years ago
Hi All,
I'm analyzing mouse samples on beadarrays, with the R beadaray
package. Using the beadarray help files, I managed to normalize my
data, and I know how to annotate them. Now I'd like to clean up the
data:
- Get rid of multiple entries per genes: my favorite option would be
to keep only the median values for genes that are represented by more
than 3 probes, and keep only the highest intensity values for genes
that are represented by 2 probes.
- Get rid of the genes that have a low expression (low intensity) and
that have a low variance across the samples (genes that are probably
not expressed).
My level in R is very limited, I can adapt scripts that I find online
to my needs, but I can't really come up with a fresh script from
scratch. If you have any links that point out to tutorials that
ressemble what I'm trying to do, or other tools that allows to do what
I'm trying to do, that would be of great help.
Thanks in advance,
Sincerely,
Laurent
(PS: it's my first post here, so please do not hesitate if the content
is inappropriate).
### sessionInfo()
R version 2.12.2 (2011-02-25)
Platform: i386-pc-mingw32/i386 (32-bit)
locale:
[1] LC_COLLATE=English_Canada.1252 LC_CTYPE=English_Canada.1252
[3] LC_MONETARY=English_Canada.1252 LC_NUMERIC=C
[5] LC_TIME=English_Canada.1252
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] beadarray_2.0.6 Biobase_2.10.0
loaded via a namespace (and not attached):
[1] limma_3.6.9
[[alternative HTML version deleted]]