-- Pergunta reformulada --
Eu tenho o seguinte código
mdl = LinearModel.fit(X_train,y_train);
pred = predict(mdl,X_test);
em que X_train
é uma matriz [m x n] com m
exemplos e n
features (valores passados da time series). y_train
é um vector coluna com a resposta (valores no instante seguinte). De notar que a minha time series é do tipo univariate, ou seja, só tenho o registo de uma variável ao longo do tempo.
Questões:
- Isto está correcto? Pela resposta do Luiz, está.
- O
y_train
pode conter o valor do sinal P passos à frente? Ex: f(t+p)
- Podem-me explicar o output do modelo?
- A função linear é aproximada com recurso ao método dos mínimos quadrados, certo?
- Existirá alguma razão para que este método dê melhores resultados quando comparado com o algoritmo SVR tanto com o kernel linear como o RBF?
EDIT1:
clc; close all;
clear all;
T = 20;
% Pode inicial aqui fazendo RV = SuaVariavel
RV = 1:T;
RV = RV';
%plot(RV)
%# Construindo as variaveis com 3 defasagens
N = 3;
first = 1;
for i=N:-1:1
if first == 1
X = RV(i:end+(i-1)-N);
first = 0;
else
X = horzcat(X,RV(i:end+(i-1)-N) );
end
end
HP = 1; % <-----------------------------------------------------
X = horzcat(ones(length(X),1),X);
y = RV(N+HP:end,end);
lim = length(X);
X_train = X(round(1:lim*0.7),:);
y_train = y(round(1:lim*0.7));
X_test = X(round(lim*0.7+0.5):end-(HP-1),:);
y_test = y(round(lim*0.7+0.5):end);
%# MQO
mdl = regress(y_train, X_train);
mdl
% Prever
for i=1:length(X_test)
pred(i)=sum(mdl'.*X_test(i,:));
end
figure()
hold on;
plot(pred,'r*-')
plot(y_test,'go-')
hold off