Converts time series data sorted within nominal classes in decreasing chronological order to a useful format for processing by a neural network.
#include <imsls.h>
float
*imsls_f_time_series_class_filter (int
n_patterns,
int n_lags,
int
n_classes,
int i_class[],
float x[], …,0)
The type double function is imsls_d_time_series_class_filter.
int
n_patterns (Input)
Number of observations. The
number of observations must be greater than max_lags.
int n_lags
(Input)
The number of lags. The number of lags must be one or
greater.
int n_classes
(Input)
The number of classes associated with these data. The number of
classes must be one or greater.
int i_class[]
(Input)
An array of length n_patterns. The
i-th element in i_class is equal to
the class associated with the i-th element of x. The classes must be
numbered from 1 to n_classes.
float x[]
(Input)
A sorted array of length n_patterns. This
array is assumed to be sorted first by class designations and then descending by
chronological order, i.e., most recent observations appear first within a
class.
A pointer to an internally allocated array of size n_patterns by n_lags columns. If errors are encountered, then NULL is returned.
#include <imsls.h>
float
*imsls_f_time_series_class_filter (int
n_patterns,
int
n_lags,
int
n_classes, int
i_class[],
float
x[],
IMSLS_LAGS, int
lag[],
IMSLS_RETURN_USER, float
z[],
0)
The type double function is imsls_d_time_series_class_filter.
IMSLS_LAGS,
int
lag[] (Input)
An array of length n_lags. The
i-th element in lag is equal to the
lag requested for the i-th column of z. Every lag
must be non-negative.
Default: lag[i]=i
IMSLS_RETURN_USER, float z[]
(Output)
A user-supplied array of size n_patterns by n_lags. The
i-th column contains the lagged values of x for a lag equal to
the number of lags in lag[i].
The function imsls_f_time_series_class_filter accepts a data array, x[], and returns a new data array, z[], containing n_lags columns, each containing a lagged version of x.
The output data array, z, can be represented symbolically as:
z = |x(0) : x(1) : x(2) : … : x(n_lags-1)|,
where x(i) is the i-th lagged column of the incoming data array, x. Notice that n_lags is the number of lags and not the maximum lag. The maximum number of lags is max_lag= n_lags-1, unless the optional input lag[] is given, the highest lag is max_lags. If n_lags =2 and the optional input lag[] is not given, then the output array contains the lags 0, 1.
Consider, an example in which n_patterns=10, n_lags =2 and
.
If and
.
then, n_classes=1 and z would contain 2 columns and 10 rows:
.
Note that since lagT = [0,1], the first column of z is formed using a lag of zero and the second is formed using a lag of two. A zero lag corresponds to no lag, which is why the first column of z in this example is equal to the original data in x.
On the other hand, if the data were organized into two classes with
,
then z is still a 2 by 10 matrix, but with the following values:
The first 5 rows of z are the lagged columns for the first class, and the last five are the lagged columns for the second class.
Suppose that the training data to the neural network is represented by the following data matrix consisting of a single nominal variable coded into two binary columns and a single time series variable:
In this case, n_patterns=8 and n_classes=2. If we wanted to lag the 3rd column by 2 time lags, i.e., n_lags=2,
,
, and
.
The resulting data matrix would have 8 rows and 2 columns:
.
int main(){
#define N_PATTERNS 8
#define N_LAGS 2
float x[N_PATTERNS] = {2.1, 2.3, 2.4, 2.5, 1.1, 1.2, 1.3, 1.4};
float *z;
int n_classes = 2;
int i_class[] = {1,1,1,1,2,2,2,2};
z = imsls_f_time_series_class_filter(N_PATTERNS, N_LAGS, n_classes,
i_class, x,
0);
imsls_f_write_matrix("z", N_PATTERNS, N_LAGS, (float*)z, 0);
}
z
1 2
1 2.1 2.3
2 2.3 2.4
3 2.4 2.5
4 2.5 ...........
5 1.1 1.2
6 1.2 1.3
7 1.3 1.4
8 1.4 ...........