olk - 経済統計練習帳

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();
 
 
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] = x[s][2];
}
 
for (s = 1; s <number+1;s++ ){
x_ols[s] = x[s][1];
}
 
 
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] = 100;
 
 
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];
}
 
$("#memo").html(sol[2]);
 
 
}