2014-10-21 2 views
0

나는 희박한 그룹 인 올가미를 사용하고 있습니다. 이것은 페널티가있는 회귀입니다. 내가 사용하고있는 패키지는 SGL이다. 내 R의 예제를 실행하려고하고, 코드는 내가 cvFit의 회귀 계수를 추출하려고불이익 회귀에서 계수 추출

set.seed(1) 
n = 50; p = 100; size.groups = 10 
index <- ceiling(1:p/size.groups) 
X = matrix(rnorm(n * p), ncol = p, nrow = n) 
beta = (-2:2) 
y = X[,1:5] %*% beta + 0.1*rnorm(n) 
data = list(x = X, y = y) 
cvFit = cvSGL(data, index, type = "linear") 

이하로 제공되지만

coef(cvFit) 
NULL 

사람이 말할 수로 밝혀 나 뭐가 잘못 됐니? 미리 감사드립니다.

+1

'str (cfFit)'을 참조하고 계수의 정확한 위치를 결정하십시오. –

+0

@ Roman Luštrik 계수의 정확한 위치를 결정한다는 것은 무엇을 의미합니까? – Xst67

+0

이 함수에 익숙하지 않고'cvFit'을 타이핑하면 어떤 값을 추출 할 수 있는지 알 수 있습니다. 즉'cvFit $ fit $ beta' 또는'cvFit $ lambdas' – user20650

답변

0

최소 람다 값으로 모델에서 계수를 추출합니다.

coef(fit,s=cvfit$lambda.min)