Triceps and subscapular skinfold thicknesses percentiles and cut-offs for overweight and obesity in a population-based sample of schoolchildren in Bogota, Colombia

OBJECTIVES: The aims of this study were to establish a Colombian smoothed centile charts and LMS tables for tríceps, subscapular and sum tríceps+subscapular skinfolds; appropriate cut-offs were selected using receiver operating characteristic analysis based in a populationbased sample of schoolchildren in Bogota, Colombia and to compare them with international studies. METHODS: A total of 9 618 children and adolescents attending public schools in Bogota, Colombia (55.7% girls; age range of 9–17.9 years). Height, weight, body mass index (BMI), waist circumference, triceps and subscapular skinfold measurements were obtained using standardized methods. We have calculated tríceps+subscapular skinfold (T+SS) sum. Smoothed percentile curves for triceps and subscapular skinfold thickness were derived by the LMS method. Receiver operating characteristics curve (ROC) analyses were used to evaluate the optimal cut-off point of tríceps, subscapular and sum tríceps+subscapular skinfolds for overweight and obesity based on the International Obesity Task Force (IOTF) definitions. Data were compared with international studies. RESULTS: Subscapular, triceps skinfolds and T+SS were significantly higher in girls than in boys (P <0.001). The median values for triceps, subscapular as well as T+SS skinfold thickness increased in a sex-specific pattern with age. The ROC analysis showed that subscapular, triceps skinfolds and T+SS have a high discrimination power in the identification of overweight and obesity in the sample population in this study. Based on the raw non-adjusted data, we found that Colombian boys and girls had high triceps and subscapular skinfolds values than their counterparts from Spain, UK, German and US. CONCLUSIONS: Our results provide sex- and age-specific normative reference standards for the triceps and subscapular skinfold thickness values in a large, population-based sample of 3 schoolchildren and adolescents from an Latin-American population. By providing LMS tables for Latin-American people based on Colombian reference data, we hope to provide quantitative tools for the study of obesity and its complications.

Valores de referencia de la fuerza prensil en escolares colombianos. Estudio FUPRECOL

