From a9427a09abdf7ce9aaeae7c7190f90c8f2e5ef52 Mon Sep 17 00:00:00 2001 From: Johannes Ranke Date: Mon, 15 Feb 2021 14:08:13 +0100 Subject: Improve README, introductory vignette and some other docs Also bump version to 1.0.3. --- docs/reference/confint.mkinfit.html | 44 ++++++++++++++++++------------------- 1 file changed, 22 insertions(+), 22 deletions(-) (limited to 'docs/reference/confint.mkinfit.html') diff --git a/docs/reference/confint.mkinfit.html b/docs/reference/confint.mkinfit.html index b71e039a..06e78459 100644 --- a/docs/reference/confint.mkinfit.html +++ b/docs/reference/confint.mkinfit.html @@ -79,7 +79,7 @@ method of Venzon and Moolgavkar (1988)." /> mkin - 1.0.0 + 1.0.3 @@ -285,13 +285,13 @@ Profile-Likelihood Based Confidence Intervals, Applied Statistics, 37, f_d_1 <- mkinfit(SFO_SFO, subset(FOCUS_2006_D, value != 0), quiet = TRUE) system.time(ci_profile <- confint(f_d_1, method = "profile", cores = 1, quiet = TRUE))
#> user system elapsed -#> 4.324 0.980 3.960
# Using more cores does not save much time here, as parent_0 takes up most of the time +#> 4.255 1.029 3.937
# Using more cores does not save much time here, as parent_0 takes up most of the time # If we additionally exclude parent_0 (the confidence of which is often of # minor interest), we get a nice performance improvement if we use at least 4 cores system.time(ci_profile_no_parent_0 <- confint(f_d_1, method = "profile", c("k_parent_sink", "k_parent_m1", "k_m1_sink", "sigma"), cores = n_cores))
#> Profiling the likelihood
#> user system elapsed -#> 1.480 0.109 0.905
ci_profile +#> 1.469 0.092 0.911
ci_profile
#> 2.5% 97.5% #> parent_0 96.456003640 1.027703e+02 #> k_parent_sink 0.040762501 5.549764e-02 @@ -349,14 +349,14 @@ Profile-Likelihood Based Confidence Intervals, Applied Statistics, 37, #> sigma 2.535612399 3.985263e+00
ci_quadratic_transformed_ff <- confint(f_d_2, method = "quadratic") ci_quadratic_transformed_ff
#> 2.5% 97.5% -#> parent_0 96.403833585 102.79311650 +#> parent_0 96.403833578 102.79311649 #> k_parent 0.090823771 0.10725430 #> k_m1 0.004012219 0.00689755 #> f_parent_to_m1 0.469118824 0.55959615 #> sigma 2.396089689 3.85491806
ci_quadratic_untransformed_ff <- confint(f_d_2, method = "quadratic", transformed = FALSE) ci_quadratic_untransformed_ff
#> 2.5% 97.5% -#> parent_0 96.403833589 1.027931e+02 +#> parent_0 96.403833583 1.027931e+02 #> k_parent 0.090491913 1.069035e-01 #> k_m1 0.003835485 6.685823e-03 #> f_parent_to_m1 0.469113477 5.598387e-01 @@ -374,15 +374,15 @@ Profile-Likelihood Based Confidence Intervals, Applied Statistics, 37, #> f_parent_to_m1 TRUE FALSE #> sigma TRUE FALSE
rel_diffs_transformed_ff
#> 2.5% 97.5% -#> parent_0 0.0005408689 0.0002217234 +#> parent_0 0.0005408690 0.0002217233 #> k_parent 0.0009598532 0.0009001864 -#> k_m1 0.0307283044 0.0290588365 -#> f_parent_to_m1 0.0046881768 0.0027780063 +#> k_m1 0.0307283041 0.0290588361 +#> f_parent_to_m1 0.0046881769 0.0027780063 #> sigma 0.0550252516 0.0327066836
rel_diffs_untransformed_ff
#> 2.5% 97.5% -#> parent_0 0.0005408689 0.0002217233 -#> k_parent 0.0046102155 0.0023732281 -#> k_m1 0.0146740688 0.0025291817 +#> parent_0 0.0005408689 0.0002217232 +#> k_parent 0.0046102156 0.0023732281 +#> k_m1 0.0146740690 0.0025291820 #> f_parent_to_m1 0.0046995211 0.0023457712 #> sigma 0.0550252516 0.0327066836
# The profiling for the following fit does not finish in a reasonable time, @@ -396,18 +396,18 @@ Profile-Likelihood Based Confidence Intervals, Applied Statistics, 37, error_model_algorithm = "direct", quiet = TRUE) confint(f_tc_2, method = "quadratic")
#> 2.5% 97.5% -#> parent_0 94.596126334 106.19944007 -#> k_M1 0.037605408 0.04490759 -#> k_M2 0.008568739 0.01087675 -#> f_parent_to_M1 0.021463787 0.62023881 -#> f_parent_to_M2 0.015166531 0.37975349 -#> k1 0.273897467 0.33388084 -#> k2 0.018614555 0.02250379 -#> g 0.671943606 0.73583278 -#> sigma_low 0.251283766 0.83992113 -#> rsd_high 0.040411014 0.07662005
confint(f_tc_2, "parent_0", method = "quadratic") +#> parent_0 94.596039609 106.19954892 +#> k_M1 0.037605368 0.04490762 +#> k_M2 0.008568731 0.01087676 +#> f_parent_to_M1 0.021462489 0.62023882 +#> f_parent_to_M2 0.015165617 0.37975348 +#> k1 0.273897348 0.33388101 +#> k2 0.018614554 0.02250378 +#> g 0.671943411 0.73583305 +#> sigma_low 0.251283495 0.83992077 +#> rsd_high 0.040411024 0.07662008
confint(f_tc_2, "parent_0", method = "quadratic")
#> 2.5% 97.5% -#> parent_0 94.59613 106.1994
# } +#> parent_0 94.59604 106.1995
# }