重回帰８ - 経済統計練習帳

var z1=new Array();
var z2=new Array();
var z3=new Array();
 
function sample(){
 
 var data;
 
  data = new XMLHttpRequest();
  data.open("GET", "data.txt", false);
  data.send(null);
 
  var str=data.responseText;
 
var x=str.split("\n");
 
number=x.length;
 
for(s=0;s<number;s++){ 
var x1=x[s].split(",");
z1[s+1]=x1[0];
z2[s+1]=x1[1];    
z3[s+1]=x1[2];  
}
 
 
 
 
number_ols=3;
datanumber_ols=number;
 
for(s=1;s<number_ols+1;s++){
x_ols[s]=new Array();
}
 
y_ols=z1;
for(s=1;s<datanumber_ols+1;s++){
x_ols[1][s]=1;
x_ols[2][s]=z2[s];
x_ols[3][s]=z3[s];   
}
 
ols();
 
 
} 
 
 
 
var a_ols=new Array();
var b_ols=new Array();
var sol_ols=new Array();
var x_ols = new Array();
var y_ols = new Array();
var number_ols,datanumber_ols;
 
function ols(){
 
for(s=1;s<number_ols+1;s++){
a_ols[s]=new Array();
}
 
 
for(s=1;s<number_ols+1;s++){    
y1 = 0;
for (sx = 1; sx <datanumber_ols;sx++ ){
y1=y1+x_ols[s][sx]*y_ols[sx];    
}
b_ols[s]=y1;
}
 
 
for(s1=1;s1<number_ols+1;s1++){    
for(s2=1;s2<number_ols+1;s2++){ 
y1 = 0;
for (sx = 1; sx <datanumber_ols;sx++ ){
y1=y1+x_ols[s1][sx]*x_ols[s2][sx];    
}
a_ols[s1][s2]=y1;
}} 
 
 
for(s=1;s<number_ols+1;s++){
for(sx=1;sx<number_ols+1;sx++){
z=a_ols[sx][s]/a_ols[s][s];
if(s==sx)z=0;
for(tx=1;tx<number_ols+1;tx++){
a_ols[sx][tx]=a_ols[sx][tx]-z*a_ols[s][tx];
}
b_ols[sx]=b_ols[sx]-z*b_ols[s];
}
}
 
 
for(s=1;s<number_ols+1;s++){
sol_ols[s]=b_ols[s]/a_ols[s][s];
}
 
$("#memo").html(sol_ols[1]);
 
 
}