A note on parallelizing the parameterized expectations algorithm

The parameterized expectations algorithm (PEA) involves a long simulation and a nonlinear least squares (NLS) fit, both embedded in a loop. Both steps are natural candidates for parallelization. This note shows that parallelization can lead to important speedups for the PEA. I provide example code for a simple model that can serve as a template for parallelization of more interesting models, as well as a download link for an image of a bootable CD that allows creation of a cluster and execution of the example code in minutes, with no need to install any software.

Auctions for government securities: a laboratory comparison of uniform, discriminatory and spanish designs

The Bank of Spain uses a unique auction format to sell government bonds, which can be seen as a hybrid of a uniform and a discriminatory auction. For winning bids above the average winning bid, buyers are charged the average winning bid, otherwise they pay their respective bids. We report on an experiment that compares this auction format to the discriminatory format, used in most other countries, and to the uniform format. Our design is based on a common value model with multi-unit supply and two-unit demand. The results show significantly higher revenue with the Spanish and the uniform formats than with the discriminatory one, while volatility of prices over time is significantly lower in the discriminatory format than in the Spanish and uniform cases. Actual price dispersion is significantly larger in the discriminatory than in the Spanish. Our data also exhibit the use of bid-spreading strategies in all three designs.

How universal is behavior? A four country comparison of spite, cooperation and errors in voluntary contribution mechanisms

This paper studies behavior in experiments with a linear voluntary contributions mechanism for public goods conducted in Japan, the Netherlands, Spain and the USA. The same experimental design was used in the four countries. Our 'contribution function' design allows us to obtain a view of subjects' behavior from two complementary points of view. If yields information about situations where, in purely pecuniary terms, it is a dominant strategy to contribute all the endowment and about situations where it is a dominant strategy to contribute nothing. Our results show, first, that differences in behavior across countries are minor. We find that when people play "the same game" they behave similarly. Second, for all four countries our data are inconsistent with the explanation that subjects contribute only out of confusion. A common cooperative motivation is needed to explain the date.

A Comparison of polarization measures

This paper provides a systematic classification of the different measures of polarization based on their properties. Together with the axioms proposed in Duclos, Esteban and Ray (2004) and in Wang and Tsui (2000) we consider three additional properties. We examine which properties are common to all indices and which set them apart.

Simulación Monte Carlo de las dosis en tumores de pulmón en tratamientos de radioterapia estereotáxica extracraneal

Estudio elaborado a partir de una estancia en el Karolinska University Hospital, Suecia, entre marzo y junio del 2006. En la radioterapia estereotáxica extracraneal (SBRT) de tumores de pulmón existen principalmente dos problemas en el cálculo de la dosis con los sistemas de planificación disponibles: la precisión limitada de los algoritmos de cálculo en presencia de tejidos con densidades muy diferentes y los movimientos debidos a la respiración del paciente durante el tratamiento. El objetivo de este trabajo ha sido llevar a cabo la simulación con el código Monte Carlo PENELOPE de la distribución de dosis en tumores de pulmón en casos representativos de tratamientos con SBRT teniendo en cuenta los movimientos respiratorios y su comparación con los resultados de varios planificadores. Se han estudiado casos representativos de tratamientos de SBRT en el Karolinska University Hospital. Los haces de radiación se han simulado mediante el código PENELOPE y se han usado para la obtención de los resultados MC de perfiles de dosis. Los resultados obtenidos para el caso estático (sin movimiento respiratorio ) ponen de manifiesto que, en comparación con la MC, la dosis (Gy/MU) calculada por los planificadores en el tumor tiene una precisión del 2-3%. En la zona de interfase entre tumor y tejido pulmonar los planificadores basados en el algoritmo PB sobrestiman la dosis en un 10%, mientras que el algoritmo CC la subestima en un 3-4%. Los resultados de la simulación mediante MC de los movimientos respiratorios indican que los resultados de los planificadores son suficientemente precisos en el tumor, aunque en la interfase hay una mayor subestimación de la dosis en comparación con el caso estático. Estos resultados son compatibles con la experiencia clínica adquirida durante 15 años en el Karolinska University Hospital. Los resultados se han publicado en la revista Acta Oncologica.

Calibración TLR y LRM en MATLAB para medidas con guías de onda

