##生成模拟数据集##生成回归观察y <-mass :: mvrnorm(n = 1,mu = ant + nat,sigma = cov)##生成强迫响应mruns <-c(1,1,1,1,1) mass :: mvrnorm(n = 1,mu = nat,sigma = cov / mruns [2]))##控制运行ctlruns <-mass <-mass :: mvrnorm(100,mu = rep = rep(0,nrow(cov)),sigma = cov cov ctlruns.sigma for Point和ctlar estions.s.s.s.s. ctlruns.sigma <-ctlruns.bhvar <-ctlruns ##位置数s <-25 ##年度步骤t <-10