Hierarchical statistical techniques are necessary to draw reliable conclusions from analysis of isolated cardiomyocyte studies

Abstract

AimsIt is generally accepted that post-MI heart failure (HF) changes a variety of aspects of sarcoplasmic reticular Ca2+ fluxes but for some aspects there is disagreement over whether there is an increase or decrease. The commonest statistical approach is to treat data collected from each cell as independent, even though they are really clustered with multiple likely similar cells from each heart. In this study, we test whether this statistical assumption of independence can lead the investigator to draw conclusions that would be considered erroneous if the analysis handled clustering with specific statistical techniques (hierarchical tests).Methods and resultsCa2+ transients were recorded in cells loaded with Fura-2AM and sparks were recorded in cells loaded with Fluo-4AM. Data were analysed twice, once with the common statistical approach (assumption of independence) and once with hierarchical statistical methodologies designed to allow for any clustering. The statistical tests found that there was significant hierarchical clustering. This caused the common statistical approach to underestimate the standard error and report artificially small P values. For example, this would have led to the erroneous conclusion that time to 50% peak transient amplitude was significantly prolonged in HF.Spark analysis showed clustering, both within each cell and also within each rat, for morphological variables. This means that a three-level hierarchical model is sometimes required for such measures. Standard statistical methodologies, if used instead, erroneously suggest that spark amplitude is significantly greater in HF and spark duration is reduced in HF.ConclusionCa2+ fluxes in isolated cardiomyocytes show so much clustering that the common statistical approach that assumes independence of each data point will frequently give the false appearance of statistically significant changes. Hierarchical statistical methodologies need a little more effort, but are necessary for reliable conclusions. We present cost-free simple tools for performing these analyses.

Late sodium current inhibitors to treat exercise‐induced obstruction in hypertrophic cardiomyopathy: an in vitro study in human myocardium

“…Statistical analysis was performed as previously described using linear mixed models (Coppini et al, ), taking into account non‐Gaussian distribution, inequality of variances and within‐subject correlation. In brief, in order to reduce the risk of type I errors resulting from the closer interrelationship between cells/trabeculae isolated from the same patient sample, we used hierarchical statistics including two nested levels (patients and cells) (Sikkel et al, ); a third hierarchical level was added when drugs were tested in a cell/trabecula, in order to allow paired comparisons. This approach was implemented using linear mixed models in Stata 12.0 (StataCorp LLC, USA).…”

Profibrotic, Electrical, and Calcium-Handling Remodeling of the Atria in Heart Failure Patients With and Without Atrial Fibrillation

“…Categorical data were analyzed using a Chi-squared test with Bonferroni correction for multiple comparisons when comparing more than two groups. For patch-clamp and RyR2 single-channel recordings, in which each patient may contribute multiple data points, hierarchical statistics were employed according to recently published methods ( Sikkel et al, 2017 ). Logarithmic transformations were applied to non-normal-distributed data (RyR2 properties and frequency of spontaneous SR Ca 2+ -release events) before applying hierarchical statistical analyses.…”

Investigation of cardiac fibroblasts using myocardial slices

“…When multiple slices from multiple animals were considered, statistical analysis was performed using hierarchical statistical techniques (random intercept mixed model) as per Sikkel et al . 33 P values *< 0.05, **<0.01, and ***<0.001 were considered statistically significant.…”