In [1]:

    

E_vals     = 2:10  # range to test of system dimensionality
tau_s_vals = 1:1 # range for lag length for manifold reconstruction
tau_p_vals = 0:15  # range to test for time lag of causal effect

ds   = readdlm("vr.raw_fixed.txt", '\t', Float64, header=false)  # read in data
para = ds[:, 2]
didi = ds[:, 3]


    

    
Out[1]:





71-element Array{Float64,1}:
   5.76
   9.05
  17.26
  41.97
  55.97
  74.91
  62.52
  27.04
  18.77
  31.11
  58.31
  73.13
  63.21
   ⋮   
  55.5 
  93.41
 117.3 
  95.02
  85.92
  82.6 
  66.08
  63.58
  37.99
  25.6 
  23.1 
  37.09

In [4]:

    

para = para[1:2:end]
didi = didi[1:2:end]


    

    
Out[4]:





36-element Array{Float64,1}:
   5.76
  17.26
  55.97
  62.52
  18.77
  58.31
  63.21
  40.07
  26.0 
  21.0 
  64.67
 103.4 
  65.4 
   ⋮   
  17.82
  65.61
  96.07
  54.79
  32.48
  35.73
  93.41
  95.02
  82.6 
  63.58
  25.6 
  37.09

In [7]:

    

using CauseMap
nboots=20
makeoptimizationplots(para, didi, 8, 36, E_vals, 
tau_s_vals, tau_p_vals, 36, 1, 'Para', 'Didi',
            nboots=nboots)


    

    




`makeoptimizationplots` has no method matching makeoptimizationplots(::Array{Float64,1}, ::Array{Float64,1}, ::Int64, ::Int64, ::UnitRange{Int64}, ::UnitRange{Int64}, ::UnitRange{Int64}, ::Int64, ::Int64)
while loading In[7], in expression starting on line 3

In [ ]: