Does the choice of neighbourhood supermarket access measure influence associations with individual-level fruit and vegetable consumption? A case study from Glasgow

Background
Previous studies have provided mixed evidence with regards to associations between food store access and dietary outcomes. This study examines the most commonly applied measures of locational access to assess whether associations between supermarket access and fruit and vegetable consumption are affected by the choice of access measure and scale.

Method
Supermarket location data from Glasgow, UK (n = 119), and fruit and vegetable intake data from the 'Health and Well-Being' Survey (n = 1041) were used to compare various measures of locational access. These exposure variables included proximity estimates (with different points-of-origin used to vary levels of aggregation) and density measures using three approaches (Euclidean and road network buffers and Kernel density estimation) at distances ranging from 0.4 km to 5 km. Further analysis was conducted to assess the impact of using smaller buffer sizes for individuals who did not own a car. Associations between these multiple access measures and fruit and vegetable consumption were estimated using linear regression models.

Results
Levels of spatial aggregation did not impact on the proximity estimates. Counts of supermarkets within Euclidean buffers were associated with fruit and vegetable consumption at 1 km, 2 km and 3 km, and for our road network buffers at 2 km, 3 km, and 4 km. Kernel density estimates provided the strongest associations and were significant at a distance of 2 km, 3 km, 4 km and 5 km. Presence of a supermarket within 0.4 km of road network distance from where people lived was positively associated with fruit consumption amongst those without a car (coef. 0.657; s.e. 0.247; p0.008).

Conclusions
The associations between locational access to supermarkets and individual-level dietary behaviour are sensitive to the method by which the food environment variable is captured. Care needs to be taken to ensure robust and conceptually appropriate measures of access are used and these should be grounded in a clear a priori reasoning.

Fidelity to foraging sites, consistency of migration routes and habitat modulation of home range by sea turtles

Aim  Resources can shape patterns of habitat utilization. Recently a broad foraging dichotomy between oceanic and coastal sites has been revealed for loggerhead sea turtles (Caretta caretta). Since oceanic and coastal foraging sites differ in prey availability, we might expect a gross difference in home-range size across these habitats. We tested this hypothesis by equipping nine adult male loggerhead sea turtles with GPS tracking devices. Location  National Marine Park of Zakynthos (NMPZ) Greece, central and eastern Mediterranean (Adriatic, Ionian and Aegean seas). Methods  In 2007, 2008 and 2009, Fastloc GPS-Argos transmitters were attached to nine male loggerheads. In addition, a Sirtrack PTT unit was attached to one male in 2007. Four of the turtles were tracked on successive years. We filtered the GPS data to ensure comparable data volumes. Route consistency between breeding and foraging sites of the four re-tracked turtles was conducted. Foraging site home range areas and within site movement patterns were investigated by the fixed kernel density method. Results  Foraging home range size ranged between circa 10 km2 at neritic habitats (coastal and open-sea on the continental shelf) to circa 1000 km2 at oceanic sites (using 90% kernel estimates), the latter most probably reflecting sparsely distributed oceanic prey. Across different years individuals did not follow exactly the same migration routes, but did show fidelity to their previous foraging sites, whether oceanic or neritic, with accurate homing in the final stages of migration. Main conclusions  The broad distribution and diverse life-history strategies of this population could complicate the identification of priority marine protected areas beyond the core breeding site.

Visualization of HIV-1 interactions with penile and foreskin epithelia: clues for female-to-male HIV transmission

To gain insight into female-to-male HIV sexual transmission and how male circumcision protects against this mode of transmission, we visualized HIV-1 interactions with foreskin and penile tissues in ex vivo tissue culture and in vivo rhesus macaque models utilizing epifluorescent microscopy. 12 foreskin and 14 cadaveric penile specimens were cultured with R5-tropic photoactivatable (PA)-GFP HIV-1 for 4 or 24 hours. Tissue cryosections were immunofluorescently imaged for epithelial and immune cell markers. Images were analyzed for total virions, proportion of penetrators, depth of virion penetration, as well as immune cell counts and depths in the tissue. We visualized individual PA virions breaching penile epithelial surfaces in the explant and macaque model. Using kernel density estimated probabilities of localizing a virion or immune cell at certain tissue depths revealed that interactions between virions and cells were more likely to occur in the inner foreskin or glans penis (from local or cadaveric donors, respectively). Using statistical models to account for repeated measures and zero-inflated datasets, we found no difference in total virions visualized at 4 hours between inner and outer foreskins from local donors. At 24 hours, there were more virions in inner as compared to outer foreskin (0.0495 +/- 0.0154 and 0.0171 +/- 0.0038 virions/image, p = 0.001). In the cadaveric specimens, we observed more virions in inner foreskin (0.0507 +/- 0.0079 virions/image) than glans tissue (0.0167 +/- 0.0033 virions/image, p<0.001), but a greater proportion was seen penetrating uncircumcised glans tissue (0.0458 +/- 0.0188 vs. 0.0151 +/- 0.0100 virions/image, p = 0.099) and to significantly greater mean depths (29.162 +/- 3.908 vs. 12.466 +/- 2.985 μm). Our in vivo macaque model confirmed that virions can breach penile squamous epithelia in a living model. In summary, these results suggest that the inner foreskin and glans epithelia may be important sites for HIV transmission in uncircumcised men.

