`rms::ols()`：如何在没有截距的情况下拟合模型

Question

我想使用rms 包中的ols()（普通最小二乘）函数进行多元线性回归，但我不希望它计算截距。使用lm() 的语法如下：

model <- lm(formula = z ~ 0 + x + y, data = myData)

0 阻止它计算截距，并且只返回两个系数，一个用于x，另一个用于y。使用ols() 时如何执行此操作？ 尝试

model <- ols(formula = z ~ 0 + x + y, data = myData)

没有用，它仍然为x 和y 分别返回一个截距和一个系数。

Here is a link to a csv file

它有五列。对于本例，只能使用前三列：

model <- ols(formula = CorrEn ~ intEn_anti_ncp + intEn_par_ncp, data = ccd)

谢谢！

Answer 1

rms::ols 使用 rms:::Design 而不是 model.frame.default。 Design 是用默认的intercept = 1 调用的，所以没有（明显的）方法来指定没有拦截。我认为这是有充分理由的，但您可以尝试使用 trace 更改 ols。