saitei6 - 経済統計練習帳

package r1;
 
class point{
 
int datanumber,number;   
int[] y=new int[201];
int[][] e=new int[201][4];
double sig;
double[] delta=new double[3];
int[][] es=new int[201][4];
 
    void makedata(){
 
        int s,sx;
 
    for(s=1;s<datanumber+1;s++){
    for(sx=1;sx<number+1;sx++){
    double z=sig;    
    int s1;
    for(s1=1;s1<number+1;s1++){
    int e1=e[s][sx];
    int h=0;
    if(s1==sx)e1=0;
    z=z+delta[sx]*e1;    
    }    
    double zz=1/(1+Math.exp(-z));
    double z1=0; 
    if(y[s]>0.5)z1=(zz-1)*(zz-1);
    if(y[s]<0.5)z1=zz*zz;
 
    z=sig;    
    for(s1=1;s1<number+1;s1++){
    int e1=e[s][sx];
    int h=0;
    if(s1==sx)e1=1;
    z=z+delta[sx]*e1;    
    }    
    zz=1/(1+Math.exp(-z));
    double z2=0; 
    if(y[s]>0.5)z2=(zz-1)*(zz-1);
    if(y[s]<0.5)z2=zz*zz;
 
    es[s][sx]=0;
    if(z1>z2)es[s][sx]=1;
    }
 
 
 
 
    }
 
}}