From c73b2f30ec836c949885784ab576e814eb8070a9 Mon Sep 17 00:00:00 2001 From: Johannes Ranke Date: Tue, 9 Mar 2021 17:35:47 +0100 Subject: Some improvements for borderline cases - fit_with_errors for saem() - test_log_parms for mean_degparms() and saem() --- docs/dev/reference/Rplot001.png | Bin 13995 -> 1011 bytes docs/dev/reference/Rplot002.png | Bin 13648 -> 16859 bytes docs/dev/reference/Rplot003.png | Bin 28745 -> 28844 bytes docs/dev/reference/Rplot004.png | Bin 49269 -> 49360 bytes docs/dev/reference/Rplot005.png | Bin 59143 -> 59216 bytes docs/dev/reference/endpoints.html | 2 +- docs/dev/reference/index.html | 2 +- docs/dev/reference/nlme-1.png | Bin 70133 -> 68233 bytes docs/dev/reference/nlme-2.png | Bin 94031 -> 90552 bytes docs/dev/reference/nlme.html | 33 ++- docs/dev/reference/nlme.mmkin-1.png | Bin 124677 -> 124827 bytes docs/dev/reference/nlme.mmkin-2.png | Bin 169523 -> 169698 bytes docs/dev/reference/nlme.mmkin-3.png | Bin 172692 -> 172809 bytes docs/dev/reference/nlme.mmkin.html | 2 +- docs/dev/reference/plot.mixed.mmkin-1.png | Bin 84734 -> 85433 bytes docs/dev/reference/plot.mixed.mmkin-2.png | Bin 173916 -> 174061 bytes docs/dev/reference/plot.mixed.mmkin-3.png | Bin 172396 -> 172540 bytes docs/dev/reference/plot.mixed.mmkin-4.png | Bin 175502 -> 175594 bytes docs/dev/reference/plot.mixed.mmkin.html | 6 +- docs/dev/reference/saem-1.png | Bin 47315 -> 47342 bytes docs/dev/reference/saem-2.png | Bin 48720 -> 48819 bytes docs/dev/reference/saem-3.png | Bin 82107 -> 82202 bytes docs/dev/reference/saem-4.png | Bin 128231 -> 128213 bytes docs/dev/reference/saem-5.png | Bin 173288 -> 173665 bytes docs/dev/reference/saem.html | 72 +++-- docs/dev/reference/summary.saem.mmkin.html | 422 ++++++++++++++--------------- 26 files changed, 288 insertions(+), 251 deletions(-) (limited to 'docs/dev/reference') diff --git a/docs/dev/reference/Rplot001.png b/docs/dev/reference/Rplot001.png index 7f498242..17a35806 100644 Binary files a/docs/dev/reference/Rplot001.png and b/docs/dev/reference/Rplot001.png differ diff --git a/docs/dev/reference/Rplot002.png b/docs/dev/reference/Rplot002.png index 54c31a3f..a9a972e5 100644 Binary files a/docs/dev/reference/Rplot002.png and b/docs/dev/reference/Rplot002.png differ diff --git a/docs/dev/reference/Rplot003.png b/docs/dev/reference/Rplot003.png index 2b011ec1..d077f01c 100644 Binary files a/docs/dev/reference/Rplot003.png and b/docs/dev/reference/Rplot003.png differ diff --git a/docs/dev/reference/Rplot004.png b/docs/dev/reference/Rplot004.png index 98dd019e..ffcd2d96 100644 Binary files a/docs/dev/reference/Rplot004.png and b/docs/dev/reference/Rplot004.png differ diff --git a/docs/dev/reference/Rplot005.png b/docs/dev/reference/Rplot005.png index 8c91d61e..dfb5965b 100644 Binary files a/docs/dev/reference/Rplot005.png and b/docs/dev/reference/Rplot005.png differ diff --git a/docs/dev/reference/endpoints.html b/docs/dev/reference/endpoints.html index c9912f9c..63bec6a8 100644 --- a/docs/dev/reference/endpoints.html +++ b/docs/dev/reference/endpoints.html @@ -78,7 +78,7 @@ advantage that the SFORB model can also be used for metabolites." /> mkin - 1.0.3.9000 + 1.0.4.9000 diff --git a/docs/dev/reference/index.html b/docs/dev/reference/index.html index 03a21517..5533a01f 100644 --- a/docs/dev/reference/index.html +++ b/docs/dev/reference/index.html @@ -71,7 +71,7 @@ mkin - 1.0.3.9000 + 1.0.4.9000 diff --git a/docs/dev/reference/nlme-1.png b/docs/dev/reference/nlme-1.png index 728cc557..fd68ae43 100644 Binary files a/docs/dev/reference/nlme-1.png and b/docs/dev/reference/nlme-1.png differ diff --git a/docs/dev/reference/nlme-2.png b/docs/dev/reference/nlme-2.png index e8167455..853cae40 100644 Binary files a/docs/dev/reference/nlme-2.png and b/docs/dev/reference/nlme-2.png differ diff --git a/docs/dev/reference/nlme.html b/docs/dev/reference/nlme.html index b850eb3d..78d132e9 100644 --- a/docs/dev/reference/nlme.html +++ b/docs/dev/reference/nlme.html @@ -75,7 +75,7 @@ datasets. They are used internally by the nlme.mmkin() method." /> mkin - 1.0.3.9000 + 1.0.4.9000 @@ -155,7 +155,7 @@ datasets. They are used internally by the nlme.m 
nlme_function(object)
 
-mean_degparms(object, random = FALSE)
+mean_degparms(object, random = FALSE, test_log_parms = FALSE, conf.level = 0.6)
 
 nlme_data(object)
@@ -170,6 +170,17 @@ datasets. They are used internally by the nlme.m random

Should a list with fixed and random effects be returned?

+ + test_log_parms +

If TRUE, log parameters are only considered in +the mean calculations if their untransformed counterparts (most likely +rate constants) pass the t-test for significant difference from zero.

+ + + conf.level +

Possibility to adjust the required confidence level +for parameter that are tested if requested by 'test_log_parms'.

+

Value

@@ -211,7 +222,7 @@ nlme for the case of a single grouping variable ds.

nlme_f <- nlme_function(f) # These assignments are necessary for these objects to be # visible to nlme and augPred when evaluation is done by -# pkgdown to generated the html docs. +# pkgdown to generate the html docs. assign("nlme_f", nlme_f, globalenv()) assign("grouped_data", grouped_data, globalenv()) @@ -226,28 +237,28 @@ nlme for the case of a single grouping variable ds.

