r - 如何创建另一个向量的随机向量?
问题描述
我正在使用常量和向量(近似长度 = 100)执行计算,我需要为此模拟正态分布 N(使用 rnorm)。对于常数(K,标准偏差 = KU),我以标准方式使用 rnorm():
K <- 2
KU <- 0.2
set.seed(123)
KN <- rnorm(n = 3, mean = K, sd = KU)
什么提供长度为 3 (KN) 的向量:
[1] 1.887905 1.953965 2.311742
现在,我需要对向量(V,标准差 VU)做同样的事情。我的第一个猜测是使用:
V <- c(1, 2, 3)
VU <- 0.1 * V
set.seed(123)
VN <- rnorm(3, V, VU)
但只生成一个包含 3 个元素的向量,每个向量元素一个:
[1] 0.9439524 1.9539645 3.4676125
这实际上是向量的第一次模拟,但是我需要这个向量的3倍。一种解决方案是创建 9 个数字,但 VN 是 9 个元素的向量:
[1] 0.9439524 1.9539645 3.4676125 1.0070508 2.0258575 3.5145195 1.0460916 1.7469878 2.7939441
不是 3 个元素的 3 个向量。我想要的是VN =
[1] 0.9439524 1.0070508 1.0460916
[2] 1.9539645 2.0258575 1.7469878
[3] 3.4676125 3.5145195 2.7939441
因此,VN 是 3 个向量,我随后可以在其他计算中使用它们,例如 KN * VN。我找到的解决方案是:
set.seed(123)
VN <- as.data.frame(t(matrix(rnorm(3 * length(V), V, VU), nrow = length(V))))
但在我看来,这是一个相当麻烦的表达式(我需要在不同的地方重复多次,变量名很长)。基础 R 中是否有更简单的方法来生成随机向量?我想看到类似的东西:
VN <- rnorm.vector(3, V, VU)
解决方案
除了@akrun 的绝佳选择之外,您还可以使用比您的方法稍微简单的方法:
matrix(rnorm(n * length(V), V, VU), nrow = n, byrow = TRUE)
# [,1] [,2] [,3]
# [1,] 0.9439524 1.953965 3.467612
# [2,] 1.0070508 2.025858 3.514519
# [3,] 1.0460916 1.746988 2.793944
或者还有从多元正态分布中抽样的MASS
包:mvrnorm
library(MASS)
mvrnorm(n, VU, diag(VU))
# [,1] [,2] [,3]
# [1,] 0.6650715 0.37923044 0.05590089
# [2,] 0.2574341 0.24949882 0.97045721
# [3,] -0.5218990 -0.04857971 0.49707815
在哪里
diag(VU)
# [,1] [,2] [,3]
# [1,] 0.1 0.0 0.0
# [2,] 0.0 0.2 0.0
# [3,] 0.0 0.0 0.3
如果您希望方差-协方差矩阵不是对角线,则后一种选择是可行的方法。
推荐阅读
- reactjs - 将 NextJS 静态站点迁移到 squarespace/cms
- javascript - Alpine.js 无法将 x-data 绑定到外部 js 文件中的 functin
- mysql - 如何在选择查询中使用案例别名
- flutter - 导入的图片库出错
- python - Windows / Linux SQLAlchemy 连接器转换类型错误
- java - 从 JNI_OnLoad 返回的错误 JNI 版本,Android 11 崩溃
- python - 'pip install 和有什么区别
' 和 'python -m pip install '? - elasticsearch - 如何在 ElasticSearch 中的索引期间从现有字段创建新字段?
- javascript - 将数组中的先前文本移动到最新文本下方而不是被覆盖
- javascript - localStorage 检查器在看起来不错时不起作用?