How to coding and plotting

Hi i have this Page hinkley test

sumPH=0
  for t=1:Nt
       mew2=(inv(Nt)*sum((y2-xhat2)));
       vari2=(inv(Nt-1)*sum((xhat2-mew2).^(2)));
       sumPH= sumPH + vv/vari2*(rg12 - mean(mew2) -vv/2);

       pp=max(0,sumPH);

    end

Now i want to plot (t,pp) using either 1 or zero. How do i coding if i want to make maximum of pp become 1 instead its own value of sumPH ?

Answer 1

[Mischa Kim]

You mean normalize the plot?

 sumPH=0
 for t=1:Nt
    mew2  = (inv(Nt)*sum((y2-xhat2)));
    vari2 = (inv(Nt-1)*sum((xhat2-mew2).^(2)));
    sumPH = sumPH + vv/vari2*(rg12 - mean(mew2) -vv/2);
    pp(t) = max(0,sumPH);         % save each data point
 end

 plot(t,pp/max(pp))               % plot normalized data

If you want to plot only 1s and 0s, then you would use

    pp(t) = sumPH>0;              % save each data point
    ...
 plot(t,pp)                       % plot data

[dab483]

eg: now the result for pp=

pp= 0 0 0 0 0.3 0.4 0.5 0.5 0.5 0.5

i want to make where pp>0 it will become

new pp = 0 0 0 0 1 1 1 1 1 1  %%(and plot this)

using

if pp>0
    new pp=1;
else
    new pp=0;
end

but all i got is all zeros. any advice? Thanks

[Mischa Kim]

Yep, as outlined above, use

 sumPH=0
 for t=1:Nt
    mew2  = (inv(Nt)*sum((y2-xhat2)));
    vari2 = (inv(Nt-1)*sum((xhat2-mew2).^(2)));
    sumPH = sumPH + vv/vari2*(rg12 - mean(mew2) -vv/2);
    pp(t) = sumPH>0;     % save each data point
 end

 plot(t,pp,'*')          % plot data

...and make sure that t and pp are both row or both column vectors. Note, with '*' parameter the data points are plotted individually using a star as the marker symbol.

[dab483]

ok. i think i got it! thanks for your help.