skrmdb.all: Function to run all three methods for determining ED50.

In ABS-dev/skrmdb: Package to estimate ED50

Function to run all three methods for determining ED50.

Description

This function is used primarily by the CVB statistics section to run all three methods included in the skrmdb package for determining ED50.

Usage

skrmdb.all(formula, data, autosort = TRUE)

Arguments

formula	A formula of the form y + n ~ x or y + n ~ x | v1 + ... + vK. y is the number responding at each dilution level, n is the number tested at each dilution level, and x is the dilution levels. v1, ..., vK are the grouping variables. y, n, x, v1, ..., vK must all be distinct.
data	A data.frame containing the titration data. Formatted as specified in the CVB Data Guide.
autosort	Default TRUE. If TRUE will sort the data according to either sort(x) or sort(-x) so that y / n appears to be increasing with the index. This is how the three methods assume the data to be ordered.

Value

A data.frame containing columns for each of the subset variables, the ED50 as computed by DragBehr, ReedMuench, the ED50 and variance as computed by SpearKarb, and columns of logical values which indicate if the data are increasing or decreasing with log_dil, has an even dilution scheme, are monotonic, and brackets the midpoint.

Examples

data(titration)
titration$log_dil = -log10(titration$dil)
skrmdb.all(positive + total ~ log_dil | Vial + Operator, titration)

titration$dil = NULL
skrmdb.all(positive + total ~ log_dil | ., titration)

ABS-dev/skrmdb documentation built on Sept. 19, 2024, 10:32 a.m.