私の問題に対する十分な答えが見つかりませんでした。おそらく誰かがここで助けてくれますか? (私はRの初心者です)
シーケンス分析を行います。状態空間は n = 10、時間空間は t = 168 (月) です。8 つのクラスターを持つクラスター ソリューションのシーケンス頻度プロットを描画しました。ただし、単一のプロットがあまりにもクリンチされているか小さすぎるため、プロットは解釈に対して実際には開かれていません。（下のグラフを参照）
これまでに次の手順を実行しました (TraMineR ヘルプ ドキュメントの手順に非常に近い)。
dist.om1 <- seqdist(neu.seq, method = "OM", indel = 1, sm = submat)
clusterward1 <- agnes(dist.om1, diss = TRUE, method = "ward")
cluster8 <- cutree(clusterward1, k = 8)
cluster8 <- factor(cluster8, labels = c("Typ 1", "Typ 2", "Typ 3", "Typ 4", "Typ 5", "Typ 6", "Typ 7", "Typ 8"))
seqfplot(neu.seq, group = cluster8, pbarw = T, withlegend = T)
マージンを再構成しようとしましたが、結果は常に同じプロットでした (添付のプロットはデフォルト設定で行われました)。そこで、代わりに、8 クラスター ソリューション内の単一クラスターのシーケンス頻度プロットを描画できるのではないかと考えました。（Stataコードでは、単一のシーケンスインデックスプロットのようなものを書きますsqindexplot if cluster8 == 4）しかし、Rでこれがどのように行われるかはわかりません。誰かがよりきれいなシーケンス頻度プロットを取得する方法を知っているなら、私はとてもありがたい！ありがとうございました！オリバー