[Rcpp-commits] r3523 - in pkg/RcppSMC: . R

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

Author: adamj
Date: 2012-03-18 12:30:21 +0100 (Sun, 18 Mar 2012)
New Revision: 3523

Modified:
   pkg/RcppSMC/ChangeLog
   pkg/RcppSMC/R/bspfGaussianOptimal.R
Log:
Plotting tweak.


Modified: pkg/RcppSMC/ChangeLog
===================================================================
--- pkg/RcppSMC/ChangeLog	2012-03-18 11:16:38 UTC (rev 3522)
+++ pkg/RcppSMC/ChangeLog	2012-03-18 11:30:21 UTC (rev 3523)
@@ -1,3 +1,8 @@
+2012-03-18 Adam Johansen <a.m.johansen at warwick.ac.uk>
+	* src/bspfGaussianOptimal.R tweaked plotting.
+	* src/pfNonlinBS.R fixed plotting bug.
+	* deprecated moved all rareEvents files here.
+
 2012-03-17  Adam Johansen <a.m.johansen at warwick.ac.uk>
 
 	* man/pfEx.Rd updated to reflect previous changes.

Modified: pkg/RcppSMC/R/bspfGaussianOptimal.R
===================================================================
--- pkg/RcppSMC/R/bspfGaussianOptimal.R	2012-03-18 11:16:38 UTC (rev 3522)
+++ pkg/RcppSMC/R/bspfGaussianOptimal.R	2012-03-18 11:30:21 UTC (rev 3523)
@@ -13,7 +13,7 @@
 	sqvect = t(res$weight) %*% res$values^2 / sum(res$weight);
 	sdvect = sqrt(sqvect - mvect^2);
  
-       plot(time, mvect, col='dark red', 'l', lty = 1, xlab = 'Iteration',
+       plot(time, mvect, col='dark red', 'l', lty = 1, lwd=3, xlab = 'Iteration',
        ylab='State', main='Mean and 1, 2 standard deviation credible
        intervals with observations', xlim = c(0,length(data)), ylim=c(min(mvect - 2.1
        * (sdvect)), max(mvect+2.1*sdvect))
@@ -22,7 +22,7 @@
        lines(time, mvect - sdvect, lty=3, col='dark blue')
        lines(time, mvect + 2 * sdvect, lty=2, col='dark blue')
        lines(time, mvect - 2 * sdvect, lty=2, col='dark blue')
-       points(time, data, col = 'dark green')
+       points(time, data, col = 'dark green', cex=0.5)
     }
 
     invisible(res)