oposSOM C binding function error
0
0
Entering edit mode
@arrigonialberto86-7608
Last seen 8.7 years ago
Italy

Hello, I am getting an error while trying to execute the example code for the oposSOM package as reported by

?opossom.run

My code is the following:

# Example with artificial data
env <- opossom.new(list(dataset.name="Example",
                        error.model="all.samples",
                        dim.1stLvlSom=20,
                        dim.2ndLvlSom=20,
                        geneset.analysis=TRUE))

env$indata <- matrix(runif(1000), 100, 10)

opossom.run(env)

And this is the output:

[2015-11-23 13:41:45][INFO] Started: Mon 23 Nov 2015 01:41:45 PM
[2015-11-23 13:41:45][INFO] Setting: Example
[2015-11-23 13:41:45][INFO] 1SOM Dim: 20
[2015-11-23 13:41:45][INFO] 2SOM Dim: 20
[2015-11-23 13:41:45][WARN] No rownames found. Set them to 1,2,3,4...
[2015-11-23 13:41:45][WARN] No colnames found. Set them to 1,2,3,4...
[2015-11-23 13:41:46][INFO] Load Annotation Data
Space required after the Public Identifier
SystemLiteral " or ' expected
SYSTEM or PUBLIC, the URI is missing
Opening and ending tag mismatch: hr line 7 and body
Opening and ending tag mismatch: body line 4 and html
Premature end of data in tag html line 2
[2015-11-23 13:41:46][WARN] biomaRt seems to be down.
[2015-11-23 13:41:46][WARN] Disabling geneset analysis.
[2015-11-23 13:41:46][INFO] Processing SOM. This may take several time until next notification.
[2015-11-23 13:41:47][INFO] Remaining ~ 1 min ~ 0 h
[2015-11-23 13:41:48][INFO] Saving environment image: Example pre.RData
[2015-11-23 13:41:48][INFO] Processing Differential Expression
[2015-11-23 13:41:48][INFO] Processing Single Genes
[2015-11-23 13:41:48][INFO] [                                                ][2015-11-23 13:41:48][INFO] [#############################                   ][2015-11-23 13:41:48][INFO] [###############################                 ][2015-11-23 13:41:48][INFO] [#################################               ][2015-11-23 13:41:48][INFO] [###################################             ][2015-11-23 13:41:48][INFO] [####################################            ][2015-11-23 13:41:48][INFO] [######################################          ][2015-11-23 13:41:48][INFO] [########################################        ][2015-11-23 13:41:48][INFO] [##########################################      ][2015-11-23 13:41:48][INFO] [############################################    ][2015-11-23 13:41:48][INFO] [##############################################  ][2015-11-23 13:41:48][INFO] [################################################][2015-11-23 13:41:48][INFO] [################################################]
[2015-11-23 13:41:48][INFO] Writing: Example - Results/LPE/all_samples.bmp
[2015-11-23 13:41:48][INFO] Processing Metagenes
[2015-11-23 13:41:48][INFO] [                                                ][2015-11-23 13:41:48][INFO] [###################                             ][2015-11-23 13:41:48][INFO] [######################                          ][2015-11-23 13:41:48][INFO] [#########################                       ][2015-11-23 13:41:48][INFO] [############################                    ][2015-11-23 13:41:48][INFO] [###############################                 ][2015-11-23 13:41:48][INFO] [##################################              ][2015-11-23 13:41:48][INFO] [####################################            ][2015-11-23 13:41:48][INFO] [#######################################         ][2015-11-23 13:41:48][INFO] [##########################################      ][2015-11-23 13:41:48][INFO] [#############################################   ][2015-11-23 13:41:48][INFO] [################################################][2015-11-23 13:41:48][INFO] [################################################]
[2015-11-23 13:41:48][INFO] Writing: Example - Results/LPE/Sigma_LPE.pdf
[2015-11-23 13:41:48][INFO] Detecting Spots
[2015-11-23 13:41:49][INFO] Plotting Sample Portraits
[2015-11-23 13:41:49][INFO] Writing: Example - Results/Expression Portraits.pdf
[2015-11-23 13:41:49][INFO] Writing: Example - Results/Expression Portraits - alternative scales.pdf
[2015-11-23 13:41:49][INFO] Writing: Example - Results/Rank Maps.pdf
[2015-11-23 13:41:49][INFO] Processing Supporting Information
[2015-11-23 13:41:49][INFO] Writing: Example - Results/Supporting Maps&Profiles/Supporting Maps.pdf
[2015-11-23 13:41:50][INFO] Writing: Example - Results/Supporting Maps&Profiles/Entropy Profiles.pdf
[2015-11-23 13:41:50][INFO] Writing: Example - Results/Supporting Maps&Profiles/Topology Profiles.pdf
[2015-11-23 13:41:50][INFO] Processing 2nd level Metagene Analysis
[2015-11-23 13:41:50][INFO] Writing: Example - Results/2nd lvl Metagene Analysis/Similarity Analysis.pdf
[2015-11-23 13:41:53][INFO] Writing: Example - Results/2nd lvl Metagene Analysis/Correlation Analysis.pdf
Error in .Call("R_igraph_layout_kamada_kawai", graph, coords, maxiter,  : 
  At structural_properties.c:5235 : cannot run Bellman-Ford algorithm, Negative loop detected while calculating shortest paths

Apart from the initial BiomaRt error, which has been reported by others (things have changed recently with the Ensembl Biomart API), the last error prevents the function from producing the other results.

 

>sessionInfo()
R version 3.2.0 (2015-04-16)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: Ubuntu precise (12.04.5 LTS)

locale:
[1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C               LC_TIME=en_US.UTF-8        LC_COLLATE=en_US.UTF-8   
[5] LC_MONETARY=en_US.UTF-8    LC_MESSAGES=en_US.UTF-8    LC_PAPER=en_US.UTF-8       LC_NAME=C                
[9] LC_ADDRESS=C               LC_TELEPHONE=C             LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C       

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
[1] oposSOM_1.4.0

loaded via a namespace (and not attached):
[1] lattice_0.20-33      ape_3.3              IRanges_2.2.9        XML_3.98-1.3         bitops_1.0-6         grid_3.2.0         
[7] GenomeInfoDb_1.4.3   nlme_3.1-122         DBI_0.3.1            stats4_3.2.0         magrittr_1.5         RSQLite_1.0.0      
[13] KernSmooth_2.23-15   som_0.3-5            pixmap_0.4-11        fdrtool_1.2.15       scatterplot3d_0.3-36 S4Vectors_0.6.6    
[19] fastICA_1.2-0        tools_3.2.0          Biobase_2.28.0       biomaRt_2.24.1       RCurl_1.95-4.7       igraph_1.0.1       
[25] parallel_3.2.0       BiocGenerics_0.14.0  AnnotationDbi_1.30.1

oposSOM Bellman-Ford bioconductor • 1.4k views
ADD COMMENT
0
Entering edit mode

Hello,

Did you manage to debug the error? I found the same problem when I tried to use this package on my data.

Thank you.

ADD REPLY
0
Entering edit mode

Is there a solution for this, as I have the same problem!

ADD REPLY

Login before adding your answer.

Traffic: 725 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