如何循环使用glm

我想对ridge和lasso进行100次循环以获得100个mse和mspe。我的最终目标是绘制一个箱线图来比较这些100个值。我已经建立了一个回归模型,但不知道如何重复这个模型。我该如何获取这些值并绘制箱线图呢?


回答:

你可以尝试以下方法:

ntimes <- 100res <- replicate(ntimes, {  cv.rr <- cv.glmnet(x=as.matrix(train[,-1]),y=as.numeric(train[,1]),alpha=0,nfolds=10,nlambda=100, intercept=FALSE)    lambda.rr=cv.rr$lambda.min        mse.rr <- mean((coef(cv.rr)[-1] - betas.true)^2)      yhat.rr <- predict(cv.rr,s="lambda.min",newx=as.matrix(test[,-1]))  mspe.rr <- mean((test[,1]-yhat.rr)^2)    list(mse=mse.rr, mspe=mspe.rr)  })library(tidyverse)res_df <- as.data.frame(apply(res, 1, function(x) unlist(x)))names(res_df) <- c('mse', 'mspe')res_df %>% gather(key='metric', value='value') %>% ggplot(aes(value, fill=metric)) + geom_boxplot()

以获得如下所示的可视化效果:

enter image description here

Related Posts

使用LSTM在Python中预测未来值

这段代码可以预测指定股票的当前日期之前的值,但不能预测…

如何在gensim的word2vec模型中查找双词组的相似性

我有一个word2vec模型,假设我使用的是googl…

dask_xgboost.predict 可以工作但无法显示 – 数据必须是一维的

我试图使用 XGBoost 创建模型。 看起来我成功地…

ML Tuning – Cross Validation in Spark

我在https://spark.apache.org/…

如何在React JS中使用fetch从REST API获取预测

我正在开发一个应用程序,其中Flask REST AP…

如何分析ML.NET中多类分类预测得分数组?

我在ML.NET中创建了一个多类分类项目。该项目可以对…

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注