multeq.diff {MultEq}R Documentation

Equivalence for differences of means of multiple endpoints

Description

Performs equivalence tests and related confidence intervals for differences of two normal means of multiple endpoints.

Usage

multeq.diff(data, grp, resp = NULL, base = 1, margin.lo = NULL, margin.up = NULL,
            method = "single.step", var.equal = FALSE, FWER = 0.05)

Arguments

data a data frame containing response variables (endpoints) and the group variable as columns, the data must have exactly two treatment groups
grp the name of the group variable in " "
resp a vector of names of the response variables (endpoints) in " "
base a single integer specifying the base/control group
margin.lo a vector of absolute lower margins under the null hypotheses relating to the endpoints
margin.up a vector of absolute upper margins under the null hypotheses relating to the endpoints
method a character string:
  • "step.up": method of Quan et al. (2001),
  • "single.step": Bonferroni-adjusted single-step procedure
  • var.equal a logical indicating homogeneous or heterogeneous variances of the data
    FWER a single numeric value specifying the familywise error rate to be controlled

    Details

    The object is to show equivalence for two treatment groups on multiple primary, normally distributed response variables (endpoints). If margin.up is not given, one-sided tests are done for the alternative hypothesis that the differences (to the base group) of the means is larger than margin.lo. Analogously, same vice versa. Only if both margin.lo and margin.up are given, a two-sided equivalence tests for differences is done. Bonferroni adjusted "two one-sided t-tests" (TOST) and related simultaneous confidence intervals are used for method "single.step"; the method of Quan et al. (2001) is applied for "step.up". Welch t-tests and related confidence intervals are used for var.equal=FALSE.

    Value

    An object of class multeq.diff containing:

    estimate a (named) vector of estimated differences
    test.stat a (named) vector of the calculated test statistics
    degr.fr either a single degree of freedom (var.equal=TRUE) or a (named) vector of degrees of freedom (var.equal=FALSE)
    p.value a (named) vector of p-values adjusted for multiplicity
    lower a (named) vector of lower confidence limits
    upper a (named) vector of upper confidence limits

    Note

    Because related to the TOST method, the two-sided confidence intervals for method="single.step" have simultaneous coverage probability (1-2alpha). The intervals for method="step.up" are stepwise adjusted and only applicable for test decisions, not for a simultaneous parameter estimation or comparing among each other.

    Author(s)

    Mario Hasler

    References

    Quan et al. (2001): Assessmant of equivalence on multiple endpoints, Statistics in Medicine 20, 3159-3173

    See Also

    multeq.rat

    Examples

    data(clinic)
    
    comp <- multeq.diff(data=clinic,grp="fact",method="step.up",margin.up=rep(0.6,5),
                        margin.lo=-rep(0.6,5))
    summary(comp)
    

    [Package MultEq version 2.0 Index]