ｃ５６６ - 経済統計練習帳

package y6;
 
 
public class pro {
 
 int datanumber,number;   
 int[] y=new int[201];
 int[][] e=new int[201][4];
 double sig;
 double[] delta=new double[3];
 
 
    public static void main(String[] args) {
        pro t=new pro();
    }
 
    pro(){
 
 
    datanumber=100;
    number=2;
 
    int s;
    for(s=1;s<datanumber+1;s++){
    e[s][1]=0;    
    if(Math.random()>0.5)e[s][1]=1;    
    e[s][2]=0;    
    if(Math.random()>0.5)e[s][2]=1;  
    int h=0;
    if(e[s][1]>0.5)h=h+1;
    if(e[s][2]>0.5)h=h+1;
    if(Math.random()>0.2)h=h+1;
    y[s]=0;
    if(h==3)y[s]=1;
    }
 
    sig=0.1;
    delta[1]=0.7;
    delta[2]=0.5;
 
 
    ex_fin sub=new ex_fin();
      sub.delta=delta;
      sub.sig=sig;
      sub.y=y;
      sub.e=e;
      sub.datanumber=datanumber;
      sub.number=number;
      double u1=sub.makedata();
 
      System.out.println(u1);
 
 
 
    }
 
 
 
 
}
 
 
package y6;
 
 
public class ex_fin {
 
 int datanumber,number;   
 int[] y=new int[201];
 int[][] e=new int[201][4];
 double sig;
 double[] delta=new double[3];
 
 
    double makedata(){
 
        int s,sx;
    double zz=0;
    for(s=1;s<datanumber+1;s++){
    double z=sig;    
    for(sx=1;sx<number+1;sx++){
    z=z+delta[sx]*e[s][sx];    
    }    
    double z1=1/(1+Math.exp(-z));
    if(y[s]>0)zz=zz+(z1-1)*(z1-1);
    if(y[s]<0)zz=zz+z1*z1;
    }
 
 
    return zz;
 
 
    }
 
}