Time-in-area represents foraging activity in a wide-ranging pelagic forager

Successful Marine Spatial Planning depends upon the identification of areas with high importance for particular species, ecosystems or processes. For seabirds, advancements in biologging devices have enabled us to identify these areas through the detailed study of at-sea behaviour. However, in many cases, only positional data are available and the presence of local biological productivity and hence seabird foraging behaviour is inferred from these data alone, under the untested assumption that foraging activity is more likely to occur in areas where seabirds spend more time. We fitted GPS devices and accelerometers to northern gannets Morus bassanus and categorised the behaviour of individuals outside the breeding colony as plunge diving, surface foraging, floating and flying. We then used the locations of foraging events to test the efficiency of 2 approaches: time-in-area and kernel density (KD) analyses, which are widely employed to detect highly-used areas and interpret foraging behaviour from positional data. For KD analyses, the smoothing parameter (h) was calculated using the ad hoc method (KDad hoc), and KDh=9.1, where h = 9.1 km, to designate core foraging areas from location data. A high proportion of foraging events occurred in core foraging areas designated using KDad hoc, KDh=9.1, and time-in-area. Our findings demonstrate that foraging activity occurs in areas where seabirds spend more time, and that both KD analysis and the time-in-area approach are equally efficient methods for this type of analysis. However, the time-in-area approach is advantageous in its simplicity, and in its ability to provide the shapes commonly used in planning. Therefore, the time-in-area approach can be used as a simple way of using seabirds to identify ecologically important locations from both tracking and survey data.

Análise de técnicas de distribuição espacial com padrões pontuais e aplicação a dados de acidentes de trânsito e a dados da dengue de Rio Claro-SP

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Epidemiologia e prevalência de infecções pelo Vírus da Diarreia Viral Bovina (BVDV) em suínos de criações não tecnificadas

Pós-graduação em Medicina Veterinária - FCAV

Morphological Homogeneity of Neurons: Searching for Outlier Neuronal Cells

We report a morphology-based approach for the automatic identification of outlier neurons, as well as its application to the NeuroMorpho.org database, with more than 5,000 neurons. Each neuron in a given analysis is represented by a feature vector composed of 20 measurements, which are then projected into a two-dimensional space by applying principal component analysis. Bivariate kernel density estimation is then used to obtain the probability distribution for the group of cells, so that the cells with highest probabilities are understood as archetypes while those with the smallest probabilities are classified as outliers. The potential of the methodology is illustrated in several cases involving uniform cell types as well as cell types for specific animal species. The results provide insights regarding the distribution of cells, yielding single and multi-variate clusters, and they suggest that outlier cells tend to be more planar and tortuous. The proposed methodology can be used in several situations involving one or more categories of cells, as well as for detection of new categories and possible artifacts.

Spatial Disparity and Dynamic Trajectory of Convergence in China: Construction and Application of a Composite Index of Regional Development

China is a large country characterized by remarkable growth and distinct regional diversity. Spatial disparity has always been a hot issue since China has been struggling to follow a balanced growth path but still confronting with unprecedented pressures and challenges. To better understand the inequality level benchmarking spatial distributions of Chinese provinces and municipalities and estimate dynamic trajectory of sustainable development in China, I constructed the Composite Index of Regional Development (CIRD) with five sub pillars/dimensions involving Macroeconomic Index (MEI), Science and Innovation Index (SCI), Environmental Sustainability Index (ESI), Human Capital Index (HCI) and Public Facilities Index (PFI), endeavoring to cover various fields of regional socioeconomic development. Ranking reports on the five sub dimensions and aggregated CIRD were provided in order to better measure the developmental degrees of 31 or 30 Chinese provinces and municipalities over 13 years from 1998 to 2010 as the time interval of three “Five-year Plans”. Further empirical applications of this CIRD focused on clustering and convergence estimation, attempting to fill up the gap in quantifying the developmental levels of regional comprehensive socioeconomics and estimating the dynamic convergence trajectory of regional sustainable development in a long run. Four clusters were benchmarked geographically-oriented in the map on the basis of cluster analysis, and club-convergence was observed in the Chinese provinces and municipalities based on stochastic kernel density estimation.

