ｖ５９９９ｐｘ - 経済統計練習帳

package lda;
 
 
public class pro {
 
  int[][] v=new int[7][7];
  double[][] p=new double[3][7];
   double[][] q=new double[7][3];
 
     double[][] qs=new double[7][3];
    double[][] ps=new double[3][7];
 
    public static void main(String[] args) {
        pro t=new pro();
    }
 
    pro(){
 
        v[1][1]=1;
        v[1][2]=1;
        v[2][2]=1;
        v[2][3]=1;
        v[3][1]=1;
        v[3][3]=1;
 
 
        v[4][4]=1;
        v[4][5]=1;
        v[5][5]=1;
        v[5][6]=1;
        v[6][4]=1;
        v[6][6]=1;
        int s;
        for(s=1;s<7;s++){
        p[1][s]=Math.random();    
        p[2][s]=Math.random();   
        q[s][1]=Math.random();
        q[s][2]=Math.random();
        }
 
        int sx,n;
        for(n=1;n<7;n++){
        double p1=1;
        double p2=1;
      for(s=1;s<7;s++){
      if(v[n][s]==1)p1=p1*p[1][s]*q[n][1]; 
      if(v[n][s]==0)p1=p1*0.01*q[1][1]; 
      if(v[n][s]==2)p2=p2*p[2][s]*q[n][2]; 
      if(v[n][s]==0)p2=p2*0.01*q[n][2]; 
      }
      qs[n][1]=p1/(p1+p2);
      qs[n][2]=p2/(p1+p2);
      System.out.println(qs[n][1]+","+qs[n][2]);  
        }
 
 
      for(n=1;n<7;n++){
      double p1=1;
      double p2=1;
      for(s=1;s<7;s++){
      if(v[s][n]==1)p1=p1*p[1][n]*q[s][1]; 
      if(v[s][n]==0)p1=p1*0.01*q[s][1]; 
      if(v[s][n]==2)p2=p2*p[2][n]*q[s][2]; 
      if(v[s][n]==0)p2=p2*0.01*q[s][2]; 
      }
      ps[1][n]=p1/(p1+p2);
      ps[2][n]=p2/(p1+p2);
      System.out.println(ps[1][n]+","+ps[2][n]);  
        }
 
 
 
 
 
 
        }
 
 
 
 
 
 
 
    }