[解決済み] 線形回帰からp値およびr二乗を取り出す
2022-02-07 01:01:46
質問
単純な線形回帰モデルから、p値(単一の説明変数の係数が0でないことの有意性)とR2乗の値を取り出すにはどうしたらよいでしょうか。例えば...
x = cumsum(c(0, runif(100, -1, +1)))
y = cumsum(c(0, runif(100, -1, +1)))
fit = lm(y ~ x)
summary(fit)
知っている
summary(fit)
ディスプレイ
のように、p値やR2乗の値を他の変数に貼り付けられるようにしたい。
どのように解決するのですか?
r二乗
: 要約オブジェクトから直接 r 二乗値を返すことができます。
summary(fit)$r.squared
. 参照
names(summary(fit))
は、直接抽出できるすべての項目のリストです。
モデルのp値。 回帰モデル全体のp値を取得したい場合。 このブログの記事 は、p値を返す関数を概説している。
lmp <- function (modelobject) {
if (class(modelobject) != "lm") stop("Not an object of class 'lm' ")
f <- summary(modelobject)$fstatistic
p <- pf(f[1],f[2],f[3],lower.tail=F)
attributes(p) <- NULL
return(p)
}
> lmp(fit)
[1] 1.622665e-05
予測変数が1つの単回帰の場合、モデルのp値と係数のp値は同じになります。
係数のp値。 もし、複数の予測変数がある場合、上記はモデルのp値を返し、係数のp値はそれを使って抽出することができます。
summary(fit)$coefficients[,4]
あるいは、係数の p 値を
anova(fit)
オブジェクトを、上記の summary オブジェクトと同様の方法で作成します。
関連
-
[解決済み】LMEモデルのレベル0、ブロック1でのバックソルブにおける特異性
-
[解決済み】 colMeans(x, na.rm = TRUE) のエラー : KNN分類では 'x' は数値でなければならない
-
[解決済み】ggplot boxplotでPosition-dodge警告?
-
[解決済み】match.fun(FUN)でのエラーについて)
-
[解決済み】RでAIC中に行数が変化するのはなぜですか?そうならないようにするにはどうしたらいいですか?
-
[解決済み】R: predict() の数値 'envir' 引数が長さ1でない。)
-
[解決済み] 因子を日付形式に変換するにはどうすればいいですか?
-
[解決済み] [Solved] read.csv warning 'EOF within quoted string' prevents complete reading of file.
-
[解決済み】ggplot2でのプロット:「Error: カテゴリ軸のY軸に "Discrete value supplied to continuous scale "と表示される。
-
[解決済み】Rで、Error: ggplot2 doesn't know how to handle of data of class numericに対処する。
最新
-
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 実装 サイバーパンク風ボタン
おすすめ
-
[解決済み】xtsオブジェクトでエラー: "antempt to set 'colnames' on the object with less than two dimension "を克服する方法
-
[解決済み】model.frame.defaultでのエラー:変数の長さが異なる
-
[解決済み】「'dimnames' [1]の長さが配列の範囲と等しくない」とはどういう意味ですか?
-
[解決済み】match.fun(FUN)でのエラーについて)
-
[解決済み】ロジスティック回帰 - eval(family$initialize) : y 値は 0 <= y <= 1 である必要があります。
-
[解決済み】R 置換する項目数が置換長の倍数でない/しかし結果は正しい
-
[解決済み】Rでデータフレームのサブセットをプロットする方法は?
-
[解決済み】Rで、Error: ggplot2 doesn't know how to handle of data of class numericに対処する。
-
[解決済み】各グループの上位3要素を計算する際に「数値型、論理型、複合型でのみ可能な演算」を回避する方法
-
[解決済み】起動時にパーをデフォルト値にリセットする方法