I suggest posting this on Statalist (http://www.statalist.org/); somersd
is a user-written package and the author, Roger Newson, appears to
follow the Statalist discussions - you could also email him directly -
his email address is at the bottom of the help file
Rich
On 5/10/16 6:00 AM, Andrew Salmon wrote:
> Hi all
> am just investigating the finer details of the Mann Whitney test and
> trying to calculate confidence intervals for the porder statistic. I
> have the following dummy data with var2 as the grouping variable
>
> using STATA's ranksum with porder option gives me a value of 0.594
> (for group = 0 over group = 1). Doing the same calculation using the
> somersd package with transf(c) option gives me a confidence interval.
> ( I used vreverse to generate rgroup)
>
> | Jackknife
> rgroup | Coef. Std. Err. z P>|z| [95% Conf. Interval]
> -------------+----------------------------------------------------------------
> var1 | .59375 .1594069 3.72 0.000 .2813182 .9061818
> -----------------------------------------------------------------------------
> but with a p value that seems rather low and doesn't seem to make
> sense with the confidence interval, since I would have thought the
> null of probability/proportion would be 0.5
>
> Any thoughts?
>
> thanks
> Andy
>
> var1 var2
> 1 0
> 5 0
> 9 0
> 14 0
> 21 0
> 29 0
> 39 0
> 49 0
> 3 1
> 4.5 1
> 7 1
> 12 1
> 19 1
> 20.5 1
> 26 1
> 30 1
You may leave the list at any time by sending the command
SIGNOFF allstat
to [log in to unmask], leaving the subject line blank.
|