5396 1

[问答] 在R语言中使用ncvreg包解决SCAD惩罚回归,代码问题如下: [推广有奖]

  • 0关注
  • 0粉丝

博士生

25%

还不是VIP/贵宾

-

威望
0
论坛币
5806 个
通用积分
13.0050
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
108 点
帖子
57
精华
0
在线时间
480 小时
注册时间
2016-8-6
最后登录
2024-11-8

100论坛币
其中的X是T*N*P,Y是T*N,beta是P*1.library(ncvreg)
X=array(runif(N*T*P,-2,2),dim=c(T,N,P))
beta=c(4,5,0,0,2,0,0,-2,0,0)
betaa=array(beta,c(T,N,P))
u=apply(X*betaa,c(1,2),sum)
e=matrix(rnorm(T*N,0,1),nrow=T,ncol=N)
y=u+e
> scad_cv<-cv.ncvreg(X,y,penalty='SCAD')
Error in cv.ncvreg(X, y, penalty = "SCAD") :
X must be a matrix or able to be coerced to a matrix


scad_cv<-cv.ncvreg(X,y,penalty='SCAD',family='mgaussian')
Error in cv.ncvreg(X, y, penalty = "SCAD", family = "mgaussian") :
X must be a matrix or able to be coerced to a matrix



沙发
AMREDA 发表于 2019-5-23 10:31:08 |只看作者 |坛友微信交流群
cv.ncvreg(X[,,1],y[,1],penalty='SCAD') 是不是可以

使用道具

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-11-10 09:42