wilcox.test returning 'NA' p-value

[Ambar Amarelo]

Hi folks, sorry for this beginner question but what means a p-value = NA on
a menn-whitney test?

v1 <- c(0.022, 0.022, 0.022, 0.022, 0.022, 0.022)
v2 <- c(0.022, 0.022, 0.022, 0.022, 0.022)
wilcox.test(v1, v2, PAIRED=FALSE)

W = 15, p-value = NA

I know that there's no statistical difference between v1 and v2, so why my
p-value is not ONE ??

Can I consider p-value= NA  as ONE ??
what i do when get this type of value on wilcox.test ???


PS: sorry for my bad english.

        [[alternative HTML version deleted]]

______________________________________________
R-help@... 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.

[Ista Zahn]

Hi Ambar,

v1 and v2 are constants. they need to be variables.

try e.g.,

v1 <- c(0.021, 0.022, 0.022, 0.022, 0.022, 0.022)
v2 <- c(0.021, 0.022, 0.022, 0.022, 0.022)
wilcox.test(v1, v2, PAIRED=FALSE)

-Ista
On Thu, Nov 5, 2009 at 7:23 PM, Ambar Amarelo <ambar.amarelo@...> wrote:

> Hi folks, sorry for this beginner question but what means a p-value = NA on
> a menn-whitney test?
>
> v1 <- c(0.022, 0.022, 0.022, 0.022, 0.022, 0.022)
> v2 <- c(0.022, 0.022, 0.022, 0.022, 0.022)
> wilcox.test(v1, v2, PAIRED=FALSE)
>
> W = 15, p-value = NA
>
> I know that there's no statistical difference between v1 and v2, so why my
> p-value is not ONE ??
>
> Can I consider p-value= NA  as ONE ??
> what i do when get this type of value on wilcox.test ???
>
>
> PS: sorry for my bad english.
>
>        [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help@... 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.
>

--
Ista Zahn
Graduate student
University of Rochester
Department of Clinical and Social Psychology
http://yourpsyche.org

______________________________________________
R-help@... 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.

[Peter Ehlers]

Ambar Amarelo wrote:

> Hi folks, sorry for this beginner question but what means a p-value = NA on
> a menn-whitney test?
>
> v1 <- c(0.022, 0.022, 0.022, 0.022, 0.022, 0.022)
> v2 <- c(0.022, 0.022, 0.022, 0.022, 0.022)
> wilcox.test(v1, v2, PAIRED=FALSE)
>
> W = 15, p-value = NA
>
> I know that there's no statistical difference between v1 and v2, so why my
> p-value is not ONE ??
>
> Can I consider p-value= NA  as ONE ??
> what i do when get this type of value on wilcox.test ???
>
>
> PS: sorry for my bad english.
>
> [[alternative HTML version deleted]]
>

Did you read the help page? If there are ties, then a
Normal approximation is used. So how would you do a
Normal approximation? Try it "by hand" - you have a
small enough amount of data.

Try also a t.test() and see what message it gives you.

BTW, there is no PAIRED argument to wilcox.test().

--
Peter Ehlers
University of Calgary

______________________________________________
R-help@... 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.

[Ambar Amarelo]

> Did you read the help page?

Hi  Ista Zahn and Peter, thanks to response  me.
Well, some times on my data I have two samples with exactly values.
(not small samples)

If the mann-whitney test (wilcox.test) return a p-value = NA
Can I consider that my two sample have a p-value = 1 ???

Or there's another's cases that this test will return NA and the samples are
not the same ???

BTW, on the help ?wilcox.test

    "...Otherwise, if both 'x' and 'y' are given and 'paired' is 'FALSE',
     a Wilcoxon rank sum test (equivalent to the Mann-Whitney test: see
     the Note) is carried out..."

This is why i set the parameter PAIRED=FALSE on wicox.test.











2009/11/5 Peter Ehlers <ehlers@...>

> Ambar Amarelo wrote:
>
>> Hi folks, sorry for this beginner question but what means a p-value = NA
>> on
>> a menn-whitney test?
>>
>> v1 <- c(0.022, 0.022, 0.022, 0.022, 0.022, 0.022)
>> v2 <- c(0.022, 0.022, 0.022, 0.022, 0.022)
>> wilcox.test(v1, v2, PAIRED=FALSE)
>>
>> W = 15, p-value = NA
>>
>> I know that there's no statistical difference between v1 and v2, so why my
>> p-value is not ONE ??
>>
>> Can I consider p-value= NA  as ONE ??
>> what i do when get this type of value on wilcox.test ???
>>
>>
>> PS: sorry for my bad english.
>>
>>        [[alternative HTML version deleted]]
>>
>>
> Did you read the help page? If there are ties, then a
> Normal approximation is used. So how would you do a
> Normal approximation? Try it "by hand" - you have a
> small enough amount of data.
>
> Try also a t.test() and see what message it gives you.
>
> BTW, there is no PAIRED argument to wilcox.test().
>
> --
> Peter Ehlers
> University of Calgary
>

--
Antigos Virus do mau,
Transformem este Nick decadente, em
Le_AmBArrrrrrr, o AmBAr Amareloooooo !

        [[alternative HTML version deleted]]

______________________________________________
R-help@... 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.

[Peter Ehlers]

Ambar Amarelo wrote:
>> Did you read the help page?
>
> Hi  Ista Zahn and Peter, thanks to response  me.
> Well, some times on my data I have two samples with exactly values.
> (not small samples)
>
> If the mann-whitney test (wilcox.test) return a p-value = NA
> Can I consider that my two sample have a p-value = 1 ???
>
You probably can, but so what? I realize that p-value obsession
is rampant, but sometimes it's actually best to study the data.

You can always check for the sort of pathological case you're
interested in before passing data to wilcox.test(), using, say,
diff() or rle() on your vectors.

You might also find the error messages in coin::wilcox_test()
more informative.

> Or there's another's cases that this test will return NA and the samples are
> not the same ???
>
> BTW, on the help ?wilcox.test
>
>     "...Otherwise, if both 'x' and 'y' are given and 'paired' is 'FALSE',
>      a Wilcoxon rank sum test (equivalent to the Mann-Whitney test: see
>      the Note) is carried out..."
>
> This is why i set the parameter PAIRED=FALSE on wicox.test.
>

Yes, but my point was that R is case-sensitive;
"paired" != "PAIRED"

  -Peter Ehlers

>
>
>
>
>
>
>
>
>
>
> 2009/11/5 Peter Ehlers <ehlers@...>
>
>
>> Ambar Amarelo wrote:
>>
>>> Hi folks, sorry for this beginner question but what means a p-value = NA
>>> on
>>> a menn-whitney test?
>>>
>>> v1 <- c(0.022, 0.022, 0.022, 0.022, 0.022, 0.022)
>>> v2 <- c(0.022, 0.022, 0.022, 0.022, 0.022)
>>> wilcox.test(v1, v2, PAIRED=FALSE)
>>>
>>> W = 15, p-value = NA
>>>
>>> I know that there's no statistical difference between v1 and v2, so why my
>>> p-value is not ONE ??
>>>
>>> Can I consider p-value= NA  as ONE ??
>>> what i do when get this type of value on wilcox.test ???
>>>
>>>
>>> PS: sorry for my bad english.
>>>
>>>        [[alternative HTML version deleted]]
>>>
>>>
>> Did you read the help page? If there are ties, then a
>> Normal approximation is used. So how would you do a
>> Normal approximation? Try it "by hand" - you have a
>> small enough amount of data.
>>
>> Try also a t.test() and see what message it gives you.
>>
>> BTW, there is no PAIRED argument to wilcox.test().
>>
>> --
>> Peter Ehlers
>> University of Calgary
>>
>
>
>

______________________________________________
R-help@... 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.