线性回归

线性回归,多么熟悉的名字,我一开始真的以为就选变量回归,看显著性得出结论就好,我真的没想到线性回归要求这么多......

在线性回归的开始应该先总结一下假设检验,但我觉得可能百度更好一点。在这里说一下分类:T检验,Z检验 ,F检验,卡方检验

再说一下模型结果里面包括的成分:

  • estimators(可以适用均值,方差)
  • p值(自己定显著度alpha)什么都不懂的我都会用p值看哪个显著哪个不显著....有时候p值不一定够,sample size增加p也增加,此时发生第一类error。所以得有点东西和p值一起看,例如CI(confidence interval), Bayes factor......
  • 置信区间(CI)因为p值可以判断显著性,却无法判断组间差异,特别是在逻辑回归中一定要多考虑CI与p值的关系

进入正题,线性回归

linear regression
包括一个continuous dependent variable(Y), independent variable(X), and disturbance term or error variable ε ,线性回归的目标就是确定βi 使得ε 最小

介绍三个概念:

  1. Heteroscadasticity(异方差性):若随机误差项ε 的方差非常数,则随机误差项有异方差性。此时我们需要考虑是否遗漏了解释变量,是否将非线性的关系强价为了线性,样本数据间是否存在误差。
  2. Multicollinearity(多重共线性):即X之间存在线性关系。在R中用vif去检验,越接近1存在多重共线性的可能越小,大于5有可能寻在多重共线性,大于10寻在多重共线性可能性极高。所以在选择解释变量的时候要仔细考虑。
  3. Homogeneity(内生性): 内生性是在线性回归中不想看到的事情。内生性的根源有四点。1⃣️X与Y之间互为因果,x变动y也变动,y变动同时对x产生影响,永无止境...2⃣️X变化对X以外的因素也改变(这些因素未纳入解释变量而存在于ε中),同时这些其他因素也在影响着Y。线性回归要求X严格外生,要防止内生性的产生。但是现在我判断内生性的就是好好看自己的数据。

我们都知道线性回归适用最小二乘法回归使得预测值与观察值的区别最小,但是最小二乘法使用时有一些条件

  • ε满足独立同分布,均值为零,方差不变(即不出现异方差性)的正态分布(正态分布可以忽略)四个条件
    p.s 当 Cov[εi,εj] = 0 则误差项满足统计独立条件
  • 严格外生
  • 解释变量(X)与预测变量(Y)符合正态分布
    我们做线性回归的时候最好满足这些条件
R里面一般是怎么检查这些条件的呢?
image.png
  1. Residuals and fitted

    我们希望在这张图上看到,红线表现出无线性(non-linear trend)并且与虚线相接近,同时residual随机散布在线两侧。若红线呈现出二次函数线性,则证明我们的回归函数中需要一个二次项调整。 若出现下图
    Heteroscadasticity.png
    则证明不满足同方差性假设
  2. Normal Q-Q
    检验样本与样本期望的正态性对比,两分布若相似,点就会分布在虚线(diagonal line)上,此时ε满足正态性假设

3.Scale-location
这个图也是用来验证同方差性假设的。若点沿线随机分布在两侧以及前后则满足不变方差假设。

  1. Residual and Leverage
    这个图主要是看outliers,想要进一步研究outliers可以使用cook‘s distance

5.方差分析
检查方差的内生性可以使用bartlett test, levene test, fingner-killeen test, brown-forsyth test.

  1. 相关系数检验
    相关系数检验可以让我们在选择解释变量时剃掉与y不相关的x,在线性回归中我们常用person collrelation coefficient

线性回归结果解释

要看模型拟合的好不好,大家都知道看R-square, 越大越好。当选择的解释变量数量有变是看adjusted r-square, 若增加解释变量数目但不会改善模型效果的项会增加惩罚项。
依旧,最后放代码

LRM3 <- lm(TTime ~ Age + AgeSq + NumVeh + LOB , data = LM_data3)
#show results plots in 2*2 format
par(mfrow=c(2,2))
plot(LRM3)
summary(LRM3)
# show confidence intervals of variables
confint(LRM3) 
# Calculates the variation inflation factors of all predictors in regression models
vif(LRM3)
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 214,444评论 6 496
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 91,421评论 3 389
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 160,036评论 0 349
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 57,363评论 1 288
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 66,460评论 6 386
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 50,502评论 1 292
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 39,511评论 3 412
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 38,280评论 0 270
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 44,736评论 1 307
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 37,014评论 2 328
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 39,190评论 1 342
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 34,848评论 5 338
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 40,531评论 3 322
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 31,159评论 0 21
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,411评论 1 268
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 47,067评论 2 365
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 44,078评论 2 352

推荐阅读更多精彩内容