ILS-ESP An efficient, simple, and parameter-free algorithm for solving the permutation flow-shop problem

From a managerial point of view, the more effcient, simple, and parameter-free (ESP) an algorithm is, the more likely it will be used in practice for solving real-life problems. Following this principle, an ESP algorithm for solving the Permutation Flowshop Sequencing Problem (PFSP) is proposed in this article. Using an Iterated Local Search (ILS) framework, the so-called ILS-ESP algorithm is able to compete in performance with other well-known ILS-based approaches, which are considered among the most effcient algorithms for the PFSP. However, while other similar approaches still employ several parameters that can affect their performance if not properly chosen, our algorithm does not require any particular fine-tuning process since it uses basic "common sense" rules for the local search, perturbation, and acceptance criterion stages of the ILS metaheuristic. Our approach defines a new operator for the ILS perturbation process, a new acceptance criterion based on extremely simple and transparent rules, and a biased randomization process of the initial solution to randomly generate different alternative initial solutions of similar quality -which is attained by applying a biased randomization to a classical PFSP heuristic. This diversification of the initial solution aims at avoiding poorly designed starting points and, thus, allows the methodology to take advantage of current trends in parallel and distributed computing. A set of extensive tests, based on literature benchmarks, has been carried out in order to validate our algorithm and compare it against other approaches. These tests show that our parameter-free algorithm is able to compete with state-of-the-art metaheuristics for the PFSP. Also, the experiments show that, when using parallel computing, it is possible to improve the top ILS-based metaheuristic by just incorporating to it our biased randomization process with a high-quality pseudo-random number generator.

Aitchison Geometry for Probability and Likelihood as a new approach to mathematical statistics

The Aitchison vector space structure for the simplex is generalized to a Hilbert space structure A2(P) for distributions and likelihoods on arbitrary spaces. Centralnotations of statistics, such as Information or Likelihood, can be identified in the algebraical structure of A2(P) and their corresponding notions in compositional data analysis, such as Aitchison distance or centered log ratio transform.In this way very elaborated aspects of mathematical statistics can be understoodeasily in the light of a simple vector space structure and of compositional data analysis. E.g. combination of statistical information such as Bayesian updating,combination of likelihood and robust M-estimation functions are simple additions/perturbations in A2(Pprior). Weighting observations corresponds to a weightedaddition of the corresponding evidence.Likelihood based statistics for general exponential families turns out to have aparticularly easy interpretation in terms of A2(P). Regular exponential families formfinite dimensional linear subspaces of A2(P) and they correspond to finite dimensionalsubspaces formed by their posterior in the dual information space A2(Pprior).The Aitchison norm can identified with mean Fisher information. The closing constant itself is identified with a generalization of the cummulant function and shown to be Kullback Leiblers directed information. Fisher information is the local geometry of the manifold induced by the A2(P) derivative of the Kullback Leibler information and the space A2(P) can therefore be seen as the tangential geometry of statistical inference at the distribution P.The discussion of A2(P) valued random variables, such as estimation functionsor likelihoods, give a further interpretation of Fisher information as the expected squared norm of evidence and a scale free understanding of unbiased reasoning

A polynomial algorithm for special case of the one-machine scheduling problem with time-lags

The standard one-machine scheduling problem consists in schedulinga set of jobs in one machine which can handle only one job at atime, minimizing the maximum lateness. Each job is available forprocessing at its release date, requires a known processing timeand after finishing the processing, it is delivery after a certaintime. There also can exists precedence constraints between pairsof jobs, requiring that the first jobs must be completed beforethe second job can start. An extension of this problem consistsin assigning a time interval between the processing of the jobsassociated with the precedence constrains, known by finish-starttime-lags. In presence of this constraints, the problem is NP-hardeven if preemption is allowed. In this work, we consider a specialcase of the one-machine preemption scheduling problem with time-lags, where the time-lags have a chain form, and propose apolynomial algorithm to solve it. The algorithm consist in apolynomial number of calls of the preemption version of the LongestTail Heuristic. One of the applicability of the method is to obtainlower bounds for NP-hard one-machine and job-shop schedulingproblems. We present some computational results of thisapplication, followed by some conclusions.

On some geometry and equivalence classes of normal form games

