Variance reduction methods for simulation of densities on Wiener space

We develop a general error analysis framework for the Monte Carlo simulationof densities for functionals in Wiener space. We also study variancereduction methods with the help of Malliavin derivatives. For this, wegive some general heuristic principles which are applied to diffusionprocesses. A comparison with kernel density estimates is made.

Tax simulation in the SHP

This paper describes a simulation package designed to estimate the annual income taxes paid by respondents of the Swiss Household Panel (SHP). In Switzerland, the 26 cantons have their own tax system. Additionally, tax levels vary between the over 2000 municipalities and over time. The simulation package takes account of this complexity by building on existing tables on tax levels which are provided by the Swiss Federal Tax Administration Office. Because these are limited to a few types of households and only 812 municipalities, they have to be extended to cover all households and municipalities. A further drawback of these tables is that they neglect several deductions. The tax simulation package fills this gap by taking additionally account of deductions for children, double-earner couples, third pillar and support for dependent persons according to cantonal legislation. The resulting variable on direct taxes not only serves to calculate household income net of taxes, but can also be a variable for analysis by its own account.

2008 Annual Report of the Iowa Railway Finance Authority

The Iowa Railway Finance Authority (IRFA) was created in 1980 by the 68th General Assembly to provide for the financing of rail facilities and to enhance and continue the operation of essential rail facilities. IRFA was authorized to offer financial assistance for the acquisition, rehabilitation, construction, refinancing, extension, replacement, maintenance, repair, or leasing of any rail facility. Rail Revolving Loan and Grant Program The 2005 legislative session amended Iowa Code 327H.20 by assigning all repayments of IRFA and other Iowa Department of Transportation (Iowa DOT) rail assistance loans to the Rail Revolving Loan and Grant Program (RRLGP). In addition, the following annual appropriations were allocated to the fund: • FY 2007: $235,000 • FY 2008: $2 million • FY 2009: $2 million Since the creation of the RRLGP in 2005, IRFA has awarded funding to 23 projects totaling over $7.2 million in grants and loans in its three rounds of competitive funding. These projects have pledged to create 1,672 jobs within two years of project completion and 1,361 jobs have been retained. Funded projects are associated with over $2 billion in total private capital investment. In 2008, the IRFA Board directed all available funds be used to help repair railroads devastated by summer flooding. On July 31, $3.9 million in deferred payment loans were offered to seven Iowa based railroads to allow for immediate repair of rail beds, including the cost of materials, labor and equipment.

An experimental test of rainfall as a control agent of Glycaspis brimblecombei Moore (Hemiptera, Psyllidae) on seedlings of Eucalyptus camaldulensis Dehn (Myrtaceae)

An experimental test of rainfall as a control agent of Glycaspis brimblecombei Moore (Hemiptera, Psyllidae) on seedlings of Eucalyptus camaldulensis Dehn (Myrtaceae). Glycaspis brimblecombei is one the greatest threats to eucalyptus plantations in Brazil. The effects of rainfall to reduce the abundance of lerp of Glycaspis brimblecombei on experimentally infested seedlings of Eucalyptus camaldulensis were assessed. The number of lerps on the adaxial and abaxial surfaces of every leaf of 60 seedlings was recorded, before and after submission to the following treatments: "artificial rain", "leaf wetting" and control. A drastic reduction in lerp abundance per plant was observed after the treatments "leaf wetting" and artificial rain (F = 53.630; p < 0.001), whereas lerp abundance remained roughly constant in the control treatment along the experiment (F = 1.450; p = 0.232). At the end of the experiment, lerp abundance was significantly lower in both the "artificial rain" and "leaf wetting" than in the control treatment. Two days of rainfall simulation were sufficient to decrease more than 50% of the lerp population, with almost 100% of effectiveness after 5 days of experiment. Our results indicate that lerp solubilization and mechanical removal by water are potential tools to the population regulation of G. brimblecombei on E. camaldulensis seedlings.

PENELOPE-2008: A Code System for Monte Carlo Simulation of Electron and Photon Transport

