Generalized eMC implementation for Monte Carlo dose calculation of electron beams from different machine types

The electron Monte Carlo (eMC) dose calculation algorithm available in the Eclipse treatment planning system (Varian Medical Systems) is based on the macro MC method and uses a beam model applicable to Varian linear accelerators. This leads to limitations in accuracy if eMC is applied to non-Varian machines. In this work eMC is generalized to also allow accurate dose calculations for electron beams from Elekta and Siemens accelerators. First, changes made in the previous study to use eMC for low electron beam energies of Varian accelerators are applied. Then, a generalized beam model is developed using a main electron source and a main photon source representing electrons and photons from the scattering foil, respectively, an edge source of electrons, a transmission source of photons and a line source of electrons and photons representing the particles from the scrapers or inserts and head scatter radiation. Regarding the macro MC dose calculation algorithm, the transport code of the secondary particles is improved. The macro MC dose calculations are validated with corresponding dose calculations using EGSnrc in homogeneous and inhomogeneous phantoms. The validation of the generalized eMC is carried out by comparing calculated and measured dose distributions in water for Varian, Elekta and Siemens machines for a variety of beam energies, applicator sizes and SSDs. The comparisons are performed in units of cGy per MU. Overall, a general agreement between calculated and measured dose distributions for all machine types and all combinations of parameters investigated is found to be within 2% or 2 mm. The results of the dose comparisons suggest that the generalized eMC is now suitable to calculate dose distributions for Varian, Elekta and Siemens linear accelerators with sufficient accuracy in the range of the investigated combinations of beam energies, applicator sizes and SSDs.

Monte Carlo source model for photon beam radiotherapy: photon source characteristics

A major barrier to widespread clinical implementation of Monte Carlo dose calculation is the difficulty in characterizing the radiation source within a generalized source model. This work aims to develop a generalized three-component source model (target, primary collimator, flattening filter) for 6- and 18-MV photon beams that match full phase-space data (PSD). Subsource by subsource comparison of dose distributions, using either source PSD or the source model as input, allows accurate source characterization and has the potential to ease the commissioning procedure, since it is possible to obtain information about which subsource needs to be tuned. This source model is unique in that, compared to previous source models, it retains additional correlations among PS variables, which improves accuracy at nonstandard source-to-surface distances (SSDs). In our study, three-dimensional (3D) dose calculations were performed for SSDs ranging from 50 to 200 cm and for field sizes from 1 x 1 to 30 x 30 cm2 as well as a 10 x 10 cm2 field 5 cm off axis in each direction. The 3D dose distributions, using either full PSD or the source model as input, were compared in terms of dose-difference and distance-to-agreement. With this model, over 99% of the voxels agreed within +/-1% or 1 mm for the target, within 2% or 2 mm for the primary collimator, and within +/-2.5% or 2 mm for the flattening filter in all cases studied. For the dose distributions, 99% of the dose voxels agreed within 1% or 1 mm when the combined source model-including a charged particle source and the full PSD as input-was used. The accurate and general characterization of each photon source and knowledge of the subsource dose distributions should facilitate source model commissioning procedures by allowing scaling the histogram distributions representing the subsources to be tuned.

Growth Codes: Intermediate Performance Analysis and Application to Video

Growth codes are a subclass of Rateless codes that have found interesting applications in data dissemination problems. Compared to other Rateless and conventional channel codes, Growth codes show improved intermediate performance which is particularly useful in applications where partial data presents some utility. In this paper, we investigate the asymptotic performance of Growth codes using the Wormald method, which was proposed for studying the Peeling Decoder of LDPC and LDGM codes. Compared to previous works, the Wormald differential equations are set on nodes' perspective which enables a numerical solution to the computation of the expected asymptotic decoding performance of Growth codes. Our framework is appropriate for any class of Rateless codes that does not include a precoding step. We further study the performance of Growth codes with moderate and large size codeblocks through simulations and we use the generalized logistic function to model the decoding probability. We then exploit the decoding probability model in an illustrative application of Growth codes to error resilient video transmission. The video transmission problem is cast as a joint source and channel rate allocation problem that is shown to be convex with respect to the channel rate. This illustrative application permits to highlight the main advantage of Growth codes, namely improved performance in the intermediate loss region.

Comparison of the effect of a CIDR-Select synch versus a long-term CIDR based AI protocol on reproductive performance in multiparous dairy cows in Swiss dairy farms

