Problem with Bioconductor ArrayExpressHTS package

[ilya.kuzovkin]

Hi,

I'm trying to make use of ArrayExpressHTS R package, but I hit the wall in the very beginning.

I used this tutorial An Introduction to ArrayExpressHTS.

And this is my code:

# Set up external tools
options('ArrayExpressHTS.fasta_formatter' = '/Applications/Bioinformatics/fastx_toolkit')
options('ArrayExpressHTS.cufflinks' = '/Applications/Bioinformatics/cufflinks')
options('ArrayExpressHTS.samtools' = '/opt/local/bin')
options('ArrayExpressHTS.bwa' = '/Applications/Bioinformatics/bwa')
options('ArrayExpressHTS.mmseq' = '/Applications/Bioinformatics/mmseq')
options('ArrayExpressHTS.bam2hits' = '/Applications/Bioinformatics/mmseq')
options('ArrayExpressHTS.bowtie' = '/opt/local/bin')
options('ArrayExpressHTS.tophat' = '/opt/local/bin')

# Load libraries
library("ArrayExpressHTS")

# Get data
aehts <- ArrayExpressHTS("E-GEOD-16190")

The problem is that I get some kind of unknown error from ArrayExpressHTS function:

> aehts <- ArrayExpressHTS("E-GEOD-16190", dir = expfolder)
Fri Jan 18 18:20:21 2013 [AEHTS]  Not Found.
Error in ArrayExpressHTS("E-GEOD-16190") :

Worst part is that it does not say what went wrong, just "Error in ArrayExpressHTS ..." and nothing else.

[ilya.kuzovkin]

Turned out to be stupid question.
The answer is that I did not specify reference directory.
This works fine:

aehts <- ArrayExpressHTS("E-GEOD-16190", refdir='reference')