[Yuima-commits] r381 - in pkg/yuima: . R

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

Author: lorenzo
Date: 2015-04-26 18:18:40 +0200 (Sun, 26 Apr 2015)
New Revision: 381

Modified:
   pkg/yuima/DESCRIPTION
   pkg/yuima/R/qmle.R
Log:
Add qmle.R

Modified: pkg/yuima/DESCRIPTION
===================================================================
--- pkg/yuima/DESCRIPTION	2015-04-26 13:28:25 UTC (rev 380)
+++ pkg/yuima/DESCRIPTION	2015-04-26 16:18:40 UTC (rev 381)
@@ -1,7 +1,7 @@
 Package: yuima
 Type: Package
 Title: The YUIMA Project Package for SDEs
-Version: 1.0.67
+Version: 1.0.68
 Date: 2015-04-26
 Depends: methods, zoo, stats4, utils, expm, cubature, mvtnorm
 Author: YUIMA Project Team

Modified: pkg/yuima/R/qmle.R
===================================================================
--- pkg/yuima/R/qmle.R	2015-04-26 13:28:25 UTC (rev 380)
+++ pkg/yuima/R/qmle.R	2015-04-26 16:18:40 UTC (rev 381)
@@ -1803,31 +1803,31 @@
 }
 
 
-yuima.Vinfinity<-function(elForVInf,v){
-  # We find the infinity stationary variance-covariance matrix
-  A<-elForVInf$A
-  sigma<-elForVInf$sigma
-#   #p<-dim(A)[1]
-#   p<-elForVInf$p
-  ATrans<-elForVInf$ATrans  
-  matrixV<-elForVInf$matrixV
-  matrixV[upper.tri(matrixV,diag=TRUE)]<-v
-  matrixV<-as.matrix(forceSymmetric(matrixV))
-#matrixV[lower.tri(matrixV)]<-matrixV[upper.tri(matrixV)]
-#  l<-rbind(matrix(rep(0,p-1),p-1,1),1)
-#  matrixV<-matrix(v,p,p)
+# yuima.Vinfinity<-function(elForVInf,v){
+#   # We find the infinity stationary variance-covariance matrix
+#   A<-elForVInf$A
+#   sigma<-elForVInf$sigma
+# #   #p<-dim(A)[1]
+# #   p<-elForVInf$p
+#   ATrans<-elForVInf$ATrans  
+#   matrixV<-elForVInf$matrixV
+#   matrixV[upper.tri(matrixV,diag=TRUE)]<-v
+#   matrixV<-as.matrix(forceSymmetric(matrixV))
+# #matrixV[lower.tri(matrixV)]<-matrixV[upper.tri(matrixV)]
+# #  l<-rbind(matrix(rep(0,p-1),p-1,1),1)
+# #  matrixV<-matrix(v,p,p)
+# 
+#   lTrans<-elForVInf$lTrans
+#   l<-elForVInf$l
+#   
+#   
+#   RigSid<-l%*%elForVInf$lTrans
+#   Matrixobj<-A%*%matrixV+matrixV%*%ATrans+sigma^2*RigSid
+#   obj<-sum(Matrixobj^2)
+#   obj
+# }
 
-  lTrans<-elForVInf$lTrans
-  l<-elForVInf$l
-  
-  
-  RigSid<-l%*%elForVInf$lTrans
-  Matrixobj<-A%*%matrixV+matrixV%*%ATrans+sigma^2*RigSid
-  obj<-sum(Matrixobj^2)
-  obj
-}
 
-
 #carma.kalman<-function(y, tt, p, q, a,bvector, sigma){
 carma.kalman<-function(y, u, p, q, a,bvector, sigma, times.obs, V_inf0){  
   #new Code