1

这个问题实际上来自我问过自己的这个问题(在序列中查找特定模式),但我相信这是一个单独的问题。

根据 Gilbert 的响应,我尝试从状态序列创建事件序列,但遇到了问题。

建议是使用

seqe=seqecreate(comp.seq,tevent="state")

然后使用

seqefsub(seqe,strsubseq="(a)-(d)")

但是当我尝试使用 seqecreate() 时,出现以下错误:

Error in `seqelength<-`(`*tmp*`, value = c(64, 64, 64, 64, 61, 62, 61,  : (...) 
s and len should be of the same size.

如果我尝试使用以下方法将其转换为事件序列,也会发生同样的情况:

seqe=seqecreate(comp.seq,tevent="transition")

尝试使用行的子集并确定哪些行导致了问题,我发现有问题的行都处于恒定状态,这意味着它们实际上没有转换,它们永久保持在同一个状态(例如 AAAAAA)。

所以我的问题是:

  1. 我可以设置任何标志或任何东西来进行转换吗?
  2. 如果不是,我如何删除这些行,因为它们具有不同的长度和缺失值。例如,我可能有如下序列:

    缺失-缺失-AAAA AA-缺失-缺失-缺失-缺失-缺失

非常感谢提前!

提供我的数据样本:

comp.seq <- seqdef(comp,NULL,states=comp.scodes,labels=comp.labels, alphabet=comp.alphabet,missing="Z") comp.seq[1:7,] 1 *-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-B-B-B-B-B-B-B-B-D-D-D-D-D-A-A-A-A-A-A-A-A-A 2 *-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-C-C-C-C-C-C-C-C-C-C-C-C-C-C-*-B-B-B-B-B-B-B-B-B-B-B-B-B-A-A-A-A-A-A 3 *-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-A-A-A-C-C-A-A-A-A-A-A-A-D-D-A-A-A-A-A-A-A-A 4 *-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-B-B-B-B-B-B-A-A-A-A-A-A-A-A-A-A-A-A-A-A-A-A-A-A-A-A 5 *-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-C-C-D-D-D-D-D-D-D-D-D-D-A-A-A-A-A 6 *-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-B-B-B-B-B-B-B-B-B-B-B-B-B-D-D-D-D-D-D-D-D-A-A-A-A 7 *-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-A-A-A-A-A-A-A-A-A-A-A-A

第 7 行是一个有问题的行。如果我尝试使用seqecreate(comp.seq[1:6,])它就可以了

4

3 回答 3

3

当缺少状态并且序列长度不同时,会发生错误。一种解决方法是right="NA"seqdef调用中设置。

这是一个最小的例子:

x1 <- "*-*-A-B"
x2 <- "*-A-A"
dat.str <- data.frame(string=rbind(x1,x2))
dat <- seqdecomp(dat.str, sep="-", miss="*")

## creating state sequence object with and without right="NA"
dat.seq.NA <- seqdef(dat, right="NA")
dat.seq.void <- seqdef(dat)

## next command works without error
dat.eseq <- seqecreate(dat.seq.NA, tevent="state")

## while this one produces the error
dat.eseq <- seqecreate(dat.seq.void, tevent="state")
于 2015-01-24T07:16:28.693 回答
1

因此,在我使用的序列中,我为缺失值设置了一个代码,即函数missing="Z"中的选项seqdef()

missing我设法通过不设置选项并创建“虚拟”状态 Z来使其工作,我将其添加到字母表和标签“Z-missing”中。我还设置了选项left="Z"right="Z".

不过对我来说仍然像一个错误。

于 2015-01-23T19:38:31.020 回答
0

我不确定这是否是答案,但在这个 Cran 页面http://cran.r-project.org/web/packages/TraMineR/NEWS上,他们谈论了 TraMineR 的开发版本 1.9.8 他们参考一个错误:

错误修复: - seqformat():从 STS 转换为 TSE 时,如果 tevent 矩阵有空字符串(即“”),则会引发错误。现在,这被视为没有事件。

情况并非完全如此,因为序列本身不是空的,但它可能是同一个问题。我相信他们在内部使用 seqformat() 函数,因此该错误可能与我有某种关联。

我将下载开发版本并将其发布给她。

于 2015-01-23T18:00:52.727 回答