R中二维正态分布的问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

R中一维正态分布的概率密度函数和累计函数的命令分别为dnorm和pnorm
请问，如果是二维正态分布，使用什么命令来计算概率密度函数和累计函数呢？
谢谢！


PS：随便问题一个小问题，二维正态概率密度函数中有一个rho，说是用来衡量两个随机变量相关程度的指标，那么在实践中可不可以等同为相关系数呢？

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：二维正态分布 正态分布 概率密度函数 pnorm dnorm 正太 正态分布

Package "fMultivar"
norm2d Bivariate Normal Distribution
  pnorm2d(x, y = x, rho = 0)
  dnorm2d(x, y = x, rho = 0)
  rnorm2d(n, rho = 0)

rho: the correlation parameter, a numeric value ranging between minus one and one,
      by default zero.

epoh 发表于 2012-12-23 21:06
Package "Multivar"
norm2d Bivariate Normal Distribution
  pnorm2d(x, y = x, rho = 0)

问题并没有被解决！尽管回答是非常好的，正是我所需要的结果。但是
我在CRAN没有发现这个软件包Multivar，那里有的下？能具体一些吗？我找不到这个软件包。


此外，我发现了Package ‘mvtnorm’，似乎可以做，但是使用比较讨厌，
dmvnorm(x, mean, sigma, log=FALSE)
sigma Covariance matrix, default is diag(ncol(x)).
我看不懂sigma的意思，是指协方差吗？如果已经知道两个随机变量的均值和方差，如何录入?
比如：
mean.values <- c(-20, 150)
sd.values     <- c(2.5, 10)
比如x是二维数值（既然是两维随机变量的话）
当然，如果能下到Multivar，这些疑问我将不管了。
谢谢！

是不是sigma=sd.values * diag(2)?
不惜赐教！

peijianshi 发表于 2012-12-24 16:02
问题并没有被解决！尽管回答是非常好的，正是我所需要的结果。但是
我在CRAN没有发现这个软件包Multivar ...

哈哈!第一个字母被我误删了
是package “fMultivar”

我感觉被您忽悠了。
下了根本没法用！需要接连下载若干个附带的包。下到fBasics后，出现如下错误信息：
> local({pkg <- select.list(sort(.packages(all.available = TRUE)),graphics=TRUE)
+ if(nchar(pkg)) library(pkg, character.only=TRUE)})
Error in loadNamespace(i, c(lib.loc, .libPaths())) :
  不存在叫‘stabledist’这个名字的程辑包
此外: 警告信息：
程辑包‘fBasics’是用R版本2.15.2 来建造的
错误: ‘fBasics’程辑包/名字空间载入失败，

peijianshi 发表于 2012-12-25 17:54
我感觉被您忽悠了。
下了根本没法用！需要接连下载若干个附带的包。下到fBasics后，出现如下错误信息：
> ...

想太多了吧,谁有那么多时间忽悠你.
library(fMultivar)
## Bivariate Normal Density:
x = (-40:40)/10
X = grid2d(x)
z = dnorm2d(X$x, X$y, rho = 0.5)
Z = list(x = x, y = x, z = matrix(z, ncol = length(x)))
persp(Z, theta = -40, phi = 30, col = "steelblue")

Package 'fMultivar' deals mainly with multivariate aspects of
economic and financial time series analysis.

epoh 发表于 2012-12-25 19:26
想太多了吧,谁有那么多时间忽悠你.
library(fMultivar)
## Bivariate Normal Density:

加载不了！

> library(fMultivar)
载入需要的程辑包：fBasics
Error in loadNamespace(i, c(lib.loc, .libPaths())) :
  不存在叫‘stabledist’这个名字的程辑包
此外: 警告信息：
1: 程辑包‘fMultivar’是用R版本2.15.2 来建造的
2: 程辑包‘fBasics’是用R版本2.15.2 来建造的
错误: 无法载入程辑包‘fBasics’

peijianshi 发表于 2012-12-26 13:20
加载不了！

> library(fMultivar)

我试了一下，可以安装那些包，没有任何问题。可能是你选的镜像有问题，或者是你的R不是最新版本的

利用epoh老师的语句，R2.15.2可以正常运行


> library(fMultivar)
载入需要的程辑包：sn
载入需要的程辑包：mnormt
Package 'sn', 0.4-17 (2011-07-13). Type 'help(SN)' for summary information
载入需要的程辑包：timeDate
载入需要的程辑包：timeSeries
载入需要的程辑包：fBasics
载入需要的程辑包：MASS

载入程辑包：‘fBasics’

The following object(s) are masked from ‘package:base’: