2

biofam データセットの使用

library(TraMineR)
data(biofam)
lab <- c("P","L","M","LM","C","LC","LMC","D")
biofam.seq <- seqdef(biofam[,10:25], states=lab)
head(biofam.seq)

 Sequence                                    
1167 P-P-P-P-P-P-P-P-P-LM-LMC-LMC-LMC-LMC-LMC-LMC
514  P-L-L-L-L-L-L-L-L-L-L-LM-LMC-LMC-LMC-LMC    
1013 P-P-P-P-P-P-P-L-L-L-L-L-LM-LMC-LMC-LMC      
275  P-P-P-P-P-L-L-L-L-L-L-L-L-L-L-L             
2580 P-P-P-P-P-L-L-L-L-L-L-L-L-LMC-LMC-LMC       
773  P-P-P-P-P-P-P-P-P-P-P-P-P-P-P-P

回帰木を適合させて表示できます。

seqt <- seqtree(biofam.seq~sex + birthyr, data=biofam)

seqtreedisplay(seqt, type="I", border=NA, withlegend= TRUE, legend.fontsize=2, legendtext = "Biofam Regression Tree")

次に、リーフ メンバーシップを識別できます。

seqt$fitted[,1]

しかし、これは私が混乱するところです。プロット内のどの葉に対応する葉の番号を知るにはどうすればよいですか? グラフには表示されないようで、実行print(seqt)してもリーフ番号が表示されないようです。

私が達成したいのは、各リーフのシーケンスを分離して、各リーフで個別に記述を実行できるようにすることです。どうすればこれを達成できますか?

4

2 に答える 2

2

実際には、ツリーによって定義されたルールを探しています。木を見るとそれらがわかります。

たとえば、例の左端のブランチはseqtルールを定義します。

birthyr <= 1940 & birthyr <= 1928

そして一番左下の葉はによって定義されます

birthyr <= 1940 & birthyr > 1928 & sex == "man"

しかし、私は恐れています。あなたが正しいこと。(your ) によって返されるdisstreeオブジェクトには、現在、その情報が明示的に含まれていません。おそらく、さらなるバージョンで。TraMineRseqt

于 2014-10-23T13:03:30.763 に答える