R version 3.1.2 (2014-10-31) -- "Pumpkin Helmet" Copyright (C) 2014 The R Foundation for Statistical Computing Platform: x86_64-pc-linux-gnu (64-bit) 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. > # $Id: test_FOCUS_2006_A.R 59 2010-07-28 12:29:15Z jranke $ > > # Copyright (C) 2008-2010 Johannes Ranke > # Contact: mkin-devel@lists.berlios.de > > # This file is part of the R package kinfit > > # kinfit is free software: you can redistribute it and/or modify it under the > # terms of the GNU General Public License as published by the Free Software > # Foundation, either version 3 of the License, or (at your option) any later > # version. > > # This program is distributed in the hope that it will be useful, but WITHOUT > # ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS > # FOR A PARTICULAR PURPOSE. See the GNU General Public License for more > # details. > > # You should have received a copy of the GNU General Public License along with > # this program. If not, see > > library(kinfit) > data(FOCUS_2006_A) > fits <- kinfit(FOCUS_2006_A, kinmodels = c("SFO", "HS")) > print(kinresults(fits)$results, digits=5) DT50 DT90 SFO 18.624 61.868 > print(kinresults(fits)$stats, digits=5) n.times df mean.means RSS err.min R2 SFO 8 6 50.054 221.81 0.083848 0.9845 > > data(FOCUS_2006_B) > fits <- kinfit(FOCUS_2006_B, kinmodels = c("SFO", "FOMC", "DFOP")) > print(kinresults(fits)$results, digits=5) DT50 DT90 SFO 8.8686 29.461 FOMC 8.6834 30.754 DFOP 8.6829 30.789 > print(kinresults(fits)$stats, digits=5) n.times df mean.means RSS err.min R2 SFO 8 6 35.015 30.656 0.044555 0.99713 FOMC 8 5 35.015 28.583 0.045886 0.99733 DFOP 8 4 35.015 28.550 0.049527 0.99733 > > data(FOCUS_2006_C) > fits <- kinfit(FOCUS_2006_C, kinmodels = c("SFO", "FOMC", "DFOP")) > print(kinresults(fits)$results, digits=5) DT50 DT90 SFO 2.2647 7.5233 FOMC 1.7852 15.1479 DFOP 1.8869 21.2507 > print(kinresults(fits)$stats, digits=5) n.times df mean.means RSS err.min R2 SFO 9 7 23.589 196.5334 0.158440 0.97136 FOMC 9 6 23.589 31.0509 0.066568 0.99548 DFOP 9 5 23.589 4.3627 0.026621 0.99936 > > data(FOCUS_2006_D) > fits <- kinfit(FOCUS_2006_D, kinmodels = c("SFO", "FOMC")) > print(kinresults(fits)$results, digits=5) DT50 DT90 SFO 7.0776 23.511 FOMC 6.9350 24.044 > print(kinresults(fits)$stats, digits=3) n.times df mean.means RSS err.min R2 SFO 9 7 39.5 208 0.0645 0.992 FOMC 9 6 39.5 205 0.0678 0.992 > > data(FOCUS_2006_E) > fits <- kinfit(FOCUS_2006_E, kinmodels = c("SFO", "FOMC", "DFOP")) > print(kinresults(fits)$results, digits=5) DT50 DT90 SFO 2.2647 7.5233 FOMC 1.7852 15.1479 DFOP 1.8869 21.2507 > print(kinresults(fits)$stats, digits=5) n.times df mean.means RSS err.min R2 SFO 9 7 23.589 196.5334 0.158440 0.97136 FOMC 9 6 23.589 31.0509 0.066568 0.99548 DFOP 9 5 23.589 4.3627 0.026621 0.99936 > > proc.time() user system elapsed 0.992 0.064 1.113