全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
24054 14
2012-12-23
R中一维正态分布的概率密度函数和累计函数的命令分别为dnorm和pnorm
请问,如果是二维正态分布,使用什么命令来计算概率密度函数和累计函数呢?
谢谢!


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

扫码加我 拉你入群

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

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

全部回复
2012-12-23 21:06:15
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.
二维码

扫码加我 拉你入群

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

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

2012-12-24 16:02:43
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)?
不惜赐教!


二维码

扫码加我 拉你入群

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

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

2012-12-24 19:05:45
peijianshi 发表于 2012-12-24 16:02
问题并没有被解决!尽管回答是非常好的,正是我所需要的结果。但是
我在CRAN没有发现这个软件包Multivar ...
哈哈!第一个字母被我误删了
是package “fMultivar”
二维码

扫码加我 拉你入群

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

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

2012-12-25 17:54:46
我感觉被您忽悠了。
下了根本没法用!需要接连下载若干个附带的包。下到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’程辑包/名字空间载入失败,

二维码

扫码加我 拉你入群

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

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

2012-12-25 19:26:05
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")
Bivariate Normal Density.jpeg
二维码

扫码加我 拉你入群

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

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

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群