Count number of distinct values in a vector

2020-08-09 05:16发布

问题:

I have a vector of scalar values of which I'm trying to get: "How many different values there are".

For instance in group <- c(1,2,3,1,2,3,4,6) unique values are 1,2,3,4,6 so I want to get 5.

I came up with:

length(unique(group))

But I'm not sure it's the most efficient way to do it. Isn't there a better way to do this?

Note: My case is more complex than the example, consisting of around 1000 numbers with at most 25 different values.

回答1:

Here are a few ideas, all points towards your solution already being very fast. length(unique(x)) is what I would have used as well:

x <- sample.int(25, 1000, TRUE)

library(microbenchmark)
microbenchmark(length(unique(x)),
               nlevels(factor(x)),
               length(table(x)),
               sum(!duplicated(x)))
# Unit: microseconds
#                 expr     min       lq   median       uq      max neval
#    length(unique(x))  24.810  25.9005  27.1350  28.8605   48.854   100
#   nlevels(factor(x)) 367.646 371.6185 380.2025 411.8625 1347.343   100
#     length(table(x)) 505.035 511.3080 530.9490 575.0880 1685.454   100
#  sum(!duplicated(x))  24.030  25.7955  27.4275  30.0295   70.446   100


回答2:

I have used this function

length(unique(array))

and it works fine, and doesn't require external libraries.



回答3:

You can use rle from base package

  x<-c(1,2,3,1,2,3,4,6)
  length(rle(sort(x))$values)

rle produces two vectors (lengths and values ). The length of values vector gives you the number of unique values.



回答4:

uniqueN function from data.table is equivalent to length(unique(group)). It is also several times faster on larger datasets, but not so much on your example.

library(data.table)
library(microbenchmark)

xSmall <- sample.int(25, 1000, TRUE)
xBig <- sample.int(2500, 100000, TRUE)
microbenchmark(length(unique(xSmall)), uniqueN(xSmall), 
               length(unique(xBig)), uniqueN(xBig))

#Unit: microseconds
#                    expr      min        lq       mean    median        uq      max neval cld
#1 length(unique(xSmall))   17.742   24.1200   34.15156   29.3520   41.1435  104.789   100 a  
#2        uniqueN(xSmall)   12.359   16.1985   27.09922   19.5870   29.1455   97.103   100 a  
#3   length(unique(xBig)) 1611.127 1790.3065 2024.14570 1873.7450 2096.5360 3702.082   100 c
#4          uniqueN(xBig)  790.576  854.2180  941.90352  896.1205  974.6425 1714.020   100 b 


回答5:

If one wants to get number of unique elements in a matrix or data frame or list, the following code would do:

  if( typeof(Y)=="list"){ # Y is a list or data frame
    # data frame to matrix
    numUniqueElems <- length( na.exclude( unique(unlist(Y)) ) )
  } else if ( is.null(dim(Y)) ){ # Y is a vector
    numUniqueElems <- length( na.exclude( unique(Y) ) )
  } else { # length(dim(Y))==2, Yis a matrix
    numUniqueElems <- length( na.exclude( unique(c(Y)) ) )
  }


回答6:

We can use n_distinct from dplyr

dplyr::n_distinct(group)
#[1] 5