r软件中如何设计模拟实验运行次数

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

n =100 #产生100个观测值

p =25 #产生25个自变量
beta =rep.int(c(1, 0), c(5, p-5)) # 定义自变量系数beta可以改变
sigma =0.5 # 定义残差系数（可以改变）
Sigma =0.5^t(sapply(1:p, function(i, j) abs(i-j), 1:p))#定义产生x矩阵时x之间的协方差阵系数0.5
x=rmvnorm(n, sigma=Sigma) # 生成x矩阵

e=rnorm(n) # 生成残差矩阵
y=c(x %*% beta + sigma * e) # 生成y矩阵

各位大神，哪位知道有什么软件包哪个命令可以将这个程序运行100次啊，这是一次的程序

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：r软件 function Sigma apply beta 如何 软件

用list 一个循环就可以了
s<-list()
for(i in 1:100)
{
n =100 #产生100个观测值
p =25 #产生25个自变量
beta =rep.int(c(1,0),c(5,p-5)) # 定义自变量系数beta可以改变
sigma =0.5 # 定义残差系数（可以改变）
Sigma =0.5^t(sapply(1:p, function(i, j) abs(i-j), 1:p))#定义产生x矩阵时x之间的协方差阵系数0.5
x=rmvnorm(n, sigma=Sigma) # 生成x矩阵
e=rnorm(n) # 生成残差矩阵
y=c(x %*% beta + sigma * e) # 生成y矩阵
s[[i]]<-y
}