[Pomp-commits] r56 - pkg/tests

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

Author: kingaa
Date: 2009-01-17 02:04:24 +0100 (Sat, 17 Jan 2009)
New Revision: 56

Added:
   pkg/tests/ou2-nlf.Rout.save
Modified:
   pkg/tests/ou2-nlf.R
Log:


Modified: pkg/tests/ou2-nlf.R
===================================================================
--- pkg/tests/ou2-nlf.R	2009-01-16 22:31:02 UTC (rev 55)
+++ pkg/tests/ou2-nlf.R	2009-01-17 01:04:24 UTC (rev 56)
@@ -1,4 +1,4 @@
-library(pomp.devel)
+library(pomp)
 
 data(ou2)
 

Added: pkg/tests/ou2-nlf.Rout.save
===================================================================
--- pkg/tests/ou2-nlf.Rout.save	                        (rev 0)
+++ pkg/tests/ou2-nlf.Rout.save	2009-01-17 01:04:24 UTC (rev 56)
@@ -0,0 +1,71 @@
+
+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’ 
+>