Background Synchronization programs have become standard in the dairy industry in many countries. In Switzerland, these programs are not routinely used for groups of cows, but predominantly as a therapy for individual problem cows. The objective of this study was to compare the effect of a CIDR-Select Synch and a 12-d CIDR protocol on the pregnancy rate in healthy, multiparous dairy cows in Swiss dairy farms. Methods Cows (N = 508) were randomly assigned to CIDR-Select Synch (N = 262) or 12-d CIDR (N = 246) protocols. Cows in the CIDR-Select Synch group received a CIDR and 2.5 ml of buserelin i.m. on d 0. On d 7, the CIDR insert was removed and 5 ml of dinoprost was administered i.m.. Cows in the 12-d CIDR group received the CIDR on d 0 and it was removed on d 12 (the routine CIDR protocol in Swiss dairies). On d 0 a milk sample for progesterone analysis was taken. Cows were inseminated upon observed estrus. Pregnancy was determined at or more than 35 days after artificial insemination. As a first step, the two groups were compared as to indication for treatment, breed, stud book, stall, pasture, and farmer's business using chi square tests or Fisher's exact test. Furthermore, groups were compared as to age, DIM, number of AI's, number of cows per farm, and yearly milk yield per cow using nonparametric ANOVA. A multiple logistic model was used to relate the success of the protocols to all of the available factors; in particular treatment (CIDR-Select Synch/12-d CIDR), milk progesterone value, age, DIM, previous treatment of the uterus, previous gynecological treatment, and number of preceding inseminations. Results The pregnancy rate was higher in cows following the CIDR-Select Synch compared to the 12-d CIDR protocol (50.4% vs. 22.4%; P < 0.0001). Conclusion The CIDR-Select Synch protocol may be highly recommended for multiparous dairy cows. The reduced time span of the progesterone insert decreased the number of days open, improved the pregnancy rate compared to the 12-d CIDR protocol and the cows did not to have to be handled more often.

Impact of weather and climate on the incidence of acute coronary syndromes

BACKGROUND: First investigations of the interactions between weather and the incidence of acute myocardial infarctions date back to 1938. The early observation of a higher incidence of myocardial infarctions in the cold season could be confirmed in very different geographical regions and cohorts. While the influence of seasonal variations on the incidence of myocardial infarctions has been extensively documented, the impact of individual meteorological parameters on the disease has so far not been investigated systematically. Hence the present study intended to assess the impact of the essential variables of weather and climate on the incidence of myocardial infarctions. METHODS: The daily incidence of myocardial infarctions was calculated from a national hospitalization survey. The hourly weather and climate data were provided by the database of the national weather forecast. The epidemiological and meteorological data were correlated by multivariate analysis based on a generalized linear model assuming a log-link-function and a Poisson distribution. RESULTS: High ambient pressure, high pressure gradients, and heavy wind activity were associated with an increase in the incidence of the totally 6560 hospitalizations for myocardial infarction irrespective of the geographical region. Snow- and rainfall had inconsistent effects. Temperature, Foehn, and lightning showed no statistically significant impact. CONCLUSIONS: Ambient pressure, pressure gradient, and wind activity had a statistical impact on the incidence of myocardial infarctions in Switzerland from 1990 to 1994. To establish a cause-and-effect relationship more data are needed on the interaction between the pathophysiological mechanisms of the acute coronary syndrome and weather and climate variables.

Circumcision and HIV infection among men who have sex with men in Britain: the insertive sexual role

The objective was to examine the association between circumcision status and self-reported HIV infection among men who have sex with men (MSM) in Britain who predominantly or exclusively engaged in insertive anal intercourse. In 2007-2008, a convenience sample of MSM living in Britain was recruited through websites, in sexual health clinics, bars, clubs, and other venues. Men completed an online survey which included questions on circumcision status, HIV testing, HIV status, sexual risk behavior, and sexual role for anal sex. The analysis was restricted to 1,521 white British MSM who reported unprotected anal intercourse in the previous 3 months and who said they only or mostly took the insertive role during anal sex. Of these men, 254 (16.7 %) were circumcised. Among men who had had a previous HIV test (n = 1,097), self-reported HIV seropositivity was 8.6 % for circumcised men (17/197) and 8.9 % for uncircumcised men (80/900) (unadjusted odds ratio [OR], 0.97; 95 % confidence interval [95 % CI], 0.56, 1.67). In a multivariable logistic model adjusted for known risk factors for HIV infection, there was no evidence of an association between HIV seropositivity and circumcision status (adjusted OR, 0.79; 95 % CI, 0.43, 1.44), even among the 400 MSM who engaged exclusively in insertive anal sex (adjusted OR, 0.84; 95 % CI, 0.25, 2.81). Our study provides further evidence that circumcision is unlikely to be an effective strategy for HIV prevention among MSM in Britain.

