modified.tstat: Kim et al.'s modified t-statistic

In samvanderpoel/partial-pairs: Statistics for Partially Matched Samples

Description

modified.tstat uses Kim et al.'s modified t-statistic to obtain a p-value for a partially matched pairs test.

Usage

modified.tstat(x, y, alternative = c("two.sided", "greater", "less"))

Arguments

x	a non-empty numeric vector containing some NA values
y	a non-empty numeric vector containing some NA values
alternative	specification of the alternative hypothesis. Takes values: "two.sided", "greater", or "less".

Details

Kim et al.’s modified t-statistic follows an approximately standard Gaussian distribution under the null hypothesis. Mathematical details are provided in [Kuan & Huang, 2013].

If proper sample size conditions are not met, then modified.tstat may exit or perform a paired or unpaired two-sample t.test, depending on the nature of the sample size issue.

If the variance of input data is close to zero, modified.tstat will return an error message.

Value

p-value associated with the hypothesis test

References

Kuan, Pei Fen, and Bo Huang. "A simple and robust method for partially matched samples using the p‐values pooling approach." Statistics in medicine 32.19 (2013): 3247-3259.

Examples

In the following, the true means are not equal:

x = rnorm(400, 0, 1)
x[sample(1:400, size=75, replace=FALSE)] = NA
y = rnorm(400, 0.4, 3)
y[sample(1:400, size=75, replace=FALSE)] = NA
modified.tstat(x, y, alternative = 'two.sided')

samvanderpoel/partial-pairs documentation built on Dec. 22, 2021, 10:14 p.m.