En este proyecto, lo primero que hemos hecho ha sido desarrollar un algoritmo en Matlab que implementara el método de calibración TRL, el funcionamiento del cual hemos comprobado en primera instancia mediante simulaciones y, posteriormente, mediante un ejemplo real. Posteriormente, hemos desarrollado otro algoritmo en Matlab para implementar el método de calibración LRM. Este algoritmo sólo lo hemos podido comprobar a nivel de simulación. A continuación, mediante los dos algoritmos, hemos realizado una comparación entre ambos sistemas de calibración a través de simulaciones. Finalmente, analizando los resultados de varias simulaciones calibradas con nuestro programa del método TRL, hemos buscado cuáles pueden ser los motivos para la aparición de picos indeseados y hemos encontrado uno de ellos.

RS virus diagnosis: comparison of isolation, immunofluorescence and enzyme immunoassay

Two techniques for rapid diagnosis, immunofluorescence (IFAT) and enzyme immunoassay (EIA), have been compared with virus isolaion in tissue culture for the detection of respiratory syncytial virus (RSV) in specimens of nasopharyngeal secretions. The specimens were obtained from children under five years of age suffering from acute respiratory iliness, during a period of six months from January to June 1982. Of 471 specimens examined 54 (11.5%) were positive by virus isolation and 180 (38.2%) were positive by immunofluorescence. The bacterial contamination of inoculated tissue cultures unfortunately prevented the isolation of virus from many samples. Specimens from 216 children were tested to compare enzyme immunoassay and immunofluorescence. Of these 60 (27%) were positive by EIA and 121 (56%) were positive by IFAT. Our results suggest that the EIA technique although highly specific is rather insensitive. This may be because by the time these tests were done the originl nasopharyngeal secretions were considerably diluted and contained more mucus fragments than the call suspension used for IFAT. Of the three techniques, IFAT gives the best results although EIA may be useful where IFAT is not possible.

A comparison of simian rotavirus SA11 preparations maintained in different laboratories

Preparations of simian SA11 maintained in different laboratories were compared with each other by polyacrylamide gel electrophoresis of genomic RNA. Differences in the migration of genome segments 4,5 and 7 allowed the classification of eight virus preparations into four electrophoretic types.

Octopus Standard Automated Perimetry, Pulsar Perimetry, Moorfields Motion Displacement Test and Heidelberg Retinal Tomography in Glaucoma Detection - A Comparison of Diagnostic Accuracy

Purpose: To investigate the accuracy of 4 clinical instruments in the detection of glaucomatous damage. Methods: 102 eyes of 55 test subjects (Age mean = 66.5yrs, range = [39; 89]) underwent Heidelberg Retinal Tomography (HRTIII), (disc area<2.43); and standard automated perimetry (SAP) using Octopus (Dynamic); Pulsar (TOP); and Moorfields Motion Displacement Test (MDT) (ESTA strategy). Eyes were separated into three groups 1) Healthy (H): IOP<21mmHg and healthy discs (clinical examination), 39 subjects, 78 eyes; 2) Glaucoma suspect (GS): Suspicious discs (clinical examination), 12 subjects, 15 eyes; 3) Glaucoma (G): progressive structural or functional loss, 14 subjects, 20 eyes. Clinical diagnostic precision was examined using the cut-off associated with the p<5% normative limit of MD (Octopus/Pulsar), PTD (MDT) and MRA (HRT) analysis. The sensitivity, specificity and accuracy were calculated for each instrument. Results: See table Conclusions: Despite the advantage of defining glaucoma suspects using clinical optic disc examination, the HRT did not yield significantly higher accuracy than functional measures. HRT, MDT and Octopus SAP yielded higher accuracy than Pulsar perimetry, although results did not reach statistical significance. Further studies are required to investigate the structure-function correlations between these instruments.

Results of computer assisted navigation in Total Knee Arthroplasty (TKA): comparison with the conventional method