Equivalence classes of normal form games are defined using the geometryof correspondences of standard equilibiurm concepts like correlated, Nash,and robust equilibrium or risk dominance and rationalizability. Resultingequivalence classes are fully characterized and compared across differentequilibrium concepts for 2 x 2 games. It is argued that the procedure canlead to broad and game-theoretically meaningful distinctions of games aswell as to alternative ways of viewing and testing equilibrium concepts.Larger games are also briefly considered.

On the efficiency and sensitivity of a pyramidal classification algorithm

In this paper we propose a Pyramidal Classification Algorithm,which together with an appropriate aggregation index producesan indexed pseudo-hierarchy (in the strict sense) withoutinversions nor crossings. The computer implementation of thealgorithm makes it possible to carry out some simulation testsby Monte Carlo methods in order to study the efficiency andsensitivity of the pyramidal methods of the Maximum, Minimumand UPGMA. The results shown in this paper may help to choosebetween the three classification methods proposed, in order toobtain the classification that best fits the original structureof the population, provided we have an a priori informationconcerning this structure.

A simple randomized algorithm for consistent sequential prediction of ergodic time series

We present a simple randomized procedure for the prediction of a binary sequence. The algorithm uses ideas from recent developments of the theory of the prediction of individual sequences. We show that if thesequence is a realization of a stationary and ergodic random process then the average number of mistakes converges, almost surely, to that of the optimum, given by the Bayes predictor.

Geometry, petrology and growth of a shallow crustal laccolith : the Torres del Paine Mafic Complex (Patagonia)