Expanding the cone location and magnitude index to include corneal thickness and posterior surface information for the detection of keratoconus

PURPOSE To extend the capabilities of the Cone Location and Magnitude Index algorithm to include a combination of topographic information from the anterior and posterior corneal surfaces and corneal thickness measurements to further improve our ability to correctly identify keratoconus using this new index: ConeLocationMagnitudeIndex_X. DESIGN Retrospective case-control study. METHODS Three independent data sets were analyzed: 1 development and 2 validation. The AnteriorCornealPower index was calculated to stratify the keratoconus data from mild to severe. The ConeLocationMagnitudeIndex algorithm was applied to all tomography data collected using a dual Scheimpflug-Placido-based tomographer. The ConeLocationMagnitudeIndex_X formula, resulting from analysis of the Development set, was used to determine the logistic regression model that best separates keratoconus from normal and was applied to all data sets to calculate PercentProbabilityKeratoconus_X. The sensitivity/specificity of PercentProbabilityKeratoconus_X was compared with the original PercentProbabilityKeratoconus, which only uses anterior axial data. RESULTS The AnteriorCornealPower severity distribution for the combined data sets are 136 mild, 12 moderate, and 7 severe. The logistic regression model generated for ConeLocationMagnitudeIndex_X produces complete separation for the Development set. Validation Set 1 has 1 false-negative and Validation Set 2 has 1 false-positive. The overall sensitivity/specificity results for the logistic model produced using the ConeLocationMagnitudeIndex_X algorithm are 99.4% and 99.6%, respectively. The overall sensitivity/specificity results for using the original ConeLocationMagnitudeIndex algorithm are 89.2% and 98.8%, respectively. CONCLUSIONS ConeLocationMagnitudeIndex_X provides a robust index that can detect the presence or absence of a keratoconic pattern in corneal tomography maps with improved sensitivity/specificity from the original anterior surface-only ConeLocationMagnitudeIndex algorithm.

The Effect Of Misleading Prior Information On The Subjective Visual Vertical

When tilted sideways participants misperceive the visual vertical assessed by means of a luminous line in otherwise complete dark- ness. A recent modeling approach (De Vrijer et al., 2009) claimed that these typical patterns of errors (known as A- and E-effects) could be explained by as- suming that participants behave in a Bayes optimal manner. In this study, we experimentally manipulate participants’ prior information about body-in-space orientation and measure the effect of this manipulation on the subjective visual vertical (SVV). Specifically, we explore the effects of veridical and misleading instructions about body tilt orientations on the SVV. We used a psychophys- ical 2AFC SVV task at roll tilt angles of 0 degrees, 16 degrees and 4 degrees CW and CCW. Participants were tilted to 4 degrees under different instruction conditions: in one condition, participants received veridical instructions as to their tilt angle, whereas in another condition, participants received the mis- leading instruction that their body position was perfectly upright. Our results indicate systematic differences between the instruction conditions at 4 degrees CW and CCW. Participants did not simply use an ego-centric reference frame in the misleading condition; instead, participants’ estimates of the SVV seem to lie between their head’s Z-axis and the estimate of the SVV as measured in the veridical condition. All participants displayed A-effects at roll tilt an- gles of 16 degrees CW and CCW. We discuss our results in the context of the Bayesian model by De Vrijer et al. (2009), and claim that this pattern of re- sults is consistent with a manipulation of precision of a prior distribution over body-in-space orientations. Furthermore, we introduce a Bayesian Generalized Linear Model for estimating parameters of participants’ psychometric function, which allows us to jointly estimate group level and individual level parameters under all experimental conditions simultaneously, rather than relying on the traditional two-step approach to obtaining group level parameter estimates.

gems: An R Package for Simulating from Disease Progression Models

