弹性网惩罚项的可视化分析(r语言)(代码片段)

统计学小王子 统计学小王子     2023-01-10     226

关键词:

0引言

在博客《变量选择——lasso、SCAD、MCP的实现(R语言》中介绍了R语言做变量选择常用的函数和包,其中msaenet和glmnet是做的弹性网,弹性网是岭回归和lasso类惩罚的回归。下面以二维为例,介绍岭回归、lasso、弹性网惩罚项的区别,因为我们知道惩罚回归是似然函数或者最小二乘求条件极值,具体的形式以限制条件给出,给定常数C有:
岭回归,
x 2 + y 2 = C , x^2 + y^2=C, x2+y2=C,
lasso:
∣ x ∣ + ∣ y ∣ = C , |x|+|y|=C, x+y=C,
弹性网,给定0-1的 α \\alpha α有:
α ( ∣ x ∣ + ∣ y ∣ ) + ( 1 − α ) ( x 2 + y 2 ) = C . \\alpha (|x|+|y|)+(1-\\alpha)(x^2+y^2)=C. α(x+y)+(1α)(x2+y2)=C.
可以看出当 α = 1 \\alpha=1 α=1时退化为lasso,当 α = 0 \\alpha=0 α=0时退化为岭回归。
下面取C=1,对上述三个可行域可视化。

1、函数构建

1.1代码

# 2.2函数调用
O <- par(mfrow = c(2, 3))
netplot(C = 1, alpha = 0)
text(0,0.8,paste0("alpha = ",0),col = "blue", cex = 4)
netplot(C = 1, alpha = 0.1)
text(0,0.8,paste0("alpha = ",0.1),col = "blue", cex = 4)
netplot(C = 1, alpha = 0.2)
text(0,0.8,paste0("alpha = ",0.2),col = "blue", cex = 4)
netplot(C = 1, alpha = 0.5)
text(0,0.8,paste0("alpha = ",0.5),col = "blue", cex = 4)
netplot(C = 1, alpha = 0.8)
text(0,0.8,paste0("alpha = ",0.8),col = "blue", cex = 4)
netplot(C = 1, alpha = 1)
text(0,0.8,paste0("alpha = ",1),col = "blue", cex = 4)
par(O)

1.2图片

1.3动图

2、总结

最后希望可以帮助大家学习R语言。水平有限发现错误还望及时评论区指正,您的意见和批评是我不断前进的动力。当然如果有更好的改进方案欢迎评论区交流。
本文的代码已整理好,实在刚需移步下面链接下载研究。

高维数据惩罚回归方法:主成分回归pcr岭回归lasso弹性网络elasticnet分析基因数据|附代码数据(代码片段)

...进行主成分回归(PCR)。使用glmnet()进行岭回归、lasso和弹性网elasticnet对这些预测模型进行评估1.1 数据集在本文中,我们将使用基因表达数据。这个数据集包含120个样本的200个基因的基因表达数据。这些数据来源于哺乳动物... 查看详情

r语言时间序列数据指数平滑法分析交互式动态可视化(代码片段)

 ? 原文链接:http://tecdat.cn/?p=13971R语言提供了丰富的功能,可用于绘制R中的时间序列数据。包括:自动绘制 xts 时间序列对象(或任何可转换为xts的对象)的图。高度可配置的轴和系列显示(包括可选的第二个Y轴... 查看详情

数据可视化应用数据统计分析的显著性标注(附python和r语言代码)(代码片段)

Python版本Python-Seaborn自定义函数绘制我们可以通过自定义绘图函数的方式在统计图表中添加显著性标注,这里我们直接使用Seaborn自带的iris数据集进行绘制,具体内容如下:自定义P值和星号对应关系由于是完全的自定义,这里需... 查看详情

《r语言实战》阅读笔记--r语言介绍(代码片段)

...专门的数据仓库R是免费的R是一个全面的统计研究平台R的可视化牛……1.2R的获取和安装R可以在CRAN上(htt 查看详情

基于弹性堆栈(elk堆栈)的日志分析存储及展示(代码片段)

...“存储”。Kibana允许用户使用Elasticsearch中的图表和图形可视化数据。ELK原理在所有需要收集日 查看详情

数据可视化应用绘制森林图(附python和r语言代码)(代码片段)

...表示出来,是一种简单直观地展示单一研究和汇总研究的可视化图表。下面,小编再附上森林图各个绘图元素间关系的解释图:KeyElementsoftheForestplot下面就介绍下如何使用R和Python绘制森林图(For 查看详情

r语言实战应用精讲50篇(三十一)-r语言入门系列-tidyverse数据分析流程(代码片段)

...verse?tidyverse是一个清晰的R包集合,在数据操作、探索和可视化领域提供统一的数据科学解决方案,这些解决方案具有共同的设计理念。它是由RStudio背后的首席科学家HadleyWickham创建的。tidyverse中的R包旨在提高统计学家和数据科... 查看详情

r获取股票数据并进行进行可视化分析(代码片段)

R获取股票数据并进行进行可视化分析#加载依赖的包library(quantmod)library(ggplot2)library(magrittr)library(broom)#设置计划分析或者处理的时间范围start=as.Date("2020-02-01")end=as.Date("2020-03-31")#通过股票代码获取时间范围内的数据,可以设置多个... 查看详情

r语言实战应用精讲50篇(二十八)-r语言时空数据分析实战案例-数据处理及可视化(代码片段)

1.数据整理时空建模和预测通常涉及大量数据,这些数据以各种形式提供给用户,但通常以CSV文件或文本文件中的表格形式提供。通常会花费大量时间来加载数据并对其进行预处理,以便将它们转化为适合分析的形式。R中有几个... 查看详情

r语言实战应用精讲50篇(二十八)-r语言时空数据分析实战案例-数据处理及可视化(代码片段)

1.数据整理时空建模和预测通常涉及大量数据,这些数据以各种形式提供给用户,但通常以CSV文件或文本文件中的表格形式提供。通常会花费大量时间来加载数据并对其进行预处理,以便将它们转化为适合分析的形式。R中有几个... 查看详情

r语言中多分类问题multicalssclassification的性能测量(代码片段)

...。本文主要关注LDA,并探讨其在理论和实践中作为分类和可视化技术的用途。由于QDA和RDA是相关技术,我不久将描述它们的主要属性以及如何在R中使用它们。线性判别分析LDA是一种分类和降 查看详情

r语言无监督学习:pca主成分分析可视化(代码片段)

...里,我们重点介绍两种技术…主成分分析:用于数据可视化或在其他监督学习方法之前进行预处理的工具。聚类:发现数据中未知组的方法。无监督学习的挑战通常,无监督学习比主观学习更具挑战性,因为它更具主观性。... 查看详情

r语言r语言数据可视化——东北大学大数据班r实训第三次作业(代码片段)

数据可视化知识点回顾基础知识点回顾:条形图、饼图、直方图、核密度图、箱线图和点图。中级知识点回顾:散点图、气泡图、折线图、相关图、马赛克图。title:“R实训第三次作业”output:html_notebook1.通过读取文件deathr... 查看详情

r语言实战应用精讲50篇(二十七)-时空数据分析-经验空间/时间均值(latex公式+r代码绘图)(代码片段)

前言数据可视化无疑是探索性数据分析的一个重要且必要的组成部分。此外,我们经常希望根据一阶(first-order)和二阶(second-order)特征的总结来探索时空数据。在这里,我们考虑经验均值(empiricalmeans)和经验协方差(empiricalcovariances)... 查看详情

r语言实战应用精讲50篇(二十七)-时空数据分析-经验空间/时间均值(latex公式+r代码绘图)(代码片段)

前言数据可视化无疑是探索性数据分析的一个重要且必要的组成部分。此外,我们经常希望根据一阶(first-order)和二阶(second-order)特征的总结来探索时空数据。在这里,我们考虑经验均值(empiricalmeans)和经验协方差(empiricalcovariances)... 查看详情

r语言方差分析(代码片段)

...d=“mcd”)#稳健单因素MANOVAvegan::adonis()#非参数MANOVA初级可视化(均值、方差、分布)group_by()%>%plot()boxpl 查看详情

r语言独立成分分析fastica谱聚类支持向量回归svr模型预测商店销量时间序列可视化(代码片段)

全文链接:http://tecdat.cn/?p=31948原文出处:拓端数据部落公众号本文利用R语言的独立成分分析(ICA)、谱聚类(CS)和支持向量回归SVR模型帮助客户对商店销量进行预测。首先,分别对商店销量的历史数据进行了独立成分分析,... 查看详情

2022最新r语言与rstudio安装配置教程(代码片段)

...装R软件是一个自由、开源软件平台,具有统计分析、可视化和科学计算的强大功能。R语言官网:https://cran.r-project.org/,进入R语言官网,DownloadeRforWindows。可以看到,下载下来,目前最新版本为 查看详情