本帖最后由 editor 于 2013-3-6 21:01 编辑
来自R语言QQ群的问题:
童鞋A: 19:18:56
想问个问题,逐步回归step()函数的参数具体解释里有个scope该如何设定?
童鞋A: 19:20:18
step(object, scope, scale = 0,
direction = c("both", "backward", "forward"),
trace = 1, keep = NULL, steps = 1000, k = 2, ...)
童鞋A: 19:27:25
matlab采用的是是stepwisefit()函数,按照R的step()定义,会给出 AIC 检验后的最优参数,我找了个教科书里的例子,用stepwisefit()选出的则是x1,x2,而用step()函数,则选出的是x2,x3,x4,我就在想应该是在参数设置上有些问题,想请教各位,如何设定参数,才能得到一样的结果呢?
童鞋B: 19:27:46
default
一般都不设
童鞋A: 19:29:01
我用step()函数,就是默认,但他选出的是x2,x3,x4,而教科书上是用stepwisefit()函数,他选出的x1,x2,我就想搞清楚,如何得到和matlab一样的结果呢?
童鞋B: 19:29:32
both,forward,backward你设好没?
童鞋A: 19:30:42
stepwisefit(X,Y,'PENTER',0.08)
而R语言是:
sales<-data.frame(y, x1, x2,x3,x4)
lm.reg<-lm(y~x1+x2+x3+x4, data=sales)
lm.step<-step(lm.reg)
summary(lm.step)
童鞋B: 19:31:16
要选 both 还是forward还是backward哎
这个也要和书上一致
童鞋A: 19:31:33
我用的是默认的,其他没有设定,用help找不到scope的设定方法,所以请教各位大人
书上用的是matlab计算的,
我想用R语言计算一下,
我在baidu一下
|