Mareen wrote:
Hi all,
I wondered whether anyone has some advice on a stats-related 'sanity check',
as I ran a nonparametric multivariate test (mulrank function as decribed by
R. Wilcox, 2005) on both systems, but got different results (please see
below for the system-specific outputs)! The functions I used are attached as
well. Any advice would be much appreciated! Thanks in advance for getting
back to me!
Perhaps you should show us a summary(data03$x) on both systems. It could
be that you are not working with the same data...
-p
Best wishes,
Mareen
------------------------------------
Mac:
data03<-selby2(data02, c(1,2), 3)
mulrank(3,6,data03$x)
$test.stat
[1] 0.9331133
$nu1
[1] 11.46300
$p.value
[,1]
[1,] 0.509296
$N
[1] 233
$q.hat
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 0.4940071 0.5256726 0.5176384 0.5476290 0.4690935 0.5265100
[2,] 0.5170627 0.4791950 0.5026431 0.4867843 0.4778865 0.5033497
[3,] 0.4680729 0.4944258 0.4889563 0.4505391 0.5311420 0.4726002
Win:
mulrank(3,6, data03$x)
$test.stat
[1] 1.114665
$nu1
[1] 8.155991
$p.value
[,1]
[1,] 0.3491221
$N
[1] 233
$q.hat
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 0.4940071 0.5406417 0.5236837 0.5656338 0.4771799 0.5324505
[2,] 0.5162776 0.4801895 0.5022244 0.4960745 0.4854234 0.4820737
[3,] 0.5013608 0.4920967 0.4810269 0.4482885 0.5326861 0.4871506
http://www.nabble.com/file/p23595008/Rallfun-v92.txt Rallfun-v92.txt
--
O__ ---- Peter Dalgaard Øster Farimagsgade 5, Entr.B
c/ /'_ --- Dept. of Biostatistics PO Box 2099, 1014 Cph. K
(*) \(*) -- University of Copenhagen Denmark Ph: (+45) 35327918
~~~~~~~~~~ - (p.dalga...@biostat.ku.dk) FAX: (+45) 35327907
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
