Power Calculations for Preclinical Studies Using a K-Sample Rank Test and the Lehmann Alternative Hypothesis

Power calculations in a small sample comparative study, with a continuous outcome measure, are typically undertaken using the asymptotic distribution of the test statistic. When the sample size is small, this asymptotic result can be a poor approximation. An alternative approach, using a rank based test statistic, is an exact power calculation. When the number of groups is greater than two, the number of calculations required to perform an exact power calculation is prohibitive. To reduce the computational burden, a Monte Carlo resampling procedure is used to approximate the exact power function of a k-sample rank test statistic under the family of Lehmann alternative hypotheses. The motivating example for this approach is the design of animal studies, where the number of animals per group is typically small.

Arterio-venous gradients of free energy change for assessment of systemic and splanchnic perfusion in cardiac surgery patients

OBJECTIVE: Adequacy of organ perfusion depends on sufficient oxygen supply in relation to the metabolic needs. The aim of this study was to evaluate the relationship between gradients of free energy change, and the more commonly used parameter for the evaluation of the adequacy of organ perfusion, such as oxygen-extraction in patients undergoing valve replacement surgery using normothermic cardiopulmonary bypass (CPB). METHODS: In 43 cardiac patients, arterial, mixed venous, and hepato-venous blood samples were taken synchronously after induction of anaesthesia (preCPB), during CPB, and 2 and 7 h after admission to the intensive care unit (ICU+2, ICU+7). Blood gas analysis, cardiac output, and hepato-splanchnic blood flow were measured. Free energy change gradients between mixed venous and arterial (-deltadeltaG(v - a)) and hepato-venous and arterial (-deltadeltaG(hv - a)) compartments were calculated. MEASUREMENTS AND RESULTS: Cardiac index (CI) increased from 1.9 (0.7) to 2.8 (1.3) L/min/m (median, inter-quartile range) (p = 0.001), and hepato-splanchnic blood flow index (HBFI) from 0.6 (0.22) to 0.8 (0.53) L/min/m (p = 0.001). Despite increasing flow, systemic oxygen extraction increased after CPB from 24 (10)% to 35 (10)% at ICU+2 (p = 0.002), and splanchnic oxygen extraction increased during CPB from 37 (19)% to 52 (14)% (p = 0.001), and remained high thereafter. After CPB, high splanchnic and systemic gradients of free energy change gradients were associated with high splanchnic and systemic oxygen extraction, respectively (p = 0.001, 0.033, respectively). CONCLUSION: Gradients of free energy change may be helpful in characterising adequacy of perfusion in cardiac surgery patients independently from measurements or calculations of data from oxygen transport.

Self-Reported Memory Symptoms with Coronary Artery Disease: A Prospective of CABG Patients and Nonsurgical Controls

Background. Subjective memory complaints are common after coronary artery bypass grafting (CABG), but previous studies have concluded that such symptoms are more closely associated with depressed mood than objective cognitive dysfunction. We compared the incidence of self-reported memory symptoms at 3 and 12 months after CABG with that of a control group of patients with comparable risk factors for coronary artery disease but without surgery. Methods. Patients undergoing CABG (n = 140) and a demographically similar nonsurgical control group with coronary artery disease (n = 92) were followed prospectively at 3 and 12 months. At each follow-up time, participants were asked about changes since the previous evaluation in areas of memory, calculations, reading, and personality. A Functional Status Questionnaire (FSQ) and self-report measure of symptoms of depression (CES-D) were also completed. Results. The frequency of self-reported changes in memory, personality, and reading at 3 months was significantly higher among CABG patients than among nonsurgical controls. By contrast, there were no differences in the frequency of self-reported symptoms relating to calculations or overall rating of functional status. After adjusting for a measure of depression (CES-D rating score), the risk for self-reported memory changes remained nearly 5 times higher among the CABG patients than control subjects. The relative risk of developing new self-reported memory symptoms between 3 and 12 months was 2.5 times higher among CABG patients than among nonsurgical controls (CI 1.24 – 5.02), and the overall prevalence of memory symptoms at 12 months was also higher among CABG patients (39%) than controls (14%). Conclusions. The frequency of self-reported memory symptoms 3 and 12 months after baseline is significantly higher among CABG patients than control patients with comparable risk factors for coronary and cerebrovascular disease. These differences could not be accounted for by symptoms of depression. The self-reported cognitive symptoms appear to be relatively specific for memory, and may reflect aspects of memory functioning that are not captured by traditional measures of new verbal learning and memory. The etiology of these self-reported memory symptoms remains unclear, but our findings as well as those of others, may implicate factors other than cardiopulmonary bypass itself.

