如何在depmixS4中使用standerror()

Question

Summary(fittedmodel)给出了过渡矩阵的协方差系数。我也很好奇如何获得站立误差值。它看起来像standerror(fittedmodel)可以工作。但我对结果的理解有问题。例如，我的模型的summary()结果和调用standardError(fm2a)的结果，但是我不理解standardror()函数的结果。

代码:

mod1 <- depmix(list(rt~1,corr~1),data=speed,transition=~Pacc,nstates=2,
family=list(gaussian(),multinomial("identity")),ntimes=c(168,134,137))
set.seed(3)
fmod1 <- fit(mod1)
standardError(fmod1)

Answer 1

depmixS4::standardError函数返回一个数据帧，其行中参数的顺序与调用getpars(fmod1)时相同。在第一列(par)中，您将找到参数估计，在第二列constr中，您将找到反映估计是否在参数空间的边界上(bnd)，固定参数(fix)，或者参数是否包含在估计的标准误差中(inc)的代码。第三列(se)提供了所包含参数的估计标准误差。bnd或fix参数的标准误差无法估计，因此该值为NA。