计算器 r-adj r squared (adj)

苹果/安卓/wp
积分 2564, 距离下一级还需 1036 积分
权限: 自定义头衔, 签名中使用图片, 设置帖子权限, 隐身, 设置回复可见
道具: 彩虹炫, 涂鸦板, 雷达卡, 热点灯, 金钱卡, 显身卡, 匿名卡, 抢沙发, 提升卡, 沉默卡下一级可获得
道具: 千斤顶
购买后可立即获得
权限: 隐身
道具: 金钱卡, 彩虹炫, 雷达卡, 热点灯, 涂鸦板
上火签到天数: 5 天连续签到: 1 天[LV.2]偶尔看看I
本人用无截距模型,作一元回归。同时用自己的函数和R中的函数lm运行结果对比,发现参数估计、残差等数据都完全吻合,但是R^2不同。我的是负数(因为此处用无截距方法,出现负数并不意外),下面是我计算的方法,由于这里不方便编辑公式,大家就领会精神吧,用到的都是最常见的那几个量。我运算的方法是1-RSS/SYY. 其中SYY=sum((y-ybar)^2), RSS=sum((y-yhat)^2)。希望各位大人解释一下,谢谢。
载入中......
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不同
>modelsummary(model)["r.squared"]
>summary(model)[c("r.squared", "adj.r.squared")]
本帖被以下文库推荐
& |主题: 828, 订阅: 12
友情帮顶。
最好的医生是自己,最好的药物是时间……
&model&-lm(y~x-1, data=data)
&summary(model)[&r.squared&]
&summary(model)[c(&r.squared&, &adj.r.squared&)]
simplystat
谢谢楼上的回答,可是麻烦看清问题。我不是不会显示lm中的R-square,是不理解我计算和程序中值的差异。
Please post your code here
SYY=t(y-ybar)%*%(y-ybar)
RSS=t(y-yfit)%*%(y-yfit)
R-Square=1-RSS/SYY
ybar是y的均值,yfit是拟合值,这些量都和lm结果相同,只有R-square不同
自己顶一下
帮顶一下。。
如果该贴对您有些许帮助,希望你能回复一下或者评一下热心指数!谢谢!
要不看下这个函数的代码呢,不就找到它的算法了吗
胜人者有力,自胜者强!
初级热心勋章
初级热心勋章
初级信用勋章
初级信用勋章
无限扩大经管职场人脉圈!每天抽选10位免费名额,现在就扫& 论坛VIP& 贵宾会员& 可免费加入
加入我们,立即就学扫码下载「就学」app& Join us!& JoinLearn&
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
如有投资本站或合作意向,请联系(010-);
邮箱:service@pinggu.org
投诉或不良信息处理:(010-)
京ICP证090565号
京公网安备号
论坛法律顾问:王进律师

我要回帖

更多关于 r square 多大 有意义 的文章

 

随机推荐