2013-12-18 45 views
-2

有人可以幫助我使用persp函數得到R中兩個變量正態分佈的3d圖嗎?我正在使用mvtnorm軟件包,它有點令人困惑...兩個變量正態分佈

任何幫助將不勝感激。

+0

這是非常模糊的。請指定您想要繪製的內容。 –

+0

使用mvtnorm的**樣本**兩個變量正態分佈的圖。 – AHHP

回答

2

基於您的評論,看來你想是這樣的:

show.2dnorm <- function(n, mean = rep(0, nrow(sigma)), sigma = diag(length(mean))) { 
    require(mvtnorm) 
    require(ggplot2) 
    norm2d <- as.data.frame(rmvnorm(n, mean, sigma)) 
    colnames(norm2d) <- c('x', 'y') 
    ggplot(norm2d, aes(x,y)) + geom_point() 
} 

# standard normal 
show.2dnorm(1e4, c(0, 0)) 
# 0.6 correlation 
show.2dnorm(1e4, sigma = matrix(c(1, 0.6, 0.6, 1), 2)) 

enter image description here enter image description here

+0

您能否幫我使用persp函數來創建3d圖? – AHHP

+0

這需要一點不同的方法。將這個問題保持原樣可能是一個好主意,而您可以通過詳細的要求打開另一個問題。 – tonytonov