[Depmix-commits] r202 - pkg/R trunk/R

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

Author: maarten
Date: 2008-07-01 10:59:49 +0200 (Tue, 01 Jul 2008)
New Revision: 202

Modified:
   pkg/R/responseMVN.R
   trunk/R/responseMVN.R
Log:
fixed remaining bug in responseMVN related to invsigma

Modified: pkg/R/responseMVN.R
===================================================================
--- pkg/R/responseMVN.R	2008-06-30 10:59:02 UTC (rev 201)
+++ pkg/R/responseMVN.R	2008-07-01 08:59:49 UTC (rev 202)
@@ -50,7 +50,7 @@
       if (length(mean) != NROW(invSigma)) {
           stop("mean and invSigma have non-conforming size")
       }
-      distval <- mahalanobis(x, center = mean, cov = invsigma, inverted=TRUE) 
+      distval <- mahalanobis(x, center = mean, cov = invSigma, inverted=TRUE)
     } 
     if(missing(logdet)) logdet <- sum(log(eigen(sigma, symmetric = TRUE, only.values = TRUE)$values))
     logretval <- -(ncol(x) * log(2 * pi) + logdet + distval)/2
Modified: trunk/R/responseMVN.R
===================================================================
--- trunk/R/responseMVN.R	2008-06-30 10:59:02 UTC (rev 201)
+++ trunk/R/responseMVN.R	2008-07-01 08:59:49 UTC (rev 202)
@@ -50,7 +50,7 @@
       if (length(mean) != NROW(invSigma)) {
           stop("mean and invSigma have non-conforming size")
       }
-      distval <- mahalanobis(x, center = mean, cov = invsigma, inverted=TRUE) 
+      distval <- mahalanobis(x, center = mean, cov = invSigma, inverted=TRUE)
     } 
     if(missing(logdet)) logdet <- sum(log(eigen(sigma, symmetric = TRUE, only.values = TRUE)$values))
     logretval <- -(ncol(x) * log(2 * pi) + logdet + distval)/2