The aim of this retrospective study was to compare the clinical and radiographic results after TKA (PFC, DePuy), performed either by computer assisted navigation (CAS, Brainlab, Johnson&Johnson) or by conventional means. Material and methods: Between May and December 2006 we reviewed 36 conventional TKA performed between 2002 and 2003 (group A) and 37 navigated TKA performed between 2005 and 2006 (group B) by the same experienced surgeon. The mean age in group A was 74 years (range 62-90) and 73 (range 58-85) in group B with a similar age distribution. The preoperative mechanical axes in group A ranged from -13° varus to +13° valgus (mean absolute deviation 6.83°, SD 3.86), in group B from -13° to +16° (mean absolute deviation 5.35, SD 4.29). Patients with a previous tibial osteotomy or revision arthroplasty were excluded from the study. Examination was done by an experienced orthopedic resident independent of the surgeon. All patients had pre- and postoperative long standing radiographs. The IKSS and the WOMAC were utilized to determine the clinical outcome. Patient's degree of satisfaction was assessed on a visual analogous scale (VAS). Results: 32 of the 37 navigated TKAs (86,5%) showed a postoperative mechanical axis within the limits of 3 degrees of valgus or varus deviation compared to only 24 (66%) of the 36 standard TKAs. This difference was significant (p = 0.045). The mean absolute deviation from neutral axis was 3.00° (range -5° to +9°, SD: 1.75) in group A in comparison to 1.54° (range -5° to +4°, SD: 1.41) in group B with a highly significant difference (p = 0.000). Furthermore, both groups showed a significant postoperative improvement of their mean IKSS-values (group A: 89 preoperative to 169 postoperative, group B 88 to 176) without a significant difference between the two groups. Neither the WOMAC nor the patient's degree of satisfaction - as assessed by VAS - showed significant differences. Operation time was significantly higher in group B (mean 119.9 min.) than in group A (mean 99.6 min., p <0.000). Conclusion: Our study showed consistent significant improvement of postoperative frontal alignment in TKA by computer assisted navigation (CAS) compared to standard methods, even in the hands of a surgeon well experienced in standard TKA implantation. However, the follow-up time of this study was not long enough to judge differences in clinical outcome. Thus, the relevance of computer navigation for clinical outcome and survival of TKA remains to be proved in long term studies to justify the longer operation time. References 1 Stulberg SD. Clin Orth Rel Res. 2003;(416):177-84. 2 Chauhan SK. JBJS Br. 2004;86(3):372-7. 3 Bäthis H, et al. Orthopäde. 2006;35(10):1056-65.

Albuterol delivery in an in vitro pediatric ventilator lung model: comparison of jet, ultrasonic, and mesh nebulizers.

OBJECTIVE: : To determine the influence of nebulizer types and nebulization modes on bronchodilator delivery in a mechanically ventilated pediatric lung model. DESIGN: : In vitro, laboratory study. SETTING: : Research laboratory of a university hospital. INTERVENTIONS: : Using albuterol as a marker, three nebulizer types (jet nebulizer, ultrasonic nebulizer, and vibrating-mesh nebulizer) were tested in three nebulization modes in a nonhumidified bench model mimicking the ventilatory pattern of a 10-kg infant. The amounts of albuterol deposited on the inspiratory filters (inhaled drug) at the end of the endotracheal tube, on the expiratory filters, and remaining in the nebulizers or in the ventilator circuit were determined. Particle size distribution of the nebulizers was also measured. MEASUREMENTS AND MAIN RESULTS: : The inhaled drug was 2.8% ± 0.5% for the jet nebulizer, 10.5% ± 2.3% for the ultrasonic nebulizer, and 5.4% ± 2.7% for the vibrating-mesh nebulizer in intermittent nebulization during the inspiratory phase (p < 0.01). The most efficient nebulizer was the vibrating-mesh nebulizer in continuous nebulization (13.3% ± 4.6%, p < 0.01). Depending on the nebulizers, a variable but important part of albuterol was observed as remaining in the nebulizers (jet and ultrasonic nebulizers), or being expired or lost in the ventilator circuit (all nebulizers). Only small particles (range 2.39-2.70 µm) reached the end of the endotracheal tube. CONCLUSIONS: : Important differences between nebulizer types and nebulization modes were seen for albuterol deposition at the end of the endotracheal tube in an in vitro pediatric ventilator-lung model. New aerosol devices, such as ultrasonic and vibrating-mesh nebulizers, were more efficient than the jet nebulizer.

Comparison of pathologic changes in mammalian hosts infected with Schistosoma mansoni, S. japonicum and S. haematobium

