아래에서 resulttitled 출력의 마지막 부분을 어떻게 추출 할 수 있는지 궁금합니다 Likelihood Ratio Test.
library(weightr)
y <- dat.bangertdrowns2004$yi
v <- dat.bangertdrowns2004$vi
(result <- weightfunct(y, v))
# .
# .
# .
Likelihood Ratio Test: ## This part
X^2(df = 1) = 1.163544, p-val = 0.28073 ## This part
또 다른 옵션은 출력을 생성하는 함수를 보는 것입니다. 이 경우 weightr:::print.weightfunct, 관련 라인은
<snip>
cat("\n")
cat("Likelihood Ratio Test:")
cat("\n")
df <- length(x[[2]]$par) - length(x[[1]]$par)
lrchisq <- 2 * (abs(x[[1]]$value - x[[2]]$value))
pvalue <- 1 - pchisq(lrchisq, df)
cat("X^2(df = ", df, ") = ", lrchisq, ", p-val = ", format.pval(pvalue),
sep = "")
</snip>
따라서 값을
x = result
df <- length(x[[2]]$par) - length(x[[1]]$par)
lrchisq <- 2 * (abs(x[[1]]$value - x[[2]]$value))
pvalue <- 1 - pchisq(lrchisq, df)
> lrchisq
#[1] 1.163544
> pvalue
#[1] 0.2807317
이 기사는 인터넷에서 수집됩니다. 재 인쇄 할 때 출처를 알려주십시오.
침해가 발생한 경우 연락 주시기 바랍니다[email protected] 삭제
몇 마디 만하겠습니다