должение этого номераСравнение критерия Манна-Уитни между группамиЯ решил создать новую тему.

Решение Руи Баррадаса помогло мне рассчитать Манна-Уитни для групп 1-2 и 1-3.

lst <- split(mydat, mydat$group)
lapply(lst[-1], function(DF) wilcox.test(DF$var, lst[[1]]$var, exact = FALSE))

Так что теперь я хочу получить описательную статистику. я используюlibrary:psych

describeBy(mydat$var,mydat$group)

Итак, я получаю следующий вывод

group: 1
   vars n mean   sd median trimmed  mad min max range skew kurtosis   se
X1    1 4 23.5 0.58   23.5    23.5 0.74  23  24     1    0    -2.44 0.29
-------------------------------------------------------------------------------------- 
group: 2
   vars n mean   sd median trimmed  mad min max range skew kurtosis   se
X1    1 4 23.5 0.58   23.5    23.5 0.74  23  24     1    0    -2.44 0.29
-------------------------------------------------------------------------------------- 
group: 3
   vars n mean   sd median trimmed  mad min max range skew kurtosis   se
X1    1 4 23.5 0.58   23.5    23.5 0.74  23  24     1    0    -2.44 0.29

Это неудобно. Мне нужно только для каждой группы среднее значение, sd, медиана и p-значение wilcox.test.

И.Е. я хочу эти выходные

       mean    sd     median      p-value
group1  23,5    0,58    23,5    -
group2  23,5    0,58    23,5    1
group3  23,5    0,58    23,5    1

Как я могу выполнить это?

редактировать

structure(list(`1` = structure(list(vars = 1, n = 4, mean = 23.5, 
    sd = 0.577350269189626, median = 23.5, trimmed = 23.5, mad = 0.7413, 
    min = 23, max = 24, range = 1, skew = 0, kurtosis = -2.4375, 
    se = 0.288675134594813), .Names = c("vars", "n", "mean", 
"sd", "median", "trimmed", "mad", "min", "max", "range", "skew", 
"kurtosis", "se"), row.names = "X1", class = c("psych", "describe", 
"data.frame")), `2` = structure(list(vars = 1, n = 4, mean = 23.5, 
    sd = 0.577350269189626, median = 23.5, trimmed = 23.5, mad = 0.7413, 
    min = 23, max = 24, range = 1, skew = 0, kurtosis = -2.4375, 
    se = 0.288675134594813), .Names = c("vars", "n", "mean", 
"sd", "median", "trimmed", "mad", "min", "max", "range", "skew", 
"kurtosis", "se"), row.names = "X1", class = c("psych", "describe", 
"data.frame")), `3` = structure(list(vars = 1, n = 4, mean = 23.5, 
    sd = 0.577350269189626, median = 23.5, trimmed = 23.5, mad = 0.7413, 
    min = 23, max = 24, range = 1, skew = 0, kurtosis = -2.4375, 
    se = 0.288675134594813), .Names = c("vars", "n", "mean", 
"sd", "median", "trimmed", "mad", "min", "max", "range", "skew", 
"kurtosis", "se"), row.names = "X1", class = c("psych", "describe", 
"data.frame"))), .Dim = 3L, .Dimnames = structure(list(group = c("1", 
"2", "3")), .Names = "group"), call = by.default(data = x, INDICES = group, 
    FUN = describe, type = type), class = c("psych", "describeBy"
))