Bonjour,

Je suis débutant sur Matlab. Je voulais notamment tester le script d'un test KPSS pour tester la stationnarité d'une série et mieux comprendre comment coder sur ce sujet là.

J'ai donc téléchargé le script en .m :

Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
function test = kpss(x,w,opt)
T = size(x,1);
if nargin == 2,
   opt = 0;
end
if opt == 0,
   e = center(x);
   else
   [par, e] = ols(x,[ones(T,1) (1:T)']);
end
prod = zeros(w,1);
for j = 1:w
   prod(j) = e(j+1:T)'*e(1:T-j);
end
s2 = e'*e + 2*(1-(1:w)/(w+1))*prod;
S = cumsum(e);
test = T^(-1)*(S'*S)/s2;
Problème comment dois je l'appliquer à une série de données ?

Merci de votre aide
Bon Dimanche