ロジスティック回帰をやるときに完全分離状態っていうのに陥るときがあって、例えば以下のようなデータセット Y <- c(1,1,1,1,0,0,0,0) X <- c(8,7,3,1,0,-2,-1,-9) dat <- data.frame(Y,X) dat > dat Y X 1 1 8 2 1 7 3 1 3 4 1 1 5 0 0 6 0 -2 7 0 -1 8 0 …

修正：2020/12/01 標準化残差の式が間違っていた(いやまあそういう定義もあるんだけど、R内での定義とは異なっていたことに気づかんかった)ので修正しました。申し訳ないです。 外れ値を含むデータ 以下の仮想データを見てみましょう set.seed(1) X <- c( rn…