我正在使用纵向数据，其中结果是每分钟的步数。我的LMM尺寸看起来像
lme(step ~ predictor, random = ~1|person, data = df, na.action = na.omit, method = "ML")
我的项目的另一个组件是从仿真数据进行功率仿真，我对如何计算从mvrnorm函数计算“ Sigma”所需的步长数据的方差和协方差感到困惑。
我的sigma函数当前看起来像
matrix(c(varSteps,rep(covSteps,(nRep))),nRep,nRep,byrow=T)
因此，我正在尝试计算构建将在我的mvrnorm函数中使用的sigma矩阵所需的“ varSteps”和“ covSteps”。
如何使用R或SAS计算这两个？对你的帮助表示感谢。