作者andrew43 (讨厌有好心推文後删文者)
看板R_Language
标题Re: [问题] 取出Residual再跑一次Regression
时间Wed Dec 2 01:53:56 2015
当你建立 lm 物件时，该物件会保留有参与的资料内容。
y <- c(NA,1,2,3, 4,5, 6)
x1 <- c( 1,2,2,2, 2,3,NA)
x2 <- c( 3,5,3,4,NA,3 ,2)
m <- lm(y ~ x1 + x2)
m$model
y x1 x2
2 1 2 5
3 2 2 3
4 3 2 4
6 5 3 3
可以看出实际使用的资料只有 lm物件$model 的内容。
这样你即可以检查了。
※ 引述《developme223 (allen)》之铭言：
: [问题类型]:
: 程式谘询(我想用R 做某件事情，但是我不知道要怎麽用R 写出来)
: [软体熟悉度]:
: 入门(写过其他程式，只是对语法不熟悉)
: [问题叙述]:
: 请简略描述你所要做的事情，或是这个程式的目的
: 程式码：
: model.1 <- lm(TAC ~ 0 +TA +REC +PPE ,data=merge.data.3)
: model.2 <- lm(merge.data.3$RIR ~ model.1$residual)
: 其中，
: merge.data.3 大概资料型态：
: Code RIR TAC TA REC PPE
: 1 20 200 0.5 . .
: 2 30 . . . .
: .
: .
: .
: 533 . . . . .
: 先跑一个简单线性回归模型(model.1)，取出此模型之残差，
: 再跑一次简单线性回归(model.2)，但此时残差(model.1$residual)
: 的排列与原先的data.frame不同(无法执行，跑出Variable lengths differ)
: 有尝试用cbind与data.frame重组资料，但好像也无法解决，想请教各位高手
: 。
: [关键字]:
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 122.117.37.172
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/R_Language/M.1448992439.A.9FA.html