PairedData (version 0.5)

rpaired.contaminated: Simulate paired samples

Description

Simulate paired data with a given correlation (Kendall's tau=(2/pi)arcsine(r)) and marginals being contaminated normal distributions: (1-e)F(x)+eF(x/K) where F is the cumulative standard normal distribution, e the percentage of contamination and K a scale parameter.

Usage

rpaired.contaminated(n,d1=c(0.1,10,1),d2=c(0.1,10,1),r=0.5)

Arguments

n

Sample size

d1

Vector of 3 parameters for the first contaminated normal distribution

d2

Vector of 3 parameters for the second contaminated normal distribution

r

Correlation

Value

• A list with two components:
• x
  First sample
• y
  Second sample

References

P.M. Grambsch. Simple robust tests for scale differences in paired data. Biometrika, 81, 359-372, 1994.

See Also

rpaired.gld

Examples

rpaired.contaminated(n=30,r=0.25)