使用“`refund::prf`”时出现错误，没有模型拦截

Question

我试图在R中使用函数标量和标量回归来运行函数数据分析模型。但是通过删除拦截，我得到了以下错误(示例基于this discussion)，

library(refund)
data(DTI)
DTI1 <- DTI[DTI$visit==1 & complete.cases(DTI),]
par(mfrow=c(1,2))

fit_af <- pfr(pasat ~ -1 +sex + case + af(cca, k=c(5, 8), bs="ps"), data = DTI1)

#Error in str2lang(x) : <text>:1:9: unexpected symbol
#1: pasat~0 sex
#            ^
#Calls: pfr -> formula -> formula.character -> str2lang

如何从pfr模型中删除拦截？

Answer 1

这是个窃听器。在调试模式下运行prf()可以识别问题。当没有拦截时，它会：

if (!attr(tf, "intercept")) {
  newfrml <- paste(newfrml, "0", sep = "")
}

我认为这应该是固定的

if (!attr(tf, "intercept")) {
  newfrml <- paste(newfrml, "0 +", sep = "")
}

考虑报告这个堆栈溢出线程，即https://stackoverflow.com/q/72856108来打包维护人员Julia : julia.wrobel@cuanschutz.edu

注：使用最新的refund_0.1~26进行测试(2022-04-16发布给CRAN )。