The computer code system PENELOPE (version 2008) performs Monte Carlo simulation of coupledelectron-photon transport in arbitrary materials for a wide energy range, from a few hundred eV toabout 1 GeV. Photon transport is simulated by means of the standard, detailed simulation scheme.Electron and positron histories are generated on the basis of a mixed procedure, which combinesdetailed simulation of hard events with condensed simulation of soft interactions. A geometry packagecalled PENGEOM permits the generation of random electron-photon showers in material systemsconsisting of homogeneous bodies limited by quadric surfaces, i.e., planes, spheres, cylinders, etc. Thisreport is intended not only to serve as a manual of the PENELOPE code system, but also to provide theuser with the necessary information to understand the details of the Monte Carlo algorithm.

Informing the judgments of fingerprint analysts using quality metric and statistical assessment tools

The aim of this research was to evaluate how fingerprint analysts would incorporate information from newly developed tools into their decision making processes. Specifically, we assessed effects using the following: (1) a quality tool to aid in the assessment of the clarity of the friction ridge details, (2) a statistical tool to provide likelihood ratios representing the strength of the corresponding features between compared fingerprints, and (3) consensus information from a group of trained fingerprint experts. The measured variables for the effect on examiner performance were the accuracy and reproducibility of the conclusions against the ground truth (including the impact on error rates) and the analyst accuracy and variation for feature selection and comparison.¦The results showed that participants using the consensus information from other fingerprint experts demonstrated more consistency and accuracy in minutiae selection. They also demonstrated higher accuracy, sensitivity, and specificity in the decisions reported. The quality tool also affected minutiae selection (which, in turn, had limited influence on the reported decisions); the statistical tool did not appear to influence the reported decisions.

Direct simulation of interface dynamics: linking capillary pressure, interfacial area and surface energy

We perform direct numerical simulations of drainage by solving Navier- Stokes equations in the pore space and employing the Volume Of Fluid (VOF) method to track the evolution of the fluid-fluid interface. After demonstrating that the method is able to deal with large viscosity contrasts and to model the transition from stable flow to viscous fingering, we focus on the definition of macroscopic capillary pressure. When the fluids are at rest, the difference between inlet and outlet pressures and the difference between the intrinsic phase average pressure coincide with the capillary pressure. However, when the fluids are in motion these quantities are dominated by viscous forces. In this case, only a definition based on the variation of the interfacial energy provides an accurate measure of the macroscopic capillary pressure and allows separating the viscous from the capillary pressure components.

Nuclear well logging in soils. Calibration of slim hole nuclear tools gamma-gamma and neutron-thermal neutron.

The calculation of elasticity parameters by sonic and ultra sonic wave propagation in saturated soils using Biot's theory needs the following variables : forpiation density and porosity (p, ø), compressional and shear wave velocities (Vp, Vs), fluid density, viscosity and compressibility (Pfi Ilfi Ki), matrix density and compressibility (p" K), The first four parameters can be determined in situ using logging probes. Because fluid and matrix characteristics are not modified during core extraction, they can be obtained through laboratory measurements. All parameters necessitate precise calibrations in various environments and for specific range of values encountered in soils. The slim diameter of boreholes in shallow geophysics and the high cost of petroleum equipment demand the use of specific probes, which usually only give qualitative results. The measurement 'of density is done with a gamma-gamma probe and the measurement of hydrogen index, in relation to porosity, by a neutron probe. The first step of this work has been carried out in synthetic formations in the laboratory using homogeneous media of known density and porosity. To establish borehole corrections different casings have been used. Finally a comparison between laboratory and in situ data in cored holes of known geometry and casing has been performed.

Brownian dynamics simulation of DNA condensation.

DNA condensation observed in vitro with the addition of polyvalent counterions is due to intermolecular attractive forces. We introduce a quantitative model of these forces in a Brownian dynamics simulation in addition to a standard mean-field Poisson-Boltzmann repulsion. The comparison of a theoretical value of the effective diameter calculated from the second virial coefficient in cylindrical geometry with some experimental results allows a quantitative evaluation of the one-parameter attractive potential. We show afterward that with a sufficient concentration of divalent salt (typically approximately 20 mM MgCl(2)), supercoiled DNA adopts a collapsed form where opposing segments of interwound regions present zones of lateral contact. However, under the same conditions the same plasmid without torsional stress does not collapse. The condensed molecules present coexisting open and collapsed plectonemic regions. Furthermore, simulations show that circular DNA in 50% methanol solutions with 20 mM MgCl(2) aggregates without the requirement of torsional energy. This confirms known experimental results. Finally, a simulated DNA molecule confined in a box of variable size also presents some local collapsed zones in 20 mM MgCl(2) above a critical concentration of the DNA. Conformational entropy reduction obtained either by supercoiling or by confinement seems thus to play a crucial role in all forms of condensation of DNA.

