[R] NA in wilcox.test

From: Michelle Emily DePrenger-Levin <michelledpl_at_gmail.com>
Date: Sat, 05 May 2007 17:01:37 -0600


Hello,

I'm trying to compare the allozyme data from two taxa. I have several columns of data (19 loci) for each species. I want to do a Mann-Whitney U-test or the wilcox.test (two sample Wilcoxon). When I try to run my code (the first two columns are 1:name of the species, 2:name of individual) I get the error message:

            "Error in wilcox.test.default(CaScSc, CaScCo, alternative = " two.sided", :

        'x' must be numeric"

I do have several NAs in the data which is all I can figure that is non-numeric.

Any suggestions as to the problem? Is it a problem with having several columns/sets of info for each individual?

Thanks for any help anyone can give (I've also used Arlequin and GDA but want non-parametric tests)

Michelle DePrenger-Levin

My code is:

scirconv2 = read.csv("CaScSc070420_2.csv", na.strings="?")

CaScSc = scirconv2[1:250,3:21]
CaScCo = scirconv2[251:475,3:21]
ScCoMWU = wilcox.test(CaScSc, CaScCo, alternative = "two.sided", mu = 0.5, paired = FALSE, exact = NULL, correct = TRUE,

            conf.int = TRUE, conf.level = 0.95)

some of the data:

    ADH TPI1 TPI2 SOD DIA1 MNR1DIA2 MNR2DIA3 ME AAT1 AAT2 G3PDH SDH SDH2 PGI2 PGD PGM2 MDH1 MDH3 IDH2
251 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 4 1 1 4
252 NA NA NA NA NA NA 1 NA NA 2 1 1 NA 2 NA NA NA NA NA
253 1 1 1 1 1 1 2 1 1 2 1 1 1 2 1 4 1 2 4
254 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 4 1 1 4
255 1 1 1 1 1 1 2 1 1 2 1 1 1 2 1 4 1 2 4

        [[alternative HTML version deleted]]



R-help_at_stat.math.ethz.ch 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. Received on Sat 05 May 2007 - 23:06:29 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Sun 06 May 2007 - 09:31:39 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.