Learn R Programming

statpsych (version 1.7.0)

sim.ci.spear: Simulates confidence interval coverage probability for a Spearman correlation

Description

Performs a computer simulation of confidence interval performance for a Spearman correlation. Sample data can be generated from bivariate population distributions with five different marginal distributions. All distributions are scaled to have standard deviations of 1.0. Bivariate random data with specified marginal skewness and kurtosis are generated using the unonr function in the mnonr package.

Usage

sim.ci.spear(alpha, n, cor, dist1, dist2, rep)

Value

Returns a 1-row matrix. The columns are:

  • Coverage - probability of confidence interval including population correlation

  • Lower Error - probability of lower limit greater than population correlation

  • Upper Error - probability of upper limit less than population correlation

  • Ave CI Width - average confidence interval width

Arguments

alpha

alpha level for 1-alpha confidence

n

sample size

cor

population Spearman correlation

dist1

type of distribution for variable 1 (1, 2, 3, 4, or 5)

dist2

type of distribution for variable 2 (1, 2, 3, 4, or 5)

  • 1 = Gaussian (skewness = 0 and excess kurtosis = 0)

  • 2 = platykurtic (skewness = 0 and excess kurtosis = -1.2)

  • 3 = leptokurtic (skewness = 0 and excess kurtosis = 6)

  • 4 = moderate skew (skewness = 1 and excess kurtosis = 1.5)

  • 5 = large skew (skewness = 2 and excess kurtosis = 6)

rep

number of Monte Carlo samples

Examples

Run this code
sim.ci.spear(.05, 30, .7, 4, 5, 1000)

# Should return (within sampling error):
#      Coverage Lower Error Upper Error Ave CI Width
# [1,]  0.96235     0.01255      0.0251    0.4257299


Run the code above in your browser using DataLab