Подсчет уникальных значений
допустим, у меня есть:
v = rep(c(1,2, 2, 2), 25)
теперь, я хочу, чтобы подсчитать, сколько раз каждое уникальное значение. unique(v) возвращает уникальные значения, а не сколько их.
> unique(v)
[1] 1 2
Я хочу что-то, что дает мне
length(v[v==1])
[1] 25
length(v[v==2])
[1] 75
но как более общий однострочный:) что-то близкое (но не совсем) вроде этого:
#<doesn't work right> length(v[v==unique(v)])
11 ответов:
возможно, таблица-это то, что вам нужно?
dummyData = rep(c(1,2, 2, 2), 25) table(dummyData) # dummyData # 1 2 # 25 75 ## or another presentation of the same data as.data.frame(table(dummyData)) # dummyData Freq # 1 1 25 # 2 2 75
это однострочный подход с помощью
aggregate.> aggregate(data.frame(count = v), list(value = v), length) value count 1 1 25 2 2 75
если у вас есть несколько факторов (= многомерного фрейма данных), вы можете использовать
dplyrпакета для подсчета уникальных значений в каждой комбинации факторов:library("dplyr") data %>% group_by(factor1, factor2) %>% summarize(count=n())Он использует оператор
%>%для цепных вызовов метода на фрейме данныхdata.
функция table () - это хороший способ пойти, как Чейз предложил. Если вы анализируете большой набор данных, альтернативным способом является использование .N функция в пакете datatable.
убедитесь, что вы установили пакет данных таблицы с помощью
install.packages("data.table")код:
# Import the data.table package library(data.table) # Generate a data table object, which draws a number 10^7 times # from 1 to 10 with replacement DT<-data.table(x=sample(1:10,1E7,TRUE)) # Count Frequency of each factor level DT[,.N,by=x]
чтобы получить безразмерный целочисленный вектор, содержащий количество уникальных значений, используйте
c().dummyData = rep(c(1, 2, 2, 2), 25) # Chase's reproducible data c(table(dummyData)) # get un-dimensioned integer vector 1 2 25 75 str(c(table(dummyData)) ) # confirm structure Named int [1:2] 25 75 - attr(*, "names")= chr [1:2] "1" "2"Это может быть полезно, если вам нужно передать количество уникальных значений в другую функцию, и это короче и более идиоматично, чем
t(as.data.frame(table(dummyData))[,2]опубликовано в комментарии к ответу Чейза. Спасибо Рикардо сапорта, который указал мне на это здесь.
Если вам нужно иметь количество уникальных значений в качестве дополнительного столбца во фрейме данных, содержащем ваши значения (столбец, который может представлять размер выборки, например), plyr предоставляет аккуратный способ:
data_frame <- data.frame(v = rep(c(1,2, 2, 2), 25)) library("plyr") data_frame <- ddply(data_frame, .(v), transform, n = length(v))
это работает для меня. Возьмите свой вектор
v
length(summary(as.factor(v),maxsum=50000))комментарий: установите maxsum, чтобы быть достаточно большим, чтобы захватить количество уникальных значений
или
magrittrпакета
v %>% as.factor %>% summary(maxsum=50000) %>% length
также делая значения категориальными и вызывая
summary()будет работать.> v = rep(as.factor(c(1,2, 2, 2)), 25) > summary(v) 1 2 25 75
вы можете попробовать также a
tidyverselibrary(tidyverse) dummyData %>% as.tibble() %>% count(value) # A tibble: 2 x 2 value n <dbl> <int> 1 1 25 2 2 75
Если вы хотите запустить unique для данных.рама (например, поезд.данные), а также получить подсчеты (которые могут быть использованы в качестве веса в классификаторах), вы можете сделать следующее:
unique.count = function(train.data, all.numeric=FALSE) { # first convert each row in the data.frame to a string train.data.str = apply(train.data, 1, function(x) paste(x, collapse=',')) # use table to index and count the strings train.data.str.t = table(train.data.str) # get the unique data string from the row.names train.data.str.uniq = row.names(train.data.str.t) weight = as.numeric(train.data.str.t) # convert the unique data string to data.frame if (all.numeric) { train.data.uniq = as.data.frame(t(apply(cbind(train.data.str.uniq), 1, function(x) as.numeric(unlist(strsplit(x, split=",")))))) } else { train.data.uniq = as.data.frame(t(apply(cbind(train.data.str.uniq), 1, function(x) unlist(strsplit(x, split=","))))) } names(train.data.uniq) = names(train.data) list(data=train.data.uniq, weight=weight) }
count_unique_words <-function(wlist) { ucountlist = list() unamelist = c() for (i in wlist) { if (is.element(i, unamelist)) ucountlist[[i]] <- ucountlist[[i]] +1 else { listlen <- length(ucountlist) ucountlist[[i]] <- 1 unamelist <- c(unamelist, i) } } ucountlist } expt_counts <- count_unique_words(population) for(i in names(expt_counts)) cat(i, expt_counts[[i]], "\n")
Comments