yuend: Paired samples robust t-tests.

Description Usage Arguments Details Value References See Also Examples

View source: R/yuend.R

Description

The function yuend performs Yuen's test on trimmed means for dependent samples. Dqcomhd compares the quantiles of the marginal distributions associated with two dependent groups via hd estimator. Tied values are allowed.

Usage

1
2
yuend(x, y, tr = 0.2)
Dqcomhd(x, y, q = c(1:9)/10, nboot = 1000, na.rm = TRUE)

Arguments

x

an numeric vector of data values (e.g. for time 1).

y

an numeric vector of data values (e.g. for time 2).

tr

trim level for the means.

q

quantiles to be compared.

nboot

number of bootstrap samples.

na.rm

whether missing values should be removed.

Details

The test statistic is a paired samples generalization of Yuen's independent samples t-test on trimmed means.

Value

yuend returns an object of class "yuen" containing:

test

value of the test statistic (t-statistic)

p.value

p-value

conf.int

confidence interval

df

degress of freedom

diff

trimmed mean difference

call

function call

Dqcomhd returns an object of class "robtab" containing:

partable

parameter table

References

Wilcox, R. (2012). Introduction to Robust Estimation and Hypothesis Testing (3rd ed.). Elsevier.

See Also

yuen, qcomhd

Examples

1
2
3
4
5
6
7
## Cholesterol data from Wilcox (2012, p. 197)
before <- c(190, 210, 300,240, 280, 170, 280, 250, 240, 220)
after <- c(210, 210, 340, 190, 260, 180, 200, 220, 230, 200)
yuend(before, after)

set.seed(123)
Dqcomhd(before, after, nboot = 200, q = c(0.25, 0.5, 0.75))

Example output

Call:
yuend(x = before, y = after)

Test statistic: 1.9886 (df = 5), p-value = 0.10343

Trimmed mean difference:  28.33333 
95 percent confidence interval:
-8.2918     64.9585 

Explanatory measure of effect size: 0.52 

Call:
Dqcomhd(x = before, y = after, q = c(0.25, 0.5, 0.75), nboot = 200)

Parameter table: 
     q n1 n2     est1     est2 est1-est.2   ci.low   ci.up p.crit p.value
1 0.25 10 10 205.9676 196.4354     9.5322 -10.2331 39.1004 0.0500    0.53
2 0.50 10 10 238.8929 211.2158    27.6770  -2.2910 56.0509 0.0167    0.07
3 0.75 10 10 271.0781 245.3351    25.7431 -20.9334 51.7463 0.0250    0.34

WRS2 documentation built on May 31, 2017, 2:07 a.m.