文献若只如初见-W9

5073d182da4e8a1f53d418528374c50a.png

1) PNAS; 2021 April; anatomical imbalance mapping (AIM); neural development; 

2) BP; 2020 Nov; Cross-disorder; ENIGMA

3) PNAS; 2021 April; MZ-twin studies; adolescent cannabis use

4) Brain Informatics; 2021 May; MRI quality; Freesurfer; paper with code

1

d6fdb97707fa5d266f625908aa060d05.png

PNAS的文章,亮点是提出了一个新的指标anatomical imbalance mapping (AIM),实质就是individual deviations the norm。随后进行了一些包括发展,个体差异在内的分析。

AIM计算方法

5170f314020c6fffc891beb7724c96b8.png

A. cortical thickness,回归掉sex/age/data quality

B. 两个region之间做相关,提取每个点的orthogonal residuals,此处反应的是被试到回归线的垂直距离。

C. 对所有区域进行同样的操作,便得到了k*k的对称矩阵,这里每个被试都能得到一个matrix。

D. 被试内平均得到每个被试的AIM平局值

E. 组间平均得到区域的anatomical imbalance

这个方法初看会以为是一个类似于brain age做残差,但是如果想象一堆身高体重不匀称的月半子在坐标轴的左上方,便会豁然开朗,anatomical imbalance是一个不错的idea。

fab106189467d168f111eeb9af53e138.png

在坐标轴左上还是右下,可能是有意义的,文章只考虑了垂直距离。增加imbalance的方向很简单,因为它和残差的正负一致,但会这有可能大大增加研究的复杂程度。SCN结构协变网络一般做的都是组间的网络,也有一些研究做个体水平的结构网络,AIM指标为个体差异的研究提供了新的思路

2

cd11a1b332107538b1da87ca3f99e9d6.png

其实是一个简单的meta-analysis,把ENIGMA的结果,从paper中整理一下,做了一个EFA。

f37ce0edcf093c432d0bbe862b07d5de.png

对于左右进行了平均,原因是某几个研究只进行了左右平均的分析,把subcortial和cortical的数据放在一起分析,在补充材料中有分开的结果。

1. 做correlation 【连label的颜色都没改!】

fb340321674453e40f72b7454653d825.png

2. 用SPSS做了EFA找overlap的pattern,因子提取方法用是PCA。Rotation用的是Oblique。可以看到development相关的disorder和F1相关较低。

2650ab2c066aad2e6b6d29f5561fe92d.png

3. 使用的几个样本,年龄不一致,为了排除这个可能的混淆因素,作者使用了两个独立的样本HCP和BiDirect,分别计算了age的effect size,并且分别加入2中的模型。结果表明age的pattern和F1不相关,所以找到的F1并不是反映aging的因素【这一点存疑】。

b222bb84a06f1cd3be1389a3cc9bd60a.png

4. 既然找到了共同的latent factor,那么就计算一下每种disorder的deviation吧。其实就是用F1预测原始的effect size,其residual就是deviation。这个值越高说明无法被共同的F1解释。

34d8ae9c166d22f2f520f6e618ff5a5c.png

文章的思路不错,但是有一些局限值得注意:

  • 文章ROI数量为34+7=41,N的数量为9+1。思路和结果看起来不错,但是这相当于用EFA分析了一个有41个items的问卷,而被试数量只是10。EFA推荐每个item至少需要10个被试,那么要410个sample结果才可靠。

  • 文章把subcortial和cortial的结果同时放在一个模型中分析是有问题的,因为之前的研究对于这两种指标的分析有可能用的不是同样的模型。

  • 把之前meta和mega的结果放在一个模型中有可能是问题的。

  • 因为某个文章只做了左右平均的分析,将剩下研究的effect size按照左右进行平均放入模型也是有问题的。

3

141134a0ef0fc01f7b2f6f503fd6cd43.png

