0

I'm an R newbie. I'm using "clogitL1" to run a regularized conditional logistic regression for a matched case-control study with 1021 independent variables (metabolites). I'm not able to extract the regression coefficients. I've tried summary(x), coef(x), coefficient(x), x$beta - none of them work. I'm able to run it OK, and if I follow it with "cv.clogitL1" I can extract the cross-validated estimated coefficients, but not the estimated coefficients for the original model. Here's some of my code:

strata=sort(data.meta$MATCHED_NEW)
condlog <- clogitL1(y=data.meta$BCR, x=data.meta$ln_metab[, data.features ], strata,
  numLambda=100, minLambdaRatio=0.000001, alpha = 1.0)

"strata" is a vector indicating pairing of cases & controls. "data.meta$BCR" is a vector inicating case or control status "data.meta$ln_metab" is a matrix with observations as rows and metabolite levels as columns "data.features" is a vector indicating which metabolites passed several dimension reduction filters.

Appreciate any suggestions.

MrFlick
  • 195,160
  • 17
  • 277
  • 295
bruce_t
  • 1
  • 1
    It's easier to help you if you include a simple [reproducible example](https://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example) with sample input and desired output that can be used to test and verify possible solutions. – MrFlick Apr 06 '21 at 19:37
  • Are you literally typing summary(x) instead of summary(condlog )? – Elin Apr 06 '21 at 20:31
  • No - I'm typing summary(condlog) - should have clarified that. – bruce_t Apr 09 '21 at 19:33

0 Answers0