r - seqistatd 命令是否计算每个状态所花费时间的标准误差？

Question

我正在使用 R 包 TraMineR 来计算和分析事件状态序列。我的字母表由 7 个状态组成。我使用该seqistatd命令计算了我感兴趣的一些亚群（例如具有不同教育水平的女性）在每个州花费的平均时间。但是，我想添加有关标准误差或在每个状态中花费的平均时间的标准偏差的信息。是否可以seqistatd在行命令中做到这一点，还是我必须以某种方式自己计算估计？谢谢你。

Answer 1

不，该函数seqistatd确实为每个单独的序列返回每个状态的时间。然而，使用该函数很容易按组计算平均时间和方差apply。我在下面用一个使用biofam提供的数据的例子来说明TraMineR。

library(TraMineR)
data(biofam)
bf.seq <- seqdef(biofam[,10:25])

bf.ldist <- seqistatd(bf.seq)
sel <- (biofam$sex == "man")

## Mean and variance of time in each state for men
(meant <- apply(bf.ldist[sel,], 2, mean, na.rm=TRUE))
(vart <- apply(bf.ldist[sel,], 2, var, na.rm=TRUE))
## and for women
(meant <- apply(bf.ldist[!sel,], 2, mean, na.rm=TRUE))
(vart <- apply(bf.ldist[!sel,], 2, var, na.rm=TRUE))