kunbb - 経済統計練習帳

package var;
import tool.*;
 
public class pro {
 
    int[] point=new int[31];
      int[] pointx=new int[31];
    String[] data=new String[500];
   int s,number,datanumber;
   double[] mu=new double[11];
  double[] sig=new double[11];
  double[] p=new double[11];
   double[] mux=new double[11];
  double[] sigx=new double[11];
 
 
  double[] mus=new double[11];
  double[] sigs=new double[11];
  double[] ps=new double[11];
 
  double[] dmu=new double[11];
  double[] dsig=new double[11];
 
    public static void main(String[] args) {
        pro t=new pro();
    }
 
    pro(){
 
 
 
 
         readfile sub2=new readfile();
        sub2.makedata("data.csv","UTF-8");
        data=sub2.data;
        datanumber=sub2.datanumber;     
 
        for(s=2;s<datanumber+1;s++){
        data[s]=rev_csv(data[s]);
        String[] x=data[s].split(",");
        point[s-1]=num_int(x[1]);
        pointx[s-1]=num_int(x[2]);
       }
 
 
        datanumber=datanumber-1;
 
        System.out.println(datanumber);
 
        number=10;
 
 
        for(s=1;s<number+1;s++){
        mu[s]=2*s;
        sig[s]=2;
        p[s]=0.1;
        mux[s]=2*s;
        sigx[s]=2;
        }
 
 
        cal sub=new cal();
        sub.point=point;
        sub.pointx=pointx;
        sub.number=number;
        sub.datanumber=datanumber;
        sub.mu=mu;
        sub.sig=sig;
        sub.mux=mux;
        sub.sigx=sigx;
        sub.p=p;
        double u1=sub.makedata();
 
        System.out.println(u1);
 
 
 
 
 
    }
 
       int num_int(String str){
 
    int nx;
 
    nx=0;
 
 try{
 
    nx=Integer.parseInt(str);
 
 }catch (Exception e) {}
 
 return nx;
 
 
 
}    
 
    String rev_csv(String str){
 
        String str5=str;
        int px=str5.indexOf("\"");
 
 
 
        while(px>-1){
       str5=rev(str5);
       px=str5.indexOf("\"");   
        }
        return str5;
    }
 
    String rev(String str){
 
        int p1=str.indexOf("\"");
        int p2=str.indexOf("\"",p1+1);
 
     String str1=str.substring(0,p1);   
        String str2=str.substring(p1+1,p2);  
        String str3=str.substring(p2+1,str.length());
        str2=str2.replace(",","");
 
        return str1+str2+str3;
 
    }
 
 
 
 
}
 
 
package var;
 
public class cal{
 
   int datanumber;
    int[] point=new int[50];
    double[] mux=new double[11];
  double[] sigx=new double[11];
    double[] mu=new double[11];
  double[] sig=new double[11];
  double[] p=new double[11];
  int number;
     int[] pointx=new int[50];
 
    double makedata(){
 
 
        int s,sx;
 
 
 
 
 
 
    double v1=0;
    for(s=1;s<number+1;s++){
    for(sx=1;sx<datanumber+1;sx++){
    v1=v1+p[s]*point[sx]*f(sx,mu[s],sig[s]);
    }
    }
 
      double v2=0;
    for(s=1;s<number+1;s++){
    for(sx=1;sx<datanumber+1;sx++){
    v2=v2+p[s]*pointx[sx]*f(sx,mux[s],sigx[s]);
    }
    }  
 
    v1=v1*v2;
 
    return v1;
 
 
 
 
 
 
 
    }
 
        double f(double x,double mux,double sigx){
double pi,f1,f2;
 
        pi=Math.PI;
 
        f2=Math.sqrt(2*pi*sigx*sigx);
        f1=-(x-mux)*(x-mux)/(2*sigx*sigx);
 
        return Math.exp(f1)/f2;
}
 
 
 
 
}