median_test: Monte-Carlo test of differences between medians

In DIDE-MRes/pkg1:

Description

This is a test we developed with DIDE's MRes.

Usage

median_test(x, fac, n = 999, ...)

Arguments

x	A vector of numeric values
fac	A factor defining 2 groups
n	The number of observations

Examples

## CREATE DATA ##
set.seed(1)
titer <- c(rexp(100, rate = 0.002), rexp(20, rate = 0.0005), # control
        rexp(55, rate = 0.006), rexp(15, rate = 0.0005)) # treatment

type <- factor(rep(c("control","treatment"), c(120,70)))

## make a data frame
VL <- data.frame(titer = titer, type = type)

DIDE-MRes/pkg1 documentation built on May 25, 2019, 8:29 p.m.