CerdoFunk / Manipular-DM-Data

Geek Repo:Geek Repo

Github PK Tool:Github PK Tool

Manipular-MD-Data

Partiendo de una trayectoria de dinámica molecular de una membrana modelo podemos usar up_or_down.py para seleccionar la monocapa que queremos analizar. Esta función devuelve un AtomGroup (MDAnalysis), con todos los átomos en la monocapa superior o inferior.

Podemos usar lo que devuelva la función up_or_down para generar un archivo .csv, ya sea con la función cartesian_to_polar o acomodar_data que devuelven un archivo.csv con las columnas siendo la selección (residuo, átomo, etc) que hayamos seleccionado. Es decir, si selección = átomoP, entonces columnas serán átomoP_1, átomoP_2,..., átomoP_N y los renglones son los frames de la simulación. El valor de cada celda en este arreglo depende de la función que usemos.

Podemos usar dicho arreglo para calcular la matriz de correlación de Pearson, tal que podamos gráficar los resultados usando un mapa de calor. Usando lo que devuelva la función (cartesian_to_polar o acomodar_data) podemos utilizar normalized_pears_corr para gráficar un mapa de valor sólo con los puntos significativos.

About


Languages

Language:Python 100.0%