kaikiXX - 経済統計練習帳

var y = new Array();
var x = new Array();
var sol = new Array();
 
function sample(){
 
    datanumber=200;
    numberx = 3;
 
    for (s = 1; s <datanumber+1;s++ ){
        x[s] = new Array();
    }
 
     for (s = 1; s <datanumber+1;s++ ){
        x[s][1] =1 ;
        x[s][2] = s;
        x[s][3] = Math.random();
        y[s] = 5 + 25 * x[s][2] +x[s][3]+ Math.random();
    }
 
    var sub = new ols();
    sub.x = x;
    sub.y = y;
    sub.numberx = numberx;
    sub.datanumber =datanumber;
    sub.mat();
    sol= sub.sol;
 
    alert(sol[2]);
 
 
}
 
 
class ols{
 
x:number[][]=new Array();
y:number[]=new Array();
datanumber:number;
numberx:number;
sol:number[]=new Array();
a:number[][]=new Array();
b:number[]=new Array();
 
mat():void{
 
var s:number;
 
for(s=1;s<this.numberx+1;s++){
this.a[s]=new Array();
}
 
var s1:number;
var s2:number;
var d1:number;
 
for(s1=1;s1<this.numberx+1;s1++){
for(s2=1;s2<this.numberx+1;s2++){
d1=0;
for(s=1;s<this.datanumber+1;s++){
d1=d1+this.x[s][s1]*this.x[s][s2];
}
this.a[s1][s2]=d1;
}}
 
for(s1=1;s1<this.numberx+1;s1++){
d1=0;
for(s=1;s<this.datanumber+1;s++){
d1=d1+this.x[s][s1]*this.y[s];
}
this.b[s1]=d1;
}
 
var z:number;
 
for(s1=1;s1<this.numberx+1;s1++){
for(s2=1;s2<this.numberx+1;s2++){
z=this.a[s2][s1]/this.a[s1][s1];
if(s1==s2)z=0;
for(s=1;s<this.numberx+1;s++){
this.a[s2][s]=this.a[s2][s]-z*this.a[s1][s];
}
this.b[s2]=this.b[s2]-z*this.b[s1]
}
}
 
 
 
 
for(s=1;s<this.numberx+1;s++){
this.sol[s]=this.b[s]/this.a[s][s];
}
 
 
 
 
}
 
}