b47 - 経済統計練習帳

var number_ols;
var a=new Array();
var b=new Array();
var sol=new Array();
var x_ols = new Array();
var y_ols = new Array();
 
var x=new Array();
var y=new Array();
var alpha,beta;
 
function sample(){
 
number=5;
 
for(s=1;s<6;s++){
x[s]=new Array();    
 
}
 
x[1][2]=1;
x[2][2]=2;
x[3][2]=3;
x[4][2]=4;
x[5][2]=5;
x[1][1]=1;
x[2][1]=1;
x[3][1]=1;
x[4][1]=1;
x[5][1]=1;
 
y[1]=0;
y[2]=0;
y[3]=0;
y[4]=1;
y[5]=1;
 
ols();
 
s=4;
alert(sol[1]*x[s][1]+sol[2]*x[s][2]);
 
 
 
 
 
 
}
 
function ols(){
 
 
number_ols = 2;    
 
for(s=1;s<number_ols+1;s++){
a[s]=new Array();
}
 
 
for (s = 1; s <number+1;s++ ){
y_ols[s] = y[s];
}
 
for (s = 1; s <number+1;s++ ){
x_ols[s] = x[s][2];
}
 
 
y1 = 0;
for (s = 1; s <number+1;s++ ){
y1=y1+y_ols[s];    
}
 
b[1] = y1;
 
y1 = 0;
for (s = 1; s < number+1;s++ ){
y1=y1+y_ols[s]*x_ols[s];    
}
 
b[2]=y1;
 
a[1][1] = number;
 
 
y1 = 0;
for (s = 1; s <number+1;s++ ){
y1=y1+x_ols[s]*x_ols[s];    
}
 
a[2][2] = y1;
 
y1 = 0;
for (s = 1; s <number+1;s++ ){
y1=y1+x_ols[s];    
}
 
a[1][2] = y1;
a[2][1] = y1;
 
for(s=1;s<number_ols+1;s++){
for(sx=1;sx<number_ols+1;sx++){
z=a[sx][s]/a[s][s];
if(s==sx)z=0;
for(tx=1;tx<number_ols+1;tx++){
a[sx][tx]=a[sx][tx]-z*a[s][tx];
}
b[sx]=b[sx]-z*b[s];
}
}
 
 
for(s=1;s<number_ols+1;s++){
sol[s]=b[s]/a[s][s];
}
 
 
 
}