[Pomp-commits] r57 - pkg/tests

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

Author: kingaa
Date: 2009-01-17 02:05:33 +0100 (Sat, 17 Jan 2009)
New Revision: 57

Removed:
   pkg/tests/ou2-nlf.Rout
Modified:
   pkg/tests/ou2-mif.Rout.save
Log:


Modified: pkg/tests/ou2-mif.Rout.save
===================================================================
--- pkg/tests/ou2-mif.Rout.save	2009-01-17 01:04:24 UTC (rev 56)
+++ pkg/tests/ou2-mif.Rout.save	2009-01-17 01:05:33 UTC (rev 57)
@@ -1,5 +1,5 @@
 
-R version 2.7.1 (2008-06-23)
+R version 2.7.2 (2008-08-25)
 Copyright (C) 2008 The R Foundation for Statistical Computing
 ISBN 3-900051-07-0
 
@@ -68,11 +68,11 @@
 > mif.fit <- continue(mif.fit,Nmif=70,max.fail=100)
 > toc <- Sys.time()
 > print(toc-tic)
-Time difference of 1.694885 mins
+Time difference of 1.222864 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
@@ -83,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)

Deleted: pkg/tests/ou2-nlf.Rout
===================================================================
--- pkg/tests/ou2-nlf.Rout	2009-01-17 01:04:24 UTC (rev 56)
+++ pkg/tests/ou2-nlf.Rout	2009-01-17 01:05:33 UTC (rev 57)
@@ -1,71 +0,0 @@
-
-R version 2.7.2 (2008-08-25)
-Copyright (C) 2008 The R Foundation for Statistical Computing
-ISBN 3-900051-07-0
-
-R is free software and comes with ABSOLUTELY NO WARRANTY.
-You are welcome to redistribute it under certain conditions.
-Type 'license()' or 'licence()' for distribution details.
-
-R is a collaborative project with many contributors.
-Type 'contributors()' for more information and
-'citation()' on how to cite R or R packages in publications.
-
-Type 'demo()' for some demos, 'help()' for on-line help, or
-'help.start()' for an HTML browser interface to help.
-Type 'q()' to quit R.
-
-> library(pomp.devel)
-Loading required package: subplex
-Loading required package: odesolve
-Loading required package: plotrix
-Loading required package: lattice
-> 
-> data(ou2)
-> 
-> set.seed(1066)
-> 
-> po <- ou2
-> coef(po,c("x1.0","x2.0","alpha.1","alpha.4")) <- c(0,0,0.1,0.2)
-> po <- simulate(po,times=(1:10000))[[1]]
-> p.truth <- coef(po)
-> guess <- p.truth
-> ## guess[c('x1.0','x2.0','alpha.1','alpha.4')] <- c(45,-60,0.8,0.9)
-> 
-> m1 <- nlf(
-+           object=po,
-+           start=guess,
-+           est=c("alpha.1","alpha.4"),
-+           lags=c(1,2),
-+           nconverge=100,
-+           nasymp=2500,
-+           method="Nelder-Mead",
-+           maxit=500,
-+           trace=1,
-+           verbose=TRUE,
-+           lql.frac = 0.025
-+           )
-  Nelder-Mead direct search function minimizer
-function value for initial parameters = 40059.679966
-  Scaled convergence tolerance is 0.000596936
-Stepsize computed as 0.020000
-BUILD              3 40059.679966 40059.013186
-REFLECTION         5 40059.024962 40058.355658
-HI-REDUCTION       7 40059.013186 40058.355658
-HI-REDUCTION       9 40058.480960 40058.355658
-REFLECTION        11 40058.428912 40058.287812
-HI-REDUCTION      13 40058.355658 40058.283261
-HI-REDUCTION      15 40058.287812 40058.270086
-HI-REDUCTION      17 40058.283261 40058.262895
-HI-REDUCTION      19 40058.270086 40058.262895
-HI-REDUCTION      21 40058.263735 40058.260296
-HI-REDUCTION      23 40058.262895 40058.260296
-Exiting from Nelder Mead minimizer
-    25 function evaluations used
-h in NLF =  0.1 
-epsilon in NLF = 0.4420236 0.2750263 
-Fitted param  1 -4.042142 -4.042142  up in  ‘nlf’ 
-Fitted param  1 -4.031832  down in  ‘NLF’ 
-Fitted param  2 -4.032517 -4.032517  up in  ‘nlf’ 
-Fitted param  2 -4.03145  down in  ‘NLF’ 
->