Statistics for 2 by K table
This function calculates one-to-others and maximum accumulated chi-squared statistics for a 2 by K contingency table.
s2k(y1,y2)
y1 |
a vector containing the first row of a 2 by K contingency table |
y2 |
a vector containing the second row of a 2 by K contingency table |
The returned value is a list containing:
x2a |
the one-to-other chisquare |
x2b |
the maximum accumulated chisquare |
col1 |
the column index for x2a |
col2 |
the column index for x2b |
p |
the corresponding p value |
Hirotsu C, Aoki S, Inada T, Kitao Y (2001) An exact test for the association between the disease and alleles at highly polymorphic loci with particular interest in the haplotype analysis. Biometrics 57:769-778
The lengths of y1 and y2 should be the same
Chihiro Hirotsu, Jing Hua Zhao
## Not run: # an example from Mike Neale # termed 'ugly' contingency table by Patrick Sullivan y1 <- c(2,15,16,35,132,30,25,7,12,24,10,10,0) y2 <- c(0, 6,31,49,120,27,15,8,14,25, 3, 9,3) result <- s2k(y1,y2) ## End(Not run)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.