hdm包中的rlassoEffects语法不起作用?
我遵循hdm包中的单套索和双套索程序,但不知何故,双套索对我来说不能正常工作。我将使用包提供的数据集,因为除了我自己的数据集之外,它也不适合我 如本文所述,我加载GrowthData数据集,然后尝试实现double LASSO过程:hdm包中的rlassoEffects语法不起作用?,r,R,我遵循hdm包中的单套索和双套索程序,但不知何故,双套索对我来说不能正常工作。我将使用包提供的数据集,因为除了我自己的数据集之外,它也不适合我 如本文所述,我加载GrowthData数据集,然后尝试实现double LASSO过程: library(hdm) DS=rlassoEffects(Outcome~. , I=~gdpsh465, data=GrowthData) summary(DS) 但是,我收到以下错误消息: > DS=rlassoEffects(Outcome~. ,
library(hdm)
DS=rlassoEffects(Outcome~. , I=~gdpsh465, data=GrowthData)
summary(DS)
但是,我收到以下错误消息:
> DS=rlassoEffects(Outcome~. , I=~gdpsh465, data=GrowthData)
Error in colnames(X) : object 'X' not found
> summary(DS)
Error in summary(DS) : object 'DS' not found
我知道为函数提供了两种不同的语法,所以对“object X”的引用似乎是指另一种语法
但是,为什么我使用的语法样式不起作用
谢谢 事实上,我在同一篇文章中偶然发现了同样的问题。请注意,rlassoEffects只是rlassoEffects的包装器。您可以尝试通过以下方式实现它 编辑: 或者不用包装就可以使用
rlassoEffect(data.matrix(dataset[, -c(1:2)]), unlist(dataset[, 1]), unlist(dataset[, 2]), method = "double selection")
奇怪的是,与你引用的文章相比,这给出了一个稍微不同的系数值
编辑2:
我找到了可能导致该行为的代码行:
I.c <- which(colnames(X) %in% cn[I.c])
I.c谢谢你的提问。是的,没错。指示行中有一个错误。我们已经在开发版本中修复了它,该版本可从中获得
您可以通过以下命令安装开发版本
devtools::install\u github(“MartinSpindler/hdm”)
你的例子现在应该有用了
我们将很快更新CRAN版本