Дискретные бины с помощью cut()
Я хочу отобразить данные [используя xyplot() решетки] в моем фрейме данных age.model на основе дискретных ячеек столбца StartAge.
Я использую следующий код:
# set up boundaries for intervals/bins
breaks <- c(0,3,4,5,6,8,13,15,17,18,19,20,22)
# specify interval/bin labels
labels <- c("<3", "3-4)", "4-5)","5-6)", "6-8)","8-13)", "13-15)","15-17)","17-18)","18-19)","19-20)",">=20")
# bucketing data points into bins
bins <- cut(age.model$StartAge, breaks, include.lowest = T, right=FALSE, labels=labels)
# inspect bins
summary(bins)
В первом аргументе функции cut() я указал столбец, по которому хочу дискретизировать. Однако возвращаемый коэффициент не включает весь DF. Как я могу это сделать?
Воспроизводимо с использованием dput:
structure(list(Height = c(0.207224416925809, -1.19429150954007,
0.0247585682642494, 0.023546515879641, 1.51423735121426, -1.09376538778425,
-0.125209484617016, -0.63639210765747, 0.305071992864995, -0.422021082477656
), Weight = c(-0.366133564723644, -1.06969961340686, -0.0793604259237282,
-0.708230200986797, 1.71593234004357, -0.685215310472794, -1.20353653394014,
-0.490399232488568, 0.742874184424376, -0.331519044995803), Training = c(19,
27, 27, 24, 35, 23, 15, 14, 47, 7), StartAge = c(13, 19, 20,
20, 14, 2, 8, 4, 17, 18)), row.names = c("1", "2", "3", "4",
"5", "6", "7", "8", "9", "10"), class = "data.frame")
Добавил, надеюсь поможет?
Насколько я понял, вы можете получить полный фрейм данных с помощью age.model$bins <- bins. В противном случае укажите ожидаемый результат.
Ответы 2
Чтобы добавить бины в фрейм данных, просто повлияйте на него в новом столбце:
age.model$bins <- bins
Если вы используете xyplot для проводить исследования ваших данных, рассмотрите возможность использования equal.count() или shingle() в своем коде. Получая (бестолковое) удовольствие от ваших данных, приблизительная линейная зависимость между весом и ростом, по-видимому, не выполняется для нижних StartAge ячеек, как показано в первом примере.
# Starting with data in age.model
library(lattice)
xyplot(Weight ~ Height | equal.count(StartAge), age.model, type = c("p", "r"))
Количество ячеек по умолчанию для equal.count равно 6. Его можно легко изменить, чтобы изучить другие группы:
# Create four groups of equal counts to explore
xyplot(Weight ~ Height | equal.count(StartAge, 4), age.model, type = c("p", "r"))
Функция shingle() позволяет перекрывать ячейки, как показано здесь.
# Create three groups that overlapping with each other
bins <- cbind(lower = c(0,8,16), upper = c(13,18,24))
xyplot(Weight ~ Height | shingle(StartAge, bins), age.model, type = c("p", "r"))
пожалуйста, предоставьте часть
age.model, чтобы дать нам воспроизводимый пример, используя, например,dput