基于CAT12的report数据做QC质量检查，排除离群点—

前言

这段时间被质量检查QC折磨的很惨，各种检查，要检查提取的被试数据质量是否有问题，病人组和正常组中是否出现个别质量差别较大，还有两组被试的质量方面是否有差异，会不会影响我们之后的数据分析处理。这都是需要考虑的，还要检查自己的代码提取的数据是否跟cat12计算的是不是一样的，这样便可以不需要cat12的GUI界面也能提取所需的指标。脱离GUI还有个好处是可用服务器直接调用，便于操作。由于看到一篇文章，在做结构网络分析时，采用欧拉数来排除可能会影响结果的被试。根据这个，下面写了一个小函数来提取IQR和欧拉数，来排除一些离群点的被试，综合考虑排除被试。仅供参考！

指标

指标提取的基于Matlab相关代码（仅供参考）

注：代码仅供参考，有的被改动过，目的希望读者能自己学会，而不是直接copy，当然这对于代码大神们就小菜一碟喽，也不需要小编这样的笔记。所以此篇写出只是给自己存档，看自己的笔记也可以重新捋一捋思路。

catNC_folder = dir([catsub_path, '\cat_NC*.mat']); for NCsi = 1:length(catNC_folder)cat_name = catNC_folder(NCsi).name;load(['',catsub_path,'\',cat_name,'']);for li=1:length(S.catlog)if ~isempty(strfind(S.catlog{li,1},'IQR'))==1IQR_NC(NCsi,1)=str2num(S.catlog{li,1}(isstrprop(S.catlog{li,1},'digit')))*0.0001;endend  Eulnum_NC(NCsi,1)=S.qualitymeasures.SurfaceEulerNumber;endcatEP_folder = dir([catsub_path, '\cat_*.mat']);for EPsi = 1:length(catEP_folder)cat_name = catEP_folder(EPsi).name;load(['' catsub_path '\' cat_name '']);for li=1:length(S.catlog)if ~isempty(strfind(S.catlog{li,1},'IQR'))==1IQR_EP(EPsi,1)=str2num(S.catlog{li,1}(isstrprop(S.catlog{li,1},'digit')))*0.0001;endend  Eulnum_EP(EPsi,1)=S.qualitymeasures.SurfaceEulerNumber;endcat_folder = dir([catsub_path, '\cat_*.mat']); for si = 1:length(cat_folder)cat_name = cat_folder(si).name;load(['' catsub_path '\' cat_name '']);for li=1:length(S.catlog)if ~isempty(strfind(S.catlog{li,1},'IQR'))==1IQR_*(si,1)=str2num(S.catlog{li,1}(isstrprop(S.catlog{li,1},'digit')))*0.0001;endend  Eulnum_*(si,1)=S.qualitymeasures.SurfaceEulerNumber;end[h1,p1,ci1,stats1] = ttest2(IQR_NC(1:25,1),IQR_NC);[h2,p2,ci2,stats2] = ttest2(IQR_NC(25:99,1),IQR);[h3,p3,ci3,stats3] = ttest2(Eulnum_NC(1:25,1),Eulnum_*);[h4,p4,ci4,stats4] = ttest2(Eulnum_NC(25:99,1),Eulnum_*);figure('color','w')subplot(2,2,1);g = [repmat(1,length(IQR_NC(1:25,1)),1);repmat(2,length(IQR_EP),1)];IQR1 = [IQR_NC(1:25,1);IQR_EP];pos = ['NC';'EP'];boxplot(IQR1,g, 'labels', pos)subplot(2,2,2);g = [repmat(1,length(IQR_NC(26:99,1)),1);repmat(2,length(IQR_*),1)];IQR2 = [IQR_NC(26:99,1);IQR_*];boxplot(IQR2,g, 'labels', pos)subplot(2,2,3);g = [repmat(1,length(Eulnum_NC(1:25,1)),1);repmat(2,length(Eulnum_EP),1)];Eulnum1 = [Eulnum_NC(1:25,1);Eulnum_EP];boxplot(Eulnum1,g, 'labels', pos)subplot(2,2,4);g = [repmat(1,length(Eulnum_NC(26:99,1)),1);repmat(2,length(Eulnum_*),1)];Eulnum2 = [Eulnum_NC(26:99,1);Eulnum_*];boxplot(Eulnum2,g, 'labels', pos)

结果

下图为IQR的提取结果：

下图为欧拉数的提取结果：

最后的问题就是要不要去除这些离群点呢？待之后的数据分析结果考虑。
很久没写笔记了，写了一篇可能也没啥用的笔记，哈哈哈，刷刷存在感吗？下篇认真写吧，应该是写FDR和FWE的操作吧！给自己定的计划和小目标，期待一下。

基于CAT12的report数据做QC质量检查，排除离群点——IQR欧拉数相关推荐

基于WiFi的CSI数据做呼吸频率检测-python版(含代码和数据)
一.概述本Demo无需机器学习模型,Demo功能涉及的理论主要参考了硕士学位论文<基于WiFi的人体行为感知技术研究>,作者是南京邮电大学的朱XX,本人用python复现了论文中呼吸频率 ...
用fastp对转录组数据做QC
链接 fastp: 极速全能的FASTQ文件自动质控+过滤+校正+预处理软件 github地址看到介绍的时候是真的心动不已↓↓↓ fastp可以仅仅扫描 FASTQ 文件一次,就完成比FASTQC ...
金融领域如何基于大数据做风控
金融主要分为资产.负债和非息业务.现在的很多互联网金融产品也都是基于这些业务衍生出来的,例如三方支付.投资理财.P2P.消费信贷.分期.白条等等,而在这其中P2P.消费信贷.分期.白条等都属于信贷类业 ...
基于Python实现的数据质量检查
目录 1:应用场景 2:外部数据数据质量评估解决方案构思一: 2.1:评估维度--"三率" 2.2:评估维度--"三性" 2.3:评估维度--"三度 ...
SARScape中用sentinel-1数据做SBAS-InSAR完整流程（1/2）
SARScape中用sentinel-1数据做SBAS-InSAR完整流程 1 SABA-InSAR原理简述 2 数据采集和预设 2.1 SAR数据采集 2.2 DEM数据下载与放置 2.3 精密轨道 ...
RM: 基于页面结构化数据生成报表,一键导出图片,生成定制图表文末有效果图 , 开放部分代码
背景开发这个工具是因为一句抱怨故事是这样的,我们公司是一个非常重视员工健康的公司,一年前老董说让HR(后面改为ZT)督促员工多多运动,可持续地位公司创造价值.并拿出了一部预算来奖励那些积极运动的人 ...
一种基于三代PacBio测序数据的补洞方法
一种基于三代PacBio测序数据的补洞方法技术领域本发明涉及生物信息技术领域,具体涉及DNA组装的补洞方法,它使用三代PacBio 测序数据来进行基因组数据的补洞. 背景技术三代PacBio测序 ...
基于Grafana的监控数据钻取功能应用实践
互联网企业中,随着机器规模以及业务量的爆发式增长,监控数据逐渐成为一种大数据,对监控大数据的分析,包括数据采集.数据缓存.数据聚合分析.数据存储.数据展现等几个阶段.不同阶段有不同的解决方案及支撑工具 ...
一步一步教你使用AgileEAS.NET基础类库进行应用开发-基础篇-基于接口驱动的数据层...
系列回顾在前面的文章中,我用了大量的篇幅对UDA及ORM的使用进行了讲解和演示,我们已经知道并熟悉的使用UDA和ORM构建简单的应用,AgileEAS.NET在应用的纵向结构上建议使用分层结构,提出 ...
如何利用大数据做金融风控？原创 2016年11月24日 17:42:03 标签：大数据 / 金融 / 风控 1594 导语：如何通过海量数据与欺诈风险进行博弈？随着金融科技、科技金融等概念的
如何利用大数据做金融风控? 原创 2016年11月24日 17:42:03 标签: 大数据 / 金融 / 风控 1594 导语:如何通过海量数据与欺诈风险进行博弈? 随着金融科技.科技金融等概念的热起 ...

基于CAT12的report数据做QC质量检查，排除离群点——IQR欧拉数

前言

指标

指标提取的基于Matlab相关代码（仅供参考）

结果

基于CAT12的report数据做QC质量检查，排除离群点——IQR欧拉数相关推荐

最新文章

热门文章