-
Notifications
You must be signed in to change notification settings - Fork 1
/
esnet.py
642 lines (519 loc) · 25.1 KB
/
esnet.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
import numpy as np
from scipy import sparse
from scipy.io import loadmat
from scipy.sparse import linalg as slinalg
from sklearn.decomposition import KernelPCA, PCA, SparsePCA, IncrementalPCA, TruncatedSVD, FactorAnalysis
from sklearn.manifold import MDS
from sklearn.linear_model import ElasticNet, Lasso, Ridge, LinearRegression, BayesianRidge
from sklearn.gaussian_process import GaussianProcessRegressor
from sklearn.metrics import mean_squared_error
from sklearn.preprocessing import StandardScaler, RobustScaler, MinMaxScaler
from sklearn.svm import NuSVR, LinearSVR
from sklearn.model_selection import GridSearchCV
#from skbayes.rvm_ard_models import RegressionARD,ClassificationARD,RVR,RVC,vrvm, VBRegressionARD
from skbayes.linear_models import VBLinearRegression, EBLinearRegression
from sklearn import linear_model
#from wpca import WPCA, EMPCA
import RobustPCA
import FastICA
import KICA
import math
import subprocess
import tga
import warnings
import sys
if not sys.warnoptions:
warnings.simplefilter('ignore')
def NRMSE(y_true, y_pred, scaler):
y_true = scaler.inverse_transform(y_true)
y_pred = scaler.inverse_transform(y_pred)
#Normalized Root Mean Squared Error
y_std = np.std(y_true)
#return mean_squared_error(y_true, y_pred)
return np.sqrt(mean_squared_error(y_true, y_pred))/y_std
class ESN(object):
def __init__(self, n_internal_units = 100, spectral_radius = 0.9, connectivity = 0.5, input_scaling = 0.5, input_shift = 0.0,
teacher_scaling = 0.5, teacher_shift = 0.0, noise_level = 0.01):
# Initialize attributes
self._n_internal_units = n_internal_units
self._spectral_radius = spectral_radius
self._connectivity = connectivity
self._input_scaling = input_scaling
self._input_shift = input_shift
self._teacher_scaling = teacher_scaling
self._teacher_shift = teacher_shift
self._noise_level = noise_level
self._dim_output = None
# The weights will be set later, when data is provided
self._input_weights = None
# Regression method and embedding method.
# Initialized to None for now. Will be set during 'fit'.
self._regression_method = None
self._embedding_method = None
# Generate internal weights
self._internal_weights = self._initialize_internal_weights(n_internal_units, connectivity, spectral_radius)
def fit(self, Xtr, Ytr, n_drop = 100, regression_method = 'linear', regression_parameters = None, embedding = 'identity', n_dim = 3, embedding_parameters = None):
_,_ = self._fit_transform(Xtr = Xtr, Ytr = Ytr, n_drop = n_drop, regression_method = regression_method, regression_parameters = regression_parameters, embedding = embedding, n_dim = n_dim, embedding_parameters = embedding_parameters)
return
def _fit_transform(self, Xtr, Ytr, n_drop = 100, regression_method = 'linear', regression_parameters = None, embedding = 'identity', n_dim = 3, embedding_parameters = None):
n_data, dim_data = Xtr.shape
_, dim_output = Ytr.shape
self._dim_output = dim_output
# If this is the first time the network is tuned, set the input weight.
# The weights are dense and uniformly distributed in [-1.0, 1.0]
if (self._input_weights is None):
self._input_weights = 2.0*np.random.rand(self._n_internal_units, dim_data) - 1.0
# Initialize regression method
if (regression_method == 'nusvr'):
# NuSVR, RBF kernel
C, nu, gamma = regression_parameters
self._regression_method = NuSVR(C = C, nu = nu, gamma = gamma)
elif (regression_method == 'linsvr'):
# NuSVR, linear kernel
#C = regression_parameters[0]
#nu = regression_parameters[1]
C, epsilon = regression_parameters
#self._regression_method = NuSVR(C = C, nu = nu, kernel='linear')
self._regression_method = LinearSVR(C = C, epsilon = epsilon)
elif (regression_method == 'enet'):
# Elastic net
alpha, l1_ratio = regression_parameters
self._regression_method = ElasticNet(alpha = alpha, l1_ratio = l1_ratio)
elif (regression_method == 'ridge'):
# Ridge regression
self._regression_method = Ridge(alpha = regression_parameters)
elif (regression_method == 'lasso'):
# LASSO
self._regression_method = Lasso(alpha = regression_parameters)
elif (regression_method == 'bayeridge'):
lambda_1, lambda_2, alpha_1, alpha_2 = regression_parameters
self._regression_method = BayesianRidge(lambda_1=lambda_1,lambda_2=lambda_2,alpha_1=alpha_1,alpha_2=alpha_2)
elif (regression_method == 'gpr'):
self._regression_method = GaussianProcessRegressor()
elif (regression_method == 'bayelinear'):
self._regression_method = EBLinearRegression()
else:
# Use canonical linear regression
self._regression_method = LinearRegression()
# Initialize embedding method
if (embedding == 'identity'):
self._embedding_dimensions = self._n_internal_units
else:
self._embedding_dimensions = n_dim
if (embedding == 'kpca'):
# Kernel PCA with RBF kernel
self._embedding_method = KernelPCA(n_components = n_dim, kernel = 'rbf', gamma = embedding_parameters)
elif (embedding == 'pca'):
# PCA
self._embedding_method = PCA(n_components = n_dim)
elif (embedding == 'fa'):
# ICA
self._embedding_method = FactorAnalysis(n_components = n_dim)
elif (embedding == 'spca'):
# Sparse PCA
self._embedding_method = SparsePCA(n_components = n_dim, alpha = embedding_parameters)
elif (embedding == 'ipca'):
# Sparse PCA
self._embedding_method = IncrementalPCA(n_components = n_dim)
elif (embedding == 'tsvd'):
# Sparse PCA
if n_dim >= self._n_internal_units:
self._embedding_method = TruncatedSVD(n_components = self._n_internal_units-1)
else:
self._embedding_method = TruncatedSVD(n_components = n_dim)
elif (embedding == 'wpca'):
# Bayesian Probabilistic PCA
self._embedding_method = WPCA(n_components=n_dim)
elif (embedding == 'rpca'):
# Bayesian Probabilistic PCA
self._embedding_method = RobustPCA.RobustPCA()
elif (embedding == 'tga'):
# Bayesian Probabilistic PCA
self._embedding_method = tga.TGA(n_components=n_dim, random_state=1)
elif (embedding == 'empca'):
# Expectation Maximization PCA
self._embedding_method = EMPCA(n_components=n_dim)
elif (embedding == 'mds'):
# Multi-Dimensional Scaling (MDS)
self._embedding_method = MDS(n_components=n_dim)
elif (embedding == 'ica'):
# Sparse PCA
alpha = embedding_parameters
self._embedding_method = FastICA.FastICA(n_components=n_dim)
#self._embedding_method = FastICA.FastICA(n_components=n_dim, fun_args={'alpha':alpha})
#self._embedding_method = FastICA.FastICA(n_components = n_dim, algorithm = 'deflation')
elif (embedding == 'kica'):
self._embedding_method = KICA.KICA(n_components=n_dim)
else:
raise(ValueError, "Unknown embedding method")
# Calculate states/embedded states.
# Note: If the embedding is 'identity', embedded states will be equal to the states.
states, embedded_states,_ = self._compute_state_matrix(X = Xtr, Y = Ytr, n_drop = n_drop)
# Train output
self._regression_method.fit(np.concatenate((embedded_states, self._scaleshift(Xtr[n_drop:, :], self._input_scaling, self._input_shift)), axis=1),
self._scaleshift(Ytr[n_drop:, :], self._teacher_scaling,self._teacher_shift).flatten())
return states, embedded_states
def predict(self, X, Y = None, n_drop = 100, error_function = NRMSE, scaler = None):
Yhat, error, _, _ = self._predict_transform(X = X, Y = Y, n_drop = n_drop, error_function = error_function, scaler = scaler)
return Yhat, error
def _predict_transform(self, X, Y = None, n_drop = 100, error_function = NRMSE, scaler = None):
# Predict outputs
states,embedded_states,Yhat = self._compute_state_matrix(X = X, n_drop = n_drop)
# Revert scale and shift
Yhat = self._uscaleshift(Yhat, self._teacher_scaling, self._teacher_shift)
# Compute error if ground truth is provided
if (Y is not None):
error = error_function(Y[n_drop:,:], Yhat, scaler)
return Yhat, error, states, embedded_states
def _compute_state_matrix(self, X, Y = None, n_drop = 100):
n_data, _ = X.shape
# Initial values
previous_state = np.zeros((1, self._n_internal_units), dtype=float)
previous_output = np.zeros((1, self._dim_output), dtype=float)
# Storage
state_matrix = np.empty((n_data - n_drop, self._n_internal_units), dtype=float)
embedded_states = np.empty((n_data - n_drop, self._embedding_dimensions), dtype=float)
outputs = np.empty((n_data - n_drop, self._dim_output), dtype=float)
for i in range(n_data):
# Process inputs
previous_state = np.atleast_2d(previous_state)
current_input = np.atleast_2d(self._scaleshift(X[i, :], self._input_scaling, self._input_shift))
# Calculate state. Add noise and apply nonlinearity.
state_before_tanh = self._internal_weights.dot(previous_state.T) + self._input_weights.dot(
current_input.T)
state_before_tanh += np.random.rand(self._n_internal_units, 1) * self._noise_level
previous_state = np.tanh(state_before_tanh).T
# Embed data and perform regression if applicable.
if (Y is not None):
# If we are training, the previous output should be a scaled and shifted version of the ground truth.
previous_output = self._scaleshift(Y[i, :], self._teacher_scaling, self._teacher_shift)
else:
# Should the data be embedded?
if (self._embedding_method is not None):
current_embedding = self._embedding_method.transform(previous_state)
else:
current_embedding = previous_state
# Perform regression
previous_output = self._regression_method.predict(
np.concatenate((current_embedding, current_input), axis=1))
# Store everything after the dropout period
if (i > n_drop - 1):
state_matrix[i - n_drop, :] = previous_state.flatten()
# Only save embedding for test data.
# In training, we do it after computing the whole state matrix.
if (Y is None):
embedded_states[i - n_drop, :] = current_embedding.flatten()
outputs[i - n_drop, :] = previous_output.flatten()
# Now, embed the data if we are in training
if (Y is not None):
if (self._embedding_method is not None):
embedded_states = self._embedding_method.fit_transform(state_matrix)
else:
embedded_states = state_matrix
return state_matrix, embedded_states, outputs
def _scaleshift(self, x, scale, shift):
# Scales and shifts x by scale and shift
return (x*scale + shift)
def _uscaleshift(self, x, scale, shift):
# Reverts the scale and shift applied by _scaleshift
return ( (x - shift)/float(scale) )
def _initialize_internal_weights(self, n_internal_units, connectivity, spectral_radius):
# The eigs function might not converge. Attempt until it does.
convergence = False
while (not convergence):
# Generate sparse, uniformly distributed weights.
internal_weights = sparse.rand(n_internal_units, n_internal_units, density=connectivity).todense()
# Ensure that the nonzero values are uniformly distributed in [-0.5, 0.5]
internal_weights[np.where(internal_weights > 0)] -= 0.5
try:
# Get the largest eigenvalue
w,_ = slinalg.eigs(internal_weights, k=1, which='LM')
convergence = True
except:
continue
# Adjust the spectral radius.
internal_weights /= np.abs(w)/spectral_radius
return internal_weights
def sigmoid(x):
return 1 / (1 + np.exp(-x))
def run_from_config(Xtr, Ytr, Xte, Yte, config, scaler):
# Instantiate ESN object
esn = ESN(n_internal_units = config['n_internal_units'],
spectral_radius = config['spectral_radius'],
connectivity = config['connectivity'],
input_scaling = config['input_scaling'],
input_shift = config['input_shift'],
teacher_scaling = config['teacher_scaling'],
teacher_shift = config['teacher_shift'],
noise_level = config['noise_level'])
# Get parameters
n_drop = config['n_drop']
regression_method = config['regression_method']
regression_parameters = config['regression_parameters']
embedding = config['embedding']
n_dim = config['n_dim']
embedding_parameters = config['embedding_parameters']
# Fit and predict
esn.fit(Xtr, Ytr, n_drop = n_drop, regression_method = regression_method, regression_parameters = regression_parameters,
embedding = embedding, n_dim = n_dim, embedding_parameters = embedding_parameters)
Yhat,error = esn.predict(Xte, Yte, scaler = scaler)
return Yhat, error
def run_from_config_return_states(Xtr, Ytr, Xte, Yte, config, scaler):
# Instantiate ESN object
esn = ESN(n_internal_units = config['n_internal_units'],
spectral_radius = config['spectral_radius'],
connectivity = config['connectivity'],
input_scaling = config['input_scaling'],
input_shift = config['input_shift'],
teacher_scaling = config['teacher_scaling'],
teacher_shift = config['teacher_shift'],
noise_level = config['noise_level'])
# Get parameters
n_drop = config['n_drop']
regression_method = config['regression_method']
regression_parameters = config['regression_parameters']
embedding = config['embedding']
n_dim = config['n_dim']
embedding_parameters = config['embedding_parameters']
# Fit and predict
train_states, train_embedding = esn._fit_transform(Xtr, Ytr, n_drop = n_drop, regression_method = regression_method, regression_parameters = regression_parameters,
embedding = embedding, n_dim = n_dim, embedding_parameters = embedding_parameters)
Yhat, error, test_states, test_embedding = esn._predict_transform(Xte, Yte, scaler = scaler)
return Yhat, error, train_states, train_embedding, test_states, test_embedding
def format_config(n_internal_units, spectral_radius, connectivity, input_scaling, input_shift, teacher_scaling, teacher_shift, noise_level,
n_drop, regression_method, regression_parameters, embedding, n_dim, embedding_parameters):
config = dict(
n_internal_units = n_internal_units,
spectral_radius = spectral_radius,
connectivity = connectivity,
input_scaling = input_scaling,
input_shift = input_shift,
teacher_scaling = teacher_scaling,
teacher_shift = teacher_shift,
noise_level = noise_level,
n_drop = n_drop,
regression_method = regression_method,
regression_parameters = regression_parameters,
embedding = embedding,
n_dim = n_dim,
embedding_parameters = embedding_parameters
)
return config
def generate_datasets(X, Y, test_percent = 0.25, val_percent = 0.25, scaler = StandardScaler):
n_data,_ = X.shape
n_te = np.ceil(test_percent*n_data).astype(int)
n_val = np.ceil(val_percent*n_data).astype(int)
#n_te = 491
#n_val = 400
n_tr = n_data - n_te - n_val
# Split dataset
Xtr = X[:n_tr, :]
Ytr = Y[:n_tr, :]
Xval = X[n_tr:-n_te, :]
Yval = Y[n_tr:-n_te, :]
Xte = X[-n_te:, :]
Yte = Y[-n_te:, :]
# Scale
Xscaler = scaler()
Yscaler = scaler()
# Fit scaler on training set
Xtr = Xscaler.fit_transform(Xtr)
Ytr = Yscaler.fit_transform(Ytr)
# Transform the rest
Xval = Xscaler.transform(Xval)
Yval = Yscaler.transform(Yval)
Xte = Xscaler.transform(Xte)
Yte = Yscaler.transform(Yte)
return Xtr, Ytr, Xval, Yval, Xte, Yte, Yscaler
def generate_datasets_1d(path, test_percent = 0.25, val_percent = 0.25, scaler = StandardScaler):
data = np.loadtxt(path, delimiter=',')
data = data.reshape((data.shape[0], 1))
n_data,_ = data.shape
n_te = np.ceil(test_percent * n_data).astype(int)
n_val = np.ceil(val_percent * n_data).astype(int)
n_tr = n_data - n_te - n_val
data_train = data[:n_tr, :]
data_val = data[n_tr:-n_te, :]
data_test = data[-n_te:, :]
Xtr = data_train[:-1, :]
Ytr = data_train[1:, :]
Xval = data_val[:-1, :]
Yval = data_val[1:, :]
Xte = data_test[:-1, :]
Yte = data_test[1:, :]
# Scale
Xscaler = scaler()
Yscaler = scaler()
# Fit scaler on training set
Xtr = Xscaler.fit_transform(Xtr)
Ytr = Yscaler.fit_transform(Ytr)
# Transform the rest
Xval = Xscaler.transform(Xval)
Yval = Yscaler.transform(Yval)
Xte = Xscaler.transform(Xte)
Yte = Yscaler.transform(Yte)
return Xtr, Ytr, Xval, Yval, Xte, Yte, Yscaler
def construct_output(X, shift):
return X[:-shift,:], X[shift:, :]
def load_from_text(path):
data = np.loadtxt(path, delimiter=',')
return np.atleast_2d(data[:, :-1]), np.atleast_2d(data[:, -1]).T
def reconstruct_input_2d(arrays,reconstructconfig):
reconstructDim = [reconstructconfig['reconstruct_dim_x'],reconstructconfig['reconstruct_dim_y']]
reconstructDelay = [reconstructconfig['reconstruct_delay_x'], reconstructconfig['reconstruct_delay_y']]
startIndex = 0
for i in range(len(reconstructDim)):
if (reconstructDim[i]-1)*reconstructDelay[i]>startIndex:
startIndex = (reconstructDim[i]-1)*reconstructDelay[i]
returnVals = []
for array in arrays:
reconstructed = None
dataDim = array.shape
for i in range(startIndex, dataDim[0]):
#curIndex = i - startIndex
construct = None
for j in range(dataDim[1]):
subSeries = array[range(i, i-reconstructDelay[j]*(reconstructDim[j]-1)-1,
-reconstructDelay[j]),j]
if construct is None:
construct = subSeries
else:
construct = np.concatenate([construct, subSeries])
if reconstructed is None:
reconstructed = construct
else:
reconstructed = np.vstack([reconstructed, construct])
returnVals.append(reconstructed)
return returnVals
def reconstruct_output_2d(arrays,reconstructconfig):
reconstructDim = [reconstructconfig['reconstruct_dim_x'], reconstructconfig['reconstruct_dim_y']]
reconstructDelay = [reconstructconfig['reconstruct_delay_x'], reconstructconfig['reconstruct_delay_y']]
startIndex = 0
for i in range(len(reconstructDim)):
if (reconstructDim[i] - 1) * reconstructDelay[i] > startIndex:
startIndex = (reconstructDim[i] - 1) * reconstructDelay[i]
returnVals = []
for array in arrays:
dataDim = array.shape
reconstructed = None
for i in range(startIndex, dataDim[0]):
subSeries = array[i, dataDim[1] - 1]
if reconstructed is None:
reconstructed = subSeries
else:
reconstructed = np.vstack([reconstructed, subSeries])
returnVals.append(reconstructed)
return returnVals
def reconstruct_input_3d(arrays,reconstructconfig):
reconstructDim = [reconstructconfig['reconstruct_dim_x'], reconstructconfig['reconstruct_dim_y'],
reconstructconfig['reconstruct_dim_z']]
reconstructDelay = [reconstructconfig['reconstruct_delay_x'], reconstructconfig['reconstruct_delay_y'],
reconstructconfig['reconstruct_delay_z']]
startIndex = 0
for i in range(len(reconstructDim)):
if (reconstructDim[i] - 1) * reconstructDelay[i] > startIndex:
startIndex = (reconstructDim[i] - 1) * reconstructDelay[i]
returnVals = []
for array in arrays:
reconstructed = None
dataDim = array.shape
for i in range(startIndex, dataDim[0]):
# curIndex = i - startIndex
construct = None
for j in range(dataDim[1]):
subSeries = array[range(i, i - reconstructDelay[j] * (reconstructDim[j] - 1) - 1,
-reconstructDelay[j]), j]
if construct is None:
construct = subSeries
else:
construct = np.concatenate([construct, subSeries])
if reconstructed is None:
reconstructed = construct
else:
reconstructed = np.vstack([reconstructed, construct])
returnVals.append(reconstructed)
return returnVals
def reconstruct_output_3d(arrays,reconstructconfig):
reconstructDim = [reconstructconfig['reconstruct_dim_x'], reconstructconfig['reconstruct_dim_y'],
reconstructconfig['reconstruct_dim_z']]
reconstructDelay = [reconstructconfig['reconstruct_delay_x'], reconstructconfig['reconstruct_delay_y'],
reconstructconfig['reconstruct_delay_z']]
startIndex = 0
for i in range(len(reconstructDim)):
if (reconstructDim[i] - 1) * reconstructDelay[i] > startIndex:
startIndex = (reconstructDim[i] - 1) * reconstructDelay[i]
returnVals = []
for array in arrays:
dataDim = array.shape
reconstructed = None
for i in range(startIndex, dataDim[0]):
subSeries = array[i, dataDim[1] - 1]
if reconstructed is None:
reconstructed = subSeries
else:
reconstructed = np.vstack([reconstructed, subSeries])
returnVals.append(reconstructed)
return returnVals
def reconstruct_input_1d(arrays,reconstructconfig):
reconstructDim = reconstructconfig['reconstruct_dim_x']
reconstructDelay = reconstructconfig['reconstruct_delay_x']
startIndex = 0
for i in range(reconstructDim):
if (reconstructDim - 1) * reconstructDelay > startIndex:
startIndex = (reconstructDim - 1) * reconstructDelay
returnVals = []
for array in arrays:
reconstructed = None
dataDim = array.shape
for i in range(startIndex, dataDim[0]):
subSeries = array[range(i, i - reconstructDelay * (reconstructDim - 1) - 1,
-reconstructDelay), 0]
if reconstructed is None:
reconstructed = subSeries
else:
reconstructed = np.vstack([reconstructed, subSeries])
returnVals.append(reconstructed)
return returnVals
def reconstruct_output_1d(arrays,reconstructconfig):
reconstructDim = [reconstructconfig['reconstruct_dim_x']]
reconstructDelay = [reconstructconfig['reconstruct_delay_x']]
startIndex = 0
for i in range(len(reconstructDim)):
if (reconstructDim[i] - 1) * reconstructDelay[i] > startIndex:
startIndex = (reconstructDim[i] - 1) * reconstructDelay[i]
returnVals = []
for array in arrays:
dataDim = array.shape
reconstructed = None
for i in range(startIndex, dataDim[0]):
subSeries = array[i, dataDim[1] - 1]
if reconstructed is None:
reconstructed = subSeries
else:
reconstructed = np.vstack([reconstructed, subSeries])
returnVals.append(reconstructed)
return returnVals
def load_from_dir(path):
Xtr_base = np.loadtxt(path + '/Xtr')
Ytr_base = np.loadtxt(path + '/Ytr')
Xval_base = np.loadtxt(path + '/Xval')
Yval_base = np.loadtxt(path + '/Yval')
Xte_base = np.loadtxt(path + '/Xte')
Yte_base = np.loadtxt(path + '/Yte')
Xtr, Ytr, Xval, Yval, Xte, Yte = np.atleast_2d(Xtr_base, Ytr_base, Xval_base, Yval_base, Xte_base, Yte_base)
# Need axis 0 to be the samples
if Xtr.shape[0] == 1:
Xtr = Xtr.T
if Ytr.shape[0] == 1:
Ytr = Ytr.T
if Xval.shape[0] == 1:
Xval = Xval.T
if Yval.shape[0] == 1:
Yval = Yval.T
if Xte.shape[0] == 1:
Xte = Xte.T
if Yte.shape[0] == 1:
Yte = Yte.T
return Xtr, Ytr, Xval, Yval, Xte, Yte
if __name__ == "__main__":
pass