ｐ９９７６ - 経済統計練習帳

package gei;
 
 
public class pro {
 
    int[] point=new int[1001];
    int s,number;
    double[] mu=new double[11];
    double[] sig=new double[11];
     double[][] rate=new double[1001][21];
      double[] rate_topic=new double[101];
 
    public static void main(String[] args) {
        pro test=new pro();
    }
    pro(){
 
        number=100;
 
        for(s=1;s<number+1;s++){
        point[s]=s;    
        }
 
 
    for(s=1;s<11;s++){
    mu[s]=25+5*s;    
    sig[s]=10;
    }    
 
 
    for(s=1;s<11;s++){
    rate_topic[s]=0.1;
    }    
 
    int tr;
    for(tr=1;tr<11;tr++){
 
    int sx;
    for(s=1;s<number+1;s++){
    double p1=0;
    for(sx=1;sx<11;sx++){  
    p1=p1+f(point[s],mu[sx],sig[sx])*rate_topic[sx];    
    }
    for(sx=1;sx<11;sx++){  
    rate[s][sx]=f(point[s],mu[sx],sig[sx])*rate_topic[sx]/p1;    
    }
    }
 
    for(sx=1;sx<11;sx++){  
    mu[sx]=seek_mu(sx);
    sig[sx]=seek_sig(sx);
 
    }
 
    System.out.println(tr+","+mu[5]);
 
    }
 
 
 
    }
 
      double seek_mu(int sx){
 
        double m1,m2;
        m1=0;
        int s;
        for(s=1;s<number+1;s++){
        m1=m1+point[s]*rate[s][sx];    
        }
 
        m2=0;
        for(s=1;s<number+1;s++){
        m2=m2+rate[s][sx];    
        }
 
 
 
        return m1/m2;
    }
 
     double seek_sig(int sx){
 
        double m1,m2;
        m1=0;
        int s;
        for(s=1;s<number+1;s++){
        m1=m1+(point[s]-mu[sx])*(point[s]-mu[sx])*rate[s][sx];    
        }
 
        m2=0;
        for(s=1;s<number+1;s++){
        m2=m2+rate[s][sx];    
        }
 
        double m3=m1/m2;
        m3=Math.sqrt(m3);
 
        return m3;
    }
 
 
 
      double f(double x1,double mux,double sigx){
double pi,f1,f2;
 
        pi=Math.PI;
 
        f2=Math.sqrt(2*pi*sigx*sigx);
        f1=-(x1-mux)*(x1-mux)/(2*sigx*sigx);
 
        return Math.exp(f1)/f2;
}
 
 
 
 
 
 
}