#> Model: value ~ nlme_f(name, time, parent_0, log_k_parent_sink) #> Data: grouped_data #> AIC BIC logLik -#> 300.6824 310.2426 -145.3412 +#> 298.2781 307.7372 -144.1391 #> #> Random effects: #> Formula: list(parent_0 ~ 1, log_k_parent_sink ~ 1) #> Level: ds #> Structure: Diagonal #> parent_0 log_k_parent_sink Residual -#> StdDev: 1.697361 0.6801209 3.666073 +#> StdDev: 0.937473 0.7098105 3.83543 #> #> Fixed effects: parent_0 + log_k_parent_sink ~ 1 #> Value Std.Error DF t-value p-value -#> parent_0 100.99378 1.3890416 46 72.70753 0 -#> log_k_parent_sink -3.07521 0.4018589 46 -7.65246 0 +#> parent_0 101.76838 1.1445443 45 88.91607 0 +#> log_k_parent_sink -3.05444 0.4195622 45 -7.28008 0 #> Correlation: #> prnt_0 -#> log_k_parent_sink 0.027 +#> log_k_parent_sink 0.034 #> #> Standardized Within-Group Residuals: -#> Min Q1 Med Q3 Max -#> -1.9942823 -0.5622565 0.1791579 0.7165038 2.0704781 +#> Min Q1 Med Q3 Max +#> -2.61693595 -0.21853231 0.05740682 0.57209372 3.04598764 #> -#> Number of Observations: 50 +#> Number of Observations: 49 #> Number of Groups: 3
plot(augPred(m_nlme, level = 0:1), layout = c(3, 1))
# augPred does not work on fits with more than one state # variable diff --git a/docs/dev/reference/nlme.mmkin-1.png b/docs/dev/reference/nlme.mmkin-1.png index 9186c135..90ede880 100644 Binary files a/docs/dev/reference/nlme.mmkin-1.png and b/docs/dev/reference/nlme.mmkin-1.png differ diff --git a/docs/dev/reference/nlme.mmkin-2.png b/docs/dev/reference/nlme.mmkin-2.png index d395fe02..0d140fd1 100644 Binary files a/docs/dev/reference/nlme.mmkin-2.png and b/docs/dev/reference/nlme.mmkin-2.png differ diff --git a/docs/dev/reference/nlme.mmkin-3.png b/docs/dev/reference/nlme.mmkin-3.png index 40518a59..8a60b52b 100644 Binary files a/docs/dev/reference/nlme.mmkin-3.png and b/docs/dev/reference/nlme.mmkin-3.png differ diff --git a/docs/dev/reference/nlme.mmkin.html b/docs/dev/reference/nlme.mmkin.html index 925cf7cf..f308d8b7 100644 --- a/docs/dev/reference/nlme.mmkin.html +++ b/docs/dev/reference/nlme.mmkin.html @@ -74,7 +74,7 @@ have been obtained by fitting the same model to a list of datasets." /> mkin - 1.0.3.9000 + 1.0.4.9000
diff --git a/docs/dev/reference/plot.mixed.mmkin-1.png b/docs/dev/reference/plot.mixed.mmkin-1.png index 9c9a2211..2224d96e 100644 Binary files a/docs/dev/reference/plot.mixed.mmkin-1.png and b/docs/dev/reference/plot.mixed.mmkin-1.png differ diff --git a/docs/dev/reference/plot.mixed.mmkin-2.png b/docs/dev/reference/plot.mixed.mmkin-2.png index 0f66ff0f..28168495 100644 Binary files a/docs/dev/reference/plot.mixed.mmkin-2.png and b/docs/dev/reference/plot.mixed.mmkin-2.png differ diff --git a/docs/dev/reference/plot.mixed.mmkin-3.png b/docs/dev/reference/plot.mixed.mmkin-3.png index 34212f1c..d18275dd 100644 Binary files a/docs/dev/reference/plot.mixed.mmkin-3.png and b/docs/dev/reference/plot.mixed.mmkin-3.png differ diff --git a/docs/dev/reference/plot.mixed.mmkin-4.png b/docs/dev/reference/plot.mixed.mmkin-4.png index c1450d24..2fd52425 100644 Binary files a/docs/dev/reference/plot.mixed.mmkin-4.png and b/docs/dev/reference/plot.mixed.mmkin-4.png differ diff --git a/docs/dev/reference/plot.mixed.mmkin.html b/docs/dev/reference/plot.mixed.mmkin.html index 630e95a3..36796580 100644 --- a/docs/dev/reference/plot.mixed.mmkin.html +++ b/docs/dev/reference/plot.mixed.mmkin.html @@ -72,7 +72,7 @@ mkin - 1.0.3.9000 + 1.0.4.9000 @@ -283,10 +283,10 @@ corresponding model prediction lines for the different datasets.

f_saem <- saem(f, transformations = "saemix")
#> Running main SAEM algorithm -#> [1] "Mon Feb 15 17:12:17 2021" +#> [1] "Tue Mar 9 17:34:35 2021" #> .... #> Minimisation finished -#> [1] "Mon Feb 15 17:12:24 2021"
plot(f_saem) +#> [1] "Tue Mar 9 17:34:42 2021"
plot(f_saem)
# We can overlay the two variants if we generate predictions pred_nlme <- mkinpredict(dfop_sfo, diff --git a/docs/dev/reference/saem-1.png b/docs/dev/reference/saem-1.png index 2df248bb..0da31388 100644 Binary files a/docs/dev/reference/saem-1.png and b/docs/dev/reference/saem-1.png differ diff --git a/docs/dev/reference/saem-2.png b/docs/dev/reference/saem-2.png index d4a2c1be..010950ba 100644 Binary files a/docs/dev/reference/saem-2.png and b/docs/dev/reference/saem-2.png differ diff --git a/docs/dev/reference/saem-3.png b/docs/dev/reference/saem-3.png index 4474b1f1..829f22bf 100644 Binary files a/docs/dev/reference/saem-3.png and b/docs/dev/reference/saem-3.png differ diff --git a/docs/dev/reference/saem-4.png b/docs/dev/reference/saem-4.png index bf24d6b0..4e976fa2 100644 Binary files a/docs/dev/reference/saem-4.png and b/docs/dev/reference/saem-4.png differ diff --git a/docs/dev/reference/saem-5.png b/docs/dev/reference/saem-5.png index 27ed3f8f..f50969b4 100644 Binary files a/docs/dev/reference/saem-5.png and b/docs/dev/reference/saem-5.png differ diff --git a/docs/dev/reference/saem.html b/docs/dev/reference/saem.html index bdb1226e..23102df3 100644 --- a/docs/dev/reference/saem.html +++ b/docs/dev/reference/saem.html @@ -74,7 +74,7 @@ Expectation Maximisation algorithm (SAEM)." /> mkin - 1.0.3.9000 + 1.0.4.9000
@@ -158,9 +158,12 @@ Expectation Maximisation algorithm (SAEM).

object, transformations = c("mkin", "saemix"), degparms_start = numeric(), + test_log_parms = FALSE, + conf.level = 0.6, solution_type = "auto", control = list(displayProgress = FALSE, print = FALSE, save = FALSE, save.graphs = FALSE), + fail_with_errors = TRUE, verbose = FALSE, quiet = FALSE, ... @@ -174,6 +177,7 @@ Expectation Maximisation algorithm (SAEM).

solution_type = "auto", transformations = c("mkin", "saemix"), degparms_start = numeric(), + test_log_parms = FALSE, verbose = FALSE, ... ) @@ -204,6 +208,18 @@ SFO or DFOP is used for the parent and there is either no metabolite or one.

degparms_start

Parameter values given as a named numeric vector will be used to override the starting values obtained from the 'mmkin' object.

+ + + test_log_parms +

If TRUE, an attempt is made to use more robust starting +values for population parameters fitted as log parameters in mkin (like +rate constants) by only considering rate constants that pass the t-test +when calculating mean degradation parameters using mean_degparms.

+ + + conf.level +

Possibility to adjust the required confidence level +for parameter that are tested if requested by 'test_log_parms'.

solution_type @@ -214,6 +230,11 @@ automatic choice is not desired

control

Passed to saemix::saemix

+ + fail_with_errors +

Should a failure to compute standard errors +from the inverse of the Fisher Information Matrix be a failure?

+ verbose

Should we print information about created objects of @@ -261,33 +282,36 @@ using mmkin.

