[解決済み] lm()関数でleave-one-outクロスバリデーションを行う簡単なコマンドはありますか?
2022-02-03 11:52:52
質問
を使ったleave-one-outクロスバリデーションを行う簡単なコマンドはありますか?
lm()
という関数をRで使うか?
具体的には、以下のような簡単なコマンドがあるのでしょうか?
x <- rnorm(1000,3,2)
y <- 2*x + rnorm(1000)
pred_error_sq <- c(0)
for(i in 1:1000) {
x_i <- x[-i]
y_i <- y[-i]
mdl <- lm(y_i ~ x_i) # leave i'th observation out
y_pred <- predict(mdl, data.frame(x_i = x[i])) # predict i'th observation
pred_error_sq <- pred_error_sq + (y[i] - y_pred)^2 # cumulate squared prediction errors
}
y_squared <- sum((y-mean(y))^2)/100 # Variation of the data
R_squared <- 1 - (pred_error_sq/y_squared) # Measure for goodness of fit
どのように解決するのですか?
もう一つの解決策は
caret
library(caret)
data <- data.frame(x = rnorm(1000, 3, 2), y = 2*x + rnorm(1000))
train(y ~ x, method = "lm", data = data, trControl = trainControl(method = "LOOCV"))
線形回帰
1000 サンプル 1 予測器
前処理なし リサンプリング:Leave-One-Out Cross-Validation まとめ サンプルサイズ:999、999、999、999、999、・・・・。 リサンプリング の結果です。
RMSE Rsquared MAE
1.050268 0.940619 0.836808チューニングパラメーター「intercept」はTRUEで一定とした
関連
-
[解決済み】library(ggplot2)でエラー:'ggplot2'というパッケージは存在しません。
-
[解決済み】データのマージ - fix.by(by.x, x)のエラー)
-
[解決済み】match.fun(FUN)でのエラーについて)
-
[解決済み】 lm.fit(x, y, offset = offset, singular.ok = singular.ok, ...) のエラー : 'y' の NA/NaN/Inf, あらゆる方法を試したが解決しなかった。
-
[解決済み】RでAIC中に行数が変化するのはなぜですか?そうならないようにするにはどうしたらいいですか?
-
[解決済み】Rでデータフレームのサブセットをプロットする方法は?
-
[解決済み】Rはプロットするが、アブラインを描画しない
-
[解決済み】randomForestの実行予測で「NA/NaN/Inf in foreign function call (arg 7)」をなくすには?
-
[解決済み】.External.graphics Rでエラーが発生しました。
-
[解決済み】dplyr: "Error in n(): 関数は直接呼ばれるべきではありません"
最新
-
nginxです。[emerg] 0.0.0.0:80 への bind() に失敗しました (98: アドレスは既に使用中です)
-
htmlページでギリシャ文字を使うには
-
ピュアhtml+cssでの要素読み込み効果
-
純粋なhtml + cssで五輪を実現するサンプルコード
-
ナビゲーションバー・ドロップダウンメニューのHTML+CSSサンプルコード
-
タイピング効果を実現するピュアhtml+css
-
htmlの選択ボックスのプレースホルダー作成に関する質問
-
html css3 伸縮しない 画像表示効果
-
トップナビゲーションバーメニュー作成用HTML+CSS
-
html+css 実装 サイバーパンク風ボタン
おすすめ
-
[解決済み】model.frame.defaultでのエラー:変数の長さが異なる
-
[解決済み】エラー。Rの次元数が正しくない
-
[解決済み] 変数の型(リスト)が無効です
-
[解決済み】添え字付き代入でNAを使用することはできません
-
[解決済み] テスト
-
[解決済み] write.tableしようとすると、未実装の型リストが表示される。
-
[解決済み】reshape2 meltの警告メッセージ
-
[解決済み】長いオブジェクトの長さは、短いオブジェクトの長さの倍数ではない?[重複]。
-
[解決済み】R: predict() の数値 'envir' 引数が長さ1でない。)
-
[解決済み】.External.graphics Rでエラーが発生しました。