hi i want to visualize my tiling array data with HilbertVis and failed using the hilbertImage function. I get a proper plot using plotLongVector but when i call the hilbertImage function i get "Error: length(vec) >= newLength is not TRUE". all operations on random data work. so what is wrong with mine? > summary(data) Min. 1st Qu. Median Mean 3rd Qu. Max. -1.48800 -0.18810 0.08333 0.48070 1.00400 3.80400 > mode(data) [1] "numeric" > plotLongVector(data) > hMat <- hilbertImage(data) Error: length(vec) >= newLength is not TRUE > sessionInfo() R version 2.9.0 (2009-04-17) x86_64-apple-darwin9.6.0 locale: en_US.UTF-8/en_US.UTF-8/C/C/en_US.UTF-8/en_US.UTF-8 attached base packages: [1] splines grid stats graphics grDevices utils datasets methods base other attached packages: [1] simpleTile_0.1.39 RColorBrewer_1.0-2 tileHMM_1.0-2 preprocessCore_1.6.0 [5] affy_1.22.0 locfdr_1.1-6 st_1.1.3 sda_1.1.0 [9] fdrtool_1.2.5 corpcor_1.5.2 entropy_1.1.3 vsn_3.12.0 [13] geneplotter_1.22.0 annotate_1.22.0 AnnotationDbi_1.6.0 Biobase_2.4.1 [17] gplots_2.7.1 caTools_1.9 bitops_1.0-4.1 gdata_2.4.2 [21] gtools_2.6.1 HilbertVis_1.2.0 lattice_0.17-22 loaded via a namespace (and not attached): [1] affyio_1.12.0 DBI_0.2-4 limma_2.18.0 RSQLite_0.7-1 tools_2.9.0 xtable_1.5-5 > ---------------------------------------------------------------------- Tobias Straub ++4989218075439 Adolf-Butenandt-Institute, M?nchen D

Hi Tobia, Tobias Straub wrote: > i want to visualize my tiling array data with HilbertVis and failed > using the hilbertImage function. > I get a proper plot using plotLongVector but when i call the > hilbertImage function i get "Error: length(vec) >= newLength is not TRUE". > all operations on random data work. so what is wrong with mine? Most liekly, your vector is too short. 'hilbertImage' requires each pixel in the image to correspond to at least one data element. When you don't specify a level, the default 9 is used, i.e., a picture with 2^9 x 2^9 pixels is produced, requiring a vector of minimum length 4^9=262,144. If your vector is shorter, just call hilbertImage with a suitable level, e.g., img = hilbertImage( data, level=7 ) 'hilbertDisplay' does not have this limitation, by the way, and is more convenient to use anyway, unless you want to batch-process something. (Of course, the limitation should be unnecessary, and I'll try to remove it in the next version.) Best regards Simon +--- | Dr. Simon Anders, Dipl. Phys. | European Bioinformatics Institute, Hinxton, Cambridgeshire, UK | office phone +44-1223-492680, mobile phone +44-7505-841692 | preferred (permanent) e-mail: sanders at fs.tum.de