state.ini = c(parent = 100), fixed_initials = "parent", quiet = TRUE) f_saem_p0_fixed <- saem(f_mmkin_parent_p0_fixed)
#> Running main SAEM algorithm -#> [1] "Mon Feb 15 17:12:32 2021" +#> [1] "Tue Mar 9 17:34:44 2021" #> .... #> Minimisation finished -#> [1] "Mon Feb 15 17:12:34 2021"
+#> [1] "Tue Mar 9 17:34:45 2021"
f_mmkin_parent <- mmkin(c("SFO", "FOMC", "DFOP"), ds, quiet = TRUE) f_saem_sfo <- saem(f_mmkin_parent["SFO", ])
#> Running main SAEM algorithm -#> [1] "Mon Feb 15 17:12:35 2021" +#> [1] "Tue Mar 9 17:34:46 2021" #> .... #> Minimisation finished -#> [1] "Mon Feb 15 17:12:36 2021"
f_saem_fomc <- saem(f_mmkin_parent["FOMC", ]) +#> [1] "Tue Mar 9 17:34:48 2021"
f_saem_fomc <- saem(f_mmkin_parent["FOMC", ])
#> Running main SAEM algorithm -#> [1] "Mon Feb 15 17:12:36 2021" +#> [1] "Tue Mar 9 17:34:48 2021" #> .... #> Minimisation finished -#> [1] "Mon Feb 15 17:12:38 2021"
f_saem_dfop <- saem(f_mmkin_parent["DFOP", ]) +#> [1] "Tue Mar 9 17:34:50 2021"
f_saem_dfop <- saem(f_mmkin_parent["DFOP", ])
#> Running main SAEM algorithm -#> [1] "Mon Feb 15 17:12:39 2021" +#> [1] "Tue Mar 9 17:34:51 2021" #> .... #> Minimisation finished -#> [1] "Mon Feb 15 17:12:42 2021"
+#> [1] "Tue Mar 9 17:34:53 2021"
# The returned saem.mmkin object contains an SaemixObject, therefore we can use # functions from saemix library(saemix)
#> Package saemix, version 3.1.9000 -#> please direct bugs, questions and feedback to emmanuelle.comets@inserm.fr
compare.saemix(list(f_saem_sfo$so, f_saem_fomc$so, f_saem_dfop$so)) -
#> Error in compare.saemix(list(f_saem_sfo$so, f_saem_fomc$so, f_saem_dfop$so)): 'compare.saemix' requires at least two models.
plot(f_saem_fomc$so, plot.type = "convergence") +#> please direct bugs, questions and feedback to emmanuelle.comets@inserm.fr
compare.saemix(f_saem_sfo$so, f_saem_fomc$so, f_saem_dfop$so) +
#> Likelihoods calculated by importance sampling
#> AIC BIC +#> 1 624.2484 622.2956 +#> 2 467.7096 464.9757 +#> 3 495.4373 491.9222
plot(f_saem_fomc$so, plot.type = "convergence")
#> Plotting convergence plots
plot(f_saem_fomc$so, plot.type = "individual.fit")
#> Plotting individual fits
plot(f_saem_fomc$so, plot.type = "npde")
#> Simulating data using nsim = 1000 simulated datasets @@ -324,11 +348,13 @@ using mmkin.

f_mmkin_parent_tc <- update(f_mmkin_parent, error_model = "tc") f_saem_fomc_tc <- saem(f_mmkin_parent_tc["FOMC", ])
#> Running main SAEM algorithm -#> [1] "Mon Feb 15 17:12:44 2021" +#> [1] "Tue Mar 9 17:34:55 2021" #> .... #> Minimisation finished -#> [1] "Mon Feb 15 17:12:49 2021"
compare.saemix(list(f_saem_fomc$so, f_saem_fomc_tc$so)) -
#> Error in compare.saemix(list(f_saem_fomc$so, f_saem_fomc_tc$so)): 'compare.saemix' requires at least two models.
+#> [1] "Tue Mar 9 17:35:00 2021"
compare.saemix(f_saem_fomc$so, f_saem_fomc_tc$so) +
#> Likelihoods calculated by importance sampling
#> AIC BIC +#> 1 467.7096 464.9757 +#> 2 469.6831 466.5586
sfo_sfo <- mkinmod(parent = mkinsub("SFO", "A1"), A1 = mkinsub("SFO"))
#> Temporary DLL for differentials generated and loaded
fomc_sfo <- mkinmod(parent = mkinsub("FOMC", "A1"), @@ -346,15 +372,15 @@ using mmkin.

# four minutes f_saem_sfo_sfo <- saem(f_mmkin["SFO-SFO", ])
#> Running main SAEM algorithm -#> [1] "Mon Feb 15 17:12:51 2021" +#> [1] "Tue Mar 9 17:35:02 2021" #> .... #> Minimisation finished -#> [1] "Mon Feb 15 17:12:56 2021"
f_saem_dfop_sfo <- saem(f_mmkin["DFOP-SFO", ]) +#> [1] "Tue Mar 9 17:35:07 2021"
f_saem_dfop_sfo <- saem(f_mmkin["DFOP-SFO", ])
#> Running main SAEM algorithm -#> [1] "Mon Feb 15 17:12:56 2021" +#> [1] "Tue Mar 9 17:35:07 2021" #> .... #> Minimisation finished -#> [1] "Mon Feb 15 17:13:05 2021"
# We can use print, plot and summary methods to check the results +#> [1] "Tue Mar 9 17:35:15 2021"
# We can use print, plot and summary methods to check the results print(f_saem_dfop_sfo)
#> Kinetic nonlinear mixed-effects model fit by SAEM #> Structural model: @@ -395,10 +421,10 @@ using mmkin.

#> SD.g_qlogis 0.44771 -0.86417 1.7596
plot(f_saem_dfop_sfo)
summary(f_saem_dfop_sfo, data = TRUE)
#> saemix version used for fitting: 3.1.9000 -#> mkin version used for pre-fitting: 1.0.3.9000 -#> R version used for fitting: 4.0.3 -#> Date of fit: Mon Feb 15 17:13:05 2021 -#> Date of summary: Mon Feb 15 17:13:06 2021 +#> mkin version used for pre-fitting: 1.0.4.9000 +#> R version used for fitting: 4.0.4 +#> Date of fit: Tue Mar 9 17:35:16 2021 +#> Date of summary: Tue Mar 9 17:35:16 2021 #> #> Equations: #> d_parent/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 * @@ -413,7 +439,7 @@ using mmkin.

#> #> Model predictions using solution type analytical #> -#> Fitted in 8.985 s using 300, 100 iterations +#> Fitted in 8.668 s using 300, 100 iterations #> #> Variance model: Constant variance #> diff --git a/docs/dev/reference/summary.saem.mmkin.html b/docs/dev/reference/summary.saem.mmkin.html index 0d661ee9..1166abb1 100644 --- a/docs/dev/reference/summary.saem.mmkin.html +++ b/docs/dev/reference/summary.saem.mmkin.html @@ -76,7 +76,7 @@ endpoints such as formation fractions and DT50 values. Optionally mkin - 1.0.3.9000 + 1.0.4.9000
@@ -260,15 +260,15 @@ saemix authors for the parts inherited from saemix.

quiet = TRUE, error_model = "tc", cores = 5) f_saem_dfop_sfo <- saem(f_mmkin_dfop_sfo)
#> Running main SAEM algorithm -#> [1] "Mon Feb 15 17:13:15 2021" +#> [1] "Tue Mar 9 17:35:19 2021" #> .... #> Minimisation finished -#> [1] "Mon Feb 15 17:13:26 2021"
summary(f_saem_dfop_sfo, data = TRUE) +#> [1] "Tue Mar 9 17:35:30 2021"
summary(f_saem_dfop_sfo, data = TRUE)
#> saemix version used for fitting: 3.1.9000 -#> mkin version used for pre-fitting: 1.0.3.9000 -#> R version used for fitting: 4.0.3 -#> Date of fit: Mon Feb 15 17:13:27 2021 -#> Date of summary: Mon Feb 15 17:13:27 2021 +#> mkin version used for pre-fitting: 1.0.4.9000 +#> R version used for fitting: 4.0.4 +#> Date of fit: Tue Mar 9 17:35:31 2021 +#> Date of summary: Tue Mar 9 17:35:31 2021 #> #> Equations: #> d_parent/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 * @@ -283,7 +283,7 @@ saemix authors for the parts inherited from saemix.