Wavelet estimation in diffusions with periodicity

Wir betrachten einen zeitlich inhomogenen Diffusionsprozess, der durch eine stochastische Differentialgleichung gegeben wird, deren Driftterm ein deterministisches T-periodisches Signal beinhaltet, dessen Periodizität bekannt ist. Dieses Signal sei in einem Besovraum enthalten. Wir schätzen es mit Hilfe eines nichtparametrischen Waveletschätzers. Unser Schätzer ist von einem Wavelet-Dichteschätzer mit Thresholding inspiriert, der 1996 in einem klassischen iid-Modell von Donoho, Johnstone, Kerkyacharian und Picard konstruiert wurde. Unter gewissen Ergodizitätsvoraussetzungen an den Prozess können wir nichtparametrische Konvergenzraten angegeben, die bis auf einen logarithmischen Term den Raten im klassischen iid-Fall entsprechen. Diese Raten werden mit Hilfe von Orakel-Ungleichungen gezeigt, die auf Ergebnissen über Markovketten in diskreter Zeit von Clémencon, 2001, beruhen. Außerdem betrachten wir einen technisch einfacheren Spezialfall und zeigen einige Computersimulationen dieses Schätzers.

Evaluating northern hardwood management using retrospective analysis and diameter distributions

Northern hardwood management was assessed throughout the state of Michigan using data collected on recently harvested stands in 2010 and 2011. Methods of forensic estimation of diameter at breast height were compared and an ideal, localized equation form was selected for use in reconstructing pre-harvest stand structures. Comparisons showed differences in predictive ability among available equation forms which led to substantial financial differences when used to estimate the value of removed timber. Management on all stands was then compared among state, private, and corporate landowners. Comparisons of harvest intensities against a liberal interpretation of a well-established management guideline showed that approximately one third of harvests were conducted in a manner which may imply that the guideline was followed. One third showed higher levels of removals than recommended, and one third of harvests were less intensive than recommended. Multiple management guidelines and postulated objectives were then synthesized into a novel system of harvest taxonomy, against which all harvests were compared. This further comparison showed approximately the same proportions of harvests, while distinguishing sanitation cuts and the future productive potential of harvests cut more intensely than suggested by guidelines. Stand structures are commonly represented using diameter distributions. Parametric and nonparametric techniques for describing diameter distributions were employed on pre-harvest and post-harvest data. A common polynomial regression procedure was found to be highly sensitive to the method of histogram construction which provides the data points for the regression. The discriminative ability of kernel density estimation was substantially different from that of the polynomial regression technique.

Sexual dimorphism dominates divergent host plant use in stick insect trophic morphology

Background: Clear examples of ecological speciation exist, often involving divergence in trophic morphology. However, substantial variation also exists in how far the ecological speciation process proceeds, potentially linked to the number of ecological axes, traits, or genes subject to divergent selection. In addition, recent studies highlight how differentiation might occur between the sexes, rather than between populations. We examine variation in trophic morphology in two host-plant ecotypes of walking-stick insects (Timema cristinae), known to have diverged in morphological traits related to crypsis and predator avoidance, and to have reached an intermediate point in the ecological speciation process. Here we test how host plant use, sex, and rearing environment affect variation in trophic morphology in this species using traditional multivariate, novel kernel density based and Bayesian morphometric analyses. Results: Contrary to expectations, we find limited host-associated divergence in mandible shape. Instead, the main predictor of shape variation is sex, with secondary roles of population of origin and rearing environment. Conclusion: Our results show that trophic morphology does not strongly contribute to host-adapted ecotype divergence in T. cristinae and that traits can respond to complex selection regimes by diverging along different intraspecific lines, thereby impeding progress toward speciation.

Statistical characterization and development of summary measures of non-normal distributions of nuclear morphometry data from prostate cancer cells for use as prognostic indicators

