对称半正定矩阵可进行choleskey分解，使用chol()函数结合tryCatch错误异常判断，即可判断矩阵是否对称半正定。

1. 仅输出是否半正定

Sigma为一个对称矩阵，但非半正定，进行choleskey分解后报error

> chol(Sigma)
Error in chol.default(Sigma) : the leading minor of order 14 is not positive definite

结合tryCatch判断矩阵半正定

tryCatch({chol(Sigma)},error = function(e){print("Sigma非对称半正定")})

输出结果

[1] "Sigma非对称半正定"

若Sigma为对称非半正定，则输出chol(Sigma)

2. 将判断结果保存在变量中用于后续处理

flag<-TRUE
tryCatch({chol(Sigma)},error = function(e){flag<<-FALSE})

Sigma半正定，则flag为TRUE；Sigma非半正定，则flag为FALSE
注意：tryCatch中赋值需用<<-，才能改变父级中的变量

参考

R语言错误异常处理tryCatch(变量作用域)

R语言判断矩阵是否对称半正定相关推荐

R语言ggplot2可视化：水平半小提琴图（Horizontal Half Violin Plots）、去除水平半小提琴图中的填充色、ggridges包的绘制山脊图
R语言ggplot2可视化:水平半小提琴图(Horizontal Half Violin Plots).去除水平半小提琴图中的填充色.ggridges包的geom_density_ridges函数绘制 ...
R语言如何计算矩阵的列均值？
R语言如何计算矩阵的列均值? 目录 R语言如何计算矩阵的列均值? R语言是解决什么问题的? R语言如何计算矩阵的列均值? R语言是解决什么问题的? R 是一个有着统计分析功能及强大作图功能的软件系统, ...
笔记—R语言做矩阵散点图
笔记-R语言做矩阵散点图 #install.packages("corrgram") #install.packages("GGally") install.p ...
r语言中正定矩阵由于误差不正定_浅谈「正定矩阵」和「半正定矩阵」
在众多的机器学习模型中,线性代数的身影无处不在,当然,我们也会时常碰到线性代数中的正定矩阵和半正定矩阵.例如,多元正态分布的协方差矩阵要求是半正定的. --------------×--------- ...
r语言中正定矩阵由于误差不正定_R语言之数据处理(一)
在上一篇小文中,提到了关于R语言导入数据的一些方法,之后的重点就转向了数据的处理上.数据处理其实在整个数据分析项目中所占用的时间是比较多的,所以根据处理的目的不同,也有不同的处理方法.在R语言中,我通 ...
r语言中矩阵QR分解_从零开始学R语言Day4|向量、矩阵和数组
从零开始学R语言Day4|向量.矩阵和数组 1.1向量 1.1.1向量在Day2中我们提及过用和c()函数来构建向量,具体实例如下. 我们还可以采用vector("类型",长度) ...
r语言中矩阵QR分解_R语言常用的矩阵操作
R语言是一门非常方便的数据分析语言,它内置了许多处理矩阵的方法.下面列出一些常用的矩阵操作方法示例. 矩阵的生成 > mat <- matrix(1:16, ncol = 4, nrow ...
R语言使用矩阵操作回归分析
一朋友问我说: 飞哥,你知道回归分析中利用的是最小二乘法,比如最简单的单变量回归分析,得到的有回归系数和截距,但是相关的标准误是如何计算的??? 我:--竟然讲不出来内心小99 作为杠精我是不服气的 ...
C语言——判断矩阵维数（sizeof、_msize）
C语言中可以利用二维数组或者malloc开辟矩阵,本章将分别介绍这两种方法是如何判断矩阵维数的. 判断矩阵维数 sizeof 以二维数组创建矩阵时,可以利用sizeof判断内存的大小例如: doub ...

R语言判断矩阵是否对称半正定

1. 仅输出是否半正定

2. 将判断结果保存在变量中用于后续处理

参考

R语言判断矩阵是否对称半正定相关推荐

最新文章

热门文章

R语言 判断矩阵是否对称半正定

1. 仅输出是否半正定

2. 将判断结果保存在变量中用于后续处理

参考

R语言 判断矩阵是否对称半正定相关推荐

最新文章

热门文章

R语言判断矩阵是否对称半正定

R语言判断矩阵是否对称半正定相关推荐