2つの単純なループの出力の違いが混乱する

Question

なぜ2つのコードが異なる結果を返すのか混乱しています。 ループの唯一の違いは、Wage $ age.cut1対age.cut1の使用です。違いの意義は何ですか？

DATA：ISLRパッケージ、賃金データ

cv.err <- rep(NA, 10) 
for (i in 2:10){ 
    Wage$age.cut1 = cut(Wage$age, i) 
    fit = glm(wage~age.cut1, data = Wage) 
    cv.err[i] = cv.glm(Wage, fit, K = 10)$delta[2] 
} 

> cv.err 
[1]  NA 1733.646 1681.587 1636.521 1632.931 1623.112 1611.965 1600.903 1609.973 
[10] 1607.234 # these are the expected results 




cv.err <- rep(NA, 10) 
for (i in 2:10){ 
    age.cut1 = cut(Wage$age, i) 
    fit = glm(wage~age.cut1, data = Wage) 
    cv.err[i] = cv.glm(Wage, fit, K = 10)$delta[2] 
} 
> cv.err 
[1]  NA 1603.255 1608.617 1602.296 1606.265 1606.139 1602.448 1606.063 1605.100 
[10] 1606.986 

Answer 1

はい、これら二つの違いは、あなたのループロジックに大きな違いを生みます。最初のループでは、 age.cut1 は$修飾子によって証明賃金データフレームの列であり、はglm式使用されます。 2番目のループでは、age.cut1はスタンドアロンの別名のベクターであり、の式での未使用のです。数式が使用されるときはいつでも、列はデータ引数で参照されるオブジェクトから派生します。記載されているパッケージとデータ構造に不慣れ

は、最も可能性の高いage.cut1列は前に（エラーが第二のループのglmコールでの参照で発生していないので）ループにWageデータフレーム内に存在しません。ただし、最初のループで更新されるのはcut(Wage$age, i)のみです。同じ名前のオブジェクトが2番目のループに割り当てられていますが、元の列データはそのままglmに残ります。