1. ABSTRACTS - RÉSUMÉSSCIENTIFIC ABSTRACT - ENGLISH VERSIONGeometry, petrology and growth of a shallow crustal laccolith: the Torres del Paine Mafi c Complex (Patagonia)The Torres del Paine intrusive complex (TPIC) is a composite mafic-granitic intrusion, ~70km2, belonging to a chain of isolated Miocene plutons in southern Patagonia. Their position is intermediate between the Mesozoic-Cenozoic calc-alkaline subduction related Patagonian batholith in the West and the late Cenozoic alkaline basaltic back-arc related plateau lavas in the East. The Torres del Paine complex formed during an important reconfiguration of the Patagonian geodynamic setting, with a migration of magmatism from the arc to the back-arc, possibly related to the Chile ridge subductionThe complex intruded the flysch of the Cretaceous Cerro Toro and Punta Barrosa Formations during the Miocene, creating a well-defined narrow contact aureole of 200-400 m width.In its eastern part, the Torres del Paine intrusive complex is a laccolith, composed of a succession of hornblende-gabbro to diorite sills at its base, with a total thickness of ~250m, showing brittle contacts with the overlying granitic sills, that form spectacular cliffs of more than 1000m. This laccolith is connected, in the western part, to its feeding system, with vertical alternating sheets of layered gabbronorite and Hbl-gabbro, surrounded and percolated by diorites. ID-TIMS U-Pb on zircons on feeder zone (FZ) gab- bros yield 12.593±0.009Ma and 12.587±0.009Ma, which is identifcal within error to the oldest granite dated so far by Michel et al. (2008). In contrast, the laccolith mafic complex is younger than than the youngest granite (12.50±0.02Ma), and has been emplaced from 12.472±0.009Ma to 12.431 ±0.006Ma, by under-accretion beneath the youngest granite at the interface with previously emplaced mafic sills.The gabbronorite crystallization sequence in the feeder zone is dominated by olivine, plagioclase, clinopyroxene and orthopyroxene, while amphibole forms late interstitial crystals. The crystallization sequence is identical in Hornblende-gabbro from the feeder zone, with higher modal hornblende. Gabbronorite and Hornblende-gabbro both display distinct Eu and Sr positive anomalies. In the laccolith, a lower Hornblende-gabbro crystallized in sills and evolved to a high alkali shoshonitic series. The Al203, Ti02, Na20, K20, Ba and Sr composition of these gabbros is highly variable and increases up to ~50wt% Si02. The lower hornblende-gabbro is characterized by kaersutite anhedral cores with inclusions of olivine, clino- and orthopyroxene and rare apatite and An70 plagioclase. Trace element modelling indicates that hornblende and clinopyroxene are in equilibrium with a liquid whose composition is similar to late basaltic trachyandesitic dikes that cut the complex. The matrix in the lower hornblende gabbro is composed of normally zoned oligoclase, Magnesio-hornblende, biotite, ilmenite and rare quartz and potassium feldspar. This assemblage crystallized in-situ from a Ba and Sr-depleted melts. In contrast, the upper Hbl-gabbro is high-K calc-alkaline. Poikilitic pargasite cores have inclusions of euhedral An70 plagioclase inclusions, and contain occasionally clinopyroxene, olivine and orthopyroxene. The matrix composition is identical to the lower hornblende-gabbro and similar to the diorite. Diorite bulk rock compositions show the same mineralogy but different modal proportions relative to hornblende-gabbrosThe Torres del Paine Intrusive Complex isotopic composition is 87Sr/86Sr=0.704, 143Nd/144Nd=0.5127, 206Pb/204Pb=18.70 and 207Pb/204Pb=15.65. Differentiated dioritic and granitic units may be linked to the gabbroic cumulates series, with 20-50% trapped interstitial melt, through fractionation of olivine-bearing gabbronorite or hornblende-gabbro fractionation The relative homogeneity of the isotopic compositions indicate that only small amounts of assimilation occurred. Two-pyroxenes thermometry, clinopyroxene barometry and amphibole-plagioclase thermometry was used to estimate pressure and temperature conditions. The early fractionation of ultramafic cumulates occurs at mid to lower crustal conditions, at temperatures exceeding 900°C. In contrast, the TPIC emplacement conditions have been estimated to ~0.7±0.5kbar and 790±60°C.Based on field and microtextural observations and geochemical modelling, fractionation of basaltic-trachyandesitic liquids at intermediate to lower crustal levels, has led to the formation of the Torres del Paine granites. Repetitive replenishment of basaltic trachy- andesitic liquid in crustal reservoirs led to mixed magmas that will ascend via the feeder zone, and crystallize into a laccolith, in the form of successive dioritic and gabbroic sills. Dynamic fractionation during emplacement concentrated hornblende rich cumulates in the center of individual sills. Variable degrees.of post-emplacement compaction led to the expulsion of felsic liquids that preferentially concentrated at the top of the sills. Incremental sills amalgamation of the entire Torres del Paine Intrusive Complex has lasted for ~160ka.RESUME SCIENTIFIQUE - VERSION FRANÇAISEGéométrie, pétrologie et croissance d'un laccolite peu profond : Le complexe ma- fique du Torres del Paine (Patagonie)Le Complexe Intrusif du Torres del Paine (CITP) est une intrusion bimodale, d'environ 70km2, appartenant à une chaîne de plutons Miocènes isolés, dans le sud de la Patago-nie. Leur position est intermédiaire entre le batholite patagonien calco-alcalin, à l'Ouest, mis en place au Mesozoïque-Cenozoïque dans un contexte de subduction, et les basal-tes andésitiques et trachybasaltes alcalins de plateau, plus jeune, à l'Est, lié à l'ouverture d'un arrière-arc.A son extrémité Est, le CITP est une succession de sills de gabbro à Hbl et de diorite, sur une épaisseur de ~250m, avec des évidences de mélange. Les contacts avec les sills de granite au-dessus, formant des parois de plus de 1000m, sont cassants. Ce laccolite est connecté, dans sa partie Ouest, à une zone d'alimentation, avec des intrusions sub-ver- ticales de gabbronorite litée et de gabbro à Hbl, en alternance. Celles-ci sont traversées et entourées par des diorites. Les zircons des gabbros de la zone d'alimentation, datés par ID-TIMS, ont cristallisés à 12.593±0.009Ma et 12.587±0.009Ma, ce qui correspond au plus vieux granite daté à ce jour par Michel et al. (2008). A l'inverse, les roches manques du laccolite se sont mises en place entre 12.472±0.009Ma et 12.431 ±0.006Ma, par sous-plaquage successifs à l'interface avec le granite le plus jeune daté à ce jour (12.50±0.02Ma).La séquence de cristallisation des gabbronorites est dominée par Ol, Plg, Cpx et Opx, alors que la Hbl est un cristal interstitiel. Elle est identique dans les gabbros à Hbl de la zone d'alimentation, avec ~30%vol de Hbl. Les gabbros de la zone d'alimentation montrent des anomalies positives en Eu et Sr distinctes. Dans le laccolite, le gabbro à Hbl inférieur évolue le long d'une série shoshonitique, riche en éléments incompatibles. Sa concentration en Al203, Ti02, Na20, K20, Ba et Sr est très variable et augmente rapide-ment jusqu'à ~50wt% Si02. Il est caractérisé par la présence de coeurs résorbés de kaer- sutite, entourés de Bt, et contenant des inclusions d'OI, Cpx et Opx, ou alors d'Ap et de rares Plg (An70). Hbl et Cpx ont cristallisés à partir d'un liquide de composition similaire aux dykes trachy-andesite basaltique du CITP. La matrice, cristallisée in-situ à partir d'un liquide pauvre en Ba et Sr, est composée d'oligoclase zoné de façon simple, de Mg-Hbl, Bt, llm ainsi que de rares Qtz et KF. Le gabbro à Hbl supérieur, quant à lui, appartient à une suite chimique calco-alcaline riche en K. Des coeurs poecilitiques de pargasite con-tiennent de nombreuses inclusions de Plg (An70) automorphe, ainsi que des Ol, Cpx et Opx. La composition de la matrice est identique à celle des gabbros à Hbl inférieurs et toutes deux sont similaires à la minéralogie des diorites. Les analyses sur roches totales de diorites montrent la même variabilité que celles de gabbros à Hbl, mais avec une ten-eur en Si02 plus élevée.La composition isotopique des liquides primitifs du CITP a été mesurée à 87Sr/86Sr=0.704, 143Nd/144Nd=0.5127, 206Pb/204Pb=18.70 et 207Pb/204Pb=15.65. Les granites et diorites différenciés peuvent être reliés à des cumulais gabbronoritiques (F=0.74 pour les granites et F=1-0.5 pour les diorites) et gabbroïques à Hbl (fractionnement supplémentaire pour les granites, avec F=0.3). La cristallisation de 20 à 50%vol de liquide interstitiel piégé dans les gabbros du CITP explique leur signature géochimique. Seules de faibles quantités de croûte continentale ont été assimilées. La température et la pression de fractionnement ont été estimées, sur la base des thermobaromètres Opx-Cpx, Hbl-Plg et Cpx, à plus de 900°C et une profondeur correspondant à la croûte inférieure-moyenne. A l'inverse, les conditions de cristallisation de la matrice des gabbros et diorites du laccolite ont été estimées à 790±60°C et ~0.7±0.5kbar.Je propose que les liquides felsiques du CITP se soient formés par cristallisation frac-tionnée en profondeur des assemblages minéralogiques observés dans les gabbros du CITP, à partir d'un liquide trachy-andesite basaltique. La percolation de magma dans les cristaux accumulés permet la remontée du mélange à travers la zone d'alimentation, vers le laccolite, où des sills se mettent en place successivement. L'amalgamation de sills dans le CITP a duré ~160ka.Le CITP s'est formé durant une reconfiguration importante du contexte géodynamique en Patagonie, avec un changement du magmatisme d'arc vers un volcanisme d'arrière- arc. Ce changement est certainement lié à la subduction de la ride du Chili.RESUME GRAND PUBLIC - VERSION FRANÇAISEGéométrie, pétrologie et croissance d'une chambre magmatique peu profonde : Le complexe mafique du Torres del Paine (Patagonie)Le pourtour de l'Océan Pacifique est caractérisé par une zone de convergence de plaques tectoniques, appelée zone de subduction, avec le plongement de croûte océa-nique sous les Andes dans le cas de la Patagonie. De nombreux volcans y sont associés, formant la ceinture de feu. Mais seuls quelques pourcents de tout le magma traversant la croûte terrestre parviennent à la surface et la majeure partie cristallise en profondeur, dans des chambres magmatiques. Quelles est leur forme, croissance, cristallisation et durée de vie ? Le complexe magmatique du Torres del Paine représente l'un des meilleurs endroits au monde pour répondre à ces questions. Il se situe au sud de la Patagonie, formant un massif de 70km2. Des réponses peuvent être trouvées à différentes échelles, variant de la montagne à des minéraux de quelques 1000ème de millimètres.Il est possible de distinguer trois types de roches : des gabbros et des diorites sur une épaisseur de 250m, surmontées par des parois de granite de plus de 1000m. Les contacts entre ces roches sont tous horizontaux. Entre granites et gabbro-diorite, le contact est net, indiquant que le second magma s'est mis en place au contact avec un magma plus ancien, totalement solidifié. Entre gabbros et diorites, les contacts sont diffus, souvent non-linéaires, indiquant à l'inverse la mise en contact de magmas encore partiellement liquides. Dans la partie Ouest de cette chambre magmatique, les contacts entre roches sont verticaux. Il s'agit certainement du lieu de remplissage de la chambre magmatique.Lors du refroidissement d'un magma, différents cristaux vont se former. Leur stabilité et leur composition varient en fonction de la pression, de la température ou de la chimie du magma. La séquence de cristallisation peut être définie sur la base d'observations microscopiques et de la composition chimique des minéraux. Différents gabbros sont ainsi distingués : le gabbro à la base est riche en hornblende, d'une taille de ~5mm, sans inclusion de plagioclase mais avec des cristaux d'olivine, clinopyroxene et orthopyroxene inclus ; le gabbro supérieur est lui-aussi riche en hornblende (~5mm), avec les mêmes inclusions additionnées de plagioclase. Ces cristaux se sont formés à une température supérieure à 900°C et une profondeur correspondant à la croûte moyenne ou inférieure. Les minéraux plus fin, se trouvant hors des cristaux de hornblende des deux gabbros, sont similaires à ceux des diorites : plagioclase, biotite, hornblende, apatite, quartz et feldspath alcalin. Ces minéraux sont caractéristiques des granites. Ils ont cristallisé à ~790°C et ~2km de profondeur.La cristallisation des minéraux et leur extraction du magma par gravité provoque un changement progressif de la composition de ce dernier. Ainsi, après extraction d'olivine et d'orthopyroxene riches en Mg, de clinopyroxene riche en Ca, de plagioclase riche en Ca et Al et d'hornblende riche en Ca, Al et Mg, le liquide final sera appauvri en ces élé-ments. Un lien peut ainsi être proposé entre les diorites dont la composition est proche du liquide de départ, les granites dont la composition est similaire au liquide final, et les gabbros dont la minéralogie correspond aux minéraux extraits.L'utilisation de zircons, un minéral riche en U dont les atomes se transforment en Pb par décomposition radioactive au cours de millions d'années, permet de dater le refroidissement des roches qui les contiennent. Ainsi, il a été observé que les roches de la zone d'alimentation, à l'Ouest du complexe magmatique, ont cristallisés il y a 12.59±0.01 Ma, en même temps que les granites les plus vieux, se trouvant au sommet de la chambre magmatique, datés par Michel et al. (2008). Les deux roches pourraient donc avoir la même origine. A l'inverse, les gabbros et diorites de la chambre magmatique ont cristallisé entre 12.47±0.01Ma et 12.43±0.01Ma, les roches les plus vieilles étant à la base.En comparant la composition des roches du Torres del Paine avec celles d'autres en-tités géologiques de Patagonie, les causes du magmatisme peuvent être recherchées. A l'Ouest, on trouve en effet des intrusions granitiques, plus anciennes, caractéristiques de zones de convergence de plaque tectonique, alors qu'à l'Est, des laves basaltiques plus jeunes sont caractéristiques d'une dynamique d'extension. Sur la base des compositions chimiques des roches de ces différentes entités, l'évolution progressive de l'une à l'autre a pu être démontrée. Elle est certainement due à l'arrivée d'une dorsale océanique (zone d'extension crustale et de création de croûte océanique par la remontée de magma) dans la zone de subduction, le long des Andes.Je propose que, dans un premier temps, des magmas granitiques sont remontés dans la chambre magmatique, laissant d'importants volumes de cristaux dans la croûte pro-fonde. Dans un second épisode, les cristaux formés en profondeur ont été transportés à travers la croûte continentale, suite au mélange avec un nouveau magma injecté. Ces magmas chargés de cristaux ont traversé la zone d'alimentation avant de s'injecter dans la chambre magmatique. Différents puises ont été distingués, injectés dans la chambre magmatique du sommet à la base concernant les granites, puis à la base du granite le plus jeune pour les gabbros et diorites. Le complexe magmatique du Torres del Paine s'est construit sur une période totale de 160'000±20'000 ans.

