[BioC] Error using Ringos readNimblegen

Joern Toedling toedling at ebi.ac.uk
Wed May 28 11:47:41 CEST 2008


Hi Morten,
sorry for the trouble, my first guess would be that the raw data pair 
files are of some unexpected format, possibly having no column named 
"PM". I would be very grateful, if you could send me the first few lines 
of one of those pair files off-list, so that I can have a more detailed 
look and extend the function 'readNgIntensitiesTxt' to handle that 
format as well.
Btw, I advise you to update your R and Bioconductor to the current 
release version (R 2.7.0 and BioC 2.2, released 1 May). Since you are 
using Linux, you are also very welcome to use the current development 
version of Ringo, which is the quickest one to include any changes.
http://www.bioconductor.org/packages/devel/bioc/html/Ringo.html
Best regards,
Joern

Morten Hansen wrote:
>
> I have a problem with the Ringo package.
>
> I am trying to read in the results from a 1.5Kb Human Promoter Array 
> (CHP0150-HP1). I get the following error shown in below in the 
> "Output" section.
>
> The target.txt, spottype.txt and sessionInfo() is shown below the 
> "Output" section.
>
> **************
> *** Output ***
> **************
> > myDir <- c("/home/mhansen/R/ChIP-chip/39037")
> > myRg <- readNimblegen("targets.txt","spottypes.txt", path = myDir)
> Reading targets file...
> Reading raw intensities...
> Read header information
> Read /home/mhansen/R/ChIP-chip/39037/39037_532_pair.txt
> Read /home/mhansen/R/ChIP-chip/39037/39037_635_pair.txt
> Error in RG$R[, i] <- obj[, columns$f] :
>   number of items to replace is not a multiple of replacement length
> >
>
>
>
> The target.txt and spottype.txt look like this:
> ******************
> *** target.txt ***
> ******************
> SlideNumber     FileNameCy3     FileNameCy5     Species Cy3     Cy5
> HSAP0406        39037_532_pair.txt      39037_635_pair.txt      Homo 
> sapiens (human)    total   HSAP0406
>
> ********************
> *** spottype.txt ***
> ********************
> SpotType        GENE_EXPR_OPTION        PROBE_ID        Color
> Probe   BLOCK1  *       black
>
>
> ********************
> *** sessionInfo()***
> ********************
> > sessionInfo()
> R version 2.6.2 (2008-02-08)
> x86_64-unknown-linux-gnu
>
> locale:
> ...
>
> attached base packages:
> [1] tools     stats     graphics  grDevices utils     datasets
> methods   base
>
> other attached packages:
>  [1] Ringo_1.2.0          RColorBrewer_1.0-2   vsn_3.2.1 
> affy_1.16.0          preprocessCore_1.0.0 affyio_1.6.1         
> limma_2.12.0
>  [8] geneplotter_1.16.0   lattice_0.17-8       annotate_1.16.1 
> xtable_1.5-2         AnnotationDbi_1.0.6  RSQLite_0.6-8
> DBI_0.2-4         [15] Biobase_1.16.3
>
> loaded via a namespace (and not attached):
> [1] grid_2.6.2         KernSmooth_2.22-22 rcompgen_0.1-17
> >



More information about the Bioconductor mailing list