Autarchy of the Private Cave

1. Install the annotationTools R package:
   source(“http://bioconductor.org/biocLite.R”)
   biocLite(“annotationTools”)
2. Download full HomoloGene data file from ftp://ftp.ncbi.nlm.nih.gov/pub/HomoloGene/current
3. library(annotationTools)
4. homologene = read.delim(“homologene.data”, header=FALSE)
5. mygenes = read.table(“file with one entrez ID of the source organism per line.txt”)
6. getHOMOLOG(unlist(mygenes), taxonomy_ID_of_target_organism, homologene) [alternatively, wrap the call to getHOMOLOG into unlist to get a vector]

It might be easier to achieve the same results with a Perl script calling NCBI’s e-utils.

• R time series tutorial (2010, a website of the “Time Series Analysis and Its Applications: With R Examples” book)
• Statistics with R (2007)
• R for programmers PDF (2008, 104 pages, linked to from here)
• Brief R tutorial (2004)
• Statistical computing with R: a tutorial (2004)
• An introduction to R (from the official r-project website, should be always up-to-date)
• R tutorial (date unknown, definitely newer than 2005)

Sometimes there is a need to remove all the probesets, which have expression values below the minimal spike-in intensity on the Affymetrix microarray. The reasoning behind this procedure is simple: minimal-expression spike-ins represent the bottom margin of microarray sensitivity, and anything below that margin cannot be reliably quantified – which also means that both fold-change and p-value of expression variance will be unreliable for these probesets.

Here’s a simple R script to do just that. It is abundantly commented, and also contains an optional (commented out) fragment which allows the removal of more low-variance, low-intensity probesets.

Read the rest of this entry »

If you get this message when opening vignettes:

Error in openPDF(vif) :
getOption(‘pdfviewer’) is ”; please use ‘options(pdfviewer=…)’

and you are tired of running this command every time:

> options(pdfviewer=”okular”)

then you should check if your system-wide Renviron file has proper PDF viewer set:
Read the rest of this entry »

First, learn about custom CDFs and why they are needed.

The aroma.affymetrix R package google group has a how-to: create a CDF annotation file from scratch.

Also useful: how to convert CDF into an R package, which has all CDF data available (as a PDF with more details).