mesu - 経済統計練習帳

package r1;
 
 
public class pro {
 
  double[][] x=new double[700][5];
  int datanumber,number;
 
  double[] alpha=new double[4];
  double[][] beta=new double[4][4];
   int[][] e=new int[700][6];
 
  int[][] es=new int[700][6];
  int[] y=new int[600]; 
   double sig;
     double[] delta=new double[4];
 
    public static void main(String[] args) {
       pro test=new pro();
    }
 
    pro(){
 
 
        int s;
 
         datanumber=100;
        number=2;
 
        for(s=1;s<datanumber+1;s++){
        x[s][1]=5*Math.random()-1;    
        x[s][2]=2*Math.random()+1; 
        if(x[s][1]+2*x[s][2]+0.01*Math.random()>3.5)y[s]=1;
        }
 
     faststep sub=new faststep();
     sub.x=x;
     sub.datanumber=datanumber;
     sub.number=number;
     sub.makedata();
     alpha=sub.alpha;
     beta=sub.beta;   
 
 
      hidden sub2=new hidden();
  sub2.datanumber=datanumber;
  sub2.number=number;
  sub2.x=x;
  sub2.alpha=alpha;
  sub2.beta=beta; 
  sub2.makedata();
  e=sub2.e;
 
 
  laststep sub3=new laststep();
  sub3.datanumber=datanumber;
  sub3.number=number;
  sub3.e=e;
  sub3.makedata();
  sig=sub3.sig;
  delta=sub3.delta;
 
int tr;
 
for(tr=1;tr<100;tr++){
 
  point sub4=new point();
  sub4.datanumber=datanumber;
  sub4.number=number;
  sub4.y=y;
  sub4.sig=sig;
  sub4.delta=delta; 
  sub4.makedata();
  es=sub4.es;
 
  for(s=1;s<number+1;s++){
  drive sub5=new drive();
  sub5.datanumber=datanumber;
  sub5.number=number;
  sub5.es=es;
  sub5.x=x;
  sub5.alpha=alpha;
  sub5.beta=beta; 
  sub5.makedata(s);
  beta=sub5.beta;
  alpha=sub5.alpha;
  }
 
  drive_fin sub6=new drive_fin();
  sub6.datanumber=datanumber;
  sub6.number=number;
  sub6.e=e;
  sub6.y=y;
  sub6.sig=sig;
  sub6.delta=delta; 
  sub6.makedata();
  sig=sub6.sig;
  delta=sub6.delta;
 
 
 
  hidden sub7=new hidden();
  sub7.datanumber=datanumber;
  sub7.number=number;
  sub7.x=x;
  sub7.alpha=alpha;
  sub7.beta=beta; 
  sub7.makedata();
  e=sub7.e;
 
    }
 
 
 
  for(s=1;s<datanumber+1;s++){
  double y1=sig+delta[1]*e[s][1]+delta[2]*e[s][2];       
  System.out.println(y[s]+","+(1/(1+Math.exp(-y1))));
    }
 
 
 
 
 
    }
}