随着recreational cannabis use逐步合法化,它对正处于生发发育期间的青少年的影响成为了研究的热点。这篇发在PNAS上的文章,使用了双生子+longitudinal的设计探究了cannabis use和mental health,cognitive ability, academic functioning和socioeconomic status的相关。

前言简单回顾了之前该方向的研究,

  • observational studies

  • neuroimaging studies

  • animal studies

但是这些研究的不足是

  • same sample or lack of replication sample

  • cross-sectional design

  •  retrospective studies

文章最大的亮点就是使用双生子的样本探究因果关系。首先他们在普通的样本中做了相关,发现各种显著,这也是大部分研究发现的,大麻使用问题太多,立法应该需要关注这些问题。

25df68d3bab80456e90d516e166f7998.png

在只用同卵双生子的样本中,只是有socioeconomic中出现了显著

0e0230120f785a301cacc6b72f914658.png

结果表明,大麻使用并没有和青少年和成年人的心理健康和认知能力的改变有明显的因果关系,但是它对于学业成绩和社会经济状态有影响。【这个研究不表明它不会成瘾,也不表明它不会对发育中的大脑造成伤害,文章中的Cognitive都是用的WAIS测的。】

因果关系的研究思路

JAMA Psychiatry: Genetically Informed Methods

4

7415ee4f343ebe642ccbb883ddf06fab.png

首先使用CAT12计算了图像的质量分数,随后人工做QC的结果做对比,表明CAT12是一个有效的指标。CAT12的优点是计算简单,使用个人笔记本即可完成。

db20f2918734f7cba3785c09e1a159f4.png

075dfe8c335ba3305f09fc7312da392a.png

之后用CAT12图像质量分数和Freesurfer常用的指标做相关,结果发现即使人工做过visual QC,图像的质量分数依然和ROI相关。

注意⚠

这里(👇图)做的只是简单的相关,并没有控制其他的变量。而文章之前明确表明图像质量和年龄相关,年龄小的被试图像质量相对较差。但是在之后的分析中都是用的两两相关的方式,并没有控制诸如年龄、性别等因素。那么这一堆pattern会不会就是age的效应呢。

Surface Area

d141caf42b22de189745a716f35ca54f.png

Cortical Thickness

d97612a31ee2b3fabb23f3790c346e41.png

Subcortical volume

3871a60f6b52e631fa935b4061f182d5.png

上面这一批图使用ggseg做的:

R| ggseg 绘制统计结果

—END—

7d4fd6e5646555753b5865de9a8e7405.png

数据预处理、建模和可视化是数据分析过程中非常重要的环节,这些环节的成功与否直接影响到数据分析的结果与价值。作为CSDN数据分析的三大关键步骤,数据预处理、建模和可视化的作用是多方面的。 首先,数据预处理的目的是对原始数据进行清洗、处理和转换,以便更好地进行后续数据分析。数据预处理有很多种方法,如数据清洗、缺失值填充、异常值处理、数据集成、数据降维等。只有在预处理阶段对数据进行充分准确的处理,才能保证后续的建模和可视化的准确性。 其次,建模是数据分析的核心步骤,是通过数据分析来建立模型的过程,目标是对数据进行预测和决策。建模过程包括选择建模方法,配置模型参数,模型训练和模型评估等步骤。建模需要根据实际问题选择合适的建模方法,并进行适量的特征工程和模型调参,以获得最优的预测和决策效果。 最后,可视化是数据分析结果的展示方式,旨在将结果呈现给用户,在数据探索和决策时加强用户的感性理解。可视化可以采用多种图表形式,例如柱状图、折线图、散点图等,也可以通过交互式可视化系统进行交互操作表示。通过可视化,数据分析器可以更直观地理解数据分析的结果,为后续的分析和决策提供有力支持。 综上所述,CSDN数据分析中的数据预处理、建模和可视化是相互依存、相辅相成的环节,只有将它们表现得尽善尽美,我们才能解决实际问题,提高数据因信而产生的价值。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值