permb {IncomPair}R Documentation

Permutation Based Tests for Incomplete Paired Data

Description

The function Performs testing the hypothesis of equality of means for the incomplete pairs setting data. The function combines the observed mean difference for the complete pairs with the difference between the two means of the independent samples. The function implements two different nonparametric tests based on permutation tests that were proposed by Einsporn and Habtzghi (2013), and Maritz (1995). The two methods are denoted by EH and Maritz, respectively.

Usage

permb(xp, yp, xu, yu, r, mu = NULL, method = "EH",
      alternative = "two.sided", verbose = TRUE)

Arguments

xp, yp

(non-empty) numeric vectors of data values of the the complete pairs

xu

a numeric vector of data on x only

yu

a numeric vector of data on y only

r

a number indicating the correlation between the complete pairs

mu

a number indicating the true value of the mean (or difference in means if you are performing a two sample test)

method

a character string specifying the different type of methods, must be one of "EH" (default) or "Maritz"

alternative

a character string specifying the alternative hypothesis, must be one of "two.sided" (default), "greater" or "less"

verbose

if TRUE, show the test used, null and alternative hypotheses in addition to the p-value

Value

A S4 object containing the following components:

Title

a character string describing the test used

Nhypothesis

a character string describing the null hypothesis

Ahypothesis

a character string describing the alternative hypothesis

Pval

the p-value for the test

References

1 Einsporn, R. L., & Habtzghi, D. (2013). Combining paired and two-sample data using a permutation test. Journal of Data Science, 11(4), 767-779.

2 Maritz, J. S. (1995). A permutation paired test allowing for missing values. Australian Journal of Statistics, 37(2), 153-159.

3 Woolson, R., Leeper, J., Cole, J. and Clarke, W. (1976). A Monte Carlo investigation of a statistic for a bivariate missing data problem. Communications in Statistics - Theory and Methods A5, 681-688.

Examples

n=20
n1=15
n2=10
r=0.8
xp=rnorm(n)
yp=r*xp+(1-r)*(rnorm(n))
xu=rnorm(n1)
yu=rnorm(n2)
mu=0
permb(xp,yp,xu,yu,r,mu,method="Maritz",alternative="two.sided")

[Package IncomPair version 0.1.0 Index]