tiratira - 経済統計練習帳

var x=new Array(); 
var y=new Array(); 
var  beta=new Array();
var alpha;
 
function sample(){
 
 
datanumber=100;
 
for(s=1;s<datanumber+1;s++){
x[s]=new Array();    
}
 
 
for(s=1;s<datanumber+1;s++){
x[s][1]=Math.random();
x[s][2]=Math.random();
y[s]=1;
if(Math.random()>0.5)y[s]=0;
}
 
beta[1]=Math.random();
beta[2]=Math.random();
 
 
b1=100;
b2=-100;
for(tr=1;tr<100;tr++){
b3=(b1+b2)/2;
n1=f(b3);
if(n1>50)b1=b3;
if(n1<50)b2=b3;
}
alpha=b3;
 
 
paint();
 
}
 
 
function paint(){
 
str="<table>";
for(s=1;s<datanumber+1;s++){
str=str+"<tr><td>"+x[s][1]+"</td>";
str=str+"<td>"+x[s][2]+"</td>";
y1=alpha+beta[1]*x[s][1]+beta[2]*x[s][2];
u1=1/(1+Math.exp(-y1));
str=str+"<td>"+u1+"</td>";
str=str+"<td>"+y[s]+"</td></tr>";
}
str=str+"</table>";
 
$("#memo").html(str);
 
 
}
 
 
 
 
function f(a1){
 
c1=0;
for(s=1;s<datanumber+1;s++){
z=a1;
for(sx=1;sx<3;sx++){
z=z+beta[sx]*x[s][sx];
} 
z1=1/(1+Math.exp(-z));
if(z1>0.5)c1=c1+1;
}
 
return c1;
 
}