求高手解答，万分感谢（岭回归）

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

执行以下代码，出来的岭迹图很奇怪，见下图1
> data<-read.csv("F:\\yang\\yang2.csv")
> data
> dim(data)
[1] 560  16
> names(data)
[1] "y"   "x1"  "x2"  "x3"  "x4"  "x5"  "x6"  "x7"  "x8"  "x9"  "x10" "x11"
[13] "x12" "x13" "x14" "x15"
> library(MASS)
> a<-lm.ridge(y~.,data,lambda=seq(0,1,0.05))
> a
> plot(lm.ridge(y~.,data,lambda=seq(0,1,0.05)))
> select(lm.ridge(y~.,data,lambda=seq(0,1,0.005)))
modified HKB estimator is 209.2754
modified L-W estimator is 184.3846
smallest value of GCV  at 1
> a$coef[,a$lam==1]
         x1          x2          x3          x4          x5          x6
   6.615381  -12.912217 -117.296826  -24.855233  -49.437182  -48.490056
         x7          x8          x9         x10         x11         x12
-47.642887   82.532577   -7.366161  -30.356692   48.048630   24.325980
        x13         x14         x15
  62.902631   12.088613   39.599596
> a$coef[,1]
         x1          x2          x3          x4          x5          x6
   6.717755  -12.982036 -117.589705  -24.943189  -49.610816  -48.558601
         x7          x8          x9         x10         x11         x12
-47.715001   82.723190   -7.352264  -30.372013   48.117154   24.544152
        x13         x14         x15
  62.973416   12.115731   39.679528

图1


以下是我用sas做的岭迹图（注：此次变量为16个，比用R多了X7）


想请教下，我的R代码是不是哪里有问题？？？
万分感谢，我纠结了一天了~o(╯□╰)o

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：求高手解答 万分感谢 求高手 岭回归 estimator library names

有可能SAS和R对于lamda的定义不同，也有可能他们中的某一个没有standardize

ntsean 发表于 2014-5-15 04:56
有可能SAS和R对于lamda的定义不同，也有可能他们中的某一个没有standardize

R的话怎么标准化呢，SAS是经过标准化的