new;

#include c:\gauss15\src\icss.src;


fname1 = "elec_bernd_sea.xls";

rname1 = "a1:a3776"; 

vls1 = reshape(error(0),9,1);

/*NSW*/

el_LNSQNSQ = xlsReadM(fname1, rname1, 1, vls1[1]);  
el_LNSQNSQ = packr(el_LNSQNSQ);
el_LNSQNSQ = el_LNSQNSQ - meanc(el_LNSQNSQ);
cri = 0|1|4;
{cp1, nbre1} = icss(el_LNSQNSQ,2,cri);

output file = test1.out on;

	format /rd 1,4;
print "vector containing the change points = " cp1;
print "number of structural breaks =" nbre1;

/*QLD*/

rname2 = "b1:b3776"; 
el_LNSQQLD = xlsReadM(fname1, rname2, 1, vls1[1]);  
el_LNSQQLD = packr(el_LNSQQLD);
el_LNSQQLD = el_LNSQQLD - meanc(el_LNSQQLD);
cri = 0|1|4;
{cp2, nbre2} = icss(el_LNSQQLD,2,cri);

output file = test1.out on;

	format /rd 1,4;
print "vector containing the change points = " cp2;
print "number of structural breaks =" nbre2;


/*SA*/
rname3 = "c1:c3776"; 
el_LNSQSA = xlsReadM(fname1, rname3, 1, vls1[1]);  
el_LNSQSA = packr(el_LNSQSA);
el_LNSQSA = el_LNSQSA - meanc(el_LNSQSA);
cri = 0|1|4;
{cp3, nbre3} = icss(el_LNSQSA,2,cri);

output file = test1.out on;

	format /rd 1,4;
print "vector containing the change points = " cp3;
print "number of structural breaks =" nbre3;


/*TAS*/
rname4 = "d1:d3776"; 
el_LNSQTAS = xlsReadM(fname1, rname1, 1, vls1[1]);  
el_LNSQTAS = packr(el_LNSQTAS);
el_LNSQTAS = el_LNSQTAS - meanc(el_LNSQTAS);
cri = 0|1|4;
{cp4, nbre4} = icss(el_LNSQTAS,2,cri);

output file = test1.out on;

	format /rd 1,4;
print "vector containing the change points = " cp4;
print "number of structural breaks =" nbre4;


/*VIC*/
rname5 = "e1:e3776"; 
el_LNSQVIC = xlsReadM(fname1, rname5, 1, vls1[1]);  
el_LNSQVIC = packr(el_LNSQVIC);
el_LNSQVIC = el_LNSQVIC - meanc(el_LNSQVIC);
cri = 0|1|4;
{cp5, nbre5} = icss(el_LNSQVIC,2,cri);

output file = test1.out on;

	format /rd 1,4;
print "vector containing the change points = " cp5;
print "number of structural breaks =" nbre5;