Nuclear morphometry (NM) uses image analysis to measure features of the cell nucleus which are classified as: bulk properties, shape or form, and DNA distribution. Studies have used these measurements as diagnostic and prognostic indicators of disease with inconclusive results. The distributional properties of these variables have not been systematically investigated although much of the medical data exhibit nonnormal distributions. Measurements are done on several hundred cells per patient so summary measurements reflecting the underlying distribution are needed.^ Distributional characteristics of 34 NM variables from prostate cancer cells were investigated using graphical and analytical techniques. Cells per sample ranged from 52 to 458. A small sample of patients with benign prostatic hyperplasia (BPH), representing non-cancer cells, was used for general comparison with the cancer cells.^ Data transformations such as log, square root and 1/x did not yield normality as measured by the Shapiro-Wilks test for normality. A modulus transformation, used for distributions having abnormal kurtosis values, also did not produce normality.^ Kernel density histograms of the 34 variables exhibited non-normality and 18 variables also exhibited bimodality. A bimodality coefficient was calculated and 3 variables: DNA concentration, shape and elongation, showed the strongest evidence of bimodality and were studied further.^ Two analytical approaches were used to obtain a summary measure for each variable for each patient: cluster analysis to determine significant clusters and a mixture model analysis using a two component model having a Gaussian distribution with equal variances. The mixture component parameters were used to bootstrap the log likelihood ratio to determine the significant number of components, 1 or 2. These summary measures were used as predictors of disease severity in several proportional odds logistic regression models. The disease severity scale had 5 levels and was constructed of 3 components: extracapsulary penetration (ECP), lymph node involvement (LN+) and seminal vesicle involvement (SV+) which represent surrogate measures of prognosis. The summary measures were not strong predictors of disease severity. There was some indication from the mixture model results that there were changes in mean levels and proportions of the components in the lower severity levels. ^

Income Inequality in the 1990s: Re-Forging a Lost Relationship?

Using data from March Current Population Surveys we find gains from economic growth over the 1990s business cycle (1989-2000) were more equitably distributed than over the 1980s business cycle (1979-1989) using summary inequality measures as well as kernel density estimations. The entire distribution of household size-adjusted income moved upwards in the 1990s with profound improvements for African Americans, single mothers and those living in households receiving welfare. Most gains occurred over the growth period 1993-2000. Improvements in average income and income inequity over the latter period are reminiscent of gains seen in the first three decades after World War II.

(Table 1) Tagging date, total length and sex of bowhead whales (Balaena mysticetus) tagged at Point Barrow, Alaska

Working with subsistence whale hunters, we tagged 19 mostly immature bowhead whales (Balaena mysticetus) with satellite-linked transmitters between May 2006 and September 2008 and documented their movements in the Chukchi Sea from late August through December. From Point Barrow, Alaska, most whales moved west through the Chukchi Sea between 71° and 74° N latitude; nine whales crossed in six to nine days. Three whales returned to Point Barrow for 13 to 33 days, two after traveling 300 km west and one after traveling ~725 km west to Wrangel Island, Russia; two then crossed the Chukchi Sea again while the other was the only whale to travel south along the Alaskan side of the Chukchi Sea. Seven whales spent from one to 21 days near Wrangel Island before moving south to northern Chukotka. Whales spent an average of 59 days following the Chukotka coast southeastward. Kernel density analysis identified Point Barrow, Wrangel Island, and the northern coast of Chukotka as areas of greater use by bowhead whales that might be important for feeding. All whales traveled through a potential petroleum development area at least once. Most whales crossed the development area in less than a week; however, one whale remained there for 30 days.

Parallel Computer Vision Algorithms for Graphics Processing Units

