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

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

Author: hinohide
Date: 2011-05-19 14:11:48 +0200 (Thu, 19 May 2011)
New Revision: 154

Modified:
   pkg/yuima/DESCRIPTION
   pkg/yuima/R/asymptotic_term.R
Log:
a bug in asymptotic term is fixed

Modified: pkg/yuima/DESCRIPTION
===================================================================
--- pkg/yuima/DESCRIPTION	2011-05-07 08:38:02 UTC (rev 153)
+++ pkg/yuima/DESCRIPTION	2011-05-19 12:11:48 UTC (rev 154)
@@ -1,8 +1,8 @@
 Package: yuima
 Type: Package
 Title: The YUIMA Project package (unstable version)
-Version: 0.1.189
-Date: 2011-05-07
+Version: 0.1.190
+Date: 2011-05-19
 Depends: methods, zoo, stats4, utils
 Suggests: cubature, mvtnorm
 Author: YUIMA Project Team.

Modified: pkg/yuima/R/asymptotic_term.R
===================================================================
--- pkg/yuima/R/asymptotic_term.R	2011-05-07 08:38:02 UTC (rev 153)
+++ pkg/yuima/R/asymptotic_term.R	2011-05-19 12:11:48 UTC (rev 154)
@@ -390,10 +390,12 @@
   # function to calculate sigma in thesis p5
   # require: aMat
   funcsigma <- function(e=0){
+
     division <- nrow(X.t0)
     sigma <- matrix(0,k.size,k.size) #size:matrix[k.size,k.size]
     for(t in 1:division){
-      sigma <- sigma+(aMat[,,t]%*%t(aMat[,,t])) /(division-1) #calculate sigma
+##      sigma <- sigma+(aMat[,,t]%*%t(aMat[,,t])) /(division-1) #calculate sigma: old code, modified 20110519
+      sigma <- sigma+(t(aMat[,,t])%*%(aMat[,,t])) /(division-1) #calculate sigma
     }
     if(any(eigen(sigma)$value<=0.0001)){
     # Singularity check