From 5015e75f7bdc8c5761317d5fe760896eb6c265be Mon Sep 17 00:00:00 2001 From: "Mattan S. Ben-Shachar" Date: Thu, 18 Dec 2025 13:28:03 +0200 Subject: [PATCH] Update utils_ci.R --- R/utils_ci.R | 11 ++++++----- 1 file changed, 6 insertions(+), 5 deletions(-) diff --git a/R/utils_ci.R b/R/utils_ci.R index f6c3f101..949a8a0b 100644 --- a/R/utils_ci.R +++ b/R/utils_ci.R @@ -10,8 +10,9 @@ probs <- c(alpha / 2, 1 - alpha / 2) lambda <- f * df + init <- rep(ifelse(lambda < 2, 2, lambda / 2), 2) ncp <- suppressWarnings(stats::optim( - par = 1.1 * rep(lambda, 2), + par = init, fn = function(x) { quan <- stats::qf(p = probs, df, df_error, ncp = x) sum(abs(quan - f)) @@ -39,9 +40,9 @@ alpha <- 1 - conf.level probs <- c(alpha / 2, 1 - alpha / 2) - + init <- rep(ifelse(abs(t) < 2, sign(t) * 2, t / 2), 2) ncp <- suppressWarnings(stats::optim( - par = 1.1 * rep(t, 2), + par = init, fn = function(x) { quan <- stats::qt(p = probs, df = df_error, ncp = x) sum(abs(quan - t)) @@ -62,9 +63,9 @@ alpha <- 1 - conf.level probs <- c(alpha / 2, 1 - alpha / 2) - + init <- rep(ifelse(chisq < 2, 2, chisq / 2), 2) ncp <- suppressWarnings(stats::optim( - par = 1.1 * rep(chisq, 2), + par = init, fn = function(x) { quan <- stats::qchisq(p = probs, df, ncp = x) sum(abs(quan - chisq))