1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
|
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 <http://www.gnu.org/licenses/>
>
> 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
|
