Dataset Information

Improving the speed of volumetric density map generation via cubic spline interpolation.

ABSTRACT: Visualizing data generated from molecular dynamics simulations can be difficult, particularly when there can be thousands to millions of trajectory frames. The creation of a 3D grid of atomic density (i.e. a volumetric map) is one way to easily view the long-time average behavior of a system. One way to generate volumetric maps is by approximating each atom with a Gaussian function centered on that atom and spread over neighboring grid cells. However the calculation of the Gaussian function requires evaluation of the exponential function, which is computationally costly. Here we report on speeding up the calculation of volumetric maps from molecular dynamics trajectory data by replacing the expensive exponential function evaluation with an approximation using interpolating cubic splines. We also discuss the errors involved in this approximation, and recommend settings for volumetric map creation based on this.

SUBMITTER: Roe DR

PROVIDER: S-EPMC7940575 | biostudies-literature | 2021 May

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Improving the speed of volumetric density map generation via cubic spline interpolation.

Roe Daniel R DR Brooks Bernard R BR

Journal of molecular graphics & modelling 20210104

Visualizing data generated from molecular dynamics simulations can be difficult, particularly when there can be thousands to millions of trajectory frames. The creation of a 3D grid of atomic density (i.e. a volumetric map) is one way to easily view the long-time average behavior of a system. One way to generate volumetric maps is by approximating each atom with a Gaussian function centered on that atom and spread over neighboring grid cells. However the calculation of the Gaussian function requ ...[more]

PMID: 33444979

Similar Datasets

Project description:BackgroundDropout is a common problem in longitudinal clinical trials and cohort studies, and is of particular concern when dropout occurs for reasons that may be related to the outcome of interest. This paper reviews common parametric models to account for dropout and introduces a Bayesian semi-parametric varying coefficient model for exponential family longitudinal data with non-ignorable dropout.MethodsTo demonstrate these methods, we present results from a simulation study and estimate the impact of drug use on longitudinal CD4 + T cell count and viral load suppression in the Women's Interagency HIV Study. Sensitivity analyses are performed to consider the impact of model assumptions on inference. We compare results between our semi-parametric method and parametric models to account for dropout, including the conditional linear model and a parametric frailty model. We also compare results to analyses that fail to account for dropout.ResultsIn simulation studies, we show that semi-parametric methods reduce bias and mean squared error when parametric model assumptions are violated. In analyses of the Women's Interagency HIV Study data, we find important differences in estimates of changes in CD4 + T cell count over time in untreated subjects that report drug use between different models used to account for dropout. We find steeper declines over time using our semi-parametric model, which makes fewer assumptions, compared to parametric models. Failing to account for dropout or to meet parametric assumptions of models to account for dropout could lead to underestimation of the impact of hard drug use on CD4 + cell count decline in untreated subjects. In analyses of subjects that initiated highly active anti-retroviral treatment, we find that the estimated probability of viral load suppression is lower in models that account for dropout.ConclusionsNon-ignorable dropout is an important consideration when analyzing data from longitudinal clinical trials and cohort studies. While methods that account for non-ignorable dropout must make some unavoidable assumptions that cannot be verified from the observed data, many methods make additional parametric assumptions. If these assumptions are not met, inferences can be biased, making more flexible methods with minimal assumptions important.

Dataset Information

Improving the speed of volumetric density map generation via cubic spline interpolation.

Publications

Improving the speed of volumetric density map generation via cubic spline interpolation.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets