この配列は、左から右へ読む必要があります。また、プロットツリーが下から上に読まれるべきであることを示している
library(PST)
data.seq <- seqdef("A-B-C-D-E-F")
S1.test <- pstree(data.seq, ymin = 0.001, lik = FALSE, with.missing = FALSE)
print(S1.test)
--(e)-[ p=(0.2,0.2,0.2,0.2,0.2,0.2) - n=6 ]
`--(A)-[ p=(0.001,0.995,0.001,0.001,0.001,0.001) - n=1 ]--|
`--(B)-[ p=(0.001,0.001,0.995,0.001,0.001,0.001) - n=1 ]
`--(A-B)-[ p=(0.001,0.001,0.995,0.001,0.001,0.001) - n=1 ]--|
`--(C)-[ p=(0.001,0.001,0.001,0.995,0.001,0.001) - n=1 ]
`--(B-C)-[ p=(0.001,0.001,0.001,0.995,0.001,0.001) - n=1 ]
`--(A-B-C)-[ p=(0.001,0.001,0.001,0.995,0.001,0.001) - n=1 ]--|
`--(D)-[ p=(0.001,0.001,0.001,0.001,0.995,0.001) - n=1 ]
`--(C-D)-[ p=(0.001,0.001,0.001,0.001,0.995,0.001) - n=1 ]
`--(B-C-D)-[ p=(0.001,0.001,0.001,0.001,0.995,0.001) - n=1 ]
`--(A-B-C-D)-[ p=(0.001,0.001,0.001,0.001,0.995,0.001) - n=1 ]--|
`--(E)-[ p=(0.001,0.001,0.001,0.001,0.001,0.995) - n=1 ]
`--(D-E)-[ p=(0.001,0.001,0.001,0.001,0.001,0.995) - n=1 ]
`--(C-D-E)-[ p=(0.001,0.001,0.001,0.001,0.001,0.995) - n=1 ]
`--(B-C-D-E)-[ p=(0.001,0.001,0.001,0.001,0.001,0.995) - n=1 ]
`--(A-B-C-D-E)-[ p=(0.001,0.001,0.001,0.001,0.001,0.995) - n=1 ]--|
plot(S1.test)
:次のコードは、この検証を提供します。