Note that there are some explanatory texts on larger screens.

plurals
  1. POR package VLMC dies if state space size exceeds 27
    primarykey
    data
    text
    <p>I am using VLMC to fit some Markov models and it dies as soon as the alphabet size reaches 28. I thought this was due to using a single letter in the alphabet by default, but it has the same behavior with "code1char = FALSE". This is true for me on real data as well as this fake example.</p> <pre><code>library(VLMC) # works fine ins &lt;- sample(seq(1,27,1),50000,replace=T) vlmc(ins, dump = 1,threshold.gen = 2, debug = TRUE) #core dump ins &lt;- sample(seq(1,28,1),50000,replace=T) vlmc(ins, dump = 1,threshold.gen = 2, debug = TRUE) </code></pre> <p>Any ideas?</p> <p>The seg fault looks like this BTW. It looks to me like the alphabet after z is being mapped to NA which is causing an array bound issue.</p> <pre><code>library(VLMC) sc &lt;- 10 amp &lt;- 13 x &lt;- round(amp*sin(seq(0,2*sc*pi,0.01))) x &lt;- amp + x + rpois(NROW(x),1) length(table(x)) length(x) vlmc(x, dump = 1,threshold.gen = 2, debug = TRUE) vlmc: Alpha = 'abcdefghijklmnopqrstuvwxyzNANANANANA' ; |X| = 31 vlmc: ctl.dump = 4 11 vlmc: n = |data| = 6284, cutoff{prune} = 21.8865, threshold{gen} = 2 vlmc: |alphabet| = 31, alphabet = abcdefghijklmnopqrstuvwxyzNA generating... *** caught segfault *** address 0x0, cause 'memory not mapped' Traceback: 1: .C("vlmc_p", data = Data, n = n, threshold.gen = as.integer(threshold.gen), cutoff.prune = as.double(cutoff.prune), alpha.len = as.integer(alpha.len), alpha = as.character(Alpha), debug = as.integer(as.logical(debug)), dump.flags = as.integer(c(dump, ctl.dump)), size = integer(4), PACKAGE = "VLMC") 2: vlmc(x, dump = 1, threshold.gen = 2, debug = TRUE) </code></pre>
    singulars
    1. This table or related slice is empty.
    plurals
    1. This table or related slice is empty.
    1. This table or related slice is empty.
    1. This table or related slice is empty.
 

Querying!

 
Guidance

SQuiL has stopped working due to an internal error.

If you are curious you may find further information in the browser console, which is accessible through the devtools (F12).

Reload