Welcome to OGeek Q&A Community for programmer and developer-Open, Learning and Share
Welcome To Ask or Share your Answers For Others

Categories

0 votes
373 views
in Technique[技术] by (71.8m points)

incorrect calculating (R-squared) in R(wrong value)


与恶龙缠斗过久,自身亦成为恶龙;凝视深渊过久,深渊将回以凝视…
Welcome To Ask or Share your Answers For Others

1 Reply

0 votes
by (71.8m points)

It is likely that x is not the predictions but is a predictor that goes into a linear regression. Perform the regression, fm, in which case the predicted values are fitted(fm) and then get the R squared from summary or get it directly as shown in the alternatives.

fm <- lm(yield ~ x, mydat)
summary(fm)$r.squared
# [1] 0.02508245

# same
cor(mydat$yield, fitted(fm))^2
# [1] 0.02508245

# same
with(mydat, cor(yield, x)^2)
# [1] 0.02508245

# same
tss <- with(mydat, sum((yield - mean(yield))^2))
rss <- deviance(fm)
1 - rss/tss
# [1] 0.02508245

# same
tss <- with(mydat, sum((yield - mean(yield))^2))
rss <- sum(resid(fm)^2)
1 - rss/tss
# [1] 0.02508245


plot(yield ~ x, mydat)
abline(fm)

screenshot


与恶龙缠斗过久,自身亦成为恶龙;凝视深渊过久,深渊将回以凝视…
OGeek|极客中国-欢迎来到极客的世界,一个免费开放的程序员编程交流平台!开放,进步,分享!让技术改变生活,让极客改变未来! Welcome to OGeek Q&A Community for programmer and developer-Open, Learning and Share
Click Here to Ask a Question

...