Confidence intervals for Ridge regression

Ana Laura Carreiras

I can't do the confidence intervals in a ridge regression. I have this model.

model5 <- glmnet(train_x,train_y,family = "gaussian",alpha=0, lambda=0.01)

And when I do the prediction I use these command:

test_pred <- predict(model5, test_x, type = "link")

Someone knows how to do the confidence interval for the predictions?

ilanman

It turns out that glmnet doesn't offer standard errors (and therefore doesn't give you confidence intervals) as explained here and also addressed in this vignette (excerpt below):

It is a very natural question to ask for standard errors of regression coefficients or other estimated quantities. In principle such standard errors can easily be calculated, e.g. using the bootstrap.

Still, this package deliberately does not provide them. The reason for this is that standard errors are not very meaningful for strongly biased estimates such as arise from penalized estimation methods. Penalized estimation is a procedure that reduces the variance of estimators by introducing substantial bias. The bias of each estimator is therefore a major component of its mean squared error, whereas its variance may contribute only a small part.

Unfortunately, in most applications of penalized regression it is impossible to obtain a sufficiently precise estimate of the bias. Any bootstrap-based calculations can only give an assessment of the variance of the estimates. Reliable estimates of the bias are only available if reliable unbiased estimates are available, which is typically not the case in situations in which penalized estimates are used.

Reporting a standard error of a penalized estimate therefore tells only part of the story. It can give a mistaken impression of great precision, completely ignoring the inaccuracy caused by the bias. It is certainly a mistake to make confidence statements that are only based on an assessment of the variance of the estimates, such as bootstrap-based confidence intervals do.

Reliable confidence intervals around the penalized estimates can be obtained in the case of low dimensional models using the standard generalized linear model theory as implemented in lm, glm and coxph. Methods for constructing reliable confidence intervals in the high-dimensional situation are, to my knowledge, not available.

However, if you insist on confidence intervals, check out this post.

この記事はインターネットから収集されたものであり、転載の際にはソースを示してください。

侵害の場合は、連絡してください[email protected]

編集
0

コメントを追加

0

関連記事

分類Dev

ILNumerics:ILMath.ridge_regression

分類Dev

Compute a kernel ridge regression in R for model selection

分類Dev

Questions about ridge regression on python : Scaling, and interpretation

分類Dev

scikit-learn Ridge Regression UnboundLocalError

分類Dev

Plotting confidence and prediction intervals with repeated entries

分類Dev

Forest plot using risk ratio and confidence intervals

分類Dev

How to compute confidence intervals and plot them on a bar plot

分類Dev

How to plot confidence intervals for glm models (gamma family)?

分類Dev

How to understand F-test based lmfit confidence intervals

分類Dev

How to calculate the 99% confidence interval for the slope in a linear regression model in python?

分類Dev

Ridge / Lasso Regressionのh(simpleError(msg、call))のエラー

分類Dev

Logistic Regression Model using Regularization (L1 / L2) Lasso and Ridge

分類Dev

Ridge / Lasso Regressionのh(simpleError(msg、call))のエラー

分類Dev

Writing a for-loop in R that populates a matrix with the output of a correlation test (estimates and confidence intervals)?

分類Dev

change line colors of confidence intervals (geom_line, line plot) manually

分類Dev

How to plot regression transformed back on original scale with colored confidence interval bands?

分類Dev

How to aggregate confidence levels in csv?

分類Dev

Stata - Scatter plot confidence interval

分類Dev

Generate dynamic custom intervals

分類Dev

Maximize sum of array intervals

分類Dev

Sort DataFrame on column of intervals

分類Dev

merge timestamps intervals in postgres

分類Dev

De-overlapping intervals

分類Dev

Creating 10% intervals for displacement

分類Dev

for loop skipping intervals javascript

分類Dev

ggplot boxplot with mean and confidence interval by group

分類Dev

T-distribution with confidence level in Javascript

分類Dev

Plot variance and confidence interval for variance with ggplot

分類Dev

Confidence interval for all rows of a dataframe in R

Related 関連記事

  1. 1

    ILNumerics:ILMath.ridge_regression

  2. 2

    Compute a kernel ridge regression in R for model selection

  3. 3

    Questions about ridge regression on python : Scaling, and interpretation

  4. 4

    scikit-learn Ridge Regression UnboundLocalError

  5. 5

    Plotting confidence and prediction intervals with repeated entries

  6. 6

    Forest plot using risk ratio and confidence intervals

  7. 7

    How to compute confidence intervals and plot them on a bar plot

  8. 8

    How to plot confidence intervals for glm models (gamma family)?

  9. 9

    How to understand F-test based lmfit confidence intervals

  10. 10

    How to calculate the 99% confidence interval for the slope in a linear regression model in python?

  11. 11

    Ridge / Lasso Regressionのh(simpleError(msg、call))のエラー

  12. 12

    Logistic Regression Model using Regularization (L1 / L2) Lasso and Ridge

  13. 13

    Ridge / Lasso Regressionのh(simpleError(msg、call))のエラー

  14. 14

    Writing a for-loop in R that populates a matrix with the output of a correlation test (estimates and confidence intervals)?

  15. 15

    change line colors of confidence intervals (geom_line, line plot) manually

  16. 16

    How to plot regression transformed back on original scale with colored confidence interval bands?

  17. 17

    How to aggregate confidence levels in csv?

  18. 18

    Stata - Scatter plot confidence interval

  19. 19

    Generate dynamic custom intervals

  20. 20

    Maximize sum of array intervals

  21. 21

    Sort DataFrame on column of intervals

  22. 22

    merge timestamps intervals in postgres

  23. 23

    De-overlapping intervals

  24. 24

    Creating 10% intervals for displacement

  25. 25

    for loop skipping intervals javascript

  26. 26

    ggplot boxplot with mean and confidence interval by group

  27. 27

    T-distribution with confidence level in Javascript

  28. 28

    Plot variance and confidence interval for variance with ggplot

  29. 29

    Confidence interval for all rows of a dataframe in R

ホットタグ

アーカイブ