Wie erstellt man eine Matrix mit einer Normalverteilung in r?

Question

Ich möchte ein Pixelraster mit den Werten einer Normalverteilung in es wie unten

im Bild

Wer weiß, wie dies in R gefüllt zu tun haben?

Answer 1

Es ist nicht klar, was Sie tun möchten. Fragen Sie, wie Sie aus einer bivariaten Normalverteilung abtasten können? Oder wie man eine Heatmap erstellt? Oder wie man Text auf einer Heatmap überlagert? Aber ok, ich beiße ...

Hier ist ein Beispiel von Grund auf, wie man Daten erzeugt und eine Heatmap mit Zahlen erzeugt, die Zählimpulse pro Behälterkasten anzeigen.

# 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); 

Oder wenn Sie in einer Heatmap nur daran interessiert sind, können Sie ggplot2 mit stat_bin2d verwenden.

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

Oder gplots::hist2d.

Siehe here Weitere Optionen ...