subgrp_perf_pred {BioPred} | R Documentation |
Subgroup Performance Evaluation for Predictive Cases
Description
This function evaluates the performance of subgroups based on different types of response variables in predictive cases.
Usage
subgrp_perf_pred(
yvar,
censorvar = NULL,
grpvar,
grpname,
trtvar,
trtname,
xvars.adj = NULL,
data,
type,
yvar.display = yvar,
grpvar.display = grpvar,
trtvar.display = trtvar
)
Arguments
yvar |
Response variable name. |
censorvar |
Censoring variable name (0-censored, 1-event). |
grpvar |
Subgroup variable name. |
grpname |
A vector of ordered subgroup names (values in the column of grpvar). |
trtvar |
Treatment variable name. |
trtname |
A vector of ordered treatment names (values in the column of trtvar). |
xvars.adj |
Other covariates to adjust when evaluating the performance. |
data |
The dataset. |
type |
"c" for continuous; "s" for "survival", and "b" for binary. |
yvar.display |
Display name of the response variable. |
grpvar.display |
Display name of the group variable. |
trtvar.display |
Display name of the treatment variable. |
Value
A list containing the comparison results, group results, and possibly a plot.
Examples
# Load a sample dataset
data <- data.frame(
response = rnorm(100, mean = 10, sd = 5), # continuous response
survival_time = rexp(100, rate = 0.1), # survival time
status = sample(c(0, 1), 100, replace = TRUE), # censoring status
group = sample(c("Low", "Medium", "High"), 100, replace = TRUE), # subgroup variable
treatment = sample(c("A", "B"), 100, replace = TRUE) # treatment variable
)
# Subgroup performance evaluation for predictive cases - survival analysis
subgrp_perf_pred(
yvar = "survival_time",
censorvar = "status",
grpvar = "group",
grpname = c("Low", "Medium", "High"),
trtvar = "treatment",
trtname = c("A", "B"),
data = data,
type = "s",
yvar.display = "Survival Time",
grpvar.display = "Risk Group",
trtvar.display = "Treatment"
)
# Subgroup performance evaluation for predictive cases - continuous outcome
subgrp_perf_pred(
yvar = "response",
grpvar = "group",
grpname = c("Low", "Medium", "High"),
trtvar = "treatment",
trtname = c("A", "B"),
data = data,
type = "c",
yvar.display = "Response",
grpvar.display = "Risk Group",
trtvar.display = "Treatment"
)
# Subgroup performance evaluation for predictive cases - binary outcome
data$binary_response <- sample(c(0, 1), 100, replace = TRUE)
subgrp_perf_pred(
yvar = "binary_response",
grpvar = "group",
grpname = c("Low", "Medium", "High"),
trtvar = "treatment",
trtname = c("A", "B"),
data = data,
type = "b",
yvar.display = "Binary Response",
grpvar.display = "Risk Group",
trtvar.display = "Treatment"
)
[Package BioPred version 1.0.2 Index]