我使用的都是R,计算结果的确不一样,给你数据集,你自己试试就知道了。一个做出来的是0.5231816,两外一个结果是0.4466259。你给的公式只是针对前一种情况。
data<-scan("shidata.txt",list(t=0,k=0))
x<-data$t
y<-data$k
#a<-RMA(x,y)
#x<-seq(1,30)
#y<--25*x+300+rnorm(30,0,5)
SSy<-sum((y-mean(y))^2)
SSx<-sum((x-mean(x))^2)
SSxy<-sum((y-mean(y))*(x-mean(x)))
slope<--sqrt(SSy/SSx)
intercept<-mean(y)-slope*mean(x)
Rsquare<-SSxy^2/(SSx*SSy)
Rsquare
Rpf<-1-sum((intercept+slope*x-y)^2)/sum((y-mean(y))^2)
Rpf
cor(x,y)^2
附件列表