关于我花了几天时间又复现了一篇 SCI 一区 Top 的代码
【MATLAB】赫尔默特方差分量估计算法
「微信公众号由于改变了推送规则，为了每次新的推送可以在第一时间出现在您的订阅列表中，记得将本公众号设为星标或置顶哦~」
有意向获取代码，请转文末观看代码获取方式~
1 文章简介
文章 DOI: 10.1109/TGRS.2023.3265508
链接：https://ieeexplore.ieee.org/document/10097458/keywords#keywords
该篇文章讲的是一种遥感反演雪深的新方法，在对每日诸多雪深反演值进行取值的过程中，应用了赫尔默特方差分量估计算法对 GPS、Galileo、GLONASS 和 BDS 卫星系统每日反演的雪深分别进行定权，并通过权重求出多系统融合的每日雪深反演结果，是一种不同于以往的雪深求取算法。虽然相关作者没有开源赫尔默特方差分量估计算法的代码，但是其理论讲解的很清晰完善，完全可以做到复现。
2、什么是赫尔默特方差分量估计算法呢？
赫尔默特方差分量估计算法（Herzberg variance component estimation）是一种用于估计数据集中各个因素对总变异贡献的方法。它可以帮助研究人员确定数据集中不同因素的相对重要性，从而更好地理解数据。该算法首先将数据集分成若干组，然后计算每组数据的方差。接下来，通过比较每组数据之间的方差大小，可以确定哪些因素对总变异贡献最大。这些因素被称为方差分量。 具体而言，赫尔默特方差分量估计算法可以分为以下步骤：
- 首先，将数据集分成若干组。
- 对于每组数据，计算其方差。
- 计算数据集总体方差。
- 计算每个因素的方差分量，即根据每组数据的方差和组内样本数量，使用方差分析方法计算出每个因素的方差分量。
- 根据方差分量的大小排序，确定哪些因素对总变异贡献最大。
需要注意的是，赫尔默特方差分量估计算法假定每个因素的方差分量是固定的，因此在实际应用中可能会存在一定的误差。此外，该算法在数据集比较小的情况下可能会面临“稀疏性”问题，即每个因素的样本数量非常少，从而导致方差分量估计不准确。因此，在使用该算法时需要根据具体应用场景进行评估和选择。
3、其应用的学科领域有哪些呢？
赫尔默特方差分量估计算法（Herzberg variance component estimation）是一种用于估计数据集中各个因素对总变异贡献的方法。它可以帮助研究人员确定数据集中不同因素的相对重要性，从而更好地理解数据。
以下是赫尔默特方差分量估计算法的一些应用：
- 生物医学研究：在生物医学领域，赫尔默特方差分量估计算法常常用于分析基因或某种治疗方式对疾病进展的影响。通过比较不同基因或治疗方式之间的方差分量，可以确定其中哪些因素对治疗效果或疾病进展具有重要作用。
- 社会科学研究：在社会科学领域，赫尔默特方差分量估计算法常常用于分析调查问卷数据。通过比较不同问题或变量之间的方差分量，可以确定哪些问题或变量对受访者行为或态度具有显著影响。
- 工程和物理学研究：在工程和物理学领域，赫尔默特方差分量估计算法常常用于分析实验数据。通过比较不同因素（如温度、压力等）之间的方差分量，可以确定哪些因素对实验结果具有显著影响。
总之，赫尔默特方差分量估计算法适用于各种领域的数据分析。它可以帮助研究人员更好地理解数据，并找出其中的规律和趋势。
通常情况下，各个学科都是适用的。比如 20 个测量值和 16个测量值，分别是由两种不同测量方式得出的测量结果，为了得到两种测量方式组合后的测量值，需要对其进行定权（说白了就是定两者的比例），然后再对测量结果合并求出最佳测量值。
应用途径及注意事项：两组不同测量方式测雪深的 20 个和 16 个测量值，用赫尔默特方差分量估计算法求解给定两组雪深测量方式的权重（通常为 1：xx，比如：1: 0.92）。赫尔默特方差分量估计算法通常情况下是确定两组不同测量方式的权重，若要测量三种测量方式的权重，则需要两两定权（比如方式 1 和方式 2 的权重为 1: 0.8，方式 1 和方式 3 的权重为 1: 0.9，则方式 1: 方式 2: 方式 3 的权重为 1: 0.8: 0.9）。
4、复现效果
原作者用的雪深值，那这里代码我也自己制作了几个雪深值用于数据处理和分析。主要是用赫尔默特方差分量估计算法进行复现的。
假设两组测雪深的数据 snow 为
第一组
10.6 cm、9.8 cm、9.4 cm、9.5 cm、9.8 cm
10.7 cm、10.9 cm、9.4 cm、9.9 cm、10.8 cm
10.9 cm、10.2 cm、10.4 cm、10.5 cm、9.5 cm
10.3 cm、10.1 cm、9.6 cm、9.2 cm、10.6 cm
第二组
11.9 cm、9.2 cm、10.4 cm、8.9 cm、9.9 cm
10.2 cm、10.1 cm、11.4 cm、9.9 cm、11.8 cm
通过下方命令行输出内容可以发现：
赫尔默特方差分量估计定权的结果 sita_res 中给出两组数据的权重之比为 1: 0.3047
第一组雪深数据的中值为 10.15；第二组雪深数据的中值为 10.15
第一组雪深数据的平均值为 10.105；第二组雪深数据的平均值为 10.37
由两组雪深数据的中值加权平均结果为 10.15
由两组雪深数据的平均值加权平均结果为10.1669
注：赫尔默特方差分量估计定权的目的是根据两组数据确定这两组测量方式的权重，然后需要根据每组数据的中值或者平均值确定赫尔默特方差分量估计后的雪深值。
5 附视频教程
https://www.bilibili.com/video/BV1GX4y1B7Y9/
6 代码获取
MATLAB 赫尔默特方差分量估计算法的开源代码请转:
https://mbd.pub/o/bread/ZJiUlp5x
「知识是无价的，但个人的时间成本是有价的，万望海涵~」
1、感谢关注 Lwcah 的个人公众号，有关资源获取，请公众号后台发送推文末的关键词，自助获取。
2、若要添加个人微信号，请后台发送关键词：微信号。
3、若要进微信群：Lwcah 科研技巧群 3。请添加个人微信号后进群（大家沉浸式科研，广告勿扰），不定时更新科研技巧类推文。可以一起探讨科研，写作，文献，代码等诸多学术问题，我们一起进步。
记得关注公众号，并设为星标哦~谢谢啦~
「万请尊重原创成果！！！」
「声明：本公众号（Lwcah）的原创成果，在未经允许的情况下，请勿用于任何商业用途！」
「And，今后我会尽可能出一些更高质量的推文与大家共享，再一次感谢大家的关注与支持~也特别感谢大家对公众号的传播与分享，每天新增的关注都是我持续更新的动力！」
「您的每一次点赞，在看，关注和分享都是对我最大的鼓励~谢谢~」
