ばかな - 経済統計練習帳

package word;
 
public class pro {
 
    double[] alpha=new double[4];
    double[][] beta=new double[4][3];
    int number,datanumber;
    int[] y=new int[101];
    double[][] x=new double[101][3];
    int[][] e=new int[3][101];
 
     double a;
    double[] b=new double[3];
 
    public static void main(String[] args) {
        pro te=new pro();
    }
 
    pro(){
 
 
        number=2;
        datanumber=50;
 
        int s;
 
          for(s=1;s<30;s++){
        x[s][1]=2+Math.random();
        x[s][2]=5+Math.random();
        y[s]=1;
        }
 
         for(s=31;s<datanumber+1;s++){
        x[s][1]=-3+Math.random();
        x[s][2]=-7+Math.random();
        y[s]=0;
        }
 
         int h;
 
        for(s=1;s<datanumber+1;s++){
        e[1][s]=0;
        h=0;    
        if(Math.random()>0.1)h=h+1;
        if(y[s]>0.5)h=h+1;
        if(h==2)e[1][s]=1;
        }
 
        for(s=1;s<datanumber+1;s++){
        e[2][s]=0;
        h=0;    
        if(Math.random()>0.2)h=h+1;
        if(y[s]>0.5)h=h+1;
        if(h==2)e[2][s]=1;
        }
 
        int[] ex=new int[100]; 
        for(s=1;s<datanumber+1;s++){
        ex[s]=e[1][s];    
        }
 
   optimal sub2=new optimal();
         sub2.number=number;
         sub2.datanumber=datanumber;
         sub2.x=x;
         sub2.y=ex;
         sub2.makedata();
         a=sub2.alpha;
         b=sub2.beta;
 
         alpha[1]=a;
         for(s=1;s<number+1;s++){
         beta[1][s]=b[s];    
         }
 
         int sx;
         for(s=1;s<datanumber+1;s++){
         double x1=alpha[1];
         for(sx=1;sx<number+1;sx++){
         x1=x1+beta[1][sx]*x[s][sx];    
         }
         double f1=1/(1+Math.exp(-x1));
         }
 
 
        for(s=1;s<datanumber+1;s++){
        ex[s]=e[2][s];    
        }
 
 
   optimal sub3=new optimal();
         sub3.number=number;
         sub3.datanumber=datanumber;
         sub3.x=x;
         sub3.y=ex;
         sub3.makedata();
         a=sub3.alpha;
         b=sub3.beta;
 
         alpha[2]=a;
         for(s=1;s<number+1;s++){
         beta[2][s]=b[s];    
         }
 
         for(s=1;s<datanumber+1;s++){
         double x1=alpha[2];
         for(sx=1;sx<number+1;sx++){
         x1=x1+beta[2][sx]*x[s][sx];    
         }
         double f1=1/(1+Math.exp(-x1));
         }
 
 
 
         double[][] q=new double[101][3];
 
        for(s=1;s<datanumber+1;s++){
        q[s][1]=e[1][s];
        q[s][2]=e[2][s];
        }
 
   optimal sub4=new optimal();
         sub4.number=number;
         sub4.datanumber=datanumber;
         sub4.x=q;
         sub4.y=y;
         sub3.makedata();
         a=sub3.alpha;
         b=sub3.beta;
 
         alpha[3]=a;
         for(s=1;s<number+1;s++){
         beta[3][s]=b[s];    
         }
 
         for(s=1;s<datanumber+1;s++){
         double x1=alpha[3];
         for(sx=1;sx<number+1;sx++){
         x1=x1+beta[3][sx]*q[s][sx];    
         }
         double f1=1/(1+Math.exp(-x1));
         System.out.println(y[s]+","+f1);
         }
 
 
 
 
 
 
 
 
    }
 
 
 
}