RVAideMemoire (version 0.9-11)

cor.multcomp: Comparison of several Pearson's linear correlation coefficients

Description

Perform comparisons of several Pearson's linear correlation coefficients. If no difference, the function returns the common correlation coefficient, its confidence interval and test for its equality to a given value. If difference is significative, the functions performs pairwise comparisons between coefficients.

Usage

cor.multcomp(var1, var2, fact, alpha = 0.05, conf.level = 0.95, theo = 0, p.method = "fdr")

Arguments

var1

numeric vector (first variable).

var2

numeric vector (second variable).

fact

factor (groups).

alpha

significance level.

conf.level

confidence level.

theo

theoretical coefficient.

p.method

method for p-values correction. See help of the p.adjust() function.

Value

• levels
  number of levels of fact.
• conf.level
  confidence level.
• alpha
  significance level.
• coeffs
  observed correlation coefficients.
• statistic.comp
  test statistics for equality of all coefficients.
• p.value.comp
  p-value for equality of all coefficients.
• comp
  data frame of results for equality of all coefficients.
• r.comm
  common correlation coefficient and confidence interval.
• r.theo
  theoretical coefficient.
• statistic.comm
  test statistics for equality of the common coefficient to the theoretical value.
• p.value.comm
  p-value for equality of the common coefficient to the theoretical value.
• comm
  data frame for equality of the common coefficient to the theoretical value.
• p.adjust.method
  method for p-values correction.
• p.multcomp
  data frame of results of pairwise comparisons.

See Also

cor.test

Examples

var1 <- c(1:15+rnorm(15,0,4),1:15+rnorm(15,0,1),1:15+rnorm(15,0,8))
var2 <- c(-1:-15+rnorm(15,0,4),1:15+rnorm(15,0,1),1:15+rnorm(15,0,8))
fact <- gl(3,15,labels=LETTERS[1:3])
cor.multcomp(var1,var2,fact)

var3 <- c(1:15+rnorm(15,0,1),1:15+rnorm(15,0,3),1:15+rnorm(15,0,2))
cor.multcomp(var1,var3,fact)