ブートパッケージで加重回帰（lm）を使用しているときの問題R

Question

ブート機能を使用して統計情報に引数を渡す際に問題が発生しています。ここでは一例です：evalの中

rm(list=ls()) 
library(faraway) 
library(boot) 
data(strongx) 
g3 <- lm(crossx ~ energy, strongx, weights=rep(1, nrow(strongx))) 

bootReg <- function(formulatostat, data, ind, weitostat) { 
    weig <- weitostat[ind] 
    fit <- lm(formula=formulatostat, data = data[ind, ], weights = weig) 
    return(coef(fit)) 
} 
wei <- rep(1, nrow(strongx)) 

bootResults <- 
    boot(
    data = strongx, 
    statistic = bootReg, 
    R = 2000, 
    formulatostat = crossx ~ energy, 
    weitostat = wei) 

エラー（exprの、ENVIR、enclos）：オブジェクトのweigは、「もちろん

を見つけていない、私の重みは1であり、その加重部分が不要です。しかし、実際の重量を使用している間も問題は残ります。これは私が別のデータを持っている問題の実証です。ブート機能は次のように述べています。 "それ以上の引数は、引数を介して統計に渡すことができます。" FormRatostatはbootRegに渡されていますが、重みビットは渡されていません。ウェイトなしで比較する：

bootReg <- function(formulatostat, data, ind) { 
    fit <- lm(formula=formulatostat, data = data[ind, ]) 
    return(coef(fit)) 
} 
wei <- rep(1, nrow(strongx)) 

bootResults <- 
    boot(
    data = strongx, 
    statistic = bootReg, 
    R = 2000, 
    formulatostat = crossx ~ energy) 

これは機能します。体重をどのように組み込むべきかアドバイスをいただければ幸いです。

Answer 1

はdata引数に指定されたオブジェクトにweigという名前のオブジェクトを探し、そして式（crossx ~ energy）が定義されている親環境にある（lm関数内で呼び出される）stats::model.frame関数のように思えます。

dataに与えられたdata.frame内にweigという名前の列がありません。また、式crossx ~ energyをbootReg関数の外に定義したため、lm関数はweigオブジェクトを見つけることに失敗しています。

data引数に指定されたオブジェクトに重みを追加するか、lm関数を呼び出すのと同じ環境で数式オブジェクトを作成して、エラーを修正できます。

解決方法1：

rm(list=ls()) 
library(faraway) 
library(boot) 
data(strongx) 
g3 <- lm(crossx ~ energy, strongx, weights=rep(1, nrow(strongx))) 

bootReg <- function(formulatostat, data, ind, weitostat) { 
    input_to_lm <- data[ind, ] 
    input_to_lm$weig <- weitostat[ind] 
    fit <- lm(formula=formulatostat, data = input_to_lm, weights = weig) 
    return(coef(fit)) 
} 
wei <- rep(1, nrow(strongx)) 

bootResults <- 
    boot(
    data = strongx, 
    statistic = bootReg, 
    R = 2000, 
    formulatostat = crossx ~ energy, 
    weitostat = wei) 

解決策2：

rm(list=ls()) 
library(faraway) 
library(boot) 
data(strongx) 
g3 <- lm(crossx ~ energy, strongx, weights=rep(1, nrow(strongx))) 

bootReg <- function(formulatostat, data, ind, weitostat) { 
    weig <- weitostat[ind] 
    fit <- lm(formula=as.formula(formulatostat), data = data[ind, ], weights = weig) 
    return(coef(fit)) 
} 
wei <- rep(1, nrow(strongx)) 

bootResults <- 
    boot(
    data = strongx, 
    statistic = bootReg, 
    R = 2000, 
    formulatostat = 'crossx ~ energy', 
    weitostat = wei)