多維直方圖，其中R

Question

假設我們有一個這樣的數據幀：

dat <- data.frame(
    a = rnorm(1000), 
    b = 1/(rnorm(1000))^2, 
    c = 1/rnorm(1000), 
    d = as.factor(sample(c(0, 1, 2), 1000, replace=TRUE)), 
    e = as.factor(sample(c('X', 'Y'), 1000, replace=TRUE)) 
) 

我們想計算此數據的直方圖在所有維度上（即A，B，C，d，e）在指定的每個維度都會打破。顯然，因素維度意味着他們已經休息了。最終的數據應該像data.frame一樣，其中每行是所有維度（中斷的組合）以及該組合的數據發生次數的向量。 Python numpy有直方圖：Multidimension histogram in python。 R中有類似的東西嗎？在R中做這件事的最好方法是什麼？謝謝。

我結束了使用以下，其中庫計數傳遞給函數的最後一行：

dat <- data.frame(
    a = rnorm(1000), 
    b = 1/(rnorm(1000))^2, 
    c = 1/rnorm(1000), 
    d = as.factor(sample(c(0, 1, 2), 1000, replace=TRUE)), 
    e = as.factor(sample(c('X', 'Y'), 1000, replace=TRUE)) 
) 

dat[nrow(dat)+1,] <- c(10,10,10,NaN,NaN) 

histnd <- function(df) { 
    res <- lapply(df, function(x) { 
    bin_idx <- length(x) 
    if (is.factor(x) || is.character(x)) { 
     return(x[-bin_idx]) 
    } 
    # 
    x_min <- min(x[-bin_idx]) 
    x_max <- max(x[-bin_idx]) 
    breaks <- seq(x_min, x_max, (x_max - x_min)/x[bin_idx]) 
    cut(x[-bin_idx], breaks) 
    }) 
    res <- do.call(data.frame, res) 
    res$FR <- as.numeric(0) 
    res <- aggregate(FR ~ ., res, length) 
} 

h <- histnd(dat) 

Answer 1

我不知道預期的結果是什麼，但是這應該提供一個起點：

histnd <- function(DF) { 
    res <- lapply(DF, function(x) { 
    if (is.factor(x) || is.character(x)) return(x) 
    breaks <- pretty(range(x), n = nclass.Sturges(x), min.n = 1) 
    cut(x, breaks) 
    }) 
    res <- do.call(data.frame, res) 
    as.data.frame(table(res)) 
} 

h <- histnd(dat)