最尤法 - プログラム倉庫

package p01;
 
 
public class pro {
 
 double[] x=new double[1001];
 double[] y=new double[1001];
 int s;   
 double mu,sig;
 double a0,delta_mu,delta_sig;
 double delta;
 
    public static void main(String[] args) {
       pro test=new pro();
    }
 
 pro(){   
 
 
     for(s=1;s<101;s++){
     x[s]=4.5+Math.random();
     }
 
     mu=5.5;
     sig=0.7;
 
     delta=0.1;
 
   a0=logz(mu,sig);
   delta_mu=logz(mu+delta,sig)-a0;  
   System.out.println(delta_mu);
 
 
 
 }
 
 double logz(double mu,double sig){
 int s;
double[] y=new double[1001];
     for(s=1;s<101;s++){
     y[s]=f(x[s],mu,sig);
     }
 
     double z;
 
    z=0;
    for(s=1;s<101;s++){
    z=z+Math.log(y[s]);
    }
   return z;
 
 }
double f(double x1,double mu,double sig){
 
    double f1,f2,f3;
    double pi;
    pi=Math.PI;
    f1=-(x1-mu)*(x1-mu)/(2*sig*sig);   
    f2=Math.exp(f1);
    f3=Math.sqrt(2*sig*sig*pi);
 
    return f2/f3; 
 
}
}