本人用无截距模型，作一元回归。同时用自己的函数和R中的函数lm运行结果对比，发现参数估计、残差等数据都完全吻合，但是R^2不同。我的是负数（因为此处用无截距方法，出现负数并不意外），下面是我计算的方法，由于这里不方便编辑公式，大家就领会精神吧，用到的都是最常见的那几个量。我运算的方法是1－RSS/SYY. 其中SYY=sum((y-ybar)^2), RSS=sum((y-yhat)^2)。希望各位大人解释一下，谢谢。