#> #> Model predictions using solution type analytical #> -#> Fitted in 12.204 s using 300, 100 iterations +#> Fitted in 12.058 s using 300, 100 iterations #> #> Variance model: Two-component variance function #> @@ -300,231 +300,231 @@ saemix authors for the parts inherited from saemix.

#> #> Likelihood computed by importance sampling #> AIC BIC logLik -#> 829.3 823.9 -400.7 +#> 825.9 820.4 -398.9 #> #> Optimised parameters: -#> est. lower upper -#> parent_0 101.29457 97.855 104.7344 -#> log_k_m1 -4.06337 -4.182 -3.9445 -#> f_parent_qlogis -0.94546 -1.307 -0.5841 -#> log_k1 -2.98794 -3.844 -2.1321 -#> log_k2 -3.47891 -4.253 -2.7050 -#> g_qlogis -0.03211 -1.157 1.0931 +#> est. lower upper +#> parent_0 101.118986 97.368 104.8695 +#> log_k_m1 -4.057591 -4.177 -3.9379 +#> f_parent_qlogis -0.933087 -1.290 -0.5763 +#> log_k1 -2.945520 -3.833 -2.0576 +#> log_k2 -3.531954 -4.310 -2.7542 +#> g_qlogis -0.009584 -1.688 1.6687 #> #> Correlation: #> prnt_0 lg_k_1 f_prn_ log_k1 log_k2 -#> log_k_m1 -0.202 -#> f_parent_qlogis -0.145 0.195 -#> log_k1 0.094 -0.099 -0.049 -#> log_k2 -0.042 0.056 0.024 -0.097 -#> g_qlogis -0.005 0.000 0.007 -0.160 -0.113 +#> log_k_m1 -0.198 +#> f_parent_qlogis -0.153 0.184 +#> log_k1 0.080 -0.077 -0.045 +#> log_k2 0.005 0.008 -0.003 -0.019 +#> g_qlogis -0.059 0.048 0.041 -0.334 -0.253 #> #> Random effects: #> est. lower upper -#> SD.parent_0 2.70085 -0.64980 6.0515 -#> SD.log_k_m1 0.08408 -0.04023 0.2084 -#> SD.f_parent_qlogis 0.39215 0.13695 0.6473 -#> SD.log_k1 0.89280 0.27466 1.5109 -#> SD.log_k2 0.82387 0.26388 1.3838 -#> SD.g_qlogis 0.36468 -0.86978 1.5991 +#> SD.parent_0 2.97797 -0.62927 6.5852 +#> SD.log_k_m1 0.09235 -0.02448 0.2092 +#> SD.f_parent_qlogis 0.38712 0.13469 0.6396 +#> SD.log_k1 0.88671 0.27052 1.5029 +#> SD.log_k2 0.80497 0.25587 1.3541 +#> SD.g_qlogis 0.36812 -3.56188 4.2981 #> #> Variance model: #> est. lower upper -#> a.1 0.65724 0.49361 0.82086 -#> b.1 0.06434 0.05034 0.07835 +#> a.1 0.85879 0.68143 1.03615 +#> b.1 0.07787 0.06288 0.09286 #> #> Backtransformed parameters: #> est. lower upper -#> parent_0 101.29457 97.85477 104.73437 -#> k_m1 0.01719 0.01526 0.01936 -#> f_parent_to_m1 0.27980 0.21302 0.35798 -#> k1 0.05039 0.02141 0.11859 -#> k2 0.03084 0.01422 0.06687 -#> g 0.49197 0.23916 0.74896 +#> parent_0 101.11899 97.36850 104.86947 +#> k_m1 0.01729 0.01534 0.01949 +#> f_parent_to_m1 0.28230 0.21587 0.35979 +#> k1 0.05257 0.02163 0.12776 +#> k2 0.02925 0.01344 0.06366 +#> g 0.49760 0.15606 0.84140 #> #> Resulting formation fractions: #> ff -#> parent_m1 0.2798 -#> parent_sink 0.7202 +#> parent_m1 0.2823 +#> parent_sink 0.7177 #> #> Estimated disappearance times: #> DT50 DT90 DT50back DT50_k1 DT50_k2 -#> parent 17.49 61.05 18.38 13.76 22.47 -#> m1 40.32 133.94 NA NA NA +#> parent 17.47 62.31 18.76 13.18 23.7 +#> m1 40.09 133.17 NA NA NA #> #> Data: -#> ds name time observed predicted residual std standardized -#> ds 1 parent 0 89.8 9.878e+01 8.98039 6.3899 1.40541 -#> ds 1 parent 0 104.1 9.878e+01 -5.31961 6.3899 -0.83251 -#> ds 1 parent 1 88.7 9.422e+01 5.52084 6.0981 0.90533 -#> ds 1 parent 1 95.5 9.422e+01 -1.27916 6.0981 -0.20976 -#> ds 1 parent 3 81.8 8.587e+01 4.06752 5.5641 0.73103 -#> ds 1 parent 3 94.5 8.587e+01 -8.63248 5.5641 -1.55147 -#> ds 1 parent 7 71.5 7.180e+01 0.29615 4.6662 0.06347 -#> ds 1 parent 7 70.3 7.180e+01 1.49615 4.6662 0.32063 -#> ds 1 parent 14 54.2 5.360e+01 -0.59602 3.5112 -0.16975 -#> ds 1 parent 14 49.6 5.360e+01 4.00398 3.5112 1.14035 -#> ds 1 parent 28 31.5 3.213e+01 0.62529 2.1691 0.28828 -#> ds 1 parent 28 28.8 3.213e+01 3.32529 2.1691 1.53306 -#> ds 1 parent 60 12.1 1.271e+01 0.60718 1.0490 0.57879 -#> ds 1 parent 60 13.6 1.271e+01 -0.89282 1.0490 -0.85108 -#> ds 1 parent 90 6.2 6.080e+00 -0.12020 0.7649 -0.15716 -#> ds 1 parent 90 8.3 6.080e+00 -2.22020 0.7649 -2.90279 -#> ds 1 parent 120 2.2 3.011e+00 0.81059 0.6852 1.18302 -#> ds 1 parent 120 2.4 3.011e+00 0.61059 0.6852 0.89113 -#> ds 1 m1 1 0.3 1.131e+00 0.83071 0.6613 1.25628 -#> ds 1 m1 1 0.2 1.131e+00 0.93071 0.6613 1.40750 -#> ds 1 m1 3 2.2 3.147e+00 0.94691 0.6877 1.37688 -#> ds 1 m1 3 3.0 3.147e+00 0.14691 0.6877 0.21361 -#> ds 1 m1 7 6.5 6.341e+00 -0.15949 0.7736 -0.20618 -#> ds 1 m1 7 5.0 6.341e+00 1.34051 0.7736 1.73290 -#> ds 1 m1 14 10.2 9.910e+00 -0.28991 0.9157 -0.31659 -#> ds 1 m1 14 9.5 9.910e+00 0.41009 0.9157 0.44783 -#> ds 1 m1 28 12.2 1.255e+01 0.34690 1.0410 0.33323 -#> ds 1 m1 28 13.4 1.255e+01 -0.85310 1.0410 -0.81949 -#> ds 1 m1 60 11.8 1.087e+01 -0.92713 0.9599 -0.96586 -#> ds 1 m1 60 13.2 1.087e+01 -2.32713 0.9599 -2.42434 -#> ds 1 m1 90 6.6 7.813e+00 1.21254 0.8274 1.46541 -#> ds 1 m1 90 9.3 7.813e+00 -1.48746 0.8274 -1.79766 -#> ds 1 m1 120 3.5 5.295e+00 1.79489 0.7403 2.42457 -#> ds 1 m1 120 5.4 5.295e+00 -0.10511 0.7403 -0.14198 -#> ds 2 parent 0 118.0 1.074e+02 -10.63436 6.9396 -1.53242 -#> ds 2 parent 0 99.8 1.074e+02 7.56564 6.9396 1.09021 -#> ds 2 parent 1 90.2 1.012e+02 10.96486 6.5425 1.67594 -#> ds 2 parent 1 94.6 1.012e+02 6.56486 6.5425 1.00342 -#> ds 2 parent 3 96.1 9.054e+01 -5.56014 5.8627 -0.94839 -#> ds 2 parent 3 78.4 9.054e+01 12.13986 5.8627 2.07069 -#> ds 2 parent 7 77.9 7.468e+01 -3.21805 4.8501 -0.66350 -#> ds 2 parent 7 77.7 7.468e+01 -3.01805 4.8501 -0.62226 -#> ds 2 parent 14 56.0 5.748e+01 1.47774 3.7563 0.39340 -#> ds 2 parent 14 54.7 5.748e+01 2.77774 3.7563 0.73948 -#> ds 2 parent 28 36.6 3.996e+01 3.36317 2.6541 1.26717 -#> ds 2 parent 28 36.8 3.996e+01 3.16317 2.6541 1.19182 -#> ds 2 parent 60 22.1 2.132e+01 -0.78225 1.5210 -0.51430 -#> ds 2 parent 60 24.7 2.132e+01 -3.38225 1.5210 -2.22369 -#> ds 2 parent 90 12.4 1.215e+01 -0.25010 1.0213 -0.24487 -#> ds 2 parent 90 10.8 1.215e+01 1.34990 1.0213 1.32169 -#> ds 2 parent 120 6.8 6.931e+00 0.13105 0.7943 0.16500 -#> ds 2 parent 120 7.9 6.931e+00 -0.96895 0.7943 -1.21994 -#> ds 2 m1 1 1.3 1.519e+00 0.21924 0.6645 0.32995 -#> ds 2 m1 3 3.7 4.049e+00 0.34891 0.7070 0.49351 -#> ds 2 m1 3 4.7 4.049e+00 -0.65109 0.7070 -0.92094 -#> ds 2 m1 7 8.1 7.565e+00 -0.53526 0.8179 -0.65448 -#> ds 2 m1 7 7.9 7.565e+00 -0.33526 0.8179 -0.40993 -#> ds 2 m1 14 10.1 1.071e+01 0.60614 0.9521 0.63663 -#> ds 2 m1 14 10.3 1.071e+01 0.40614 0.9521 0.42657 -#> ds 2 m1 28 10.7 1.224e+01 1.54440 1.0260 1.50526 -#> ds 2 m1 28 12.2 1.224e+01 0.04440 1.0260 0.04327 -#> ds 2 m1 60 10.7 1.056e+01 -0.14005 0.9453 -0.14815 -#> ds 2 m1 60 12.5 1.056e+01 -1.94005 0.9453 -2.05226 -#> ds 2 m1 90 9.1 8.089e+00 -1.01088 0.8384 -1.20577 -#> ds 2 m1 90 7.4 8.089e+00 0.68912 0.8384 0.82197 -#> ds 2 m1 120 6.1 5.855e+00 -0.24463 0.7576 -0.32292 -#> ds 2 m1 120 4.5 5.855e+00 1.35537 0.7576 1.78911 -#> ds 3 parent 0 106.2 1.095e+02 3.30335 7.0765 0.46680 -#> ds 3 parent 0 106.9 1.095e+02 2.60335 7.0765 0.36788 -#> ds 3 parent 1 107.4 9.939e+01 -8.01282 6.4287 -1.24641 -#> ds 3 parent 1 96.1 9.939e+01 3.28718 6.4287 0.51133 -#> ds 3 parent 3 79.4 8.365e+01 4.24698 5.4222 0.78326 -#> ds 3 parent 3 82.6 8.365e+01 1.04698 5.4222 0.19309 -#> ds 3 parent 7 63.9 6.405e+01 0.14704 4.1732 0.03523 -#> ds 3 parent 7 62.4 6.405e+01 1.64704 4.1732 0.39467 -#> ds 3 parent 14 51.0 4.795e+01 -3.04985 3.1546 -0.96681 -#> ds 3 parent 14 47.1 4.795e+01 0.85015 3.1546 0.26950 -#> ds 3 parent 28 36.1 3.501e+01 -1.09227 2.3465 -0.46549 -#> ds 3 parent 28 36.6 3.501e+01 -1.59227 2.3465 -0.67858 -#> ds 3 parent 60 20.1 2.012e+01 0.02116 1.4520 0.01457 -#> ds 3 parent 60 19.8 2.012e+01 0.32116 1.4520 0.22119 -#> ds 3 parent 90 11.3 1.206e+01 0.76096 1.0170 0.74826 -#> ds 3 parent 90 10.7 1.206e+01 1.36096 1.0170 1.33825 -#> ds 3 parent 120 8.2 7.230e+00 -0.97022 0.8052 -1.20493 -#> ds 3 parent 120 7.3 7.230e+00 -0.07022 0.8052 -0.08721 -#> ds 3 m1 0 0.8 -5.684e-13 -0.80000 0.6572 -1.21722 -#> ds 3 m1 1 1.8 2.045e+00 0.24538 0.6703 0.36608 -#> ds 3 m1 1 2.3 2.045e+00 -0.25462 0.6703 -0.37987 -#> ds 3 m1 3 4.2 5.136e+00 0.93594 0.7356 1.27228 -#> ds 3 m1 3 4.1 5.136e+00 1.03594 0.7356 1.40822 -#> ds 3 m1 7 6.8 8.674e+00 1.87438 0.8623 2.17381 -#> ds 3 m1 7 10.1 8.674e+00 -1.42562 0.8623 -1.65335 -#> ds 3 m1 14 11.4 1.083e+01 -0.56746 0.9580 -0.59233 -#> ds 3 m1 14 12.8 1.083e+01 -1.96746 0.9580 -2.05369 -#> ds 3 m1 28 11.5 1.098e+01 -0.51762 0.9651 -0.53637 -#> ds 3 m1 28 10.6 1.098e+01 0.38238 0.9651 0.39623 -#> ds 3 m1 60 7.5 8.889e+00 1.38911 0.8713 1.59436 -#> ds 3 m1 60 8.6 8.889e+00 0.28911 0.8713 0.33183 -#> ds 3 m1 90 7.3 6.774e+00 -0.52608 0.7886 -0.66708 -#> ds 3 m1 90 8.1 6.774e+00 -1.32608 0.7886 -1.68150 -#> ds 3 m1 120 5.3 4.954e+00 -0.34584 0.7305 -0.47345 -#> ds 3 m1 120 3.8 4.954e+00 1.15416 0.7305 1.58004 -#> ds 4 parent 0 104.7 9.957e+01 -5.13169 6.4403 -0.79681 -#> ds 4 parent 0 88.3 9.957e+01 11.26831 6.4403 1.74966 -#> ds 4 parent 1 94.2 9.644e+01 2.23888 6.2400 0.35879 -#> ds 4 parent 1 94.6 9.644e+01 1.83888 6.2400 0.29469 -#> ds 4 parent 3 78.1 9.054e+01 12.43946 5.8627 2.12180 -#> ds 4 parent 3 96.5 9.054e+01 -5.96054 5.8627 -1.01669 -#> ds 4 parent 7 76.2 8.004e+01 3.83771 5.1918 0.73919 -#> ds 4 parent 7 77.8 8.004e+01 2.23771 5.1918 0.43101 -#> ds 4 parent 14 70.8 6.511e+01 -5.69246 4.2406 -1.34238 -#> ds 4 parent 14 67.3 6.511e+01 -2.19246 4.2406 -0.51702 -#> ds 4 parent 28 43.1 4.454e+01 1.43744 2.9401 0.48890 -#> ds 4 parent 28 45.1 4.454e+01 -0.56256 2.9401 -0.19134 -#> ds 4 parent 60 21.3 2.132e+01 0.02005 1.5211 0.01318 -#> ds 4 parent 60 23.5 2.132e+01 -2.17995 1.5211 -1.43310 -#> ds 4 parent 90 11.8 1.182e+01 0.02167 1.0053 0.02156 -#> ds 4 parent 90 12.1 1.182e+01 -0.27833 1.0053 -0.27687 -#> ds 4 parent 120 7.0 6.852e+00 -0.14780 0.7914 -0.18675 -#> ds 4 parent 120 6.2 6.852e+00 0.65220 0.7914 0.82408 -#> ds 4 m1 0 1.6 -5.684e-14 -1.60000 0.6572 -2.43444 -#> ds 4 m1 1 0.9 6.918e-01 -0.20821 0.6587 -0.31607 -#> ds 4 m1 3 3.7 1.959e+00 -1.74131 0.6692 -2.60204 -#> ds 4 m1 3 2.0 1.959e+00 -0.04131 0.6692 -0.06173 -#> ds 4 m1 7 3.6 4.076e+00 0.47590 0.7076 0.67252 -#> ds 4 m1 7 3.8 4.076e+00 0.27590 0.7076 0.38989 -#> ds 4 m1 14 7.1 6.698e+00 -0.40189 0.7859 -0.51135 -#> ds 4 m1 14 6.6 6.698e+00 0.09811 0.7859 0.12483 -#> ds 4 m1 28 9.5 9.175e+00 -0.32492 0.8835 -0.36779 -#> ds 4 m1 28 9.3 9.175e+00 -0.12492 0.8835 -0.14141 -#> ds 4 m1 60 8.3 8.818e+00 0.51810 0.8683 0.59671 -#> ds 4 m1 60 9.0 8.818e+00 -0.18190 0.8683 -0.20949 -#> ds 4 m1 90 6.6 6.645e+00 0.04480 0.7841 0.05713 -#> ds 4 m1 90 7.7 6.645e+00 -1.05520 0.7841 -1.34581 -#> ds 4 m1 120 3.7 4.648e+00 0.94805 0.7221 1.31293 -#> ds 4 m1 120 3.5 4.648e+00 1.14805 0.7221 1.58991 -#> ds 5 parent 0 110.4 1.026e+02 -7.81752 6.6333 -1.17853 -#> ds 5 parent 0 112.1 1.026e+02 -9.51752 6.6333 -1.43482 -#> ds 5 parent 1 93.5 9.560e+01 2.10274 6.1865 0.33989 -#> ds 5 parent 1 91.0 9.560e+01 4.60274 6.1865 0.74399 -#> ds 5 parent 3 71.0 8.356e+01 12.55799 5.4165 2.31846 -#> ds 5 parent 3 89.7 8.356e+01 -6.14201 5.4165 -1.13394 -#> ds 5 parent 7 60.4 6.550e+01 5.09732 4.2653 1.19506 -#> ds 5 parent 7 59.1 6.550e+01 6.39732 4.2653 1.49984 -#> ds 5 parent 14 56.5 4.641e+01 -10.09145 3.0576 -3.30044 -#> ds 5 parent 14 47.0 4.641e+01 -0.59145 3.0576 -0.19344 -#> ds 5 parent 28 30.2 2.982e+01 -0.37647 2.0284 -0.18560 -#> ds 5 parent 28 23.9 2.982e+01 5.92353 2.0284 2.92028 -#> ds 5 parent 60 17.0 1.754e+01 0.53981 1.3060 0.41332 -#> ds 5 parent 60 18.7 1.754e+01 -1.16019 1.3060 -0.88834 -#> ds 5 parent 90 11.3 1.175e+01 0.45050 1.0018 0.44969 -#> ds 5 parent 90 11.9 1.175e+01 -0.14950 1.0018 -0.14923 -#> ds 5 parent 120 9.0 7.915e+00 -1.08476 0.8315 -1.30462 -#> ds 5 parent 120 8.1 7.915e+00 -0.18476 0.8315 -0.22220 -#> ds 5 m1 0 0.7 0.000e+00 -0.70000 0.6572 -1.06507 -#> ds 5 m1 1 3.0 3.062e+00 0.06170 0.6861 0.08992 -#> ds 5 m1 1 2.6 3.062e+00 0.46170 0.6861 0.67290 -#> ds 5 m1 3 5.1 8.209e+00 3.10938 0.8432 3.68760 -#> ds 5 m1 3 7.5 8.209e+00 0.70938 0.8432 0.84130 -#> ds 5 m1 7 16.5 1.544e+01 -1.05567 1.1914 -0.88605 -#> ds 5 m1 7 19.0 1.544e+01 -3.55567 1.1914 -2.98436 -#> ds 5 m1 14 22.9 2.181e+01 -1.08765 1.5498 -0.70181 -#> ds 5 m1 14 23.2 2.181e+01 -1.38765 1.5498 -0.89539 -#> ds 5 m1 28 22.2 2.404e+01 1.83624 1.6805 1.09270 -#> ds 5 m1 28 24.4 2.404e+01 -0.36376 1.6805 -0.21647 -#> ds 5 m1 60 15.5 1.875e+01 3.25390 1.3741 2.36805 -#> ds 5 m1 60 19.8 1.875e+01 -1.04610 1.3741 -0.76131 -#> ds 5 m1 90 14.9 1.380e+01 -1.09507 1.1050 -0.99102 -#> ds 5 m1 90 14.2 1.380e+01 -0.39507 1.1050 -0.35753 -#> ds 5 m1 120 10.9 1.002e+01 -0.88429 0.9205 -0.96069 -#> ds 5 m1 120 10.4 1.002e+01 -0.38429 0.9205 -0.41749
# } +#> ds name time observed predicted residual std standardized +#> ds 1 parent 0 89.8 9.838e+01 8.584661 7.7094 1.113536 +#> ds 1 parent 0 104.1 9.838e+01 -5.715339 7.7094 -0.741350 +#> ds 1 parent 1 88.7 9.388e+01 5.182489 7.3611 0.704041 +#> ds 1 parent 1 95.5 9.388e+01 -1.617511 7.3611 -0.219739 +#> ds 1 parent 3 81.8 8.563e+01 3.825382 6.7229 0.569010 +#> ds 1 parent 3 94.5 8.563e+01 -8.874618 6.7229 -1.320062 +#> ds 1 parent 7 71.5 7.169e+01 0.188290 5.6482 0.033336 +#> ds 1 parent 7 70.3 7.169e+01 1.388290 5.6482 0.245795 +#> ds 1 parent 14 54.2 5.361e+01 -0.586595 4.2624 -0.137621 +#> ds 1 parent 14 49.6 5.361e+01 4.013405 4.2624 0.941587 +#> ds 1 parent 28 31.5 3.219e+01 0.688936 2.6496 0.260011 +#> ds 1 parent 28 28.8 3.219e+01 3.388936 2.6496 1.279016 +#> ds 1 parent 60 12.1 1.278e+01 0.678998 1.3145 0.516562 +#> ds 1 parent 60 13.6 1.278e+01 -0.821002 1.3145 -0.624595 +#> ds 1 parent 90 6.2 6.157e+00 -0.043461 0.9835 -0.044188 +#> ds 1 parent 90 8.3 6.157e+00 -2.143461 0.9835 -2.179316 +#> ds 1 parent 120 2.2 3.076e+00 0.876218 0.8916 0.982775 +#> ds 1 parent 120 2.4 3.076e+00 0.676218 0.8916 0.758453 +#> ds 1 m1 1 0.3 1.134e+00 0.833749 0.8633 0.965750 +#> ds 1 m1 1 0.2 1.134e+00 0.933749 0.8633 1.081583 +#> ds 1 m1 3 2.2 3.157e+00 0.957400 0.8933 1.071763 +#> ds 1 m1 3 3.0 3.157e+00 0.157400 0.8933 0.176202 +#> ds 1 m1 7 6.5 6.369e+00 -0.130995 0.9917 -0.132090 +#> ds 1 m1 7 5.0 6.369e+00 1.369005 0.9917 1.380438 +#> ds 1 m1 14 10.2 9.971e+00 -0.229362 1.1577 -0.198112 +#> ds 1 m1 14 9.5 9.971e+00 0.470638 1.1577 0.406513 +#> ds 1 m1 28 12.2 1.265e+01 0.447735 1.3067 0.342637 +#> ds 1 m1 28 13.4 1.265e+01 -0.752265 1.3067 -0.575683 +#> ds 1 m1 60 11.8 1.097e+01 -0.832027 1.2112 -0.686945 +#> ds 1 m1 60 13.2 1.097e+01 -2.232027 1.2112 -1.842825 +#> ds 1 m1 90 6.6 7.876e+00 1.275985 1.0553 1.209109 +#> ds 1 m1 90 9.3 7.876e+00 -1.424015 1.0553 -1.349381 +#> ds 1 m1 120 3.5 5.336e+00 1.835829 0.9540 1.924292 +#> ds 1 m1 120 5.4 5.336e+00 -0.064171 0.9540 -0.067263 +#> ds 2 parent 0 118.0 1.092e+02 -8.812058 8.5459 -1.031142 +#> ds 2 parent 0 99.8 1.092e+02 9.387942 8.5459 1.098529 +#> ds 2 parent 1 90.2 1.023e+02 12.114268 8.0135 1.511724 +#> ds 2 parent 1 94.6 1.023e+02 7.714268 8.0135 0.962654 +#> ds 2 parent 3 96.1 9.066e+01 -5.436165 7.1122 -0.764344 +#> ds 2 parent 3 78.4 9.066e+01 12.263835 7.1122 1.724339 +#> ds 2 parent 7 77.9 7.365e+01 -4.245773 5.7995 -0.732090 +#> ds 2 parent 7 77.7 7.365e+01 -4.045773 5.7995 -0.697604 +#> ds 2 parent 14 56.0 5.593e+01 -0.073803 4.4389 -0.016626 +#> ds 2 parent 14 54.7 5.593e+01 1.226197 4.4389 0.276236 +#> ds 2 parent 28 36.6 3.892e+01 2.320837 3.1502 0.736737 +#> ds 2 parent 28 36.8 3.892e+01 2.120837 3.1502 0.673248 +#> ds 2 parent 60 22.1 2.136e+01 -0.741020 1.8719 -0.395868 +#> ds 2 parent 60 24.7 2.136e+01 -3.341020 1.8719 -1.784841 +#> ds 2 parent 90 12.4 1.251e+01 0.113999 1.2989 0.087765 +#> ds 2 parent 90 10.8 1.251e+01 1.713999 1.2989 1.319575 +#> ds 2 parent 120 6.8 7.338e+00 0.537708 1.0315 0.521281 +#> ds 2 parent 120 7.9 7.338e+00 -0.562292 1.0315 -0.545113 +#> ds 2 m1 1 1.3 1.576e+00 0.276176 0.8675 0.318352 +#> ds 2 m1 3 3.7 4.177e+00 0.476741 0.9183 0.519146 +#> ds 2 m1 3 4.7 4.177e+00 -0.523259 0.9183 -0.569801 +#> ds 2 m1 7 8.1 7.724e+00 -0.376365 1.0485 -0.358970 +#> ds 2 m1 7 7.9 7.724e+00 -0.176365 1.0485 -0.168214 +#> ds 2 m1 14 10.1 1.077e+01 0.674433 1.2006 0.561738 +#> ds 2 m1 14 10.3 1.077e+01 0.474433 1.2006 0.395158 +#> ds 2 m1 28 10.7 1.212e+01 1.416179 1.2758 1.110010 +#> ds 2 m1 28 12.2 1.212e+01 -0.083821 1.2758 -0.065699 +#> ds 2 m1 60 10.7 1.041e+01 -0.294930 1.1807 -0.249793 +#> ds 2 m1 60 12.5 1.041e+01 -2.094930 1.1807 -1.774316 +#> ds 2 m1 90 9.1 8.079e+00 -1.020859 1.0646 -0.958929 +#> ds 2 m1 90 7.4 8.079e+00 0.679141 1.0646 0.637941 +#> ds 2 m1 120 6.1 5.968e+00 -0.131673 0.9765 -0.134843 +#> ds 2 m1 120 4.5 5.968e+00 1.468327 0.9765 1.503683 +#> ds 3 parent 0 106.2 1.036e+02 -2.638248 8.1101 -0.325303 +#> ds 3 parent 0 106.9 1.036e+02 -3.338248 8.1101 -0.411614 +#> ds 3 parent 1 107.4 9.580e+01 -11.600063 7.5094 -1.544743 +#> ds 3 parent 1 96.1 9.580e+01 -0.300063 7.5094 -0.039958 +#> ds 3 parent 3 79.4 8.297e+01 3.574516 6.5182 0.548391 +#> ds 3 parent 3 82.6 8.297e+01 0.374516 6.5182 0.057457 +#> ds 3 parent 7 63.9 6.517e+01 1.272397 5.1472 0.247200 +#> ds 3 parent 7 62.4 6.517e+01 2.772397 5.1472 0.538618 +#> ds 3 parent 14 51.0 4.821e+01 -2.790075 3.8512 -0.724475 +#> ds 3 parent 14 47.1 4.821e+01 1.109925 3.8512 0.288205 +#> ds 3 parent 28 36.1 3.385e+01 -2.250573 2.7723 -0.811811 +#> ds 3 parent 28 36.6 3.385e+01 -2.750573 2.7723 -0.992168 +#> ds 3 parent 60 20.1 1.964e+01 -0.455700 1.7543 -0.259760 +#> ds 3 parent 60 19.8 1.964e+01 -0.155700 1.7543 -0.088753 +#> ds 3 parent 90 11.3 1.210e+01 0.795458 1.2746 0.624068 +#> ds 3 parent 90 10.7 1.210e+01 1.395458 1.2746 1.094792 +#> ds 3 parent 120 8.2 7.451e+00 -0.749141 1.0364 -0.722816 +#> ds 3 parent 120 7.3 7.451e+00 0.150859 1.0364 0.145558 +#> ds 3 m1 0 0.8 3.695e-13 -0.800000 0.8588 -0.931542 +#> ds 3 m1 1 1.8 1.740e+00 -0.059741 0.8694 -0.068714 +#> ds 3 m1 1 2.3 1.740e+00 -0.559741 0.8694 -0.643812 +#> ds 3 m1 3 4.2 4.531e+00 0.331379 0.9285 0.356913 +#> ds 3 m1 3 4.1 4.531e+00 0.431379 0.9285 0.464618 +#> ds 3 m1 7 6.8 8.113e+00 1.312762 1.0661 1.231333 +#> ds 3 m1 7 10.1 8.113e+00 -1.987238 1.0661 -1.863971 +#> ds 3 m1 14 11.4 1.079e+01 -0.613266 1.2013 -0.510507 +#> ds 3 m1 14 12.8 1.079e+01 -2.013266 1.2013 -1.675923 +#> ds 3 m1 28 11.5 1.133e+01 -0.174252 1.2310 -0.141553 +#> ds 3 m1 28 10.6 1.133e+01 0.725748 1.2310 0.589558 +#> ds 3 m1 60 7.5 8.948e+00 1.448281 1.1059 1.309561 +#> ds 3 m1 60 8.6 8.948e+00 0.348281 1.1059 0.314922 +#> ds 3 m1 90 7.3 6.665e+00 -0.634932 1.0034 -0.632752 +#> ds 3 m1 90 8.1 6.665e+00 -1.434932 1.0034 -1.430004 +#> ds 3 m1 120 5.3 4.795e+00 -0.504936 0.9365 -0.539199 +#> ds 3 m1 120 3.8 4.795e+00 0.995064 0.9365 1.062586 +#> ds 4 parent 0 104.7 9.985e+01 -4.850494 7.8227 -0.620050 +#> ds 4 parent 0 88.3 9.985e+01 11.549506 7.8227 1.476402 +#> ds 4 parent 1 94.2 9.676e+01 2.556304 7.5834 0.337093 +#> ds 4 parent 1 94.6 9.676e+01 2.156304 7.5834 0.284346 +#> ds 4 parent 3 78.1 9.092e+01 12.817485 7.1318 1.797230 +#> ds 4 parent 3 96.5 9.092e+01 -5.582515 7.1318 -0.782764 +#> ds 4 parent 7 76.2 8.050e+01 4.297338 6.3270 0.679204 +#> ds 4 parent 7 77.8 8.050e+01 2.697338 6.3270 0.426320 +#> ds 4 parent 14 70.8 6.562e+01 -5.179989 5.1816 -0.999687 +#> ds 4 parent 14 67.3 6.562e+01 -1.679989 5.1816 -0.324222 +#> ds 4 parent 28 43.1 4.499e+01 1.886936 3.6069 0.523140 +#> ds 4 parent 28 45.1 4.499e+01 -0.113064 3.6069 -0.031346 +#> ds 4 parent 60 21.3 2.151e+01 0.214840 1.8827 0.114114 +#> ds 4 parent 60 23.5 2.151e+01 -1.985160 1.8827 -1.054433 +#> ds 4 parent 90 11.8 1.190e+01 0.098528 1.2633 0.077990 +#> ds 4 parent 90 12.1 1.190e+01 -0.201472 1.2633 -0.159475 +#> ds 4 parent 120 7.0 6.886e+00 -0.113832 1.0125 -0.112431 +#> ds 4 parent 120 6.2 6.886e+00 0.686168 1.0125 0.677724 +#> ds 4 m1 0 1.6 4.263e-14 -1.600000 0.8588 -1.863085 +#> ds 4 m1 1 0.9 7.140e-01 -0.185984 0.8606 -0.216112 +#> ds 4 m1 3 3.7 2.022e+00 -1.678243 0.8731 -1.922160 +#> ds 4 m1 3 2.0 2.022e+00 0.021757 0.8731 0.024919 +#> ds 4 m1 7 3.6 4.207e+00 0.607229 0.9192 0.660633 +#> ds 4 m1 7 3.8 4.207e+00 0.407229 0.9192 0.443044 +#> ds 4 m1 14 7.1 6.912e+00 -0.188339 1.0135 -0.185828 +#> ds 4 m1 14 6.6 6.912e+00 0.311661 1.0135 0.307506 +#> ds 4 m1 28 9.5 9.449e+00 -0.050714 1.1309 -0.044843 +#> ds 4 m1 28 9.3 9.449e+00 0.149286 1.1309 0.132004 +#> ds 4 m1 60 8.3 8.997e+00 0.697403 1.1083 0.629230 +#> ds 4 m1 60 9.0 8.997e+00 -0.002597 1.1083 -0.002343 +#> ds 4 m1 90 6.6 6.697e+00 0.096928 1.0047 0.096472 +#> ds 4 m1 90 7.7 6.697e+00 -1.003072 1.0047 -0.998348 +#> ds 4 m1 120 3.7 4.622e+00 0.921607 0.9312 0.989749 +#> ds 4 m1 120 3.5 4.622e+00 1.121607 0.9312 1.204537 +#> ds 5 parent 0 110.4 1.045e+02 -5.942426 8.1795 -0.726502 +#> ds 5 parent 0 112.1 1.045e+02 -7.642426 8.1795 -0.934338 +#> ds 5 parent 1 93.5 9.739e+01 3.893915 7.6327 0.510162 +#> ds 5 parent 1 91.0 9.739e+01 6.393915 7.6327 0.837700 +#> ds 5 parent 3 71.0 8.519e+01 14.188275 6.6891 2.121098 +#> ds 5 parent 3 89.7 8.519e+01 -4.511725 6.6891 -0.674487 +#> ds 5 parent 7 60.4 6.684e+01 6.439546 5.2753 1.220701 +#> ds 5 parent 7 59.1 6.684e+01 7.739546 5.2753 1.467133 +#> ds 5 parent 14 56.5 4.736e+01 -9.138979 3.7868 -2.413407 +#> ds 5 parent 14 47.0 4.736e+01 0.361021 3.7868 0.095338 +#> ds 5 parent 28 30.2 3.033e+01 0.131178 2.5132 0.052195 +#> ds 5 parent 28 23.9 3.033e+01 6.431178 2.5132 2.558936 +#> ds 5 parent 60 17.0 1.771e+01 0.705246 1.6243 0.434177 +#> ds 5 parent 60 18.7 1.771e+01 -0.994754 1.6243 -0.612409 +#> ds 5 parent 90 11.3 1.180e+01 0.504856 1.2580 0.401315 +#> ds 5 parent 90 11.9 1.180e+01 -0.095144 1.2580 -0.075631 +#> ds 5 parent 120 9.0 7.917e+00 -1.083499 1.0571 -1.024928 +#> ds 5 parent 120 8.1 7.917e+00 -0.183499 1.0571 -0.173579 +#> ds 5 m1 0 0.7 3.553e-15 -0.700000 0.8588 -0.815100 +#> ds 5 m1 1 3.0 3.204e+00 0.204414 0.8943 0.228572 +#> ds 5 m1 1 2.6 3.204e+00 0.604414 0.8943 0.675845 +#> ds 5 m1 3 5.1 8.586e+00 3.485889 1.0884 3.202858 +#> ds 5 m1 3 7.5 8.586e+00 1.085889 1.0884 0.997722 +#> ds 5 m1 7 16.5 1.612e+01 -0.376855 1.5211 -0.247743 +#> ds 5 m1 7 19.0 1.612e+01 -2.876855 1.5211 -1.891237 +#> ds 5 m1 14 22.9 2.267e+01 -0.228264 1.9633 -0.116267 +#> ds 5 m1 14 23.2 2.267e+01 -0.528264 1.9633 -0.269072 +#> ds 5 m1 28 22.2 2.468e+01 2.480178 2.1050 1.178211 +#> ds 5 m1 28 24.4 2.468e+01 0.280178 2.1050 0.133099 +#> ds 5 m1 60 15.5 1.860e+01 3.099615 1.6838 1.840794 +#> ds 5 m1 60 19.8 1.860e+01 -1.200385 1.6838 -0.712883 +#> ds 5 m1 90 14.9 1.326e+01 -1.636345 1.3433 -1.218195 +#> ds 5 m1 90 14.2 1.326e+01 -0.936345 1.3433 -0.697072 +#> ds 5 m1 120 10.9 9.348e+00 -1.551535 1.1258 -1.378133 +#> ds 5 m1 120 10.4 9.348e+00 -1.051535 1.1258 -0.934014
# }
-- cgit v1.2.1