Description Usage Arguments Value Examples
View source: R/model_netmeta.R
Frequentist net-meta model for MD, HR, and RR
This function is the warpper of netmeta and pairwise functions from netmeta package.
R packgae netmeta is required.
1 2 3 4 5 6 7 | model_netmeta(
long.data,
id.treatments,
treatment,
reference,
outcome = c("MD", "HR", "RR")
)
|
long.data |
data.frame to be analyzed should be formatted in long format "MD" needs treat, n, mean, sd; "RR" needs treat, event, n "HR" needs treat, event, time |
id.treatments |
data.frame to specify the id and treatments |
treatment |
treatment names or labels |
reference |
the referential id in the net-meta |
outcome |
the outcome should be MD-mean difference, HR-hazard ratio, and RR-risk ratio |
Summary list of the results
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 | LDT1 <- read.csv(system.file("extdata", "HR_SH_D.csv", package = "net.meta"))
trt1 <- read.table(system.file("extdata", "HR_SH_D.txt", package = "net.meta"),
header=TRUE,quote = '"', stringsAsFactors=FALSE)
trt1$description <- factor(trt1$description, trt1$description)
LDT1$study <- factor(LDT1$study, unique(LDT1$study))
trt1$label <- paste0(trt1$id,"-", trt1$description)
LDT1$label <- factor(LDT1$treatment, labels = trt1$label)
set.seed(1)
nmt.lab1 <- model_netmeta(
long.data=LDT1[,-2],
treatment=LDT1$label,
id.treatments=trt1,
reference = "A",
outcome="HR")
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.