[Yuima-commits] r647 - pkg/yuima/R

[noreply at r-forge.r-project.org]

Author: lorenzo
Date: 2018-05-17 10:38:01 +0200 (Thu, 17 May 2018)
New Revision: 647

Modified:
   pkg/yuima/R/AllClasses.R
   pkg/yuima/R/AuxMethodforPPR.R
   pkg/yuima/R/PointProcessClasses.R
Log:


Modified: pkg/yuima/R/AllClasses.R
===================================================================
--- pkg/yuima/R/AllClasses.R	2018-05-17 06:01:40 UTC (rev 646)
+++ pkg/yuima/R/AllClasses.R	2018-05-17 08:38:01 UTC (rev 647)
@@ -133,6 +133,8 @@
                             contains="mle"
          )
 
+
+
 setClass("yuima.qmle",representation(
 model = "yuima.model"),
 contains="mle"
@@ -179,5 +181,7 @@
 contains="summary.mle"
 )
 
+
+
 # The yuima.carma.qmle extends the S4 class "mle". It contains three slots: Estimated Levy,
 # The description of the carma model and the mle.

Modified: pkg/yuima/R/AuxMethodforPPR.R
===================================================================
--- pkg/yuima/R/AuxMethodforPPR.R	2018-05-17 06:01:40 UTC (rev 646)
+++ pkg/yuima/R/AuxMethodforPPR.R	2018-05-17 08:38:01 UTC (rev 647)
@@ -38,7 +38,7 @@
   }
   ret <- -logLik/sum(cond2,na.rm=TRUE)#+sum((param-oldpar)^2*param^2)/2
   # line 40 necessary for the development of the cod
-  # cat("\n ",logLik, param)
+  cat("\n ",logLik, param)
   
   #assign("oldpar",param,envir = my.envd1)
   

Modified: pkg/yuima/R/PointProcessClasses.R
===================================================================
--- pkg/yuima/R/PointProcessClasses.R	2018-05-17 06:01:40 UTC (rev 646)
+++ pkg/yuima/R/PointProcessClasses.R	2018-05-17 08:38:01 UTC (rev 647)
@@ -77,3 +77,30 @@
 setClass("yuima.Hawkes",
          contains="yuima.PPR"
 )
+
+# Class yuima.PPR.qmle
+
+setClass("yuima.PPR.qmle",representation(
+  model = "yuima.PPR"),
+  contains="mle"
+)
+
+
+setClass("summary.yuima.PPR.qmle",
+         representation(
+           model = "yuima.PPR"),
+         contains="summary.mle"
+)
+
+setMethod("show", "summary.yuima.PPR.qmle",
+          function (object)
+          {
+            
+            cat("Quasi-Maximum likelihood estimation\n\nCall:\n")
+            print(object at call)
+            cat("\nCoefficients:\n")
+            print(coef(object))
+            cat("\n-2 log L:", object at m2logL, "\n")
+            
+          }
+)