Skip to main content
adicionou 548 caracteres ao conteúdo
Fonte Link

EDIT 2: Neste link está uma imagem com o gráfico de autocorrelacção (cima) e o PACF (baixo). As minhas conclusões são que ele evidência que um modelo AR funciona dado que o valor de autocorrelacção vai baixando lentamente à medida que o valor do lag aumenta. Para além disso, o máximo de lag que se poderá considerar é 6. O que acontece se considerar-mos mais? Simplesmente as previsões no futuro não serão tão certas (poderão conter mais ruído devido à inclusão de mais instantes passados) ?

EDIT 2: Neste link está uma imagem com o gráfico de autocorrelacção (cima) e o PACF (baixo). As minhas conclusões são que ele evidência que um modelo AR funciona dado que o valor de autocorrelacção vai baixando lentamente à medida que o valor do lag aumenta. Para além disso, o máximo de lag que se poderá considerar é 6. O que acontece se considerar-mos mais? Simplesmente as previsões no futuro não serão tão certas (poderão conter mais ruído devido à inclusão de mais instantes passados) ?

adicionou 1106 caracteres ao conteúdo
Fonte Link

-- 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:

  1. Isto está correcto? Pela resposta do Luiz, está.
  2. O y_train pode conter o valor do sinal P passos à frente? Ex: f(t+p)
  3. Podem-me explicar o output do modelo?
  4. A função linear é aproximada com recurso ao método dos mínimos quadrados, certo?
  5. 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

-- 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:

  1. Isto está correcto? Pela resposta do Luiz, está.
  2. O y_train pode conter o valor do sinal P passos à frente? Ex: f(t+p)
  3. Podem-me explicar o output do modelo?
  4. A função linear é aproximada com recurso ao método dos mínimos quadrados, certo?
  5. 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?

-- 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:

  1. Isto está correcto? Pela resposta do Luiz, está.
  2. O y_train pode conter o valor do sinal P passos à frente? Ex: f(t+p)
  3. Podem-me explicar o output do modelo?
  4. A função linear é aproximada com recurso ao método dos mínimos quadrados, certo?
  5. 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
adicionou 69 caracteres ao conteúdo
Fonte Link

-- 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 éDe notar que a minha time series é do tipo univariate, ou seja, só tenho o registo de uma variável ao longo do tipo univariate!tempo.

Questões:

  1. Isto está correcto? Pela resposta do Luiz, está.
  2. O y_train pode conter o valor do sinal P passos à frente? Ex: f(t+p)
  3. Podem-me explicar o output do modelo?
  4. A função linear é aproximada com recurso ao método dos mínimos quadrados, certo?
  5. 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?

-- 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!

Questões:

  1. Isto está correcto? Pela resposta do Luiz, está.
  2. O y_train pode conter o valor do sinal P passos à frente? Ex: f(t+p)
  3. Podem-me explicar o output do modelo?
  4. A função linear é aproximada com recurso ao método dos mínimos quadrados, certo?
  5. 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?

-- 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:

  1. Isto está correcto? Pela resposta do Luiz, está.
  2. O y_train pode conter o valor do sinal P passos à frente? Ex: f(t+p)
  3. Podem-me explicar o output do modelo?
  4. A função linear é aproximada com recurso ao método dos mínimos quadrados, certo?
  5. 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?
+Info
Fonte Link
Carregando
Pergunta reformulada
Fonte Link
Carregando
Fonte Link
Carregando