練習問題３－１ - 経済統計練習帳

var number;
var z1=new Array();
var z2=new Array();
 
function sample(){
 
var data;
 
data = new XMLHttpRequest();
data.open("GET", "data3.csv", 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[1];
z2[s+1]=x1[2];    
}
 
number_ols=2;    
datanumber_ols=10;
 
for(s=1;s<number_ols+1;s++){ 
x_ols[s]=new Array();    
}
 
for(s=1;s<11;s++){   
x_ols[1][s]=1;
x_ols[2][s]=z1[s];    
y_ols[s]=z2[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]+","+sol_ols[2]);
 
 
}
 
1960,157,143
1961,162,146
1962,169,153
1963,176,160
1964,188,169
1965,200,180
1966,211,190
1967,220,196
1968,230,207
1969,237,215
1970,247,220
1971,256,228
1972,268,242
1973,287,253
1974,285,251
1975,290,257
1976,301,271
1977,311,283
1978,326,295
1979,335,302
1980,337,301
1981,345,305
1982,348,308
1983,358,324
1984,384,341
1985,396,357
1986,409,371
1987,415,382
1988,432,397
1989,440,406
1990,448,413
1991,449,411
1992,463,422
1993,470,434
1994,477,449
1995,491,461
1996,504,475
1997,518,491
1998,535,515
1999,577,542
2000,605,567
2001,617,581
2002,636,599
2003,651,619
2004,674,642