パッケージを使用してrpart
、分類/回帰ツリーをモデル化します。変数 x、y、s があります。ここで、x は {-1,1} で、y は [0,1] で連続しており、s$ は 3 レベルの因子です。
私が使う
fit <- rpart(x~y+s, data=data, method="class")
最終的なモデルは完全に理にかなっています. を使用してプロットできますfancyRpartPlot(fit)
. ツリーを剪定することもできます。決定事項は 1 つだけです。
fit = prune(fit,cp=0.1)
$y$ は何らかのしきい値よりも小さくする必要があると書かれています。このしきい値を抽出するにはどうすればよいですか? データセットでブートストラップを実行し、しきい値を配列に保持したいと思います。