我正在编写一些代码,我注意到一些奇怪的事情。
当我在某些面板数据的子集上运行 LM 时,它工作正常,如下所示:
library('plm')
data(Cigar)
lm(log(price) ~ log(pop) + log(ndi), data=Cigar, subset=Cigar$state==1)
Call:
lm(formula = log(price) ~ log(pop) + log(ndi), data = Cigar,
subset = Cigar$state == 1)
Coefficients:
(Intercept) log(pop) log(ndi)
-26.4919 3.2749 0.4265
但是当我尝试将其包装在函数中时,我得到:
myfunction <- function(formula, data, subset){
return(lm(formula, data, subset))
}
myfunction(formula = log(price) ~ log(pop) + log(ndi), data = Cigar, subset = Cigar$state==1)
Error in xj[i] : invalid subscript type 'closure'
我真的不明白这里发生了什么,但它破坏了我编写的其他一些代码,所以我想知道。
问题似乎不在于子集。当我更改为时,我从你的函数中得到同样的错误subset = (state == 1)
。您的函数的参数未正确传递和评估。
我认为你最好使用do.call
myfunction <- function(formula, data, subset) {
do.call("lm", as.list(match.call()[-1]))
}
myfunction(log(price) ~ log(pop) + log(ndi), Cigar, state == 1)
# Call:
# lm(formula = log(price) ~ log(pop) + log(ndi), data = Cigar,
# subset = state == 1)
#
# Coefficients:
# (Intercept) log(pop) log(ndi)
# -26.4919 3.2749 0.4265
本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)