Difference between revisions of "MD simulation analysis TSD Journal"
From Bioinformatikpedia
| Line 5: | Line 5: | ||
=Trajectories= |
=Trajectories= |
||
PNGs were created with [https://gist.github.com/3557863 createTrajectories.pml] and then combined into an animation with [https://gist.github.com/3557814 createAnim.sh]. |
PNGs were created with [https://gist.github.com/3557863 createTrajectories.pml] and then combined into an animation with [https://gist.github.com/3557814 createAnim.sh]. |
||
| + | |||
| + | =Thermodynamics= |
||
| + | Averages were calculated with R's "summary" function on grep-truncated xvg files. |
||
=Ramachandran plots= |
=Ramachandran plots= |
||
| Line 15: | Line 18: | ||
Additionally the data for the background distribution was downloaded from [http://www2.warwick.ac.uk/fac/sci/moac/people/students/peter_cock/r/ramachandran/#downloads|here].<br> |
Additionally the data for the background distribution was downloaded from [http://www2.warwick.ac.uk/fac/sci/moac/people/students/peter_cock/r/ramachandran/#downloads|here].<br> |
||
Then the scipt [https://gist.github.com/3557539|draw_rama.R] was called. |
Then the scipt [https://gist.github.com/3557539|draw_rama.R] was called. |
||
| + | =B-Factor= |
||
| + | Plots for comparing B-Factor based on all atoms and only C-alpha atoms were created with [https://gist.github.com/3557912 plotBFactors.pml]. |
||
=Significance calculation= |
=Significance calculation= |
||
Revision as of 20:32, 31 August 2012
Back to results.
Most of the base data was calculated according to the tutorial with a slightly modified version of scripts by Hemochromatosis. These are TODO.
Trajectories
PNGs were created with createTrajectories.pml and then combined into an animation with createAnim.sh.
Thermodynamics
Averages were calculated with R's "summary" function on grep-truncated xvg files.
Ramachandran plots
At first the .xvg files were parsed into the following format:
ID Psi Phi Aa 1HMP:chain0:Pro5 -92.926842 12.941312 Proline 1HMP:chain0:Gly6 65.796807 -162.229709 Glycine 1HMP:chain0:Val7 -81.132082 121.413022 General
Additionally the data for the background distribution was downloaded from [1].
Then the scipt [2] was called.
B-Factor
Plots for comparing B-Factor based on all atoms and only C-alpha atoms were created with plotBFactors.pml.
Significance calculation
The p-value was computed with the script statTD.pl.
