3

さまざまな可能な線形モデルの予測力を自動的に評価するためのスクリプトを R で作成しようとしています。モデルの予測力を評価するために、パッケージ DAAG の関数 CVlm を使用する相互検証から得られる全体的な平均平方を品質指標として使用します。私の質問は、CVlm から得られた全体的な平均平方の値を自動化された方法で取得する方法です (CVlm のテキスト出力を観察する必要はありません)。

たとえば、http://maths-people.anu.edu.au/~johnm/r-book/3edn/scripts/reg1.Rの次のコード

houseprices.lm <- lm(sale.price ~ area, data=houseprices) 
CVlm(houseprices, houseprices.lm, plotit=TRUE) 

形式の出力があります

フォールド 1 テスト セット内の観察: ...

フォールド 2 テスト セット内の観察: ...

全体ミリ秒 2023

各実行の ms (2023) の値にアクセス/保存するにはどうすればよいですか?

4

1 に答える 1

2

CVlm結果を変数に格納し、ms属性にアクセスする必要があります。

houseprices.lm <- lm(sale.price ~ area, data=houseprices) 
cv <- CVlm(houseprices, houseprices.lm, plotit=TRUE) 
attr(cv, "ms")
# [1] 3934
于 2013-07-11T12:48:02.840 に答える