Mathematical models of disease progression predict disease outcomes and are useful epidemiological tools for planners and evaluators of health interventions. The R package gems is a tool that simulates disease progression in patients and predicts the effect of different interventions on patient outcome. Disease progression is represented by a series of events (e.g., diagnosis, treatment and death), displayed in a directed acyclic graph. The vertices correspond to disease states and the directed edges represent events. The package gems allows simulations based on a generalized multistate model that can be described by a directed acyclic graph with continuous transition-specific hazard functions. The user can specify an arbitrary hazard function and its parameters. The model includes parameter uncertainty, does not need to be a Markov model, and may take the history of previous events into account. Applications are not limited to the medical field and extend to other areas where multistate simulation is of interest. We provide a technical explanation of the multistate models used by gems, explain the functions of gems and their arguments, and show a sample application.

Thrombus formation in the left ventricle after large myocardial infarction – assessment with cardiac magnetic resonance imaging

INTRODUCTION Left ventricular thrombus (LVT) formation may worsen the post-infarct outcome as a result of thromboembolic events. It also complicates the use of modern antiplatelet regimens, which are not compatible with long-term oral anticoagulation. The knowledge of the incidence of LVT may therefore be of importance to guide antiplatelet and antithrombotic therapy after acute myocardial infarction (AMI). METHODS In 177 patients with large, mainly anterior AMI, standard cardiac magnetic resonance imaging (CMR) including cine and late gadolinium enhancement (LGE) imaging was performed shortly after AMI as per protocol. CMR images were analysed at an independent core laboratory blinded to the clinical data. Transthoracic echocardiography (TTE) was not mandatory for the trial, but was performed in 64% of the cases following standard of care. In a logistic model, 3 out of 61 parameters were used in a multivariable model to predict LVT. RESULTS LVT was detected by use of CMR in 6.2% (95% confidence interval [CI] 3.1%-10.8%). LGE sequences were best to detect LVT, which may be missed in cine sequences. We identified body mass index (odds ratio 1.18; p = 0.01), baseline platelet count (odds ratio 1.01, p = 0.01) and infarct size as assessed by use of CMR (odds ratio 1.03, p = 0.02) as best predictors for LVT. The agreement between TTE and CMR for the detection of LVT is substantial (kappa = 0.70). DISCUSSION In the current analysis, the incidence of LVT shortly after AMI is relatively low, even in a patient population at high risk. An optimal modality for LVT detection is LGE-CMR but TTE has an acceptable accuracy when LGE-CMR is not available.

Age-specific and sex-specific weight gain norms to monitor antiretroviral therapy in children in low-income and middle-income countries.

BACKGROUND Viral load and CD4% are often not available in resource-limited settings for monitoring children's responses to antiretroviral therapy (ART). We aimed to construct normative curves for weight gain at 6, 12, 18, and 24 months following initiation of ART in children, and to assess the association between poor weight gain and subsequent responses to ART. DESIGN Analysis of data from HIV-infected children younger than 10 years old from African and Asian clinics participating in the International epidemiologic Databases to Evaluate AIDS. METHODS The generalized additive model for location, scale, and shape was used to construct normative percentile curves for weight gain at 6, 12, 18, and 24 months following ART initiation. Cox proportional models were used to assess the association between lower percentiles (< 50th) of weight gain distribution at the different time points and subsequent death, virological suppression, and virological failure. RESULTS Among 7173 children from five regions of the world, 45% were underweight at baseline. Weight gain below the 50th percentile at 6, 12, 18, and 24 months of ART was associated with increased risk of death, independent of baseline characteristics. Poor weight gain was not associated with increased hazards of virological suppression or virological failure. CONCLUSION Monitoring weight gain on ART using age-specific and sex-specific normative curves specifically developed for HIV-infected children on ART is a simple, rapid, sustainable tool that can aid in the identification of children who are at increased risk of death in the first year of ART.

Synchrotron X-ray interlaced microbeams suppress paroxysmal oscillations in neuronal networks initiating generalized epilepsy

