ご指摘を受け、コードのミスを修正しました。生成したYを使い忘れていた)ありがとうございます。変えて一個抜き交差検証法で評価した対数尤度が結構なめらかな形になったのでブレント法で最適なバンド幅を選んでみた。対象としたデータはこれ。変えて一個抜き交差検証法で評価した対数尤度のプロットはこちら。グラフとりあえず図とコードだけ貼ります。縦の棒はオッズ比の95%信頼区間です。メンタルヘルス精神疾患の軽症化が進んでいる。統合失調症の新規発症が激減し、入院を要するレベルの気分障害がほとんどいなくなり、病棟経営の支えである慢性の統合失調症患者は急速に高齢化が進んでいる。グラフアクセス解析の分野ではユーザーの離脱率や維持率を把握したいというニーズがあります。離脱するユーザーは単にサイトへの訪問を止めるだけで、わざわざ離脱を申告することはめったにありません。不完全データ追記:一応パッケージ化しました。生存関数のプロットは便利だけど密度関数のプロットと比べると分布の形状を把握しにくい。グラフ右打ち切りのあるデータを視覚化するのにはカプラン・マイヤープロットが便利です。層の数が多くなると判例と見比べるのが大変になる。そこである程度情報を落としてプロットしてもいいかなと思い、打ち切りデータの箱ひげ図を提案にします。不完全データ層分けが多くなると凡例と見比べるのが大変なのでなんとかしたいと思った。