使用abagen包提取AHBA基因表达数据

介绍使用Python环境下abagen包提取AHBA(Allen Human Brain Atlas)基因表达数据的基本方法。

《使用abagen包提取AHBA基因表达数据》有6个想法

  1. 还想请教一下,如果我想在abagen这里用xcp_d的内置分区保证我后面的分析匹配,比如想用schaefer400图谱,我应该怎么做呢?我尝试下载了图谱,但是一般给的分区模板的信息都和文档里说的不一样

    1. 我不太清楚你具体遇到了什么问题。因为我没有实际看过xcp_d内置的分区文件是什么情况(我xcp_d用的也不多),我印象中好像是fs_LR空间的。如果xcp_d内置的分区和abagen要求的不一样,就需要进行相应的转换。比如将分区从fs_LR空间转换到fsaverage。如果是我的话,我会直接考虑用Yeo实验室提供的fsaverage空间的分区文件,而不是xcp_d内置的,只要脑区能对应上就行。

      1. 好嘞谢谢!是我在xcp_d处理的时候没有选择别的分区模板,使用了它自己内置的模板,包括了这些4S1056Parcels’, ‘4S156Parcels’, ‘4S256Parcels’, ‘4S356Parcels’, ‘4S456Parcels’, ‘4S556Parcels’, ‘4S656Parcels’, ‘4S756Parcels’, ‘4S856Parcels’, ‘4S956Parcels’, ‘Glasser’, ‘Gordon’, ‘HCP’, ‘MIDB’, ‘MyersLabonte’, ‘Tian,所以在选择abagen的时候想要跟原本的对应,处理起来就比较麻烦

        1. xcp_d内置的分区有些是合并了多个分区模板构成的,比如4S1056Parcels,但是abagen可能没办法直接用这些分区模板。比如要得到4S1056Parcels的基因表达矩阵,我个人觉得比较简单的做法可能是,先找出组成这个分区的子分区模板有哪些,比如有的子分区是MNI空间的,有的是fsaverage空间,然后分别放入abagen进行分析,最后再把不同分区的结果合并在一起。

  2. 或者是在xcp_d处理的时候外接HCPex的图谱,但是xcp_d外接图谱得用BIDS格式,也挺麻烦,找不到做好的,还得自己整

    1. 可以根据XCP_D处理后的dense数据去提取其他分区图谱的时间序列和计算相关矩阵,这些步骤并没有必要一定要用XCP_D。

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注