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

Avaluació del procés de desintegració de papers recuperats

La desintegració és una etapa important en la recuperació de paper vell, ja que té importants conseqüències en consum d'energia i en el comportament de les etapes posteriors. Per això els objectius es centren en analitzar la desintegració des del punt de vista del temps de desintegració, els aspectes energètics, modelització de la màquina de desintegració utilitzada i anàlisi dels factors de cisallament calculats com a mesura global de les forces implicades en la desintegració. Els autors que hi han treballat donen diferents explicacions a aquestes forces. Fins avui només s'ha pogut avaluar qualitativament la influència que tenen cada un dels mecanismes en el temps necessari per a desintegrar i en el consum energètic. Les característiques reològiques de les suspensions papereres, i el seu comportament no newtonià tenen una clara influència en el consum energètic i les forces de desfibrat en el desintegrador. Els experiments de desintegració s'han realitzat en un púlper convencional, amb tres tipus de paper recuperat: paper estucat d'alta qualitat imprès offset (PQ), paper revista estucat imprès en color (PR), paper blanc imprès en impresora làsser (PF). Anàlisi del temps de desintegració Per cada un del papers estudiats (PQ, PR i PF), les fraccions màssiques des de 0.06 fins a la màxima que estat possible per cada paper (de 0.14 a 0.18), i a dues velocitats d'agitació diferents, s'ha determinat el temps de desintegració (tD) fins a aconseguir un índex de Sommerville de 0.01%. S'obté que en augmentar la fracció màssica disminueix potencialment el temps de desintegració. S'ha estudiat la velocitat de desintegració, la producció teòrica del púlper en cada cas, i la seva relació amb les forces d'impacte i de fregament que produeixen la desintegració. Aspectes energètics El consum específic d'energia (SEC), definit com l'energia consumida per a desintegrar 1 kg de paper recuperat, disminueix molt en augmentar Xm, ja que a més de disminuir l'energia consumida en cada desintegració, el contingut en paper és més elevat. Pel disseny de desintegradors, cal tenir en compte que en augmentar Xm i en augmentar la velocitat, sempre augmenta la potència consumida. Però així com els beneficis de treballar a Xm alt són de 10 vegades en termes de SEC i de producció, l'augment de potència és només de l'ordre de 2 vegades la necessària respecte de la Xm baixa. Viscositat aparent i energia de fluidització S'estudia la relació entre el temps de desintegració, les forces de fregament i els valors de viscositat aparent de la bibliografia. Per cada paper i velocitat s'ha observat que el consum específic d'energia disminueix en funció de la viscositat aparent. Reologia del púlper Utilitzant el mètode de Metzner i Otto (1957) per determinar la viscositat aparent mitjana de les suspensions papereres, modificat per Roustan, s'ha caracteritzat el pulper mitjançant el model: Np= K· Rex·Fry S'han utilitzat dissolucions de glicerina com a fluid newtonià per a calcular les constants d'ajust, i a partir d'aquí, aïllar la viscositat aparent en funció de la potència neta i els paràmetres d'agitació. La viscositat aparent, d'acord amb Fabry (1999) es substitueix pel concepte de factor de cisallament. Factor de cisallament Calculat el factor de cisallament per a cada paperot i condicions d'agitació, s'ha relacionat amb Xm, SEC, tD, consum de potència, potència instal·lada i fracció cel·lulòsica. El factor de cisallament és un paràmetre útil per a quantificar les forces globals implicades en la desintegració.

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.