Latitudinal gradient effect on the wing geometry of Auca coctei (Guérin)(Lepidoptera, Nymphalidae)

Latitudinal gradient effect on the wing geometry of Auca coctei (Guérin) (Lepidoptera, Nymphalidae). When the environmental conditions change locally, the organisms and populations may also change in response to the selection pressure, so that the development of individuals may become affected in different degrees. There have been only a few studies in which the patterns of wing morphology variation have been looked into along a latitudinal gradient by means of geometric morphometrics. The aim of this work was to assess the morphologic differentiation of wing among butterfly populations of the species Auca coctei. For this purpose, 9 sampling locations were used which are representative of the distribution range of the butterfly and cover a wide latitudinal range in Chile. The wing morphology was studied in a total of 202 specimens of A. coctei (150 males and 52 females), based on digitization of 17 morphologic landmarks. The results show variation of wing shape in both sexes; however, for the centroid size there was significant variation only in females. Females show smaller centroid size at higher latitudes, therefore in this study the Bergmann reverse rule is confirmed for females of A. coctei. Our study extends morphologic projections with latitude, suggesting that wing variation is an environmental response from diverse origins and may influence different characteristics of the life history of a butterfly.

MSISpIC: A Probabilistic Scan Matching Algorithm Using a Mechanical Scanned Imaging Sonar

