如何在r中创建一个具有正态分布的矩阵？

Question

我想要一个像素网格，其中填充正态分布的值，如下图所示

有人知道如何在 R 中做到这一点吗？

Answer 1

不清楚你想做什么。您是在问如何从双变量正态分布中采样？或者如何创建热图？或者如何将文本叠加到热图上？不过好吧，我咬一口……

这是一个从头开始的示例，如何生成数据并生成热图，其中的数字表示每个 bin 框的计数。

# Create sample data
require(MASS);
set.seed(2017);
mu <- c(5, 5);
sigma <- diag(c(2, 2));
df <- as.data.frame(mvrnorm(10000, mu = mu, Sigma = sigma));
colnames(df) <- c("x1", "x2");

# Bin breaks
nbins <- 20;
breaks_x <- seq(floor(min(df$x1)), ceiling(max(df$x1)), length.out = nbins);
breaks_y <- seq(floor(min(df$x2)), ceiling(max(df$x2)), length.out = nbins);

# Frequency table
freq <-  as.data.frame(table(
    as.numeric(cut(df$x1, breaks = breaks_x)),
    as.numeric(cut(df$x2, breaks = breaks_y))));

# Plot matrix
mat <- diag(nbins) * 0;
mat[cbind(freq[, 1], freq[, 2])] <- freq[, 3];
image(breaks_x, breaks_y, mat, col= rainbow(10));
text(breaks_x[freq$Var1], breaks_y[freq$Var2], freq$Freq, cex = 0.8);

或者，如果您只对热图感兴趣，您可以使用 ggplot2 和 stat_bin2d。

require(ggplot2);
ggplot(data.frame(df), aes(x1, x2)) + stat_bin2d(bins = 20) + scale_fill_gradientn(colours = rainbow(32));

或gplots::hist2d。

更多选项请参见here...