Package edu.gmu.seor.prognos.unbbayesplugin.cps.Jama

Examples of edu.gmu.seor.prognos.unbbayesplugin.cps.Jama.Matrix.det()

2721
2722
2723
2724
2725
2726
2727
2728
2729
2730
2731
       
        // G component
        double d = 0;
        k11 = new Matrix(getK11(ghk.getK().getkMatrix(),ghkPot));
        h1 = new Matrix(getH1(ghk.getH().gethMatrix(), ghkPot));
        double logValue = k11.det() > 0 ? Math.log(k11.det()):k11.det();
        d = ghk.getG()
            + 0.5
            * (y1.size() * Math.log(2 * Math.PI)
                - logValue + (h1.transpose().times(
                    matInverse(k11)).times(h1)).det());
View Full Code Here
3042
3043
3044
3045
3046
3047
3048
3049
3050
3051
3052
     
      // G component
      double d = 0;
      k11 = new Matrix(getK11(ghkPot.getKValue(i).getkMatrix(),ghkPot));
      h1 = new Matrix(getH1(ghkPot.getHValue(i).gethMatrix(), ghkPot));
      double logValue = k11.det() > 0 ? Math.log(k11.det()):k11.det();
      System.out.println("Calculated values are "+ghkPot.getGValue(i)+" and "+ (0.5* (y1.size() * Math.log(2 * Math.PI)- logValue + (h1.transpose().times(matInverse(k11)).times(h1)).det())));
      d = ghkPot.getGValue(i)
          + 0.5
          * (y1.size() * Math.log(2 * Math.PI)
              - logValue + (h1.transpose().times(matInverse(k11)).times(h1)).det());
View Full Code Here
3042
3043
3044
3045
3046
3047
3048
3049
3050
3051
3052
     
      // G component
      double d = 0;
      k11 = new Matrix(getK11(ghkPot.getKValue(i).getkMatrix(),ghkPot));
      h1 = new Matrix(getH1(ghkPot.getHValue(i).gethMatrix(), ghkPot));
      double logValue = k11.det() > 0 ? Math.log(k11.det()):k11.det();
      System.out.println("Calculated values are "+ghkPot.getGValue(i)+" and "+ (0.5* (y1.size() * Math.log(2 * Math.PI)- logValue + (h1.transpose().times(matInverse(k11)).times(h1)).det())));
      d = ghkPot.getGValue(i)
          + 0.5
          * (y1.size() * Math.log(2 * Math.PI)
              - logValue + (h1.transpose().times(matInverse(k11)).times(h1)).det());
View Full Code Here
3042
3043
3044
3045
3046
3047
3048
3049
3050
3051
3052
     
      // G component
      double d = 0;
      k11 = new Matrix(getK11(ghkPot.getKValue(i).getkMatrix(),ghkPot));
      h1 = new Matrix(getH1(ghkPot.getHValue(i).gethMatrix(), ghkPot));
      double logValue = k11.det() > 0 ? Math.log(k11.det()):k11.det();
      System.out.println("Calculated values are "+ghkPot.getGValue(i)+" and "+ (0.5* (y1.size() * Math.log(2 * Math.PI)- logValue + (h1.transpose().times(matInverse(k11)).times(h1)).det())));
      d = ghkPot.getGValue(i)
          + 0.5
          * (y1.size() * Math.log(2 * Math.PI)
              - logValue + (h1.transpose().times(matInverse(k11)).times(h1)).det());
View Full Code Here
3097
3098
3099
3100
3101
3102
3103
3104
3105
3106
3107
     
      // G component
      double d = 0;
      k11 = new Matrix(getK11(ghk.getK().getkMatrix(),ghkPot));
      h1 = new Matrix(getH1(ghk.getH().gethMatrix(), ghkPot));
      double logValue = k11.det() > 0 ? Math.log(k11.det()):k11.det();
      d = ghk.getG()
          + 0.5
          * (y1.size() * Math.log(2 * Math.PI)
              - logValue + (h1.transpose().times(
                  matInverse(k11)).times(h1)).det());
View Full Code Here
3097
3098
3099
3100
3101
3102
3103
3104
3105
3106
3107
     
      // G component
      double d = 0;
      k11 = new Matrix(getK11(ghk.getK().getkMatrix(),ghkPot));
      h1 = new Matrix(getH1(ghk.getH().gethMatrix(), ghkPot));
      double logValue = k11.det() > 0 ? Math.log(k11.det()):k11.det();
      d = ghk.getG()
          + 0.5
          * (y1.size() * Math.log(2 * Math.PI)
              - logValue + (h1.transpose().times(
                  matInverse(k11)).times(h1)).det());
View Full Code Here
3097
3098
3099
3100
3101
3102
3103
3104
3105
3106
3107
     
      // G component
      double d = 0;
      k11 = new Matrix(getK11(ghk.getK().getkMatrix(),ghkPot));
      h1 = new Matrix(getH1(ghk.getH().gethMatrix(), ghkPot));
      double logValue = k11.det() > 0 ? Math.log(k11.det()):k11.det();
      d = ghk.getG()
          + 0.5
          * (y1.size() * Math.log(2 * Math.PI)
              - logValue + (h1.transpose().times(
                  matInverse(k11)).times(h1)).det());
View Full Code Here
3450
3451
3452
3453
3454
3455
3456
3457
3458
3459
3460
      ghk.setH(hc);
     
      //calculate jpd
      temp = new Matrix(kc.getkMatrix());
      temp2 = new Matrix(mchar.getmChar().getMeanMatrix().getMeanMatrix());
      double logValue = temp.det() > 0 ? Math.log(temp.det()):temp.det();
      double t = Math.log(mchar.getmChar().getProbability())+ 0.5*(logValue - mchar.getContinuousNodeList().size()*Math.log(2*Math.PI) - temp2.transpose().times(temp).times(temp2).det());
      ghk.setG(t);
     
      res.setGhkPot(ghk);
      res.hasNoDescreteStates = true;
View Full Code Here
3450
3451
3452
3453
3454
3455
3456
3457
3458
3459
3460
      ghk.setH(hc);
     
      //calculate jpd
      temp = new Matrix(kc.getkMatrix());
      temp2 = new Matrix(mchar.getmChar().getMeanMatrix().getMeanMatrix());
      double logValue = temp.det() > 0 ? Math.log(temp.det()):temp.det();
      double t = Math.log(mchar.getmChar().getProbability())+ 0.5*(logValue - mchar.getContinuousNodeList().size()*Math.log(2*Math.PI) - temp2.transpose().times(temp).times(temp2).det());
      ghk.setG(t);
     
      res.setGhkPot(ghk);
      res.hasNoDescreteStates = true;
View Full Code Here
3450
3451
3452
3453
3454
3455
3456
3457
3458
3459
3460
      ghk.setH(hc);
     
      //calculate jpd
      temp = new Matrix(kc.getkMatrix());
      temp2 = new Matrix(mchar.getmChar().getMeanMatrix().getMeanMatrix());
      double logValue = temp.det() > 0 ? Math.log(temp.det()):temp.det();
      double t = Math.log(mchar.getmChar().getProbability())+ 0.5*(logValue - mchar.getContinuousNodeList().size()*Math.log(2*Math.PI) - temp2.transpose().times(temp).times(temp2).det());
      ghk.setG(t);
     
      res.setGhkPot(ghk);
      res.hasNoDescreteStates = true;
View Full Code Here
TOP
Copyright © 2018 www.massapi.com. All rights reserved.
All source code are property of their respective owners. Java is a trademark of Sun Microsystems, Inc and owned by ORACLE Inc. Contact coftware#gmail.com.