From e5d1df9a9b1f0951d7dfbaf24eee4294470b73e2 Mon Sep 17 00:00:00 2001
From: Johannes Ranke
An mkinfit
object, or an mmkin
column
object containing two fits to the same data.
Optionally, another mkinfit object fitted to the same data.
Argument to mkinfit
, passed to
update.mkinfit
for creating the alternative fitted object.
# \dontrun{
-test_data <- subset(synthetic_data_for_UBA_2014[[12]]$data, name == "parent")
-sfo_fit <- mkinfit("SFO", test_data, quiet = TRUE)
-dfop_fit <- mkinfit("DFOP", test_data, quiet = TRUE)
-lrtest(dfop_fit, sfo_fit)
+ # \dontrun{
+test_data <- subset(synthetic_data_for_UBA_2014[[12]]$data, name == "parent")
+sfo_fit <- mkinfit("SFO", test_data, quiet = TRUE)
+dfop_fit <- mkinfit("DFOP", test_data, quiet = TRUE)
+lrtest(dfop_fit, sfo_fit)
#> Likelihood ratio test
#>
#> Model 1: DFOP with error model const
@@ -139,7 +150,7 @@ lower number of fitted parameters (null hypothesis).
#> 2 3 -63.954 -2 43.002 4.594e-10 ***
#> ---
#> Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
-lrtest(sfo_fit, dfop_fit)
+lrtest(sfo_fit, dfop_fit)
#> Likelihood ratio test
#>
#> Model 1: DFOP with error model const
@@ -149,14 +160,14 @@ lower number of fitted parameters (null hypothesis).
#> 2 3 -63.954 -2 43.002 4.594e-10 ***
#> ---
#> Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
-
-# The following two examples are commented out as they fail during
-# generation of the static help pages by pkgdown
-#lrtest(dfop_fit, error_model = "tc")
-#lrtest(dfop_fit, fixed_parms = c(k2 = 0))
-
-# However, this equivalent syntax also works for static help pages
-lrtest(dfop_fit, update(dfop_fit, error_model = "tc"))
+
+# The following two examples are commented out as they fail during
+# generation of the static help pages by pkgdown
+#lrtest(dfop_fit, error_model = "tc")
+#lrtest(dfop_fit, fixed_parms = c(k2 = 0))
+
+# However, this equivalent syntax also works for static help pages
+lrtest(dfop_fit, update(dfop_fit, error_model = "tc"))
#> Likelihood ratio test
#>
#> Model 1: DFOP with error model tc
@@ -166,7 +177,7 @@ lower number of fitted parameters (null hypothesis).
#> 2 5 -42.453 -1 15.731 7.302e-05 ***
#> ---
#> Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
-lrtest(dfop_fit, update(dfop_fit, fixed_parms = c(k2 = 0)))
+lrtest(dfop_fit, update(dfop_fit, fixed_parms = c(k2 = 0)))
#> Likelihood ratio test
#>
#> Model 1: DFOP with error model const
@@ -176,7 +187,7 @@ lower number of fitted parameters (null hypothesis).
#> 2 4 -57.340 -1 29.776 4.851e-08 ***
#> ---
#> Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
-# }
+# }