しｇｌｌ - 経済統計練習帳

var x=new Array();
var x1=new Array();
var y=new Array();
var beta=new Array();
var sig=new Array();
 
function sample(){
 
y[1]=1;
y[2]=1;
y[3]=1;
y[4]=0;
y[5]=0;
 
for(s=1;s<6;s++){
x[s]=new Array();    
}
 
x[1][1]=2;
x[2][1]=3;
x[3][1]=4;
x[4][1]=1;
x[5][1]=0;
 
x[1][2]=4;
x[2][2]=5;
x[3][2]=6;
x[4][2]=2;
x[5][2]=3;
 
for(s=1;s<3;s++){
beta[s]=new Array();   
}
 
beta[1][1]=2*Math.random()-1;
beta[1][2]=2*Math.random()-1;
 
beta[2][1]=2*Math.random()-1;
beta[2][2]=2*Math.random()-1;
 
faststep(1);
faststep(2);
 
for(s=1;s<6;s++){
x1[s]=new Array();    
}
 
for(s=1;s<6;s++){
for(sx=1;sx<3;sx++){
z1=beta[sx][0]+beta[sx][1]*x[s][1]+beta[sx][2]*x[s][2];
x1[s][sx]=0;
if(z1>0)x1[s][sx]=1;
}    
}
 
sig[1]=2*Math.random()-1;
sig[2]=2*Math.random()-1;
 
 
 
 
$("#memo").html(beta[1][0]);
 
}
 
 
function faststep(m){
 
m1=100;
m2=-100;
for(k=1;k<100;k++){
m3=(m1+m2)/2;
de=0;
for(s=1;s<6;s++){
z=beta[m][1]*x[s][1]+beta[m][2]*x[s][2]+m3;
if(z>0)de=de+1;
}
if(de>2)m1=m3;
if(de<3)m2=m3;
}
 
beta[m][0]=m3;
 
 
}