This paper compares two well known scan matching algorithms: the MbICP and the pIC. As a result of the study, it is proposed the MSISpIC, a probabilistic scan matching algorithm for the localization of an Autonomous Underwater Vehicle (AUV). The technique uses range scans gathered with a Mechanical Scanning Imaging Sonar (MSIS), and the robot displacement estimated through dead-reckoning with the help of a Doppler Velocity Log (DVL) and a Motion Reference Unit (MRU). The proposed method is an extension of the pIC algorithm. Its major contribution consists in: 1) using an EKF to estimate the local path traveled by the robot while grabbing the scan as well as its uncertainty and 2) proposing a method to group into a unique scan, with a convenient uncertainty model, all the data grabbed along the path described by the robot. The algorithm has been tested on an AUV guided along a 600m path within a marina environment with satisfactory results

An Efficient Nominal Unification Algorithm

Nominal Unification is an extension of first-order unification where terms can contain binders and unification is performed modulo α equivalence. Here we prove that the existence of nominal unifiers can be decided in quadratic time. First, we linearly-reduce nominal unification problems to a sequence of freshness and equalities between atoms, modulo a permutation, using ideas as Paterson and Wegman for first-order unification. Second, we prove that solvability of these reduced problems may be checked in quadràtic time. Finally, we point out how using ideas of Brown and Tarjan for unbalanced merging, we could solve these reduced problems more efficiently

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.

