计算化学公社

 找回密码 Forget password
 注册 Register
Views: 14897|回复 Reply: 10
打印 Print 上一主题 Last thread 下一主题 Next thread

[Amber] 利用cpptraj进行主成分分析,生成每一帧的主成分三维散点图求助

[复制链接 Copy URL]

71

帖子

0

威望

694

eV
积分
765

Level 4 (黑子)

参考如下方法进行主成分分析,有两个问题想请教

1、想利用数据做出生成类似下图的每一帧的主成分三维散点图,请教该怎么做?

2、做出的前三个主成分直方图归一化图有什么意义?主成分的轨迹动画有什么意义?

谢谢!

6万

帖子

99

威望

6万

eV
积分
125150

管理员

公社社长

2#
发表于 Post on 2021-6-16 18:09:29 | 只看该作者 Only view this author
仔细看此文了解基本原理
浅谈PCA与g_covar+g_anaeig+ddtdp+sigmaplot做自由能面图的方法
http://sobereva.com/73

得到每一帧在PC1、PC2、PC3投影的数值,放到Origin里就能绘制三维散点图

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

71

帖子

0

威望

694

eV
积分
765

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2021-6-16 22:05:14 | 只看该作者 Only view this author
sobereva 发表于 2021-6-16 18:09
仔细看此文了解基本原理
浅谈PCA与g_covar+g_anaeig+ddtdp+sigmaplot做自由能面图的方法
http://sobereva ...

多谢sob老师,第一个问题解决了。

71

帖子

0

威望

694

eV
积分
765

Level 4 (黑子)

4#
 楼主 Author| 发表于 Post on 2021-6-16 22:28:03 | 只看该作者 Only view this author
sobereva 发表于 2021-6-16 18:09
仔细看此文了解基本原理
浅谈PCA与g_covar+g_anaeig+ddtdp+sigmaplot做自由能面图的方法
http://sobereva ...

还有一个问题想请教sob老师,我现在是把散点图画出来了,怎么在图中将不同的簇用不同的颜色表示那?

6万

帖子

99

威望

6万

eV
积分
125150

管理员

公社社长

5#
发表于 Post on 2021-6-17 11:49:52 | 只看该作者 Only view this author
landian666 发表于 2021-6-16 22:28
还有一个问题想请教sob老师,我现在是把散点图画出来了,怎么在图中将不同的簇用不同的颜色表示那?

把不同的簇的散点分别绘制,用不同的颜色就完了
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

71

帖子

0

威望

694

eV
积分
765

Level 4 (黑子)

6#
 楼主 Author| 发表于 Post on 2021-6-17 12:56:02 | 只看该作者 Only view this author
sobereva 发表于 2021-6-17 11:49
把不同的簇的散点分别绘制,用不同的颜色就完了

我是用amber的CPPTRAJ做的簇分析,看输出文件里只有噪音帧的序号,想请教下sob老师怎么确认每个团簇的帧序号?

6万

帖子

99

威望

6万

eV
积分
125150

管理员

公社社长

7#
发表于 Post on 2021-6-18 18:14:51 | 只看该作者 Only view this author
landian666 发表于 2021-6-17 12:56
我是用amber的CPPTRAJ做的簇分析,看输出文件里只有噪音帧的序号,想请教下sob老师怎么确认每个团簇的帧序 ...

很久不用cpptraj,忘了
gromacs的gmx cluster都会给出每个簇里对应的帧,cpptraj应该也有办法给出
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

71

帖子

0

威望

694

eV
积分
765

Level 4 (黑子)

8#
 楼主 Author| 发表于 Post on 2021-6-19 16:28:56 | 只看该作者 Only view this author
sobereva 发表于 2021-6-18 18:14
很久不用cpptraj,忘了
gromacs的gmx cluster都会给出每个簇里对应的帧,cpptraj应该也有办法给出

谢谢sob老师,问题解决了,还有两个问题想请教sob老师:
1  就是我在研究一个动力学过程中先用GROMACS做动力学计算,然后用gmx_mpi covar 和gmx_mpi anaeig做PCA,然后把轨迹和拓扑转成amber后,用cpptraj中dbscan做密度聚类的簇分析,可以这样混合使用工具做分析吗?
2   通过上面方法得出的PCA三维散点图和簇分析结果能相互自洽吗,也就是说是不是PCA散点图中相对成簇的一些点正好是簇分析某一个簇中所有帧数的集合?

6万

帖子

99

威望

6万

eV
积分
125150

管理员

公社社长

9#
发表于 Post on 2021-6-20 17:48:20 | 只看该作者 Only view this author
landian666 发表于 2021-6-19 16:28
谢谢sob老师,问题解决了,还有两个问题想请教sob老师:
1  就是我在研究一个动力学过程中先用GROMACS做动 ...

1 原理上可以
2 不一定自洽,毕竟归簇算法不同,而且PC只考虑了前三个主成分,而帧之间归簇靠RMSD考虑的是所有原子
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

71

帖子

0

威望

694

eV
积分
765

Level 4 (黑子)

10#
 楼主 Author| 发表于 Post on 2021-6-22 15:29:20 | 只看该作者 Only view this author
sobereva 发表于 2021-6-20 17:48
1 原理上可以
2 不一定自洽,毕竟归簇算法不同,而且PC只考虑了前三个主成分,而帧之间归簇靠RMSD考虑的 ...

谢谢sob老师~

71

帖子

0

威望

694

eV
积分
765

Level 4 (黑子)

11#
 楼主 Author| 发表于 Post on 2021-6-25 14:48:41 | 只看该作者 Only view this author
sobereva 发表于 2021-6-20 17:48
1 原理上可以
2 不一定自洽,毕竟归簇算法不同,而且PC只考虑了前三个主成分,而帧之间归簇靠RMSD考虑的 ...

sob老师,我想针对一个结构中一部分残基进行簇分析,用的如下方法,想请教您是不是可以:
我用cpptraj中dbscan做密度聚类做簇分析,假如输入的结构有200个残基,我在cluster的rms命令中只考虑前27个残基(命令如 rms :1-27&!@H \),是不是可以说我是针对了这27个残基做了簇分析,并没有考虑其他残基对簇分析的影响?
完整命令如下,谢谢:
cpptraj
parm dna1.prmtop   #
trajin dna1.mdcrd     #输入结构含有200个残基

cluster C3 \
        dbscan minpoints 4 epsilon 2.9 sievetoframe \
        rms :1-27&!@H \
        sieve 10 \
        random \
        out cnumvtime.dat \
        summary summary.dat \
        info info.dat \
        cpopvtime cpopvtime.agr normframe \
        clusterout trajfile \
        repout rep repfmt pdb \
        singlerepout singlerep.nc singlerepfmt netcdf \
        avgout Avg avgfmt restart
run

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-23 05:56 , Processed in 0.184538 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list