Nassou1983
Messages postés1Date d'inscriptionmardi 18 avril 2017StatutMembreDernière intervention20 avril 2017
-
20 avril 2017 à 20:46
bonjour,
qui pourrait me corriger ce code
rate = panel(:,78)
for h=[1 2 3 6 12]
for lagy=1:3
for lagx=1:2
lag1 = lagmatrix(Y_AR,1:lagy);
lag2=lagmatrix(rate(h:end),1:lagx)
X =[ones(size(Y_AR,1)-lagy,1) lag1(lagy+1:end,:) lag2(lagy+1:end,:)];
[B,~,R] = regress(Y_AR(lagy+1:end),X);% on recupere les résidus pour calculer le BIC
BIC_ADL(lagy,lagx)=log(sum(R.^2)/(size(X,1)-lagy-1))+log(size(X,1))/size(X,1)*(lagy+1)
end
end
[R,P] = find(BIC_ADL == min(min(BIC_ADL)))
lag_rp_ADL(h,1) = R
lag_rp_ADL(h,2) = P
% POUR FAIRE LA PRÉVISION, ON DÉFINIT LES VARIABLES ET ON ESTIME
% LE MODELE AVEC LES NOUVEAUX P ET R .
Y_ADL = Y_AR ;
lag1 = lagmatrix(Y_AR,1:R);
lag2=lagmatrix(rate(h:end),1:P)
X =[ones(size(Y_AR,1)-R,1) lag1(R+1:end,:) lag2(R+1:end,:)]
mdl=fitlm(X, Y_ADL(R+1:end))
Coefficient.ADL.estimat{h,:}=mdl.Coefficients(:,1);
beta=table2array(Coefficient.ADL.estimat{h,:})
%beta(2,:)=[];
Yhat_adl(h,1) = X(end,:)*beta
% Significativité des coefficients et R carré ajusté
Coefficient.ADL.p_value{h,:}=mdl.Coefficients(:,4);
Rsquared_adj.ADL(h,1)=mdl.Rsquared.Adjusted(1,1);
end
end