La evolución de los teléfonos móviles inteligentes, dotados de cámaras digitales, está provocando una creciente demanda de aplicaciones cada vez más complejas que necesitan algoritmos de visión artificial en tiempo real; puesto que el tamaño de las señales de vídeo no hace sino aumentar y en cambio el rendimiento de los procesadores de un solo núcleo se ha estancado, los nuevos algoritmos que se diseñen para visión artificial han de ser paralelos para poder ejecutarse en múltiples procesadores y ser computacionalmente escalables. Una de las clases de procesadores más interesantes en la actualidad se encuentra en las tarjetas gráficas (GPU), que son dispositivos que ofrecen un alto grado de paralelismo, un excelente rendimiento numérico y una creciente versatilidad, lo que los hace interesantes para llevar a cabo computación científica. En esta tesis se exploran dos aplicaciones de visión artificial que revisten una gran complejidad computacional y no pueden ser ejecutadas en tiempo real empleando procesadores tradicionales. En cambio, como se demuestra en esta tesis, la paralelización de las distintas subtareas y su implementación sobre una GPU arrojan los resultados deseados de ejecución con tasas de refresco interactivas. Asimismo, se propone una técnica para la evaluación rápida de funciones de complejidad arbitraria especialmente indicada para su uso en una GPU. En primer lugar se estudia la aplicación de técnicas de síntesis de imágenes virtuales a partir de únicamente dos cámaras lejanas y no paralelas—en contraste con la configuración habitual en TV 3D de cámaras cercanas y paralelas—con información de color y profundidad. Empleando filtros de mediana modificados para la elaboración de un mapa de profundidad virtual y proyecciones inversas, se comprueba que estas técnicas son adecuadas para una libre elección del punto de vista. Además, se demuestra que la codificación de la información de profundidad con respecto a un sistema de referencia global es sumamente perjudicial y debería ser evitada. Por otro lado se propone un sistema de detección de objetos móviles basado en técnicas de estimación de densidad con funciones locales. Este tipo de técnicas es muy adecuada para el modelado de escenas complejas con fondos multimodales, pero ha recibido poco uso debido a su gran complejidad computacional. El sistema propuesto, implementado en tiempo real sobre una GPU, incluye propuestas para la estimación dinámica de los anchos de banda de las funciones locales, actualización selectiva del modelo de fondo, actualización de la posición de las muestras de referencia del modelo de primer plano empleando un filtro de partículas multirregión y selección automática de regiones de interés para reducir el coste computacional. Los resultados, evaluados sobre diversas bases de datos y comparados con otros algoritmos del estado del arte, demuestran la gran versatilidad y calidad de la propuesta. Finalmente se propone un método para la aproximación de funciones arbitrarias empleando funciones continuas lineales a tramos, especialmente indicada para su implementación en una GPU mediante el uso de las unidades de filtraje de texturas, normalmente no utilizadas para cómputo numérico. La propuesta incluye un riguroso análisis matemático del error cometido en la aproximación en función del número de muestras empleadas, así como un método para la obtención de una partición cuasióptima del dominio de la función para minimizar el error. ABSTRACT The evolution of smartphones, all equipped with digital cameras, is driving a growing demand for ever more complex applications that need to rely on real-time computer vision algorithms. However, video signals are only increasing in size, whereas the performance of single-core processors has somewhat stagnated in the past few years. Consequently, new computer vision algorithms will need to be parallel to run on multiple processors and be computationally scalable. One of the most promising classes of processors nowadays can be found in graphics processing units (GPU). These are devices offering a high parallelism degree, excellent numerical performance and increasing versatility, which makes them interesting to run scientific computations. In this thesis, we explore two computer vision applications with a high computational complexity that precludes them from running in real time on traditional uniprocessors. However, we show that by parallelizing subtasks and implementing them on a GPU, both applications attain their goals of running at interactive frame rates. In addition, we propose a technique for fast evaluation of arbitrarily complex functions, specially designed for GPU implementation. First, we explore the application of depth-image–based rendering techniques to the unusual configuration of two convergent, wide baseline cameras, in contrast to the usual configuration used in 3D TV, which are narrow baseline, parallel cameras. By using a backward mapping approach with a depth inpainting scheme based on median filters, we show that these techniques are adequate for free viewpoint video applications. In addition, we show that referring depth information to a global reference system is ill-advised and should be avoided. Then, we propose a background subtraction system based on kernel density estimation techniques. These techniques are very adequate for modelling complex scenes featuring multimodal backgrounds, but have not been so popular due to their huge computational and memory complexity. The proposed system, implemented in real time on a GPU, features novel proposals for dynamic kernel bandwidth estimation for the background model, selective update of the background model, update of the position of reference samples of the foreground model using a multi-region particle filter, and automatic selection of regions of interest to reduce computational cost. The results, evaluated on several databases and compared to other state-of-the-art algorithms, demonstrate the high quality and versatility of our proposal. Finally, we propose a general method for the approximation of arbitrarily complex functions using continuous piecewise linear functions, specially formulated for GPU implementation by leveraging their texture filtering units, normally unused for numerical computation. Our proposal features a rigorous mathematical analysis of the approximation error in function of the number of samples, as well as a method to obtain a suboptimal partition of the domain of the function to minimize approximation error.