楼主: ltx5151
20512 8

R软件lm函数中R Square是如何计算的? [推广有奖]

  • 0关注
  • 26粉丝

VIP

教授

26%

还不是VIP/贵宾

-

威望
1
论坛币
48929 个
通用积分
63.8334
学术水平
100 点
热心指数
143 点
信用等级
68 点
经验
2871 点
帖子
529
精华
0
在线时间
1312 小时
注册时间
2009-5-7
最后登录
2024-5-26

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信:zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

感谢您参与论坛问题回答

经管之家送您两个论坛币!

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

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

关键词:Square ARE r软件 参数估计 无截距 软件 函数 Square

回帖推荐

ltx5151 发表于6楼  查看完整内容

5# vinkwai SYY=t(y-ybar)%*%(y-ybar) RSS=t(y-yfit)%*%(y-yfit) R-Square=1-RSS/SYY ybar是y的均值,yfit是拟合值,这些量都和lm结果相同,只有R-square不同

simplystat 发表于3楼  查看完整内容

>modelsummary(model)["r.squared"] >summary(model)[c("r.squared", "adj.r.squared")]

本帖被以下文库推荐

沙发
trier2006 发表于 2010-9-25 09:29:39 |只看作者 |坛友微信交流群
友情帮顶。
最好的医生是自己,最好的药物是时间……

使用道具

藤椅
simplystat 发表于 2010-9-25 12:34:37 |只看作者 |坛友微信交流群
>model<-lm(y~x-1, data=data)
>summary(model)["r.squared"]
>summary(model)[c("r.squared", "adj.r.squared")]

使用道具

板凳
ltx5151 发表于 2010-9-25 13:10:41 |只看作者 |坛友微信交流群
3# simplystat
谢谢楼上的回答,可是麻烦看清问题。我不是不会显示lm中的R-square,是不理解我计算和程序中值的差异。

使用道具

报纸
vinkwai 发表于 2010-9-25 13:37:12 |只看作者 |坛友微信交流群
4# ltx5151

Please post your code here

使用道具

地板
ltx5151 发表于 2010-9-25 14:35:22 |只看作者 |坛友微信交流群
5# vinkwai
SYY=t(y-ybar)%*%(y-ybar)
RSS=t(y-yfit)%*%(y-yfit)
R-Square=1-RSS/SYY

ybar是y的均值,yfit是拟合值,这些量都和lm结果相同,只有R-square不同

使用道具

7
ltx5151 发表于 2010-9-26 15:36:59 |只看作者 |坛友微信交流群
自己顶一下

使用道具

8
kk22boy 发表于 2011-4-6 21:14:48 |只看作者 |坛友微信交流群
帮顶一下。。
如果该贴对您有些许帮助,希望你能回复一下或者评一下热心指数!谢谢!

使用道具

9
shenbaiseshatan 在职认证  发表于 2011-4-7 20:34:03 |只看作者 |坛友微信交流群
要不看下这个函数的代码呢,不就找到它的算法了吗
胜人者有力,自胜者强!

使用道具

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-6-3 23:11