在`bwplot`中对`caret`模型进行排序

如何解决在`bwplot`中对`caret`模型进行排序

我正在绘制使用caret训练的年度模型的重采样准确性得分的箱线图。 这些模型以其引用的年份命名:2000、2001、2002,...,2010。 我希望模型以年份(即模型名称)为基础,以升序显示在方框图中。

基于以下代码的重采样摘要

fit.year.res <- resamples(fit.year)
summary(fit.year.res)

看起来像这样:

enter image description here

但是,箱形图中不同的年度模型未排序:

scales <- list(x=list(relation="free"),y=list(relation="free"))
bwplot(fit.year.res,scales=scales)

enter image description here

我曾尝试将重采样fit.year.res$models的models元素转换为角色特征,但这没什么区别。

解决方法

我不知道使用插入符号包中的bwplot方法的简单解决方案。也许只有一种,但是我缺乏格子技巧。我建议使用ggplot2手动绘制方框图。这样,您将可以更好地控制最终情节。

由于您没有发布带有数据的示例,因此我将使用?caret:::bwplot.resamples

中的示例之一
library(caret)
library(party)
library(RWeka)

load(url("http://topepo.github.io/caret/exampleModels.RData"))

resamps <- resamples(list(CART = rpartFit,CondInfTree = ctreeFit,MARS = earthFit))

bwplot(resamps,metric = "RMSE")

产生:

enter image description here

要使用ggplot手动制作图,您将需要一些数据操作:

library(tidyverse)
resamps$values %>% #extract the values
  select(1,ends_with("RMSE")) %>% #select the first column and all columns with a name ending with "RMSE"
  gather(model,RMSE,-1) %>% #convert to long table
  mutate(model = sub("~RMSE","",model)) %>% #leave just the model names
  ggplot()+ #call ggplot
  geom_boxplot(aes(x = RMSE,y = model)) -> p1 #and plot the box plot

p1

enter image description here

要在y轴上设置特定顺序:

p1 +
  scale_y_discrete(limits = c("MARS","CART","CondInfTree"))

enter image description here

如果您喜欢晶格

library(lattice)

resamps$values %>%
  select(1,ends_with("RMSE")) %>%
  gather(model,-1) %>%
  mutate(model = sub("~RMSE",model)) %>%
  {bwplot(model ~ RMSE,data = .)}

enter image description here

更改顺序可更改模型级别(此方法也适用于ggplot2):

resamps$values %>%
  select(1,model),model = factor(model,levels = c("MARS","CondInfTree"))) %>%
    {bwplot(model ~ RMSE,data = .)}

enter image description here

,

函数bwplot.resamples用于生成该图,如果您查看underlying code,变量将根据它们在关注指标下的平均性能进行分解。

下面有执行分解的相关代码:

bwplot.resamples <- function (x,data = NULL,models = x$models,metric = x$metric,...)
{
....
  avPerf <- ddply(subset(plotData,Metric == metric[1]),.(Model),function(x) c(Median = median(x$value,na.rm = TRUE)))
  avPerf <- avPerf[order(avPerf$Median),]

    ......
}

我想您需要做的是手动绘制情节:

data(BloodBrain)
gbmFit <- train(bbbDescr[,-3],logBBB,"gbm",tuneLength=6,trControl = trainControl(method = "cv"),verbose=FALSE)
     
