Ties, no ties, ties, no ties ... [Nonparametrics]
Thanks for your very elaborate answer!
❝ Right. But if ties are considered, you are stirring up a hornets' nest!
I see

❝ My quick and dirty solution:
❝ PKt$pdiff <- (PKt$P1 - PKt$P2)/2
❝ PKt$sequence <- relevel(PKt$sequence, ref="TR")
❝ ...
Interesting: "relevel". Was not aware of that handy tool.
But astonishing the result: A change in order gives a different lower bound!
❝ Hhm; don't use library(exactRankTests)
?
Seems a good recommendation. But it is claimed frequently that Rusers should use it in case of ties. Moreover coin has a much bigger footprint.
❝ ... they refer to Hollander and Wolfe. In my edition the word 'tie' is not even mentioned in the subject index...
Very strange. Everyone writing about rank methods is citing this.
❝ With my mid-1990 implementation in Rocky Mountain BASIC (not
❝ corrected for ties) I get:
❝ PE 90.26% CI
❝ 0.25 -0.25 +0.75
Not 'correcting' for ties to obtain the CI is exactly what Charles Geyer recommends. With respect to the Wilcoxon sign rank test he states:
"First, neither ties nor zeros should make any difference in calculating point estimators or confidence intervals.
- The point estimate is the median of the Walsh averages ...
- The end points of the confidence interval are k in from each end of the sorted Walsh averages ...
- Ties and zeros affect only hypothesis tests.
And this may be transferable to the rank-sum test

If one follows his code in the interval section of the Wilcoxon rank sum test one gets your Rocky results.
❝ Now for StatXact 6.1. using the coding according to Senn (2002) ...
What exactly did you here? Senn's Basic estimator? Then a test of sequence group differences would give an estimate of the period effect. Or am I wrong? See the expected period differences below.
❝ ... Open questions:
❝ • Signs... (Sigh!) I've checked the coding in all programs numerous times,
❝
but don't know what's going on here. Fishbrain.
I came across with: Sometimes the smaller group is used as the first, sometimes the coding order, sometimes the character order, sometimes the first occurred level.
My
![[image]](img/uploaded/image2.gif)
❝ I guess the reason might be that you are comparing period differences rather than treatment differences.
The reason for using the period differences rather then the treatment diffs comes from the paper of Hauschke. It is based on the expected values
sequence P1-P2
TR µ(T-R) + p1-p2
RT -µ(T-R) + p1-p2
Thus a test of the differences by sequence group gets rid of the period effects. This is not the case with the treatment diffs.
❝ ... I tried ties.method="mid-ranks"
and ties.method="average-scores"
and got:
❝ Warning: In independence_test.IndependenceProblem(object, teststat = "scalar", :
❝ additional arguments ties.method will be ignored
I also observed this. Although mentioned on the help page it seems that the ties.method is not available in the wilcoxon_test function and other location tests of package coin. ties.method is only accepted by the function normal_test() which performs the "Exact Normal Quantile (van der Waerden) Test".
❝ If I exclude subjects with a zero difference ...
Called by Charles Geyer "The Zero fudge"

❝ ... You gave me a nice excuse not to spend the afternoon struggling with my income tax return.
If it is that horrible as here in Germanien I did you the favor very willingly

Regards,
Detlew
Complete thread:
- tmax in case of ties: R vs. R vs. SAS d_labes 2010-09-10 12:05 [Nonparametrics]
- tmax in case of ties: StatXact, Phoenix, and... Helmut 2010-09-10 18:28
- Ties, no ties, ties, no ties ...d_labes 2010-09-13 13:21
- simulation martin 2010-09-14 08:33
- simulants of the world unite d_labes 2010-09-15 10:14
- simulants of the world unite martin 2010-09-15 14:53
- simulants of the world unite Helmut 2010-09-15 16:17
- simul ants with no ties d_labes 2010-09-16 11:19
- simul ants with no ties Helmut 2010-09-16 14:01
- two different approaches martin 2010-09-16 15:35
- simul ants with no ties Helmut 2010-09-16 14:01
- simul ants with no ties d_labes 2010-09-16 11:19
- simulants of the world unite Helmut 2010-09-15 16:17
- simulants of the world unite martin 2010-09-15 14:53
- simulants of the world unite d_labes 2010-09-15 10:14
- simulation martin 2010-09-14 08:33
- Ties, no ties, ties, no ties ...d_labes 2010-09-13 13:21
- evaluation of tmax: use of relative effects? martin 2010-09-16 18:04
- Not positive about that Helmut 2010-09-16 19:52
- Not positive about that martin 2010-09-16 21:00
- Stupidity Helmut 2010-09-17 13:10
- Not positive about that martin 2010-09-16 21:00
- Not positive about that Helmut 2010-09-16 19:52
- tmax in case of ties: R vs. R vs. SAS Jack 2010-09-20 14:03
- R packages Helmut 2010-09-20 14:20
- R packages code d_labes 2010-09-27 11:27
- R packages Jack 2010-09-27 16:14
- Is Exact exact? d_labes 2010-09-27 09:53
- R packages Helmut 2010-09-20 14:20
- tmax in case of ties: StatXact, Phoenix, and... Helmut 2010-09-10 18:28