计算化学公社

 找回密码 Forget password
 注册 Register

如何计算Rg、RMSD曲线图中的各自平均值?

查看数: 1552 | 评论数: 8 | 收藏 Add to favorites 0
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2024-6-21 21:25

正文摘要:

各位老师好,这是某篇文献中对于生物大分子体系分子动力学模拟计算部分,使用gromacs计算。图中标红部分是它提到的平均回旋半径值,想问一下各位老师这个平均回旋半径值,包括后面我还看到了平均RMSD值,这个ave ...

回复 Reply

12313 发表于 Post on 2024-6-24 10:43:24
dayu8278 发表于 2024-6-23 21:47
但是也得取平衡后的均值才具有统计意义?

谢谢大佬的提醒
dayu8278 发表于 Post on 2024-6-23 21:47:46
12313 发表于 2024-6-22 20:57
好像真的是这样,谢谢你

但是也得取平衡后的均值才具有统计意义?
12313 发表于 Post on 2024-6-22 20:57:26
dayu8278 发表于 2024-6-22 16:26
如果我没理解错,你把数据放在excel里不就能得到均值么?

好像真的是这样,谢谢你
12313 发表于 Post on 2024-6-22 20:56:38
student0618 发表于 2024-6-22 13:06
If you use gromacs to do analysis, gmx rms to calculate rmsd, gmx gyrate to calculate radius of gyra ...

I get it.Thanks  for your help
dayu8278 发表于 Post on 2024-6-22 16:26:27
如果我没理解错,你把数据放在excel里不就能得到均值么?
student0618 发表于 Post on 2024-6-22 13:06:10
本帖最后由 student0618 于 2024-6-22 14:14 编辑

If you use gromacs to do analysis, gmx rms to calculate rmsd, gmx gyrate to calculate radius of gyration. See the tutorial or help message on how to use the tool. The output file e.g. rmsd.xvg is actually a dataset with the rmsd/radius of gyration of each frame that allows you to plot the graph.

For example, here we have the contents of an rmsd.xvg I generated from Dr Justin Lemkul's gromacs tutorial 1:


On top it has some information from gromacs (lines start with '#'), then plotting commands for xmgrace/qtgrace about graph labels and title etc (lines start with '@'). The list of numbers below is the xy dataset with the following format:
time rmsd

That you can use the rmsd values to calculate average with whatever tool you want to use.
12313 发表于 Post on 2024-6-22 08:54:15
student0618 发表于 2024-6-22 01:50
If you have the value of each frame e.g. calculated by gmx analysis tools in a file, you can use man ...

sorry,but how to obtain the value of each frame?
student0618 发表于 Post on 2024-6-22 01:50:34
本帖最后由 student0618 于 2024-6-22 02:02 编辑

If you have the value of each frame e.g. calculated by gmx analysis tools in a file, you can use many statistical tools to get the average. From awk script to python or even excel.

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

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

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