hdm包中的rlassoEffects语法不起作用?

hdm包中的rlassoEffects语法不起作用?,r,R,我遵循hdm包中的单套索和双套索程序,但不知何故,双套索对我来说不能正常工作。我将使用包提供的数据集,因为除了我自己的数据集之外,它也不适合我 如本文所述,我加载GrowthData数据集,然后尝试实现double LASSO过程: library(hdm) DS=rlassoEffects(Outcome~. , I=~gdpsh465, data=GrowthData) summary(DS) 但是,我收到以下错误消息: > DS=rlassoEffects(Outcome~. ,

我遵循hdm包中的单套索和双套索程序,但不知何故,双套索对我来说不能正常工作。我将使用包提供的数据集,因为除了我自己的数据集之外,它也不适合我

如本文所述,我加载GrowthData数据集,然后尝试实现double LASSO过程:

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版本