宁波做网站有哪些公司公司,微信公众号网站开发模板,广告设计公司宣传语,张家港企业网站设计目录 一、mvrnorm()函数使用介绍
例1#xff1a;生成服从多元正态分布的数据
例2:生成一组服从多元正态分布的观测 一、mvrnorm()函数使用介绍
获取来自给定均值向量和协方差阵的多元正态分布的数据。
MASS包中的mvrnorm()函数可以让这个问题变得很容易#xff0c;其调用…目录 一、mvrnorm()函数使用介绍
例1生成服从多元正态分布的数据
例2:生成一组服从多元正态分布的观测 一、mvrnorm()函数使用介绍
获取来自给定均值向量和协方差阵的多元正态分布的数据。
MASS包中的mvrnorm()函数可以让这个问题变得很容易其调用格式为 mvrnorm(n, mean, sigma) 其中n是你想要的样本大小mean是均值向量sigma是方差-协方差矩阵或相关矩阵。 可以想一下产生的随机数据一定是一个矩阵形式有n行即样本量为n产生了n个观测有p列因为是多元随机变量每一行是一个观测相当于人的体重、身高、臂长、体脂率、年龄等数据信息。
例1生成服从多元正态分布的数据
library(MASS)
options(digits3)
set.seed(1234) #设置随机种子让结果可以重现reproducible#指定均值向量和协方差矩阵
mean - c(230.7, 146.7, 3.6)
sigma - matrix(c(15360.8, 6721.2, -47.1, 6721.2, 4700.9, -16.5, -47.1, -16.5, 0.3), nrow3, ncol3)mydata - mvrnorm(500, mean, sigma)#为了方便将结果从矩阵转换为数据框
#并为变量指定了名称
mydata - as.data.frame(mydata)
names(mydata) - c(y,x1,x2)dim(mydata)
head(mydata, n10) 结果 例2:生成一组服从多元正态分布的观测
mvrnorm(1,rep(0,10),diag(rep(0.2^2,10)))
