RVAideMemoire (version 0.9-7)

cochran.qtest: Cochran Q test

Description

Perform the Cochran Q test for unreplicated randomized block design experiments whith a binary response variable and paired data. If the p-value of the test is significative, the function performs pairwise comparisons using Wilcoxon sign test.

Usage

cochran.qtest(formula, data=NULL, alpha = 0.05, p.method = "fdr")

Arguments

formula

a formula of the form a ~ b | c, where a, b and c give the data values and corresponding groups and blocks, respectively.

data

an optional data frame containing the variables in the formula formula. By default the variables are taken from environment(formula).

alpha

significance level to compute pairwise comparisons.

p.method

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

Value

• data
  variables name.
• alpha
  significance level.
• Q
  test statistics.
• p.value
  p-value of the global test.
• tab.test
  data frame of global test result.
• p.method
  method for p-values correction.
• p.multcomp
  corrected p-values.
• multcomp
  data frame of pairwise comparisons result.

Examples

response<-c(0,1,1,0,0,1,0,1,1,1,1,1,0,0,1,1,0,1,0,1,1,0,0,1,0,1,1,0,0,1)
factor<-factor(rep(LETTERS[1:3],10))
block<-factor(rep(letters[1:10],each=3))
cochran.qtest(response~factor|block)