Nothing
mrmc_tests <- function(data, cov, design) {
tests <- data[[2]]
data <- data[-2]
formula <- formula(data)
mapply(function(test) {
is_test <- tests == test
structure(
list(
design = design,
aov = aov_mrmc(formula, data = subset(data, is_test)),
cov = cov[is_test, is_test],
data = data
),
class = "mrmc_tests"
)
}, levels(tests), SIMPLIFY = FALSE)
}
summary.mrmc_tests_rrrc <- function(object, conf.level = 0.95, ...) {
comps <- vcov_comps(object)
n <- comps$n
MS <- comps$MS
cov <- comps$cov
df <- data.frame(
Estimate = mean(object),
`MS(R)` = MS[["R"]],
Cov2 = cov[2],
StdErr = sqrt(MS[["R"]] / n[["reader"]] + max(cov[2], 0)),
df = (MS[["R"]] + n[["reader"]] * max(cov[2], 0))^2 /
(MS[["R"]]^2 / (n[["reader"]] - 1)),
check.names = FALSE
)
df$CI <- with(df, {
ci <- Estimate +
qt((1 + conf.level) / 2, df) * StdErr %o% c(Lower = -1, Upper = 1)
trunc_ci(object, ci)
})
df
}
summary.mrmc_tests_frrc <- function(object, conf.level = 0.95, ...) {
comps <- vcov_comps(object)
n <- comps$n
MS <- comps$MS
cov <- comps$cov
df <- data.frame(
Estimate = mean(object),
`Var(Error)` = comps$var,
Cov2 = cov[2],
StdErr = sqrt((comps$var + (n[["reader"]] - 1) * max(cov[2], 0)) /
n[["reader"]]),
check.names = FALSE
)
df$CI <- with(df, {
ci <- Estimate +
qnorm((1 + conf.level) / 2) * StdErr %o% c(Lower = -1, Upper = 1)
trunc_ci(object, ci)
})
df
}
summary.mrmc_tests_rrfc <- function(object, conf.level = 0.95, ...) {
n <- dim(object)
MS <- meansq(object)
df <- data.frame(
Estimate = mean(object),
`MS(R)` = MS[["R"]],
StdErr = sqrt(MS[["R"]] / n[["reader"]]),
df = n[["reader"]] - 1,
check.names = FALSE
)
df$CI <- with(df, {
ci <- Estimate +
qt((1 + conf.level) / 2, df) * StdErr %o% c(Lower = -1, Upper = 1)
trunc_ci(object, ci)
})
df
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.