[R] Mann-Whitney by group

Oxenstierna david.chertudi at gmail.com
Fri Jul 6 22:34:29 CEST 2012


Hi David,

Thank you for the insight:  I could have sworn I added a picture of the
data, but providing the actual data is worlds easier to deal with, I'm sure. 
I've never used dput(), so I entered it using the dataframe in question as
the object, and I've pasted the results below.

Essentially, I would like to run the two-sample independent test, comparing
Group 5 (CD8.14 through CD8.17) to Group 6 (CD8.18 through CD8.21). So, for
A, test Group 5 agains Group 6, for B, test Group 5 against Group 6, and so
on.  I'm not going to muddy the waters by telling you what I've tried;
suffice it to say that I'm looking for insights into how to structure R
commands to compare groups of data of this format.

Many thanks in advance,

David


structure(list(Gene = structure(c(1L, 12L, 19L, 20L, 21L, 22L, 
23L, 24L, 25L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 13L, 
14L, 15L, 16L, 17L, 18L), .Label = c("CD8.1", "CD8.10", "CD8.11", 
"CD8.12", "CD8.13", "CD8.14", "CD8.15", "CD8.16", "CD8.17", "CD8.18", 
"CD8.19", "CD8.2", "CD8.20", "CD8.21", "CD8.22", "CD8.23", "CD8.24", 
"CD8.25", "CD8.3", "CD8.4", "CD8.5", "CD8.6", "CD8.7", "CD8.8", 
"CD8.9"), class = "factor"), Group = structure(c(8L, 8L, 8L, 
9L, 9L, 9L, 10L, 10L, 10L, 11L, 11L, 11L, 11L, 12L, 12L, 12L, 
12L, 13L, 13L, 13L, 13L, 14L, 14L, 14L, 14L), .Label = c("Fabbf Ova CD40", 
"Fabbf Ova MHC2", "Fabbf Ova WT", "Fabbf WT", "Naïve CD40", "Naïve MHC II", 
"Naïve WT", "1", "2", "3", "4", "5", "6", "7"), class = "factor"), 
    A = c(19.4701946749544, 0.679440926463348, 0.69035683372563, 
    0.347105466158261, 0.435480792190284, 0.338699910286907, 
    0.651378057031152, 0.707065053752258, 0.685244609506316, 
    0.816673858871597, 0.597009097584509, 0.592331304482431, 
    0.709359033358704, 0.628406759227531, 0.78158729467231, 1.0422377526669, 
    0.61560003251142, 0.463755016733183, 0.419700860701392,
0.380946898502731, 
    0.41156961153081, 0.471790683365624, 0.552773224145722,
0.466787799928649, 
    0.767624372542755), B = c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 
    1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1), C = c(0.914649979331863, 
    1.443086801592, 0.928280641141244, 0.467498974059775, 0.668742025741347, 
    0.568417298005388, 0.778907650835673, 0.778765181169635, 
    1.01103488277517, 0.674133483128923, 0.830400022230133,
0.652687178870627, 
    0.746893950266518, 0.765498308522646, 1.01273201749333,
1.00309616129672, 
    0.63015975419947, 0.615905247119739, 0.615687625199691,
0.503136087800137, 
    0.566164026974035, 0.631519467967541, 0.822126705285366, 
    0.543873075815645, 1.00596108625425), D = c(13.3272657341526, 
    3.08914950309865, 1.71836820240434, 0.723301573710509, 1.21014411624732, 
    1.92899377364865, 1.80280408189187, 2.25057819266424, 2.23876060313374, 
    1.30849425313072, 1.58782967140617, 1.19199809794126, 1.64151140806787, 
    0.241017500596534, 0.364896032519483, 0.322953808735804, 
    0.2052110581509, 0.927601295331376, 0.808910781520832,
0.538033121081646, 
    0.655348783504307, 0.564449549672088, 0.521729926793001, 
    0.414305517285192, 0.507084483980948), E = c(56.2830291897158, 
    9.76091939190267, 4.80922410182105, 17.0056576949022, 20.851046177766, 
    17.9057247086369, 5.93332779160845, 4.73058157592946, 5.59155211460608, 
    9.67484467290805, 5.92374864612388, 7.12393623733123, 5.33576126730867, 
    10.3943422629275, 10.8732527705049, 12.4861085370674, 12.0918705721064, 
    13.3210661695018, 10.9410344557684, 15.1298307761675, 13.0708078246191, 
    9.4445293976312, 6.94340249514349, 5.07888688780375, 8.33846787814466
    ), F = c(15.0459568981729, 21.6362955612539, 9.66673955488981, 
    27.2276698483913, 18.1090094072926, 20.0952712980862, 24.9249499974856, 
    23.5540183530194, 29.6638363657906, 28.9779309040733, 42.0402820641407, 
    33.8068160394092, 51.7299064374737, 37.8306751403421, 43.1955470199259, 
    45.5125262939585, 40.3109474523637, 23.6341894633273, 23.9721353180788, 
    20.4920649252818, 24.8898447627354, 34.2686409607416, 31.3815198841165, 
    31.4947528368753, 43.2686436885025), G = c(15.9703031418086, 
    5.30495997585743, 3.07594974529074, 8.28703732907722, 10.9437825143868, 
    6.91196232523896, 2.97808148581742, 2.6386825521864, 2.2415006913088, 
    5.00747306438661, 2.65238188782831, 3.18277515130905, 3.14638620532385, 
    2.5149505923191, 2.48862112414046, 2.97170069886913, 1.91643165326171, 
    7.50682774199005, 5.39102206185423, 5.2498453524987, 7.17519969844757, 
    2.66448841457179, 2.78444235996995, 1.62286520735228, 2.48760726398266
    ), H = c(-1, -1, -1, 0.0036561124481055, 0.010723007432761, 
    0.0196616746380801, 0.0371046164124276, -1, 3.27378673314144e-05, 
    -1, -1, -1, -1, 0.00816215783906802, 21.5424904701651, 9.57900616157724, 
    0.00735942043489242, 1.50346040901698, 0.0909450037365435, 
    1.49237001404701e-05, 0.000775741472561218, -1, 17.9023582944659, 
    0.0176891314806093, 33.4326253626981)), .Names = c("Gene", 
"Group", "A", "B", "C", "D", "E", "F", "G", "H"), row.names = c(NA, 
25L), class = "data.frame")

--
View this message in context: http://r.789695.n4.nabble.com/Mann-Whitney-by-group-tp4635618p4635667.html
Sent from the R help mailing list archive at Nabble.com.



More information about the R-help mailing list