The hepatic, intestinal and cardiopulmonary lesions produced by Schistosoma mansoni, S. haematobium and S. japonicum in man and experimental animals often bear striking similarities but usually have distinctive features as well. These are often related to parasitologic differences. Thus S. japonicum and S. haematobium lay their eggs in clusters which elicit the formation of large composite granulomas. The worms of these two species also tend to be sedentary, remaining in a single location for prolonged periods, thus producing large focal lesions in the intestines or urinary tract. Worm pairs of these two species also are gregarious and many worm pairs are often found in a single lesion. The size of circumoval granulomas, and the degree of fibrosis, are T cell dependent. The modulation of granuloma size is largely T cell dependent in mice infected with S. mansoni but is mostly regulated by serum factors in S. japonicum infected mice. In spite of these differences in egg laying and immunoregulation both S. mansoni and S. japonicum produce Symmers' fibrosis in the chimpanzee while S. haematobium does not, despite the presence of numerous eggs in the liver.

Comparison of two osteoporotic fracture management pathways: experience at 1 year

Introduction: Osteoporosis presenting as low-impact fractures to traumatology units is often undiagnosed and under-treated. Results from the Osteocare study in Lausanne (a nurse based intervention, passive pathway) showed that only 19% of patients received management for osteoporosis, and in the literature [1], the rate is between 10-25%. We have evaluated a different management concept, based on the systematic assessment of patients with osteoporotic fractures during and after hospitalization (active pathway). Methods: Inpatients admitted to the Department of Musculoskeletal Medicine for a fragility fracture were identified by a nurse according to a predefined questionnaire and were then clinically evaluated by a doctor. Based on the results, a management plan was proposed to the patients. Patients could choose between follow up either by their GP or by the Centre of Bone Disease of the CHUV. For patients who chose follow-up in our Centre, we assessed their adherence to medical follow-up 1 year inclusion. The results of patients who had been evaluated in our cohort between the 1 November 2008 and the 1 December 2009 were analysed. Results: 573 inpatients received specific management of their osteoporotic fracture over 18 months. The mean age was 77 y (31-99), 81% were women (203 hip fractures, 40 pelvis fractures, 101 arm fractures, 57 vertebral fractures, 63 ankle fractures, and 25 others sites). During the study period, 303 patients received a proposition of a specific treatment. 39 (13%) chose a follow up with the GP, 19 (6%) dead and 245 (81%) preferred a follow up in our Centre. After 1 year, 166 (67%) patients are under follow up in our outpatient clinic. Conclusion: With an active clinical pathway that starts during the hospitalization, consisting on a nursing evaluation followed by a medical consultation by an expert in osteoporosis, the adherence increased from 19% to 67% in terms of follow up. These results lead us to propose a consultation with a doctor experienced in osteoporosis after all osteoporotic fractures.

Enteral versus parenteral nutrition: comparison of energy metabolism in healthy subjects.

Continuous respiratory exchange measurements were performed on 10 healthy young women for 1 h before, 3 h during, and 3 h after either parenteral (iv) or intragastric (ig) administration of a nutrient mixture (52% glucose, 18% amino acid, and 30% lipid energy) infused at twice the postabsorptive resting energy expenditure (REE). REE rose from 0.98 +/- 0.02 (iv) and 0.99 +/- 0.02 kcal/min (ig) postabsorptively to 1.13 +/- 0.03 (iv) and 1.13 +/- 0.02 kcal/min (ig), resulting in nutrient-induced thermogenesis of 10 +/- 0.6 and 9.3 +/- 0.9%, respectively, when related to the metabolizable energy. The respiratory quotient rose from preinfusion values of 0.81 +/- 0.02 (iv) and 0.80 +/- 0.01 (ig) to 0.86 +/- 0.01 (iv) and 0.85 +/- 0.01 (ig). After nutrient administration the respiratory quotient fell significantly to below the preinfusion values. Plasma glucose and insulin concentrations rose during nutrient administration but were higher during the intravenous route. It is concluded that, although the response time to intragastric administration was delayed, the thermic effects and overall substrate oxidations were comparable during intravenous or intragastric administration, albeit, at lower plasma glucose and insulin concentrations via the intragastric route.