[Pomp-commits] r47 - pkg/tests

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

Author: kingaa
Date: 2008-08-15 12:47:51 +0200 (Fri, 15 Aug 2008)
New Revision: 47

Modified:
   pkg/tests/ou2-mif.Rout.save
Log:
fix bug in Euler step

Modified: pkg/tests/ou2-mif.Rout.save
===================================================================
--- pkg/tests/ou2-mif.Rout.save	2008-08-15 10:33:57 UTC (rev 46)
+++ pkg/tests/ou2-mif.Rout.save	2008-08-15 10:47:51 UTC (rev 47)
@@ -16,6 +16,7 @@
 Type 'q()' to quit R.
 
 > library(pomp)
+Loading required package: odesolve
 > 
 > data(ou2)
 > 
@@ -67,11 +68,11 @@
 > mif.fit <- continue(mif.fit,Nmif=70,max.fail=100)
 > toc <- Sys.time()
 > print(toc-tic)
-Time difference of 1.542085 mins
+Time difference of 1.396902 mins
 > cat("PF estimated log likelihood at MIF MLE\n")
 PF estimated log likelihood at MIF MLE
 > print(pfilter(mif.fit)$loglik,digits=4)
-[1] -422.6
+[1] -421.7
 > 
 > cat("coefficients at truth\n")
 coefficients at truth
@@ -82,7 +83,7 @@
 MIF MLE
 > print(coef(mif.fit,c('x1.0','x2.0','alpha.1','alpha.4')),digits=4)
     x1.0     x2.0  alpha.1  alpha.4 
- 51.1941 -51.2739   0.8990   0.9736 
+ 51.2028 -50.0965   0.9032   0.9738 
 > 
 > plot(mif.fit)
 > compare.mif(mif.fit)