The primary aim of this study was to generate normative handgrip strength (HG) data for 10- to 17.9-year-olds. The secondary aim was to determine the relative proportion of Colombian children and adolescents that fall into established Health Benefit Zones (HBZ). This cross-sectional study is enrolling 7268 schoolchildren (boys n=3129 and girls n=4139, age 12.7 (2.4) years old. HG was measured using a hand dynamometer with an adjustable grip. Five HBZs (Needs Improvement, Fair, Good, Very Good, and Excellent) have been established that correspond to combined-HG. Centile smoothed curves, percentile and tables for the 3rd, 10th, 25th, 50th, 75th, 90th and 97th percentile were calculated using Cole’s LMS method. HG peaked in the sample at 22.2 (8.9) kg in boys and 18.5 (5.5) kg in girls. The increase in HG was greater for boys than for girls, but the peak HG was lower in girls than in boys. The HBZ data indicated that a higher overall percentage of boys than girls at each age group fell into the “Needs Improvement” zone, with differences particularly pronounced during adolescence. Our results provide, for the first time, sex- and age-specific HG reference standards for Colombian schoolchildren aged 9-17.9 years.

Normative reference of standing long jump for Colombian schoolchildren aged 9-17.9 years: The FUPRECOL study

El propósito del presente estudio era generar los valores normativos de salto largo para niños de 9-17.9 años, e investigar las diferencias de sexo y grupo de edad

SPHERA, a new Italian convection-permitting regional reanalysis: improving heavy-rainfall representation and describing hail-prone environments

Extreme weather events related to deep convection are high-impact critical phenomena whose reliable numerical simulation is still challenging. High-resolution (convection-permitting) modeling setups allow to switch off physical parameterizations accountable for substantial errors in convection representation. A new convection-permitting reanalysis over Italy (SPHERA) has been produced at ARPAE to enhance the representation and understanding of extreme weather situations. SPHERA is obtained through a dynamical downscaling of the global reanalysis ERA5 using the non-hydrostatic model COSMO at 2.2 km grid spacing over 1995-2020. This thesis aims to verify the expectations placed on SPHERA by analyzing two weather phenomena that are particularly challenging to simulate: heavy rainfall and hail. A quantitative statistical analysis over Italy during 2003-2017 for daily and hourly precipitation is presented to compare the performance of SPHERA with its driver ERA5 considering the national network of rain gauges as reference. Furthermore, two extreme precipitation events are deeply investigated. SPHERA shows a quantitative added skill over ERA5 for moderate to severe and rapid accumulations in terms of adherence to the observations, higher detailing of the spatial fields, and more precise temporal matching. These results prompted the use of SPHERA for the investigation of hailstorms, for which the combination of multiple information is crucial to reduce the substantial uncertainties permeating their understanding. A proxy for hail is developed by combining hail-favoring environmental numerical predictors with observations of ESWD hail reports and satellite overshooting top detections. The procedure is applied to the extended summer season (April-October) of 2016-2018 over the whole SPHERA spatial domain. The results indicate maximum hail likelihood over pre-Alpine regions and the northern Adriatic sea around 15 UTC in June-July, in agreement with recent European hail climatologies. The method demonstrates enhanced performance in case of severe hail occurrences and the ability to separate between ambient signatures depending on hail severity.

Maximum Likelihood Estimation of Mixture Densities for Binned and Truncated Multivariate Data

Binning and truncation of data are common in data analysis and machine learning. This paper addresses the problem of fitting mixture densities to multivariate binned and truncated data. The EM approach proposed by McLachlan and Jones (Biometrics, 44: 2, 571-578, 1988) for the univariate case is generalized to multivariate measurements. The multivariate solution requires the evaluation of multidimensional integrals over each bin at each iteration of the EM procedure. Naive implementation of the procedure can lead to computationally inefficient results. To reduce the computational cost a number of straightforward numerical techniques are proposed. Results on simulated data indicate that the proposed methods can achieve significant computational gains with no loss in the accuracy of the final parameter estimates. Furthermore, experimental results suggest that with a sufficient number of bins and data points it is possible to estimate the true underlying density almost as well as if the data were not binned. The paper concludes with a brief description of an application of this approach to diagnosis of iron deficiency anemia, in the context of binned and truncated bivariate measurements of volume and hemoglobin concentration from an individual's red blood cells.

Potential-Scour Assessments and Estimates of Maximum Scour at Selected Bridges in Iowa, HR-344, 1995

This report presents the results of potential-scour assessments at 130 bridges and estimates of maximum scour at 10 bridges, in Iowa. All of the bridges evaluated in the study are constructed bridges (not culverts) that are sites of active or discontinued streamflow-gaging stations and peak-stage measurement sites. The period of the study was from October 1991 to September 1994. The potential-scour assessments were made using a potential-scour index developed by the U.S. Geological Survey for a study in Tennessee. Higher values of the index suggest a greater likelihood of scour-related problems occurring at a bridge. The estimates of maximum scour were made using scour equations recommended by the Federal Highway Administration. In this study, the long term aggradation or degradation that occurred during the period of streamflow data collection at each site was evaluated. Although the abutment-scour equation predicted deep scour holes at many of the sites, the only significant abutment scour that was measured was erosion of the embankment at the left abutment at one bridge after a flood.

Improved maximum likelihood estimators in a heteroskedastic errors-in-variables model

This paper develops a bias correction scheme for a multivariate heteroskedastic errors-in-variables model. The applicability of this model is justified in areas such as astrophysics, epidemiology and analytical chemistry, where the variables are subject to measurement errors and the variances vary with the observations. We conduct Monte Carlo simulations to investigate the performance of the corrected estimators. The numerical results show that the bias correction scheme yields nearly unbiased estimates. We also give an application to a real data set.

Genetic parameters of milk, fat, and protein yields in the first three lactations, using an animal model and restricted maximum likelihood

Milk, fat, and protein yields of Holstein cows from the States of New York and California in the United States were used to estimate (co)variances among yields in the first three lactations, using an animal model and a derivative-free restricted maximum likelihood (REML) algorithm, and to verify if yields in different lactations are the same trait. The data were split in 20 samples, 10 from each state, with means of 5463 and 5543 cows per sample from California and New York. Mean heritability estimates for milk, fat, and protein yields for California data were, respectively, 0.34, 0.35, and 0.40 for first; 0.31, 0.33, and 0.39 for second; and 0.28, 0.31, and 0.37 for third lactations. For New York data, estimates were 0.35, 0.40, and 0.34 for first; 0.34, 0.44, and 0.38 for second; and 0.32, 0.43, and 0.38 for third lactations. Means of estimates of genetic correlations between first and second, first and third, and second and third lactations for California data were 0.86, 0.77, and 0.96 for milk; 0.89, 0.84, and 0.97 for fat; and 0.90, 0.84, and 0.97 for protein yields. Mean estimates for New York data were 0.87, 0.81, and 0.97 for milk; 0.91, 0.86, and 0.98 for fat; and 0.88, 0.82, and 0.98 for protein yields. Environmental correlations varied from 0.30 to 0.50 and were larger between second and third lactations. Phenotypic correlations were similar for both states and varied from 0.52 to 0.66 for milk, fat and protein yields. These estimates are consistent with previous estimates obtained with animal models. Yields in different lactations are not statistically the same trait but for selection programs such yields can be modelled as the same trait because of the high genetic correlations.

Maximum Likelihood Estimation of Modal Parameters in Structures Using the Expectation Maximization Algorithm

This paper presents a time-domain stochastic system identification method based on maximum likelihood estimation (MLE) with the expectation maximization (EM) algorithm. The effectiveness of this structural identification method is evaluated through numerical simulation in the context of the ASCE benchmark problem on structural health monitoring. The benchmark structure is a four-story, two-bay by two-bay steel-frame scale model structure built in the Earthquake Engineering Research Laboratory at the University of British Columbia, Canada. This paper focuses on Phase I of the analytical benchmark studies. A MATLAB-based finite element analysis code obtained from the IASC-ASCE SHM Task Group web site is used to calculate the dynamic response of the prototype structure. A number of 100 simulations have been made using this MATLAB-based finite element analysis code in order to evaluate the proposed identification method. There are several techniques to realize system identification. In this work, stochastic subspace identification (SSI)method has been used for comparison. SSI identification method is a well known method and computes accurate estimates of the modal parameters. The principles of the SSI identification method has been introduced in the paper and next the proposed MLE with EM algorithm has been explained in detail. The advantages of the proposed structural identification method can be summarized as follows: (i) the method is based on maximum likelihood, that implies minimum variance estimates; (ii) EM is a computational simpler estimation procedure than other optimization algorithms; (iii) estimate more parameters than SSI, and these estimates are accurate. On the contrary, the main disadvantages of the method are: (i) EM algorithm is an iterative procedure and it consumes time until convergence is reached; and (ii) this method needs starting values for the parameters. Modal parameters (eigenfrequencies, damping ratios and mode shapes) of the benchmark structure have been estimated using both the SSI method and the proposed MLE + EM method. The numerical results show that the proposed method identifies eigenfrequencies, damping ratios and mode shapes reasonably well even in the presence of 10% measurement noises. These modal parameters are more accurate than the SSI estimated modal parameters.

Maximum likelihood estimation of a migration matrix and effective population sizes in n subpopulations by using a coalescent approach

A maximum likelihood estimator based on the coalescent for unequal migration rates and different subpopulation sizes is developed. The method uses a Markov chain Monte Carlo approach to investigate possible genealogies with branch lengths and with migration events. Properties of the new method are shown by using simulated data from a four-population n-island model and a source–sink population model. Our estimation method as coded in migrate is tested against genetree; both programs deliver a very similar likelihood surface. The algorithm converges to the estimates fairly quickly, even when the Markov chain is started from unfavorable parameters. The method was used to estimate gene flow in the Nile valley by using mtDNA data from three human populations.

A maximum likelihood algorithm for genome mapping of cytogenetic loci from meiotic configuration data.

Frequencies of meiotic configurations in cytogenetic stocks are dependent on chiasma frequencies in segments defined by centromeres, breakpoints, and telomeres. The expectation maximization algorithm is proposed as a general method to perform maximum likelihood estimations of the chiasma frequencies in the intervals between such locations. The estimates can be translated via mapping functions into genetic maps of cytogenetic landmarks. One set of observational data was analyzed to exemplify application of these methods, results of which were largely concordant with other comparable data. The method was also tested by Monte Carlo simulation of frequencies of meiotic configurations from a monotelodisomic translocation heterozygote, assuming six different sample sizes. The estimate averages were always close to the values given initially to the parameters. The maximum likelihood estimation procedures can be extended readily to other kinds of cytogenetic stocks and allow the pooling of diverse cytogenetic data to collectively estimate lengths of segments, arms, and chromosomes.

Didanosine-loaded Chitosan Microspheres Optimized By Surface-response Methodology: A Modified Maximum Likelihood Classification" Approach Formulation For Reverse Transcriptase Inhibitors."

Didanosine-loaded chitosan microspheres were developed applying a surface-response methodology and using a modified Maximum Likelihood Classification. The operational conditions were optimized with the aim of maintaining the active form of didanosine (ddI), which is sensitive to acid pH, and to develop a modified and mucoadhesive formulation. The loading of the drug within the chitosan microspheres was carried out by ionotropic gelation technique with sodium tripolyphosphate (TPP) as cross-linking agent and magnesium hydroxide (Mg(OH)2) to assure the stability of ddI. The optimization conditions were set using a surface-response methodology and applying the Maximum Likelihood Classification, where the initial chitosan concentration, TPP and ddI concentration were set as the independent variables. The maximum ddI-loaded in microspheres (i.e. 1433mg of ddI/g chitosan), was obtained with 2% (w/v) chitosan and 10% TPP. The microspheres depicted an average diameter of 11.42μm and ddI was gradually released during 2h in simulated enteric fluid.

Genetic trend estimates of meat quality traits in a male broiler line

The present research was conducted to estimate the genetic trends for meat quality traits in a male broiler line. The traits analyzed were initial pH, pH at 6 h after slaughter, final pH, initial range of falling pH, final range of falling pH, lightness, redness, yellowness, weep loss, drip loss, shrink loss, and shear force. The number of observations varied between 618 and 2125 for each trait. Genetic values were obtained by restricted maximum likelihood, and the numerator relationship matrix had 107,154 animals. The genetic trends were estimated by regression of the broiler average genetic values with respect to unit of time (generations), and the average genetic trend was estimated by regression coefficients. Generally, for the traits analyzed, small genetic trends were obtained, except for drip loss and shear force, which were higher. The small magnitude of the trends found could be a consequence of the absence of selection for meat quality traits in the line analyzed. The estimates of genetic trends obtained were an indication of an improvement in the meat quality traits in the line analyzed, except for drip loss.

Estimates of genetic trend for carcass traits in a commercial broiler line

Data from the slaughter of 24,001 chickens that were part of a selection program for the production of commercial broilers were used to estimate genetic trend for absolute carcass (CW), breast meat (BRW), and leg (LW) weights, and relative carcass (CY), breast meat (BRY), and leg (LY) weights. The components of (co) variance and breeding values of individuals were obtained by the restricted maximum likelihood method applied to animal models. The relationship matrix was composed of 132,442 birds. The models included as random effects, maternal additive genetic and permanent environmental for CW, BRW, LW, CY, and BRY, and only maternal permanent environmental for LY, besides the direct additive genetic and residual effects, and as fixed effects, hatch week, parents' mating group and sex. The estimates of genetic trend were obtained by average regression of breeding value on generation, and the average genetic trend was estimated by regression coefficients. The genetic trends for CW (+ 6.0336 g/generation), BRW (+ 3.6723 g/generation), LW (+ 1.5846 g/generation), CY (+ 0.1195%/generation), and BRY (+ 0.1388%/generation) were positive, and they were in accordance with the objectives of the selection program for these traits. The genetic trend for LY(-0.0019%/generation) was negative, possibly due to the strong emphasis on selection for BRY and the negative correlations between these two traits.

An integrated molecular map of yellow passion fruit based on simultaneous maximum-likelihood estimation of linkage and linkage phases

The development of genetic maps for auto-incompatible species, such as the yellow passion fruit (Passiflora edulis Sims f.flavicarpa Deg.) is restricted due to the unfeasibility of obtaining traditional mapping populations based on inbred lines. For this reason, yellow passion fruit linkage maps were generally constructed using a strategy known as two-way pseudo-testeross, based on monoparental dominant markers segregating in a 1:1 fashion. Due to the lack of information from these markers in one of the parents, two individual (parental) maps were obtained. However, integration of these maps is essential, and biparental markers can be used for such an operation. The objective of our study was to construct an integrated molecular map for a full-sib population of yellow passion fruit combining different loci configuration generated from amplified fragment length polymorphisms (AFLPs) and microsatellite markers and using a novel approach based on simultaneous maximum-likelihood estimation of linkage and linkage phases, specially designed for outcrossing species. Of the total number of loci, approximate to 76%, 21%, 0.7%, and 2.3% did segregate in 1:1, 3:1, 1:2:1, and 1:1:1:1 ratios, respectively. Ten linkage groups (LGs) were established with a logarithm of the odds (LOD) score >= 5.0 assuming a recombination fraction : <= 0.35. On average, 24 markers were assigned per LG, representing a total map length of 1687 cM, with a marker density of 6.9 cM. No markers were placed as accessories on the map as was done with previously constructed individual maps.