在R语言中使用ncvreg包解决SCAD惩罚回归,代码问题如下:
1 个回复 - 5387 次查看
其中的X是T*N*P,Y是T*N,beta是P*1.library(n
cvreg)
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) ...
2017-10-31 16:46 - 忠义似君稀1 - R语言论坛