我在这个归档的 r-help 线程上使用了 Terry T. 提供的 Deming 函数。我正在比较两种方法,所以我的数据如下所示:
y x stdy stdx
1 1.2 0.23 0.67
2 1.8 0.05 0.89
4 7.5 1.13 0.44
... ... ... ...
我已经完成了我的戴明回归(也称为“总最小二乘回归”),我得到了一个斜率和截距。我想得到一个相关系数,所以我开始计算。我手动输入了公式:
R2 <- function(coef,i,x,y,sdty){
predy <- (coef*x)+i
stdyl <- sum((y-predy)^2) ### The calculated std like if it was a lm (SSres)
Reelstdy <- sum(stdy) ### the real stdy from the data (SSres real)
disty <- sum((y-mean(y))^2) ### SS tot
R2 <- 1-(stdyl/disty) ### R2 formula
R2avecstdyconnu <- 1-(Reelstdy/disty) ### R2 with the known stdy
return(data.frame(R2, R2avecstdyconnu, stdy, Reelstdy))
}
这个公式有效并给了我输出。
- 两个哪个更有意义?(我个人认为他们两个都有偏见。)
- 有没有办法从总最小二乘回归中获得相关系数?
戴明回归的输出:
Call:
deming(x = Data$DS, y = Data$DM, xstd = Data$SES, ystd = Data$SEM, dfbeta = T)
Coef se(coef) z p
Intercept 0.3874572 0.2249302 3.1004680 2.806415e-10
Slope 1.2546922 0.1140142 0.8450883 4.549709e-02
Scale= 0.7906686
>