じんｋｋ - 経済統計練習帳

var x=new Array();
var y=new Array();
var alpha=new Array();
var a=new Array();
 
function myFunction() {
 
 
for(s=1;s<101;s++){
x[s]=new Array();
}
 
 
for(s=1;s<101;s++){
x[s][1]=Math.random();
x[s][2]=Math.random();
y[s]=0;
if(x[s][1]+x[s][2]>1)y[s]=1;
}
 
 
for(sx=0;sx<3;sx++){
alpha[sx]=0;
}
 
minz=9999;
for(tr=1;tr<1001;tr++){
 
z=0;
for(sx=0;sx<3;sx++){
m1=20*Math.random()-10;
a[sx]=Math.floor(m1)+1;
}
 
 
 
 
for(s=1;s<101;s++){
a1=a[0];
for(sx=1;sx<3;sx++){
a1=a1+a[sx]*x[s][sx];
}
z1=1/(1+Math.exp(-a1));
z=z+(y[s]-z1)*(y[s]-z1);
}
 
 
 
if(z<minz){
 
for(sx=0;sx<3;sx++){
alpha[sx]=a[sx];
}
 
minz=z;
}
 
}
 
Logger.log(minz);
 
 
 var id="14QxZ4Y_0lRBpl2jMD2bMP3vcOQLLghfbmVvD7Mb2rTA"; 
 ex1 = SpreadsheetApp.openById(id);
 var sh = ex1.getSheetByName("data"); 
 
 
for(s=1;s<101;s++){ 
z=alpha[0]+alpha[1]*x[s][1]+alpha[2]*x[s][2];
sh.getRange(s+1,1).setValue(y[s]);  
sh.getRange(s+1,2).setValue(z); 
}
 
 
 
}