对 lavaan 潜在增长曲线模型运行功效分析

Question

我正在尝试使用 semTools 进行功率分析使用 lavaan 估计的潜在增长曲线模型.见下文:

library(RCurl)
library(lavaan)
library(semTools)
x <- getURL("https://gist.githubusercontent.com/aronlindberg/dfa0115f1d80b84ebd48b3ed52f9c5ac/raw/3abf0f280a948d6273a61a75415796cc103f20e7/growth_data.csv")
growth_data <- read.csv(text = x)

model_regressions <- ' i =~ 1*t1 + 1*t2 + 1*t3 + 1*t4 + 1*t5 + 1*t6 + 1*t7 + 1*t8 + 1*t9 + 1*t10 + 1*t11 + 1*t12 + 1*t13+ 1*t14 + 1*t15 + 1*t16 + 1*t17 + 1*t18 + 1*t19 + 1*t20
s =~ 0*t1 + 1*t2 + 2*t3 + 3*t4 + 4*t5 + 5*t6 + 6*t7 + 7*t8 + 8*t9 + 9*t10 + 10*t11 + 11*t12 + 12*t13 + 13*t14 + 14*t15 + 15*t16 + 16*t17 + 17*t18 + 18*t19 + 19*t20

# fixing error-variances
t8 ~~ 0.01*t8
t17 ~~ 0.01*t17
t18 ~~ 0.01*t18
# regressions
s ~ h_index
i ~ h_index'

SSpower(powerModel = model_regressions, popModel = model_regressions, n = c(87, 125), fun = "growth")

然而，这似乎不起作用。我的总体问题是:如何使用 semTools 进行功效分析对于使用 lavaan 估计的潜在增长曲线模型?更具体地说，我应该用什么来指定 powerModel和 popModel ?

Answer 1

如何使用 semTools 运行功效分析对于使用 lavaan 估计的潜在增长曲线模型
SSpower来自 semTools应该管用。

更具体地说，我应该使用什么来指定 powerModel 和 popModel？

从语法的角度来说，你的 model_regressions对象似乎是有效的熔岩对象，您可以将其传递给 SSpower如 powerModel参数，用于描述要分析的模型，也可以作为 popModel参数来指定数据生成模型。尽管如此，您还需要指定 alpha - I 类错误率和 nparam - powerModel 中无效约束的数量.此外，除非您正在分析多组模型，并且只是尝试评估两个不同样本量的功效，否则您可以运行 SSpower - 命令分别使用 n=87 和 n=125。例如，当 n=87 时，代码如下所示:

  SSpower(powerModel = model_regressions, n = 87, nparam = 1, popModel = model_regressions,  fun = "growth", alpha=0.05)
    # [1] 0.05044634

注意:如果您尝试执行多组分析，我不确定此代码是否会以正确的方式工作。对于多组分析，我可能会尝试 manual 建议的另一个选项其中，而不是使用 popModel ，人们将“..指定人口模型中的所有非零参数，...通过提交人口模型隐含的人口协方差矩阵( Sigma)和可选的平均向量( mu)。”。您可以通过说明 mu 来做到这一点。和 sigma内 SSpower函数调用:

SSpower(...Sigma = , mu = , ...)