生成多变量对数正态随机数

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

如何用R 生成对数正态分布随机数，且期望为单位向量？
本人写了一下，感觉不正确
谢谢了

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：多变量 随机数 对数正态分布 正态分布 如何用 正态分布 如何

http://artax.karlin.mff.cuni.cz/ ... ns/html/rlnorm.html

hugebear 发表于 2015-5-6 01:53
http://artax.karlin.mff.cuni.cz/r-help/library/compositions/html/rlnorm.html

谢谢楼上的回答！
其实我曾尝试网页上的这个命令
发现不是自己想要的结果
你看看下面的几行
omega=matrix(c(0.5,0.3),2)
A=matrix(c(0.4,0.15,0.1,0.3),2)            
B=matrix(c(0.3,0.1,0.2,0.4),2)
icovm<- matrix(c(
0.2,0.1,
0.1,0.2
),byrow=TRUE,nrow=2)

nTa=100
  nobs=200
k=2
  xx_mean=solve(diag(k)-A-B)%*%omega
  m=matrix(0,k,nobs+nTa)
  x=matrix(0,k,nobs+nTa)
  m[,1]=xx_mean
  x[,1]=xx_mean
  
err1 <- compositions::rnorm.aplus(nobs+nTa,rep(1,k),icovm)   
err=t(err1)
   

   for( i in 2:(nobs+nTa))
      {
       m[,i]=omega+A%*%x[,i-1] +B%*%m[,i-1]
       x[,i]=m[,i]*err[,i]
      }
   x=x[,(nTa+1):ncol(x)]

打算生成vector multiplicative error model的随机数
发现跟预先设想的差别太大了
上面生成的数据太大