Radiotherapy has shown some efficacy for epilepsies but the insufficient confinement of the radiation dose to the pathological target reduces its indications. Synchrotron-generated X-rays overcome this limitation and allow the delivery of focalized radiation doses to discrete brain volumes via interlaced arrays of microbeams (IntMRT). Here, we used IntMRT to target brain structures involved in seizure generation in a rat model of absence epilepsy (GAERS). We addressed the issue of whether and how synchrotron radiotherapeutic treatment suppresses epileptic activities in neuronal networks. IntMRT was used to target the somatosensory cortex (S1Cx), a region involved in seizure generation in the GAERS. The antiepileptic mechanisms were investigated by recording multisite local-field potentials and the intracellular activity of irradiated S1Cx pyramidal neurons in vivo. MRI and histopathological images displayed precise and sharp dose deposition and revealed no impairment of surrounding tissues. Local-field potentials from behaving animals demonstrated a quasi-total abolition of epileptiform activities within the target. The irradiated S1Cx was unable to initiate seizures, whereas neighboring non-irradiated cortical and thalamic regions could still produce pathological oscillations. In vivo intracellular recordings showed that irradiated pyramidal neurons were strongly hyperpolarized and displayed a decreased excitability and a reduction of spontaneous synaptic activities. These functional alterations explain the suppression of large-scale synchronization within irradiated cortical networks. Our work provides the first post-irradiation electrophysiological recordings of individual neurons. Altogether, our data are a critical step towards understanding how X-ray radiation impacts neuronal physiology and epileptogenic processes.

Missing C-terminal filaggrin expression, NFkappaB activation and hyperproliferation identify the dog as a putative model to study epidermal dysfunction in atopic dermatitis

Filaggrin loss-of-function mutations resulting in C-terminal protein truncations are strong predisposing factors in human atopic dermatitis (AD). To assess the possibility of similar truncations in canine AD, an exclusion strategy was designed on 16 control and 18 AD dogs of various breeds. Comparative immunofluorescence microscopy was performed with an antibody raised against the canine filaggrin C-terminus and a commercial N-terminal antibody. Concurrent with human AD-like features such as generalized NFKB activation and hyperproliferation, four distinctive filaggrin expression patterns were identified in non-lesional skin. It was found that 10/18 AD dogs exhibited an identical pattern for both antibodies with comparable (category I, 3/18) or reduced (category II, 7/18) expression to that of controls. In contrast, 4/18 dogs displayed aberrant large vesicles revealed by the C-terminal but not the N-terminal antibody (category III), while 4/18 showed a control-like N-terminal expression but lacked the C-terminal protein (category IV). The missing C-terminal filaggrin in category IV strongly points towards loss-of function mutations in 4/18 (22%) of all AD dogs analysed.

A coronary heart disease risk model for predicting the effect of potent antiretroviral therapy in HIV-1 infected men

BACKGROUND: Many HIV-infected patients on highly active antiretroviral therapy (HAART) experience metabolic complications including dyslipidaemia and insulin resistance, which may increase their coronary heart disease (CHD) risk. We developed a prognostic model for CHD tailored to the changes in risk factors observed in patients starting HAART. METHODS: Data from five cohort studies (British Regional Heart Study, Caerphilly and Speedwell Studies, Framingham Offspring Study, Whitehall II) on 13,100 men aged 40-70 and 114,443 years of follow up were used. CHD was defined as myocardial infarction or death from CHD. Model fit was assessed using the Akaike Information Criterion; generalizability across cohorts was examined using internal-external cross-validation. RESULTS: A parametric model based on the Gompertz distribution generalized best. Variables included in the model were systolic blood pressure, total cholesterol, high-density lipoprotein cholesterol, triglyceride, glucose, diabetes mellitus, body mass index and smoking status. Compared with patients not on HAART, the estimated CHD hazard ratio (HR) for patients on HAART was 1.46 (95% CI 1.15-1.86) for moderate and 2.48 (95% CI 1.76-3.51) for severe metabolic complications. CONCLUSIONS: The change in the risk of CHD in HIV-infected men starting HAART can be estimated based on typical changes in risk factors, assuming that HRs estimated using data from non-infected men are applicable to HIV-infected men. Based on this model the risk of CHD is likely to increase, but increases may often be modest, and could be offset by lifestyle changes.

On Graphically Checking Goodness-of-fit of Binary Logistic Regression Models

OBJECTIVES: This paper is concerned with checking goodness-of-fit of binary logistic regression models. For the practitioners of data analysis, the broad classes of procedures for checking goodness-of-fit available in the literature are described. The challenges of model checking in the context of binary logistic regression are reviewed. As a viable solution, a simple graphical procedure for checking goodness-of-fit is proposed. METHODS: The graphical procedure proposed relies on pieces of information available from any logistic analysis; the focus is on combining and presenting these in an informative way. RESULTS: The information gained using this approach is presented with three examples. In the discussion, the proposed method is put into context and compared with other graphical procedures for checking goodness-of-fit of binary logistic models available in the literature. CONCLUSION: A simple graphical method can significantly improve the understanding of any logistic regression analysis and help to prevent faulty conclusions.