A game theoretic simulation model for quality oriented timber supply to sawmills

Summary

Estudi i guia docent per a modelatge i simulació de sistemes mitjançant l'entron ARENA

Aquest projecte es basa en el modelatge i simulació de sistemes utilitzant un simulador digital, i pretén ser una guia docent com a eina d’ajuda per a una assignatura que, a priori, s’impartirà a la Universitat de Vic. La simulació és una tècnica que permet representar el comportament de processos (físics, productius, de serveis, etc.) sense necessitat d’accedir al sistema real. Per analitzar, estudiar i millorar el comportament d’un sistema mitjançant la tècnica de la simulació digital és necessari primer desenvolupar un model conceptual que descrigui les variables d’interès, i després implementar-lo en un simulador per poder analitzar els resultats. ARENA és el software de simulació que s’estudia en aquest projecte i es presenta com una eina que permet la descripció complerta de l’experiència que una entitat desenvolupa a l’interior del sistema mentre flueix a través d’aquest. En concret s’utilitza la versió ARENA 10.0. Pel que fa a l’estructura del projecte, primerament s’introdueixen conceptes teòrics referents a la simulació, així com avantatges i inconvenients i els camps d’aplicació de la simulació. Seguidament i ja centrats en l’Arena, s’analitza un exemple senzill per començar-ne a veure el funcionament. Posteriorment, es van estudiant varis exemples amb complexitat progressiva. Aquests exemples es desenvolupen pas a pas de manera que es puguin anar provant amb el simulador. En el transcurs del projecte es van estudiant les eines de l’Arena i les seves possibilitats, així com els resultats obtinguts i les interpretacions d’aquests. Aquest projecte pretén, doncs, donar conceptes introductoris en el camp de la simulació en general, i, en particular, descriure eines bàsiques sobre el funcionament de l’Arena.

Projectes de modelatge i simulació avançada de sistemes mitjançant l'entorn ARENA

RESUM En l’exercici de les funcions típiques de qualsevol activitat humana, l’home ha de prendre decisions d’un o altre tipus permanentment i aquesta situació porta afegit, per lo general risc i incertesa, el que compromet la qualitat i l’èxit de la decisió. Per ajudar i contrarestar aquesta situació, l’home ha desenvolupat a través del temps una diversitat d’eines que li permeten minimitzar el risc i la incertesa en la presa de decisions. La simulació és una d’aquestes eines. Mitjançant la seva aplicació no només s’aconsegueix la comesa anterior, sinó que també es minimitzen els costos involucrats en la decisió mitjançant un millor ús dels recursos, la disminució del temps emprat per a la seva realització i la minimització de les probabilitats de risc. Aquest projecte representa la continuació d’un primer estudi introductori de la simulació anomenat: Estudi i guia docent per a modelatge i simulació de sistemes mitjançant l’entorn ARENA [9]. Realitzat per l’alumna d’Enginyeria en Organització Industrial de la Universitat de Vic, Montse Carbonell Crosas, l’any 2008 i sota la codirecció del director d’aquest segon projecte, el professor Juli Ordeix Rigo. Aquest nou projecte s’inicia amb una primera part teòrica, continguda dins del primer volum, la qual reforça els conceptes teòrics referents a la simulació amb ARENA, ja vistos en l’anterior projecte. Complementant aquells considerats bàsics i els de més utilitat i finalitza introduint nous conceptes avançats. Els nous capítols de temàtica avançada, junt als primers més bàsics de la primera part són exercitats dins de la segona, formant el segon volum d’aquest projecte. El mateix requereix la participació activa de l’alumne, per tal de realitzar cadascun dels 89 exercicis pràctics que es plantegen i poder onsolidar l’aprenentatge teòric d’aquesta eina avançada de simulació fent consultes als apartats teòrics recomanats dins de cada exercici. La complexitat dels exercicis anirà augmentant gradualment i s’insisteix en seguir la metodologia presentada en el projecte per a realitzar-los tots de forma ordenada i ascendent. L’ alumne quan acabi la part pràctica, haurà consolidat tota la part teòrica i serà capaç d’exercir com analista per tal de generar els seus propis projecte de simulació.