Collapsed cone convolution and analytical anisotropic algorithm dose calculations compared to VMC++ Monte Carlo simulations in clinical cases

The purpose of this work was to study and quantify the differences in dose distributions computed with some of the newest dose calculation algorithms available in commercial planning systems. The study was done for clinical cases originally calculated with pencil beam convolution (PBC) where large density inhomogeneities were present. Three other dose algorithms were used: a pencil beam like algorithm, the anisotropic analytic algorithm (AAA), a convolution superposition algorithm, collapsed cone convolution (CCC), and a Monte Carlo program, voxel Monte Carlo (VMC++). The dose calculation algorithms were compared under static field irradiations at 6 MV and 15 MV using multileaf collimators and hard wedges where necessary. Five clinical cases were studied: three lung and two breast cases. We found that, in terms of accuracy, the CCC algorithm performed better overall than AAA compared to VMC++, but AAA remains an attractive option for routine use in the clinic due to its short computation times. Dose differences between the different algorithms and VMC++ for the median value of the planning target volume (PTV) were typically 0.4% (range: 0.0 to 1.4%) in the lung and -1.3% (range: -2.1 to -0.6%) in the breast for the few cases we analysed. As expected, PTV coverage and dose homogeneity turned out to be more critical in the lung than in the breast cases with respect to the accuracy of the dose calculation. This was observed in the dose volume histograms obtained from the Monte Carlo simulations.

Conversion of CT numbers into tissue parameters for Monte Carlo dose calculations: a multi-centre study

The conversion of computed tomography (CT) numbers into material composition and mass density data influences the accuracy of patient dose calculations in Monte Carlo treatment planning (MCTP). The aim of our work was to develop a CT conversion scheme by performing a stoichiometric CT calibration. Fourteen dosimetrically equivalent tissue subsets (bins), of which ten bone bins, were created. After validating the proposed CT conversion scheme on phantoms, it was compared to a conventional five bin scheme with only one bone bin. This resulted in dose distributions D(14) and D(5) for nine clinical patient cases in a European multi-centre study. The observed local relative differences in dose to medium were mostly smaller than 5%. The dose-volume histograms of both targets and organs at risk were comparable, although within bony structures D(14) was found to be slightly but systematically higher than D(5). Converting dose to medium to dose to water (D(14) to D(14wat) and D(5) to D(5wat)) resulted in larger local differences as D(5wat) became up to 10% higher than D(14wat). In conclusion, multiple bone bins need to be introduced when Monte Carlo (MC) calculations of patient dose distributions are converted to dose to water.

A multiple source model for 6 MV photon beam dose calculations using Monte Carlo

A multiple source model (MSM) for the 6 MV beam of a Varian Clinac 2300 C/D was developed by simulating radiation transport through the accelerator head for a set of square fields using the GEANT Monte Carlo (MC) code. The corresponding phase space (PS) data enabled the characterization of 12 sources representing the main components of the beam defining system. By parametrizing the source characteristics and by evaluating the dependence of the parameters on field size, it was possible to extend the validity of the model to arbitrary rectangular fields which include the central 3 x 3 cm2 field without additional precalculated PS data. Finally, a sampling procedure was developed in order to reproduce the PS data. To validate the MSM, the fluence, energy fluence and mean energy distributions determined from the original and the reproduced PS data were compared and showed very good agreement. In addition, the MC calculated primary energy spectrum was verified by an energy spectrum derived from transmission measurements. Comparisons of MC calculated depth dose curves and profiles, using original and PS data reproduced by the MSM, agree within 1% and 1 mm. Deviations from measured dose distributions are within 1.5% and 1 mm. However, the real beam leads to some larger deviations outside the geometrical beam area for large fields. Calculated output factors in 10 cm water depth agree within 1.5% with experimentally determined data. In conclusion, the MSM produces accurate PS data for MC photon dose calculations for the rectangular fields specified.

