HilbertVis troubles
1
0
Entering edit mode
Tobias Straub ▴ 420
@tobias-straub-2182
Last seen 6.8 years ago
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
HilbertVis HilbertVis • 702 views
ADD COMMENT
0
Entering edit mode
Simon Anders ▴ 150
@simon-anders-2626
Last seen 6.8 years ago
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
ADD COMMENT

Login before adding your answer.

Traffic: 273 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6