Inference of distributional parameters from compositional samples containing nondetects

Low concentrations of elements in geochemical analyses have the peculiarity of being compositional data and, for a given level of significance, are likely to be beyond the capabilities of laboratories to distinguish between minute concentrations and complete absence, thus preventing laboratories from reporting extremely low concentrations of the analyte. Instead, what is reported is the detection limit, which is the minimum concentration that conclusively differentiates between presence and absence of the element. A spatially distributed exhaustive sample is employed in this study to generate unbiased sub-samples, which are further censored to observe the effect that different detection limits and sample sizes have on the inference of population distributions starting from geochemical analyses having specimens below detection limit (nondetects). The isometric logratio transformation is used to convert the compositional data in the simplex to samples in real space, thus allowing the practitioner to properly borrow from the large source of statistical techniques valid only in real space. The bootstrap method is used to numerically investigate the reliability of inferring several distributional parameters employing different forms of imputation for the censored data. The case study illustrates that, in general, best results are obtained when imputations are made using the distribution best fitting the readings above detection limit and exposes the problems of other more widely used practices. When the sample is spatially correlated, it is necessary to combine the bootstrap with stochastic simulation

Prediction of glucose excursions under uncertain parameters and food intake in intensive insulin therapy for type 1 diabetes mellitus

Considering the difficulty in the insulin dosage selection and the problem of hyper- and hypoglycaemia episodes in type 1 diabetes, dosage-aid systems appear as tremendously helpful for these patients. A model-based approach to this problem must unavoidably consider uncertainty sources such as the large intra-patient variability and food intake. This work addresses the prediction of glycaemia for a given insulin therapy face to parametric and input uncertainty, by means of modal interval analysis. As result, a band containing all possible glucose excursions suffered by the patient for the given uncertainty is obtained. From it, a safer prediction of possible hyper- and hypoglycaemia episodes can be calculated

Estructura i dinàmica del sistema bentònic en llacunes costaneres dels aiguamolls de l'Empordà

S'ha estudiat l'estructura i dinàmica del sistema bentònic de llacunes costaneres de la maresma dels aiguamolls de l'Empordà (NE de la península Ibèrica) i els factors que les determinen. Amb aquesta finalitat es van prendre mostres d'organismes del bentos (mensualment), nutrients i pigments de l'aigua (setmanalment) i el sediment (mensualment), així com paràmetres físics de l'aigua (setmanalment), durant dos cicles d'inundació (1997-1998 i 1998-1999). En aquesta maresma els factors determinants per als organismes bentònics i el contingut en nutrients del sediment són principalment físics: el grau de permanència de l'aigua, el confinament i la granulometria. De manera que els factors tròfics tenen menys pes que els físics tal com es posa de manifest amb les associacions de nematodes. El règim hídric de les llacunes integra des del grau de permanència de l'aigua fins a la salinitat i el confinament, per la qual cosa a partir de la caracterització dels diferents tipus de règim hídric es poden discriminar les tres tipologies d'ambients amb comunitats bentòniques diferenciades en l'estructura, la composició i la dinàmica temporal. Així, les tres comunitats diferenciades són: (1) la d'aigües permanents caracteritzada per una riquesa i diversitat elevades i constants en el temps, i per un nombre més elevat de tàxons sense adaptacions per passar situacions adverses; (2) la d'aigües semipermanents i temporànies d'inundació pulsativa, caracteritzada per una menor riquesa i diversitat, i per la dominància d'organismes amb estructures de resistència per passar la fase seca i aguantar condicions desfavorables, i (3) la d'aigües temporànies d'inundació continuada, caracteritzada per una major variabilitat de la riquesa i diversitat al llarg del cicle i per una singularitat taxonòmica més elevada, ja que s'observen cladòcers i una major representació taxonòmica dels insectes. Durant la realització d'aquest estudi els dos cicles d'inundació van suposar entrades d'aigua en el sistema diferents, per això es va diferenciar un cicle humit (1997-1998) i un de sec (1998-1999). En relació amb la major importància de les entrades d'aigua en el sistema, es va observar un contingut en nutrients de l'aigua significativament superior durant el cicle humit. Pel que fa al sediment, la seva major inèrcia explicaria la manca de diferències significatives en el contingut en nutrients del sediment del cicle humit i del sec. Ara bé, sí que s'observa diferència en la composició d'organismes. Així, la variabilitat interanual en l'estructura de la comunitat bentònica de la maresma es pot atribuir a una situació d'estrès provocada pel cicle sec. Aquesta situació va comportar un augment de la diversitat per disminució de la dominància a totes les llacunes estudiades, determinada, en gran manera, per la disminució de l'abundància original de poblacions dominants i característiques d'ambients permanents, Corophium orientale, o temporànies, Gammarus aequicauda. No existeix un patró estacional en les comunitats estudiades, sinó que els principals canvis en el plàncton i el bentos són deguts a pertorbacions hídriques. Aquestes pertorbacions provoquen que s'estableixi una situació en què dominen uns tàxons característics en el plàncton (Synchaeta i Eutintinnus) i, en canvi, que en el bentos es perdi l'estructura pròpia de les comunitats característiques dels diferents ambients (comunitats dominades per Corophium orinetale en aigües permanents i per Chironomus salinarius en temporànies). La severitat és diferent en funció del grau de permanència de l'aigua, i és més alta en ambients permanents i més baixa en els temporanis. La major severitat en ambients permanents ve determinada per la menor adaptació del bentos d'aigües permanents a les fluctuacions ambientals i a l'efecte de la predació sobre el plàncton. A més, tant en el bentos com en el plàncton es donen fenòmens d'histèresi, segons els quals l'efecte de la pertorbació no solament depèn de la causa pertorbadora sinó també de l'estat en què es troba la comunitat en el moment de la pertorbació.

Application of modal interval analysis to the simulation of the behaviour of dynamic systems with uncertain parameters

Els models matemàtics quantitatius són simplificacions de la realitat i per tant el comportament obtingut per simulació d'aquests models difereix dels reals. L'ús de models quantitatius complexes no és una solució perquè en la majoria dels casos hi ha alguna incertesa en el sistema real que no pot ser representada amb aquests models. Una forma de representar aquesta incertesa és mitjançant models qualitatius o semiqualitatius. Un model d'aquest tipus de fet representa un conjunt de models. La simulació del comportament de models quantitatius genera una trajectòria en el temps per a cada variable de sortida. Aquest no pot ser el resultat de la simulació d'un conjunt de models. Una forma de representar el comportament en aquest cas és mitjançant envolupants. L'envolupant exacta és complete, és a dir, inclou tots els possibles comportaments del model, i correcta, és a dir, tots els punts dins de l'envolupant pertanyen a la sortida de, com a mínim, una instància del model. La generació d'una envolupant així normalment és una tasca molt dura que es pot abordar, per exemple, mitjançant algorismes d'optimització global o comprovació de consistència. Per aquesta raó, en molts casos s'obtenen aproximacions a l'envolupant exacta. Una aproximació completa però no correcta a l'envolupant exacta és una envolupant sobredimensionada, mentre que una envolupant correcta però no completa és subdimensionada. Aquestes propietats s'han estudiat per diferents simuladors per a sistemes incerts.