diff options
| author | Johannes Ranke <jranke@uni-bremen.de> | 2017-07-21 14:42:14 +0200 |
|---|---|---|
| committer | Johannes Ranke <jranke@uni-bremen.de> | 2017-07-21 14:42:14 +0200 |
| commit | 4a918da6d5f971335b74b0fc83cb08f5c3163f95 (patch) | |
| tree | ce67e7a1b2f4763342bee6f4ef9af283b74a27c0 /docs/articles/FOCUS_Z.R | |
| parent | 4f73b627fefd78405d9d996ca096cd9dd2b958c9 (diff) | |
Rename twa to max_twa_parent, update docs
Diffstat (limited to 'docs/articles/FOCUS_Z.R')
| -rw-r--r-- | docs/articles/FOCUS_Z.R | 17 |
1 files changed, 17 insertions, 0 deletions
diff --git a/docs/articles/FOCUS_Z.R b/docs/articles/FOCUS_Z.R index 0d5f0940..5c70b57e 100644 --- a/docs/articles/FOCUS_Z.R +++ b/docs/articles/FOCUS_Z.R @@ -3,6 +3,23 @@ require(knitr) opts_chunk$set(engine='R', tidy = FALSE, cache = TRUE) options(width=70) +## ----FOCUS_2006_Z_data, echo=TRUE, eval=TRUE------------------------ +require(mkin) +LOD = 0.5 +FOCUS_2006_Z = data.frame( + t = c(0, 0.04, 0.125, 0.29, 0.54, 1, 2, 3, 4, 7, 10, 14, 21, + 42, 61, 96, 124), + Z0 = c(100, 81.7, 70.4, 51.1, 41.2, 6.6, 4.6, 3.9, 4.6, 4.3, 6.8, + 2.9, 3.5, 5.3, 4.4, 1.2, 0.7), + Z1 = c(0, 18.3, 29.6, 46.3, 55.1, 65.7, 39.1, 36, 15.3, 5.6, 1.1, + 1.6, 0.6, 0.5 * LOD, NA, NA, NA), + Z2 = c(0, NA, 0.5 * LOD, 2.6, 3.8, 15.3, 37.2, 31.7, 35.6, 14.5, + 0.8, 2.1, 1.9, 0.5 * LOD, NA, NA, NA), + Z3 = c(0, NA, NA, NA, NA, 0.5 * LOD, 9.2, 13.1, 22.3, 28.4, 32.5, + 25.2, 17.2, 4.8, 4.5, 2.8, 4.4)) + +FOCUS_2006_Z_mkin <- mkin_wide_to_long(FOCUS_2006_Z) + ## ----FOCUS_2006_Z_fits_1, echo=TRUE, fig.height=6------------------- Z.2a <- mkinmod(Z0 = mkinsub("SFO", "Z1"), Z1 = mkinsub("SFO")) |