Description Usage Format Details Source References Examples
The mercury level in blood, the proportion of cells with abnormalities, and the proportion of cells with chromosome aberrations in consumers of mercurycontaminated fish and a control group.
1 
A data frame with 39 observations on 4 variables.
group
a factor with levels "control"
and "exposed"
.
mercury
mercury level in blood.
abnormal
the proportion of cells with structural abnormalities.
ccells
the proportion of C_u cells, i.e., cells with asymmetrical or incompletesymmetrical chromosome aberrations.
Control subjects ("control"
) and subjects who ate contaminated fish for
more than three years ("exposed"
) are under study.
Rosenbaum (1994) proposed a coherence criterion defining a partial ordering, i.e., an observation is smaller than another when all responses are smaller, and a score reflecting the “ranking” is attached to each observation. The corresponding partially ordered set (POSET) test can be used to test if the distribution of the scores differ between the groups. Alternatively, a multivariate test can be applied.
Skerfving, S., Hansson, K., Mangs, C., Lindsten, J. and Ryman, N. (1974). Methylmercuryinduced chromosome damage in men. Environmental Research 7(1), 83–98. doi: 10.1016/00139351(74)900784
Hothorn, T., Hornik, K., van de Wiel, M. A. and Zeileis, A. (2006). A Lego system for conditional inference. The American Statistician 60(3), 257–263. doi: 10.1198/000313006X118430
Rosenbaum, P. R. (1994). Coherence in observational studies. Biometrics 50(2), 368–374. doi: 10.2307/2533380
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45  ## Coherence criterion
coherence < function(data) {
x < as.matrix(data)
matrix(apply(x, 1, function(y)
sum(colSums(t(x) < y) == ncol(x)) 
sum(colSums(t(x) > y) == ncol(x))), ncol = 1)
}
## Asymptotic POSET test
poset < independence_test(mercury + abnormal + ccells ~ group,
data = mercuryfish, ytrafo = coherence)
## Linear statistic (T in the notation of Rosenbaum, 1994)
statistic(poset, type = "linear")
## Expectation
expectation(poset)
## Variance
## Note: typo in Rosenbaum (1994, p. 371, Sec. 2, last paragraph)
variance(poset)
## Standardized statistic
statistic(poset)
## Pvalue
pvalue(poset)
## Exact POSET test
independence_test(mercury + abnormal + ccells ~ group,
data = mercuryfish, ytrafo = coherence,
distribution = "exact")
## Asymptotic multivariate test
mvtest < independence_test(mercury + abnormal + ccells ~ group,
data = mercuryfish)
## Global pvalue
pvalue(mvtest)
## Singlestep adjusted pvalues
pvalue(mvtest, method = "singlestep")
## Stepdown adjusted pvalues
pvalue(mvtest, method = "stepdown")

Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.