[R] MAC OSX vs Win XP: Different stats test results!

2009-05-18 Thread Mareen

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!

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 

-- 
View this message in context: 
http://www.nabble.com/MAC-OSX-vs-Win-XP%3A-Different-stats-test-results%21-tp23595008p23595008.html
Sent from the R help mailing list archive at Nabble.com.

__
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.


Re: [R] MAC OSX vs Win XP: Different stats test results!

2009-05-18 Thread Peter Dalgaard

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.


Re: [R] MAC OSX vs Win XP: Different stats test results!

2009-05-18 Thread Mike Lawrence
If mulrank does any sort of random number generation or non-exhaustive
randomization, you should set the seed of the random number generator
first:

set.seed(1)
mulrank(3,6,data03$x)


On Mon, May 18, 2009 at 7:37 AM, 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!
>
> 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
>
> --
> View this message in context: 
> http://www.nabble.com/MAC-OSX-vs-Win-XP%3A-Different-stats-test-results%21-tp23595008p23595008.html
> Sent from the R help mailing list archive at Nabble.com.
>
> __
> 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.
>



-- 
Mike Lawrence
Graduate Student
Department of Psychology
Dalhousie University

Looking to arrange a meeting? Check my public calendar:
http://tr.im/mikes_public_calendar

~ Certainty is folly... I think. ~

__
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.


Re: [R] MAC OSX vs Win XP: Different stats test results!

2009-05-18 Thread Mareen Weber

Hello Peter,
thanks for your quick response - you have hinted into the right direction. I 
triple-checked my data and - how embarrassing - I actually did have different 
data matrices. Due to a faulty Win Excel macro that sorted data into groups 
after data collection, I ended up with different data in my 3 groups in the Win 
"test" ... All good now :-). Thanks again! 

Best wishes,
Mareen





__
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.


Re: [R] MAC OSX vs Win XP: Different stats test results!

2009-05-18 Thread Mike Lawrence
A good example of why manipulating data in external programs is
dangerous. See a discussion of this at
http://tolstoy.newcastle.edu.au/R/e6/help/09/05/13697.html

On Mon, May 18, 2009 at 10:56 AM, Mareen Weber  wrote:
>
> Hello Peter,
> thanks for your quick response - you have hinted into the right direction. I 
> triple-checked my data and - how embarrassing - I actually did have different 
> data matrices. Due to a faulty Win Excel macro that sorted data into groups 
> after data collection, I ended up with different data in my 3 groups in the 
> Win "test" ... All good now :-). Thanks again!
>
> Best wishes,
> Mareen
>
>
>
>
>
> __
> 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.
>



-- 
Mike Lawrence
Graduate Student
Department of Psychology
Dalhousie University

Looking to arrange a meeting? Check my public calendar:
http://tr.im/mikes_public_calendar

~ Certainty is folly... I think. ~

__
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.