我正在尝试使用gamsel R程序包来适应稀疏的广义加性模型,但似乎无法使它用于实际数据。当我按照软件包文档中的描述运行综合数据时,一切运行良好:
library(gamsel)
data=gendata(n=500,p=12,k.lin=3,k.nonlin=3,deg=8,sigma=0.5)
attach(data)
bases=pseudo.bases(X,degree=10,df=6)
gamsel.out=gamsel(X,y,bases=bases)
Run Code Online (Sandbox Code Playgroud)
但是当我在真实数据上运行时,出现以下错误:
library(gamsel)
X = as.matrix(read.csv("X.csv"),header=FALSE)
y = as.matrix(read.csv("y.csv"),header=FALSE)
gam_fit = gamsel(X,y)
Run Code Online (Sandbox Code Playgroud)
if(abs((df-current.df)/ df)<1e-04 |迭代== 1)中的错误return(list(lambda = lambda,:缺少需要TRUE / FALSE的值)
您可以在此处访问将重现此结果的示例数据文件。关于如何解决此错误的任何想法?
r ×1