我正在使用glmnet包在R中工作。这是一个我可以做脊和包的套餐。套索回归。在这个R包中,代码与fortran代码连接,这是整个计算内容发生的地方。 我的目标是在此代码中更改某些内容,然后尝试将新的fortran代码与r包连接起来(即计划)。
现在我的问题:
我找到了fortran代码glmnet5.f90,所以我得出结论,这段代码是一个免费的格式文件。当我运行此代码时,我得到以下错误:
error 399 - Unrecognised Statement, did you Forget to add \FIXED_FORMAT?
在互联网上我发现类似于将代码复制到固定格式文件中但我只会得到更多错误。
我注意到的另一件事是评论是用" C"在第一列,但fortran 90使用"!"对于评论,也许它是一个强大的77代码,但我不知道。编译器是90/95,所以这不应该是一个问题吗?
现在我正在使用silverfrost的FTN95编译器。代码有超过7000行,这就是为什么我会给你代码的链接。 CRAN-package-glmnet 你可以在glmnet \ src \ glmnet5.f90找到代码。