kaikibunsek - 経済統計練習帳

var x=new Array();
var y=new Array();
var point=new Array();
var sol=new Array(); 
 
function sample() {
 
number=100;
 
    for(s=1;s<number+1;s++){
     x[s]=s;
     y[s]=5+4*s;
    }
 
ols();
 
$("#memo").html(sol[2]);
 
}
 
function ols(){
 
var a=new Array();
var b=new Array();
 
 
 
for(s=1;s<3;s++){
a[s]=new Array();    
}
 
 
b1=0;    
for(s=1;s<number+1;s++){
b1=b1+y[s];    
}    
b[1]=b1;
 
 
 
b1=0;    
for(s=1;s<number+1;s++){
b1=b1+y[s]*x[s];    
}    
b[2]=b1;
 
 
a1=0;    
for(s=1;s<number+1;s++){
a1=a1+1*1;    
}    
a[1][1]=a1;  
 
a1=0;    
for(s=1;s<number+1;s++){
a1=a1+x[s];    
}    
a[1][2]=a1;  
a[2][1]=a1;
 
a1=0;    
for(s=1;s<number+1;s++){
a1=a1+x[s]*x[s];    
}    
a[2][2]=a1;  
 
 
for(s1=1;s1<3;s1++){
for(s2=1;s2<3;s2++){
z=a[s2][s1]/a[s1][s1];
if(s1==s2)z=0;
for(s=1;s<3;s++){
a[s2][s]=a[s2][s]-z*a[s1][s];
}
b[s2]=b[s2]-z*b[s1]
}
}
 
 
 
 
for(s=1;s<3;s++){
sol[s]=b[s]/a[s][s];
}
 
 
 
 
}