У меня есть генетические данные для SNP, которые были разделены на 5 квантилей. Я хочу найти медиану этих квантилей для каждого SNP (т.е. для каждого человека).
Я использовал эту команду для создания столбца для медианных значений:
data$median<-apply(data[,2:181],1, median, na.rm=TRUE)
Тогда я хотел посчитать, сколько случаев и контролей у меня есть для каждого из моих фенотипов, но похоже, что он неправильно вычисляет медиану. Моя команда выглядит следующим образом:
table(data$anyMI, data$median)
Вывод показывает:
1 1.5 2 2.5 3 3.5 4 4.5 5
0 2044 62 7470 221 11163 248 8389 74 1659
1 102 3 357 11 557 21 404 2 85
Я не уверен, почему я получаю половинные значения, когда они должны быть только 1-5, целые числа.
Что здесь не так и почему он показывает половинные значения?