Validation of a clinical algorithm to identify low-risk patients with pulmonary embolism.

Summary Background: We previously derived a clinical prognostic algorithm to identify patients with pulmonary embolism (PE) who are at low-risk of short-term mortality who could be safely discharged early or treated entirely in an outpatient setting. Objectives: To externally validate the clinical prognostic algorithm in an independent patient sample. Methods: We validated the algorithm in 983 consecutive patients prospectively diagnosed with PE at an emergency department of a university hospital. Patients with none of the algorithm's 10 prognostic variables (age >/= 70 years, cancer, heart failure, chronic lung disease, chronic renal disease, cerebrovascular disease, pulse >/= 110/min., systolic blood pressure < 100 mm Hg, oxygen saturation < 90%, and altered mental status) at baseline were defined as low-risk. We compared 30-day overall mortality among low-risk patients based on the algorithm between the validation and the original derivation sample. We also assessed the rate of PE-related and bleeding-related mortality among low-risk patients. Results: Overall, the algorithm classified 16.3% of patients with PE as low-risk. Mortality at 30 days was 1.9% among low-risk patients and did not differ between the validation and the original derivation sample. Among low-risk patients, only 0.6% died from definite or possible PE, and 0% died from bleeding. Conclusions: This study validates an easy-to-use, clinical prognostic algorithm for PE that accurately identifies patients with PE who are at low-risk of short-term mortality. Low-risk patients based on our algorithm are potential candidates for less costly outpatient treatment.

A fast Bayesian reconstruction algorithm for emission tomography with entropy prior converging to feasible images

The development and tests of an iterative reconstruction algorithm for emission tomography based on Bayesian statistical concepts are described. The algorithm uses the entropy of the generated image as a prior distribution, can be accelerated by the choice of an exponent, and converges uniformly to feasible images by the choice of one adjustable parameter. A feasible image has been defined as one that is consistent with the initial data (i.e. it is an image that, if truly a source of radiation in a patient, could have generated the initial data by the Poisson process that governs radioactive disintegration). The fundamental ideas of Bayesian reconstruction are discussed, along with the use of an entropy prior with an adjustable contrast parameter, the use of likelihood with data increment parameters as conditional probability, and the development of the new fast maximum a posteriori with entropy (FMAPE) Algorithm by the successive substitution method. It is shown that in the maximum likelihood estimator (MLE) and FMAPE algorithms, the only correct choice of initial image for the iterative procedure in the absence of a priori knowledge about the image configuration is a uniform field.