Simple beam models for Monte Carlo photon beam dose calculations in radiotherapy

Monte Carlo (code GEANT) produced 6 and 15 MV phase space (PS) data were used to define several simple photon beam models. For creating the PS data the energy of starting electrons hitting the target was tuned to get correct depth dose data compared to measurements. The modeling process used the full PS information within the geometrical boundaries of the beam including all scattered radiation of the accelerator head. Scattered radiation outside the boundaries was neglected. Photons and electrons were assumed to be radiated from point sources. Four different models were investigated which involved different ways to determine the energies and locations of beam particles in the output plane. Depth dose curves, profiles, and relative output factors were calculated with these models for six field sizes from 5x5 to 40x40cm2 and compared to measurements. Model 1 uses a photon energy spectrum independent of location in the PS plane and a constant photon fluence in this plane. Model 2 takes into account the spatial particle fluence distribution in the PS plane. A constant fluence is used again in model 3, but the photon energy spectrum depends upon the off axis position. Model 4, finally uses the spatial particle fluence distribution and off axis dependent photon energy spectra in the PS plane. Depth dose curves and profiles for field sizes up to 10x10cm2 were not model sensitive. Good agreement between measured and calculated depth dose curves and profiles for all field sizes was reached for model 4. However, increasing deviations were found for increasing field sizes for models 1-3. Large deviations resulted for the profiles of models 2 and 3. This is due to the fact that these models overestimate and underestimate the energy fluence at large off axis distances. Relative output factors consistent with measurements resulted only for model 4.

Tenofovir use is associated with a reduction in calculated glomerular filtration rates in the Swiss HIV Cohort Study

BACKGROUND: A growing number of case reports have described tenofovir (TDF)-related proximal renal tubulopathy and impaired calculated glomerular filtration rates (cGFR). We assessed TDF-associated changes in cGFR in a large observational HIV cohort. METHODS: We compared treatment-naive patients or patients with treatment interruptions > or = 12 months starting either a TDF-based combination antiretroviral therapy (cART) (n = 363) or a TDF-sparing regime (n = 715). The predefined primary endpoint was the time to a 10 ml/min reduction in cGFR, based on the Cockcroft-Gault equation, confirmed by a follow-up measurement at least 1 month later. In sensitivity analyses, secondary endpoints including calculations based on the modified diet in renal disease (MDRD) formula were considered. Endpoints were modelled using pre-specified covariates in a multiple Cox proportional hazards model. RESULTS: Two-year event-free probabilities were 0.65 (95% confidence interval [CI] 0.58-0.72) and 0.80 (95% CI 0.76-0.83) for patients starting TDF-containing or TDF-sparing cART, respectively. In the multiple Cox model, diabetes mellitus (hazard ratio [HR] = 2.34 [95% CI 1.24-4.42]), higher baseline cGFR (HR = 1.03 [95% CI 1.02-1.04] by 10 ml/min), TDF use (HR = 1.84 [95% CI 1.35-2.51]) and boosted protease inhibitor use (HR = 1.71 [95% CI 1.30-2.24]) significantly increased the risk for reaching the primary endpoint. Sensitivity analyses showed high consistency. CONCLUSION: There is consistent evidence for a significant reduction in cGFR associated with TDF use in HIV-infected patients. Our findings call for a strict monitoring of renal function in long-term TDF users with tests that distinguish between glomerular dysfunction and proximal renal tubulopathy, a known adverse effect of TDF.