r에 정규 분포가있는 행렬을 만드는 방법은 무엇입니까?

Question

나는 누군가가 R에서이 작업을 수행하는 방법을 알고 있나요

아래 그림처럼의 정규 분포의 값 가득 픽셀 격자를 갖고 싶어?

Answer 1

당신이하고 싶은 것이 명확하지 않습니다. 2 변수 정규 분포에서 표본 추출 방법을 묻는 중입니까? 아니면 히트 맵을 만드는 방법? 또는 히트 맵에 텍스트를 오버레이하는 방법은 무엇입니까? 하지만 괜찮아요. ...

데이터를 생성하고 빈 상자 당 개수를 나타내는 숫자가 포함 된 히트 맵을 생성하는 방법을 처음부터 보여주는 예입니다.

아니면 히트 맵에만 관심이 있다면, 당신은 stat_bin2d와 ggplot2을 사용할 수 있습니다

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

.

또는 gplots::hist2d

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

.

더 많은 옵션은 을 참조하십시오 ...