A package to estimate local statistical complexity (LSC), a measure for automated pattern discovery in spatiotemporal data using optimal predictors (see References).
This package is very tightly linked to the LICORS
package, which can be used to estimate these optimal
predictors and state space from data. The LSC
builds on a known or estimated state space; most
estimation is handled by LICORS (see
?LICORS
).
There are two ways the state space can be represented: either as a unique state label or as a vector of weights. These two are the principal arguments in the functions of this package:
weight.matrix
an N \times K matrix, where N are the samples and K are the states. That is, each row contains a vector of length K that adds up to one (the mixture weights).
states
a vector of length N with entry i being the label k = 1, …, K of PLC i
This is an early release: some function names and arguments might/will (slightly) change in the future, so regularly check with new package updates.
Georg M. Goerg gmg@stat.cmu.edu
Shalizi, C. R., R. Haslinger, J.B. Rouquier, K. L. Klinkner, and C. Moore (2006). “Automatic filters for the detection of coherent structure in spatiotemporal systems.” Physical Review E 73, 036104
Shalizi, C. R., K. L. Klinkner, and R. Haslinger (2004a). “Quantifying selforganization with optimal predictors.” Physical Review Letters 93, 118701.
The main functions in this package are
states2LSC
to estimate LSC from the state
space, and
LICORS2LSC
which is a
wrapper for estimating LSC from a "LICORS"
class
estimate.
Since pattern discovery without visualization
is only of very limited use, the plot.LSC
function shows informative plots for (1+1)D and
(2+1)D systems.
1 2 3 4 5 6 7 8 9 10 11 12  ## known predictive state space with a statevector
data(contCA00)
ll < states2LSC(states = contCA00$predictive_states  min(contCA00$predictive_states) +
1)
image2(ll, density = TRUE, legend = FALSE)
# An example using estimates from LICORS
## Not run:
example(LICORS) # this will give an object 'mod' of class 'LICORS'
image2(LICORS2LSC(mod))
## End(Not run)

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.
All documentation is copyright its authors; we didn't write any of that.