我正在尝试使用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)
但是当我运行真实数据时,我收到以下错误:
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)
if(abs((df - current.df)/ df)< 1e-04 | iterations == 1)出错 return(list(lambda = lambda,:缺少值,其中TRUE / FALSE 需要
您可以访问将重现此结果here的示例数据文件。有关如何解决此错误的任何想法?