demo/manuscript/deut_analysis.R

In cboettig/warningsignals: Detection of Early Warning Signals for Catastrophic Bifurcations

# deut_analysis.R
rm(list=ls())
require(warningsignals)


load("~/flickr/5904550426.Rdat")

###############
require(socialR)
script <- "deut_analysis.R"
gitaddr <- gitcommit(script)
tags="warningsignals, stochpop"
tweet_errors(script, tags=tags)
###############
source("analysis.R")

data(deuterium)
i <- 3 ## Which deut?


cpu <- 1
nboot <- 2
taus <- 
reformat_tau_dists(
  bootstrap_tau(m$X, m$const, m$timedep, 
                cpu=cpu, nboot=nboot, windowsize=20))