[Pomp-commits] r971 - pkg/pomp/R

noreply at r-forge.r-project.org noreply at r-forge.r-project.org
Tue Jun 10 15:21:48 CEST 2014


Author: kingaa
Date: 2014-06-10 15:21:48 +0200 (Tue, 10 Jun 2014)
New Revision: 971

Modified:
   pkg/pomp/R/abc-methods.R
   pkg/pomp/R/pmcmc-methods.R
Log:
- minor tweaks to diagnostic plotting routines

Modified: pkg/pomp/R/abc-methods.R
===================================================================
--- pkg/pomp/R/abc-methods.R	2014-06-09 16:32:00 UTC (rev 970)
+++ pkg/pomp/R/abc-methods.R	2014-06-10 13:21:48 UTC (rev 971)
@@ -135,7 +135,6 @@
     oma.multi <- c(6,0,5,0)
     xx <- z[[1]]
     estnames <- pars
-    
     ## plot abc convergence diagnostics
     other.diagnostics <- c()
     plotnames <- c(other.diagnostics,estnames)
@@ -152,7 +151,7 @@
       hi <- min(low+n.per.page-1,nplots)
       for (i in seq(from=low,to=hi,by=1)) {
         n <- i-low+1
-        dat <- sapply(z,function(po,label) conv.rec(po,label),label=plotnames[i])
+        dat <- sapply(z,conv.rec,pars=plotnames[i])
         matplot(
                 y=dat, 
                 x=iteration,

Modified: pkg/pomp/R/pmcmc-methods.R
===================================================================
--- pkg/pomp/R/pmcmc-methods.R	2014-06-09 16:32:00 UTC (rev 970)
+++ pkg/pomp/R/pmcmc-methods.R	2014-06-10 13:21:48 UTC (rev 971)
@@ -135,7 +135,7 @@
     hi <- min(low+n.per.page-1,nplots)
     for (i in seq(from=low,to=hi,by=1)) {
       n <- i-low+1
-      dat <- sapply(z,function(po,label) conv.rec(po,label),label=plotnames[i])
+      dat <- sapply(z,conv.rec,pars=plotnames[i])
       matplot(
               y=dat, 
               x=iteration,



More information about the pomp-commits mailing list