glmnetFit <- train(bbbDescr[,"glmnet",trControl = trainControl(method = "cv"))

rfFit <- train(bbbDescr[,"rf",trControl = trainControl(method = "cv"))

knnFit <- train(bbbDescr[,"knn",trControl = trainControl(method = "cv"))

resamps <- resamples(list(gbm = gbmFit,glmnet=glmnetFit,knn=knnFit,rf=rfFit))

如果进行绘图,您会看到它们根据中位数(实心点)进行了排序:

bwplot(resamps,metric="MAE")

enter image description here

您可以访问$ values下的值并制作一个函数来绘制它,如下所示:

plotMet = function(obj,metric,var_order){

mat = obj$values
mat = mat[,grep(metric,colnames(mat))]
colnames(mat) = gsub("[~][^ ]*",colnames(mat))
boxplot(mat[,var_order],horizontal=TRUE,las=2,xlab=metric)

}

plotMet(resamps,"MAE",c("rf","glmnet"))

enter image description here

用数字命名模型也不是一个好主意。尝试使用诸如model_2000,model_2001等

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。

相关推荐


依赖报错 idea导入项目后依赖报错,解决方案:https://blog.csdn.net/weixin_42420249/article/details/81191861 依赖版本报错:更换其他版本 无法下载依赖可参考:https://blog.csdn.net/weixin_42628809/a
错误1:代码生成器依赖和mybatis依赖冲突 启动项目时报错如下 2021-12-03 13:33:33.927 ERROR 7228 [ main] o.s.b.d.LoggingFailureAnalysisReporter : *************************** APPL
错误1:gradle项目控制台输出为乱码 # 解决方案:https://blog.csdn.net/weixin_43501566/article/details/112482302 # 在gradle-wrapper.properties 添加以下内容 org.gradle.jvmargs=-Df
错误还原:在查询的过程中,传入的workType为0时,该条件不起作用 &lt;select id=&quot;xxx&quot;&gt; SELECT di.id, di.name, di.work_type, di.updated... &lt;where&gt; &lt;if test=&qu
报错如下,gcc版本太低 ^ server.c:5346:31: 错误:‘struct redisServer’没有名为‘server_cpulist’的成员 redisSetCpuAffinity(server.server_cpulist); ^ server.c: 在函数‘hasActiveC
解决方案1 1、改项目中.idea/workspace.xml配置文件,增加dynamic.classpath参数 2、搜索PropertiesComponent,添加如下 &lt;property name=&quot;dynamic.classpath&quot; value=&quot;tru
删除根组件app.vue中的默认代码后报错:Module Error (from ./node_modules/eslint-loader/index.js): 解决方案:关闭ESlint代码检测,在项目根目录创建vue.config.js,在文件中添加 module.exports = { lin
查看spark默认的python版本 [root@master day27]# pyspark /home/software/spark-2.3.4-bin-hadoop2.7/conf/spark-env.sh: line 2: /usr/local/hadoop/bin/hadoop: No s
使用本地python环境可以成功执行 import pandas as pd import matplotlib.pyplot as plt # 设置字体 plt.rcParams[&#39;font.sans-serif&#39;] = [&#39;SimHei&#39;] # 能正确显示负号 p
错误1:Request method ‘DELETE‘ not supported 错误还原:controller层有一个接口,访问该接口时报错:Request method ‘DELETE‘ not supported 错误原因:没有接收到前端传入的参数,修改为如下 参考 错误2:cannot r
错误1:启动docker镜像时报错:Error response from daemon: driver failed programming external connectivity on endpoint quirky_allen 解决方法:重启docker -&gt; systemctl r
错误1:private field ‘xxx‘ is never assigned 按Altʾnter快捷键,选择第2项 参考:https://blog.csdn.net/shi_hong_fei_hei/article/details/88814070 错误2:启动时报错,不能找到主启动类 #
报错如下,通过源不能下载,最后警告pip需升级版本 Requirement already satisfied: pip in c:\users\ychen\appdata\local\programs\python\python310\lib\site-packages (22.0.4) Coll
错误1:maven打包报错 错误还原:使用maven打包项目时报错如下 [ERROR] Failed to execute goal org.apache.maven.plugins:maven-resources-plugin:3.2.0:resources (default-resources)
错误1:服务调用时报错 服务消费者模块assess通过openFeign调用服务提供者模块hires 如下为服务提供者模块hires的控制层接口 @RestController @RequestMapping(&quot;/hires&quot;) public class FeignControl
错误1:运行项目后报如下错误 解决方案 报错2:Failed to execute goal org.apache.maven.plugins:maven-compiler-plugin:3.8.1:compile (default-compile) on project sb 解决方案:在pom.
参考 错误原因 过滤器或拦截器在生效时,redisTemplate还没有注入 解决方案:在注入容器时就生效 @Component //项目运行时就注入Spring容器 public class RedisBean { @Resource private RedisTemplate&lt;String
使用vite构建项目报错 C:\Users\ychen\work&gt;npm init @vitejs/app @vitejs/create-app is deprecated, use npm init vite instead C:\Users\ychen\AppData\Local\npm-