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

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

3479
3480
3481
3482
3483
3484
3485
3486
3487
3488
3489
        resH.add(hc);
       
        //calculate jpd
        temp = new Matrix(kc.getkMatrix());
        temp2 = new Matrix(mchar.getMeanMatValue(i).getMeanMatrix());
        double logValue = temp.det() > 0 ? Math.log(temp.det()):temp.det();
        double t = Math.log(mchar.getJPDValue(i))+ 0.5*(logValue - mchar.getContinuousNodeList().size()*Math.log(2*Math.PI) - temp2.transpose().times(temp).times(temp2).det());
        resG.add(t);
        }
       
        res.setgComponents(resG);
View Full Code Here
3479
3480
3481
3482
3483
3484
3485
3486
3487
3488
3489
        resH.add(hc);
       
        //calculate jpd
        temp = new Matrix(kc.getkMatrix());
        temp2 = new Matrix(mchar.getMeanMatValue(i).getMeanMatrix());
        double logValue = temp.det() > 0 ? Math.log(temp.det()):temp.det();
        double t = Math.log(mchar.getJPDValue(i))+ 0.5*(logValue - mchar.getContinuousNodeList().size()*Math.log(2*Math.PI) - temp2.transpose().times(temp).times(temp2).det());
        resG.add(t);
        }
       
        res.setgComponents(resG);
View Full Code Here
3479
3480
3481
3482
3483
3484
3485
3486
3487
3488
3489
        resH.add(hc);
       
        //calculate jpd
        temp = new Matrix(kc.getkMatrix());
        temp2 = new Matrix(mchar.getMeanMatValue(i).getMeanMatrix());
        double logValue = temp.det() > 0 ? Math.log(temp.det()):temp.det();
        double t = Math.log(mchar.getJPDValue(i))+ 0.5*(logValue - mchar.getContinuousNodeList().size()*Math.log(2*Math.PI) - temp2.transpose().times(temp).times(temp2).det());
        resG.add(t);
        }
       
        res.setgComponents(resG);
View Full Code Here
3605
3606
3607
3608
3609
3610
3611
3612
3613
3614
3615
      double t1;
      double t2;
      try {
        temp = new Matrix(ghkPot.getKValue(i).getkMatrix());
        temp2 = new Matrix(mm.getMeanMatrix());
        double logValue = temp.det() > 0 ? Math.log(temp.det()):temp.det();
        t1 = ghkPot.getGValue(i);
        t2 = - 0.5*(logValue - ghkPot.getContinuousNodeList().size()*Math.log(2*Math.PI) - temp2.transpose().times(temp).times(temp2).det());
      } catch (Exception e) {
        t1 = ghkPot.getGValue(i);
        t2 = - ghkPot.getContinuousNodeList().size()*Math.log(2*Math.PI);
 
View Full Code Here
3605
3606
3607
3608
3609
3610
3611
3612
3613
3614
3615
      double t1;
      double t2;
      try {
        temp = new Matrix(ghkPot.getKValue(i).getkMatrix());
        temp2 = new Matrix(mm.getMeanMatrix());
        double logValue = temp.det() > 0 ? Math.log(temp.det()):temp.det();
        t1 = ghkPot.getGValue(i);
        t2 = - 0.5*(logValue - ghkPot.getContinuousNodeList().size()*Math.log(2*Math.PI) - temp2.transpose().times(temp).times(temp2).det());
      } catch (Exception e) {
        t1 = ghkPot.getGValue(i);
        t2 = - ghkPot.getContinuousNodeList().size()*Math.log(2*Math.PI);
 
View Full Code Here
3605
3606
3607
3608
3609
3610
3611
3612
3613
3614
3615
      double t1;
      double t2;
      try {
        temp = new Matrix(ghkPot.getKValue(i).getkMatrix());
        temp2 = new Matrix(mm.getMeanMatrix());
        double logValue = temp.det() > 0 ? Math.log(temp.det()):temp.det();
        t1 = ghkPot.getGValue(i);
        t2 = - 0.5*(logValue - ghkPot.getContinuousNodeList().size()*Math.log(2*Math.PI) - temp2.transpose().times(temp).times(temp2).det());
      } catch (Exception e) {
        t1 = ghkPot.getGValue(i);
        t2 = - ghkPot.getContinuousNodeList().size()*Math.log(2*Math.PI);
 
View Full Code Here
3646
3647
3648
3649
3650
3651
3652
3653
3654
3655
3656
      mc.setMeanMatrix(mm);
     
      //calculate jpd
      temp = new Matrix(ghkPot.getGhkPot().getK().getkMatrix());
      temp2 = new Matrix(mm.getMeanMatrix());
      double logValue = temp.det() > 0 ? Math.log(temp.det()):temp.det();
      double t1 = ghkPot.getGhkPot().getG();
      double t2 = - 0.5*(logValue - ghkPot.getContinuousNodeList().size()*Math.log(2*Math.PI) - temp2.transpose().times(temp).times(temp2).det());
      System.out.println("t1 is "+t1);
      System.out.println("t2 is "+t2);
      double t = Math.exp(t1 + t2);
View Full Code Here
3646
3647
3648
3649
3650
3651
3652
3653
3654
3655
3656
      mc.setMeanMatrix(mm);
     
      //calculate jpd
      temp = new Matrix(ghkPot.getGhkPot().getK().getkMatrix());
      temp2 = new Matrix(mm.getMeanMatrix());
      double logValue = temp.det() > 0 ? Math.log(temp.det()):temp.det();
      double t1 = ghkPot.getGhkPot().getG();
      double t2 = - 0.5*(logValue - ghkPot.getContinuousNodeList().size()*Math.log(2*Math.PI) - temp2.transpose().times(temp).times(temp2).det());
      System.out.println("t1 is "+t1);
      System.out.println("t2 is "+t2);
      double t = Math.exp(t1 + t2);
View Full Code Here
3646
3647
3648
3649
3650
3651
3652
3653
3654
3655
3656
      mc.setMeanMatrix(mm);
     
      //calculate jpd
      temp = new Matrix(ghkPot.getGhkPot().getK().getkMatrix());
      temp2 = new Matrix(mm.getMeanMatrix());
      double logValue = temp.det() > 0 ? Math.log(temp.det()):temp.det();
      double t1 = ghkPot.getGhkPot().getG();
      double t2 = - 0.5*(logValue - ghkPot.getContinuousNodeList().size()*Math.log(2*Math.PI) - temp2.transpose().times(temp).times(temp2).det());
      System.out.println("t1 is "+t1);
      System.out.println("t2 is "+t2);
      double t = Math.exp(t1 + t2);
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.