A general class of zero-or-one inflated beta regression models

This paper proposes a general class of regression models for continuous proportions when the data contain zeros or ones. The proposed class of models assumes that the response variable has a mixed continuous-discrete distribution with probability mass at zero or one. The beta distribution is used to describe the continuous component of the model, since its density has a wide range of different shapes depending on the values of the two parameters that index the distribution. We use a suitable parameterization of the beta law in terms of its mean and a precision parameter. The parameters of the mixture distribution are modeled as functions of regression parameters. We provide inference, diagnostic, and model selection tools for this class of models. A practical application that employs real data is presented. (C) 2011 Elsevier B.V. All rights reserved.

Bayesian analysis of random regression models using B-splines to model test-day milk yield of Holstein cattle in Brazil

The objective of this paper is to model variations in test-day milk yields of first lactations of Holstein cows by RR using B-spline functions and Bayesian inference in order to fit adequate and parsimonious models for the estimation of genetic parameters. They used 152,145 test day milk yield records from 7317 first lactations of Holstein cows. The model established in this study was additive, permanent environmental and residual random effects. In addition, contemporary group and linear and quadratic effects of the age of cow at calving were included as fixed effects. Authors modeled the average lactation curve of the population with a fourth-order orthogonal Legendre polynomial. They concluded that a cubic B-spline with seven random regression coefficients for both the additive genetic and permanent environment effects was to be the best according to residual mean square and residual variance estimates. Moreover they urged a lower order model (quadratic B-spline with seven random regression coefficients for both random effects) could be adopted because it yielded practically the same genetic parameter estimates with parsimony. (C) 2012 Elsevier B.V. All rights reserved.

Suicide rates and trends in Sao Paulo, Brazil, according to gender, age and demographic aspects: a joinpoint regression analysis

Objective: To evaluate suicide rates and trends in Sao Paulo by sex, age-strata, and methods. Methods: Data was collected from State registry from 1996 to 2009. Population was estimated using the National Census. We utilized joinpoint regression analysis to explore temporal trends. We also evaluated marital status, ethnicity, birthplace and methods for suicide. Results: In the period analyzed, 6,002 suicides were accrued with a rate of 4.6 per 100,000 (7.5 in men and 2.0 in women); the male-to-female ratio was around 3.7. Trends for men presented a significant decline of 5.3% per year from 1996 to 2002, and a significant increase of 2.5% from 2002 onwards. Women did not present significant changes. For men, the elderly (> 65 years) had a significant reduction of 2.3% per year, while younger men (25-44 years) presented a significant increase of 8.6% from 2004 onwards. Women did not present significant trend changes according to age. Leading suicide methods were hanging and poisoning for men and women, respectively. Other analyses showed an increased suicide risk ratio for singles and foreigners. Conclusions: Specific epidemiological trends for suicide in the city of Sao Paulo that warrant further investigation were identified. High-risk groups - such as immigrants - could benefit from targeted strategies of suicide prevention.

Random regression analyses using B-spline functions to model growth of Nellore cattle

The objective of this study was to estimate (co)variance components using random regression on B-spline functions to weight records obtained from birth to adulthood. A total of 82 064 weight records of 8145 females obtained from the data bank of the Nellore Breeding Program (PMGRN/Nellore Brazil) which started in 1987, were used. The models included direct additive and maternal genetic effects and animal and maternal permanent environmental effects as random. Contemporary group and dam age at calving (linear and quadratic effect) were included as fixed effects, and orthogonal Legendre polynomials of age (cubic regression) were considered as random covariate. The random effects were modeled using B-spline functions considering linear, quadratic and cubic polynomials for each individual segment. Residual variances were grouped in five age classes. Direct additive genetic and animal permanent environmental effects were modeled using up to seven knots (six segments). A single segment with two knots at the end points of the curve was used for the estimation of maternal genetic and maternal permanent environmental effects. A total of 15 models were studied, with the number of parameters ranging from 17 to 81. The models that used B-splines were compared with multi-trait analyses with nine weight traits and to a random regression model that used orthogonal Legendre polynomials. A model fitting quadratic B-splines, with four knots or three segments for direct additive genetic effect and animal permanent environmental effect and two knots for maternal additive genetic effect and maternal permanent environmental effect, was the most appropriate and parsimonious model to describe the covariance structure of the data. Selection for higher weight, such as at young ages, should be performed taking into account an increase in mature cow weight. Particularly, this is important in most of Nellore beef cattle production systems, where the cow herd is maintained on range conditions. There is limited modification of the growth curve of Nellore cattle with respect to the aim of selecting them for rapid growth at young ages while maintaining constant adult weight.

Predicting smear negative pulmonary tuberculosis with classification trees and logistic regression: a cross-sectional study

Abstract Background Smear negative pulmonary tuberculosis (SNPT) accounts for 30% of pulmonary tuberculosis cases reported yearly in Brazil. This study aimed to develop a prediction model for SNPT for outpatients in areas with scarce resources. Methods The study enrolled 551 patients with clinical-radiological suspicion of SNPT, in Rio de Janeiro, Brazil. The original data was divided into two equivalent samples for generation and validation of the prediction models. Symptoms, physical signs and chest X-rays were used for constructing logistic regression and classification and regression tree models. From the logistic regression, we generated a clinical and radiological prediction score. The area under the receiver operator characteristic curve, sensitivity, and specificity were used to evaluate the model's performance in both generation and validation samples. Results It was possible to generate predictive models for SNPT with sensitivity ranging from 64% to 71% and specificity ranging from 58% to 76%. Conclusion The results suggest that those models might be useful as screening tools for estimating the risk of SNPT, optimizing the utilization of more expensive tests, and avoiding costs of unnecessary anti-tuberculosis treatment. Those models might be cost-effective tools in a health care network with hierarchical distribution of scarce resources.

Suicide rates and trends in São Paulo, Brazil, according to gender, age and demographic aspects: a joinpoint regression analysis

OBJECTIVE: To evaluate suicide rates and trends in São Paulo by sex, age-strata, and methods. METHODS: Data was collected from State registry from 1996 to 2009. Population was estimated using the National Census. We utilized joinpoint regression analysis to explore temporal trends. We also evaluated marital status, ethnicity, birthplace and methods for suicide. RESULTS: In the period analyzed, 6,002 suicides were accrued with a rate of 4.6 per 100,000 (7.5 in men and 2.0 in women); the male-to-female ratio was around 3.7. Trends for men presented a significant decline of 5.3% per year from 1996 to 2002, and a significant increase of 2.5% from 2002 onwards. Women did not present significant changes. For men, the elderly (> 65 years) had a significant reduction of 2.3% per year, while younger men (25-44 years) presented a significant increase of 8.6% from 2004 onwards. Women did not present significant trend changes according to age. Leading suicide methods were hanging and poisoning for men and women, respectively. Other analyses showed an increased suicide risk ratio for singles and foreigners. CONCLUSIONS: Specific epidemiological trends for suicide in the city of São Paulo that warrant further investigation were identified. High-risk groups - such as immigrants - could benefit from targeted strategies of suicide prevention.

Natural antibody response to Plasmodium falciparum merozoite antigens MSP5, MSP9 and EBA175 is associated to clinical protection in the Brazilian Amazon

BACKGROUND: Antibodies have an essential role in the acquired immune response against blood stage P. falciparum infection. Although several antigens have been identified as important antibody targets, it is still elusive which antigens have to be recognized for clinical protection. Herein, we analyzed antibodies from plasmas from symptomatic or asymptomatic individuals living in the same geographic area in the Western Amazon, measuring their recognition of multiple merozoite antigens. METHODS: Specific fragments of genes encoding merozoite proteins AMA1 and members of MSP and EBL families from circulating P. falciparum field isolates present in asymptomatic and symptomatic patients were amplified by PCR. After cloning and expression of different versions of the antigens as recombinant GST-fusion peptides, we tested the reactivity of patients' plasmas by ELISA and the presence of IgG subclasses in the most reactive plasmas. RESULTS: 11 out of 24 recombinant antigens were recognized by plasmas from either symptomatic or asymptomatic infections. Antibodies to MSP9 (X2(DF=1) = 9.26/p = 0.0047) and MSP5 (X2(DF=1) = 8.29/p = 0.0069) were more prevalent in asymptomatic individuals whereas the opposite was observed for MSP1 block 2-MAD20 (X2(DF=1) = 6.41/p = 0.0206, Fisher's exact test). Plasmas from asymptomatic individuals reacted more intensely against MSP4 (U = 210.5, p < 0.03), MSP5 (U = 212, p < 0.004), MSP9 (U = 189.5, p < 0.002) and EBA175 (U = 197, p < 0.014, Mann-Whitney's U test). IgG1 and IgG3 were predominant for all antigens, but some patients also presented with IgG2 and IgG4. The recognition of MSP5 (OR = 0.112, IC95% = 0.021-0.585) and MSP9 (OR = 0.125, IC95% = 0.030-0.529, cross tab analysis) predicted 8.9 and 8 times less chances, respectively, to present symptoms. Higher antibody levels against MSP5 and EBA175 were associated by odds ratios of 9.4 (IC95% = 1.29-69.25) and 5.7 (IC95% = 1.12-29.62, logistic regression), respectively, with an asymptomatic status. CONCLUSIONS: Merozoite antigens were targets of cytophilic antibodies and antibodies against MSP5, MSP9 and EBA175 were independently associated with decreased symptoms.

Discontinuity preserving in TV-L1 optical flow methods

[EN] We analyze the discontinuity preserving problem in TV-L1 optical flow methods. This type of methods typically creates rounded effects at flow boundaries, which usually do not coincide with object contours. A simple strategy to overcome this problem consists in inhibiting the diffusion at high image gradients. In this work, we first introduce a general framework for TV regularizers in optical flow and relate it with some standard approaches. Our survey takes into account several methods that use decreasing functions for mitigating the diffusion at image contours. Consequently, this kind of strategies may produce instabilities in the estimation of the optical flows. Hence, we study the problem of instabilities and show that it actually arises from an ill-posed formulation. From this study, it is possible to come across with different schemes to solve this problem. One of these consists in separating the pure TV process from the mitigating strategy. This has been used in another work and we demonstrate here that it has a good performance. Furthermore, we propose two alternatives to avoid the instability problems: (i) we study a fully automatic approach that solves the problem based on the information of the whole image; (ii) we derive a semi-automatic approach that takes into account the image gradients in a close neighborhood adapting the parameter in each position. In the experimental results, we present a detailed study and comparison between the different alternatives. These methods provide very good results, especially for sequences with a few dominant gradients. Additionally, a surprising effect of these approaches is that they can cope with occlusions. This can be easily achieved by using strong regularizations and high penalizations at image contours.

Regularization strategies for discontinuity-preserving optical flow methods

[EN]The aim of this work is to study several strategies for the preservation of flow discontinuities in variational optical flow methods. We analyze the combination of robust functionals and diffusion tensors in the smoothness assumption. Our study includes the use of tensors based on decreasing functions, which has shown to provide good results. However, it presents several limitations and usually does not perform better than other basic approaches. It typically introduces instabilities in the computed motion fields in the form of independent \textit{blobs} of vectors with large magnitude...

Influence of habitat and climate change on European bird communities

Die Fragmentierung von Waldgebieten, der Verlust geeigneter Habitate, die Invasion exotischer Arten und globale Klimaveränderung haben auf Artengemeinschaften erhebliche Auswirkungen. Vögel dienen in vielen Fällen als Indikatorarten für Umweltveränderungen und, besonders, für Veränderungen im Zusammenhang mit globaler Erwärmung. In meiner Arbeit habe ich zuerst einen Literaturüberblick über die Auswirkungen globaler Klimaveränderung auf die Verbreitungsgebiete, den Artenreichtum und die Zusammensetzung von Vogelgemeinschaften dargestellt. Zahlreiche Untersuchungen zeigen, daß die Grenzen der Verbreitungsgebiete der meisten Vogelarten mit klimatischen Faktoren korrelieren. Verschiebungen der Verbreitungsgebiete in nördliche Richtung oder in höhere Regionen im Gebirge konnten bereits für viele temperate Vogelarten beobachtet werden. Weiterhin wurde ein zunehmender Artenreichtum besonders in nördlichen Breiten und in höheren Lagen für viele temperate Vogelgemeinschaften vorhergesagt. In trockenen Gebieten ist dagegen mit einer Abnahme des Artenreichtums zu rechnen. Im zweiten Teil meiner Arbeit habe ich untersucht, ob beobachtete Veränderungen in der Zusammensetzung europäischer Vogelgemeinschaften tatsächlich durch aktuelle Klimaveränderungen beeinflußt werden. Das Zugverhalten der Arten war dabei ein Schwerpunkt der Untersuchung, weil zu erwarten war, daß Vogelarten mit verschiedenem Zugverhalten unterschiedlich auf Klimaveränderungen reagieren. Ich habe ein Regressionsmodell genutzt, welches die räumliche Beziehung zwischen dem Anteil von Langstreckenziehern, Kurzstreckenziehern und Standvögeln in europäischen Vogelgemeinschaften und verschiedenen Klimavariablen beschreibt. Für 21 Gebiete in Europa habe ich Daten über beobachtete Veränderungen in der Struktur der Vogelgemeinschaften und isochrone Klimaveränderungen zusammengetragen. Mit Hilfe dieser Klimaveränderungen und dem räumlichen Regressionsmodell konnte ich berechnen, welche Veränderungen in den Vogelgemeinschaften aufgrund der veränderten Klimabedingungen zu erwarten wären und sie mit beobachteten Veränderungen vergleichen. Beobachtete und berechnete Veränderungen korrelierten signifikant miteinander. Die beobachteten Veränderungen konnten nicht durch räumliche Autokorrelationseffekte oder durch alternative Faktoren, wie z.B. Veränderungen in der Landnutzung, erklärt werden. Im dritten Teil der Arbeit untersuchte ich für eine mitteleuropäische Vogelgemeinschaft welchen Einfluß Habitatveränderungen, die Invasion exotischer Arten und die Klimaveränderung auf Veränderungen der Häufigkeit und Verbreitungsgröße der 159 Vogelarten am Bodensee zwischen 1980-1981 und 2000-2002 hatten. Dabei konnte gezeigt werden, daß Veränderungen in der regionalen Abundanz sowohl durch Habitatveränderungen als auch durch Klimavänderungen hervorgerufen wurden. Exotische Arten schienen in dieser Zeit keinen bedeutenden Einfluß zu haben. Besonders bei Agrarlandarten, Arten mit nördlicheren Verbreitungsgebieten und bei Langstreckenziehern konnten signifikante Abnahmen in der Abundanz beobachtet werden. Vor allem die anhaltenden negativen Bestandsveränderungen bei Langstreckenziehern und die in den letzten zehn Jahren aufgetretenen Abnahmen nördlicher verbreiteter Vogelarten deuten darauf hin, daß die Klimaveränderung aktuell als der größte Einfluß für Vögel in Europa angesehen werden muß. Insgesamt zeigen die Ergebnisse dieser Arbeit, daß sich der anhaltende Druck auf die Umwelt in erster Linie durch Habitat- und Klimaveränderungen manifestiert.

Detection and Modeling of Boundary Layer Height

This thesis tackles the problem of the automated detection of the atmospheric boundary layer (BL) height, h, from aerosol lidar/ceilometer observations. A new method, the Bayesian Selective Method (BSM), is presented. It implements a Bayesian statistical inference procedure which combines in an statistically optimal way different sources of information. Firstly atmospheric stratification boundaries are located from discontinuities in the ceilometer back-scattered signal. The BSM then identifies the discontinuity edge that has the highest probability to effectively mark the BL height. Information from the contemporaneus physical boundary layer model simulations and a climatological dataset of BL height evolution are combined in the assimilation framework to assist this choice. The BSM algorithm has been tested for four months of continuous ceilometer measurements collected during the BASE:ALFA project and is shown to realistically diagnose the BL depth evolution in many different weather conditions. Then the BASE:ALFA dataset is used to investigate the boundary layer structure in stable conditions. Functions from the Obukhov similarity theory are used as regression curves to fit observed velocity and temperature profiles in the lower half of the stable boundary layer. Surface fluxes of heat and momentum are best-fitting parameters in this exercise and are compared with what measured by a sonic anemometer. The comparison shows remarkable discrepancies, more evident in cases for which the bulk Richardson number turns out to be quite large. This analysis supports earlier results, that surface turbulent fluxes are not the appropriate scaling parameters for profiles of mean quantities in very stable conditions. One of the practical consequences is that boundary layer height diagnostic formulations which mainly rely on surface fluxes are in disagreement to what obtained by inspecting co-located radiosounding profiles.

Landwirtschaftliche Schaderregerprognose mit Hilfe von geographischen Informationssystemen

In der hier vorliegenden Arbeit wurde am Beispiel der Kraut- und Knollenfäule an Kartoffeln Phytophthora infestans und des Kartoffelkäfers Leptinotarsa decemlineata untersucht, ob durch den Einsatz von Geographischen Informationssystemen (GIS) landwirtschaftliche Schader¬reger¬prognosen für jeden beliebigen Kartoffelschlag in Deutschland erstellt werden können. Um dieses Ziel zu erreichen, wurden die Eingangsparameter (Temperatur und relative Luftfeuchte) der Prognosemodelle für die beiden Schaderreger (SIMLEP1, SIMPHYT1, SIMPHYT3 and SIMBLIGHT1) so aufbereitet, dass Wetterdaten flächendeckend für Deutschland zur Verfügung standen. Bevor jedoch interpoliert werden konnte, wurde eine Regionalisierung von Deutschland in Interpolationszonen durchgeführt und somit Naturräume geschaffen, die einen Vergleich und eine Bewertung der in ihnen liegenden Wetterstationen zulassen. Hierzu wurden die Boden-Klima-Regionen von SCHULZKE und KAULE (2000) modifiziert, an das Wetterstationsnetz angepasst und mit 5 bis 10 km breiten Pufferzonen an der Grenze der Interpolationszonen versehen, um die Wetterstationen so häufig wie möglich verwenden zu können. Für die Interpolation der Wetterdaten wurde das Verfahren der multiplen Regression gewählt, weil dieses im Vergleich zu anderen Verfahren die geringsten Abweichungen zwischen interpolierten und gemessenen Daten aufwies und den technischen Anforderungen am besten entsprach. Für 99 % aller Werte konnten bei der Temperaturberechnung Abweichungen in einem Bereich zwischen -2,5 und 2,5 °C erzielt werden. Bei der Berechnung der relativen Luftfeuchte wurden Abweichungen zwischen -12 und 10 % relativer Luftfeuchte erreicht. Die Mittelwerte der Abweichungen lagen bei der Temperatur bei 0,1 °C und bei der relativen Luftfeuchte bei -1,8 %. Zur Überprüfung der Trefferquoten der Modelle beim Betrieb mit interpolierten Wetterdaten wurden Felderhebungsdaten aus den Jahren 2000 bis 2007 zum Erstauftreten der Kraut- und Knollenfäule sowie des Kartoffelkäfers verwendet. Dabei konnten mit interpolierten Wetterdaten die gleichen und auch höhere Trefferquoten erreicht werden, als mit der bisherigen Berechnungsmethode. Beispielsweise erzielte die Berechnung des Erstauftretens von P. infestans durch das Modell SIMBLIGHT1 mit interpolierten Wetterdaten im Schnitt drei Tage geringere Abweichungen im Vergleich zu den Berechnungen ohne GIS. Um die Auswirkungen interpretieren zu können, die durch Abweichungen der Temperatur und der relativen Luftfeuchte entstanden wurde zusätzlich eine Sensitivitätsanalyse zur Temperatur und relativen Luftfeuchte der verwendeten Prognosemodelle durchgeführt. Die Temperatur hatte bei allen Modellen nur einen geringen Einfluss auf das Prognoseergebnis. Veränderungen der relativen Luftfeuchte haben sich dagegen deutlich stärker ausgewirkt. So lag bei SIMBLIGHT1 die Abweichung durch eine stündliche Veränderung der relativen Luftfeuchte (± 6 %) bei maximal 27 Tagen, wogegen stündliche Veränderungen der Temperatur (± 2 °C) eine Abweichung von maximal 10 Tagen ausmachten. Die Ergebnisse dieser Arbeit zeigen, dass durch die Verwendung von GIS mindestens die gleichen und auch höhere Trefferquoten bei Schaderregerprognosen erzielt werden als mit der bisherigen Verwendung von Daten einer nahegelegenen Wetterstation. Die Ergebnisse stellen einen wesentlichen Fortschritt für die landwirtschaftlichen Schaderregerprognosen dar. Erstmals ist es möglich, bundesweite Prognosen für jeden beliebigen Kartoffelschlag zur Bekämpfung von Schädlingen in der Landwirtschaft bereit zu stellen.

Applications of remote sensing and geographic information systems in geomorphological studies: Safaga - El Quseir area, Red Sea coast, Egypt as an example

In den letzten drei Jahrzehnten sind Fernerkundung und GIS in den Geowissenschaften zunehmend wichtiger geworden, um die konventionellen Methoden von Datensammlung und zur Herstellung von Landkarten zu verbessern. Die vorliegende Arbeit befasst sich mit der Anwendung von Fernerkundung und geographischen Informationssystemen (GIS) für geomorphologische Untersuchungen. Durch die Kombination beider Techniken ist es vor allem möglich geworden, geomorphologische Formen im Überblick und dennoch detailliert zu erfassen. Als Grundlagen werden in dieser Arbeit topographische und geologische Karten, Satellitenbilder und Klimadaten benutzt. Die Arbeit besteht aus 6 Kapiteln. Das erste Kapitel gibt einen allgemeinen Überblick über den Untersuchungsraum. Dieser umfasst folgende morphologische Einheiten, klimatischen Verhältnisse, insbesondere die Ariditätsindizes der Küsten- und Gebirgslandschaft sowie das Siedlungsmuster beschrieben. Kapitel 2 befasst sich mit der regionalen Geologie und Stratigraphie des Untersuchungsraumes. Es wird versucht, die Hauptformationen mit Hilfe von ETM-Satellitenbildern zu identifizieren. Angewandt werden hierzu folgende Methoden: Colour Band Composite, Image Rationing und die sog. überwachte Klassifikation. Kapitel 3 enthält eine Beschreibung der strukturell bedingten Oberflächenformen, um die Wechselwirkung zwischen Tektonik und geomorphologischen Prozessen aufzuklären. Es geht es um die vielfältigen Methoden, zum Beispiel das sog. Image Processing, um die im Gebirgskörper vorhandenen Lineamente einwandfrei zu deuten. Spezielle Filtermethoden werden angewandt, um die wichtigsten Lineamente zu kartieren. Kapitel 4 stellt den Versuch dar, mit Hilfe von aufbereiteten SRTM-Satellitenbildern eine automatisierte Erfassung des Gewässernetzes. Es wird ausführlich diskutiert, inwieweit bei diesen Arbeitsschritten die Qualität kleinmaßstäbiger SRTM-Satellitenbilder mit großmaßstäbigen topographischen Karten vergleichbar ist. Weiterhin werden hydrologische Parameter über eine qualitative und quantitative Analyse des Abflussregimes einzelner Wadis erfasst. Der Ursprung von Entwässerungssystemen wird auf der Basis geomorphologischer und geologischer Befunde interpretiert. Kapitel 5 befasst sich mit der Abschätzung der Gefahr episodischer Wadifluten. Die Wahrscheinlichkeit ihres jährlichen Auftretens bzw. des Auftretens starker Fluten im Abstand mehrerer Jahre wird in einer historischen Betrachtung bis 1921 zurückverfolgt. Die Bedeutung von Regentiefs, die sich über dem Roten Meer entwickeln, und die für eine Abflussbildung in Frage kommen, wird mit Hilfe der IDW-Methode (Inverse Distance Weighted) untersucht. Betrachtet werden außerdem weitere, regenbringende Wetterlagen mit Hilfe von Meteosat Infrarotbildern. Genauer betrachtet wird die Periode 1990-1997, in der kräftige, Wadifluten auslösende Regenfälle auftraten. Flutereignisse und Fluthöhe werden anhand von hydrographischen Daten (Pegelmessungen) ermittelt. Auch die Landnutzung und Siedlungsstruktur im Einzugsgebiet eines Wadis wird berücksichtigt. In Kapitel 6 geht es um die unterschiedlichen Küstenformen auf der Westseite des Roten Meeres zum Beispiel die Erosionsformen, Aufbauformen, untergetauchte Formen. Im abschließenden Teil geht es um die Stratigraphie und zeitliche Zuordnung von submarinen Terrassen auf Korallenriffen sowie den Vergleich mit anderen solcher Terrassen an der ägyptischen Rotmeerküste westlich und östlich der Sinai-Halbinsel.

Reproduction and maternal health care among young women in Kenya: geographic and socio-economic determinants

Many factors influence the propensity of young women to seek appropriate maternal healthcare, and they need to be considered when analyzing these women’s reproductive behavior. This study aimed to contribute to the analysis concerning Kenyan young women’s determinants on maternal healthcare-seeking behavior for the 5 years preceding the 2008/9 Kenya Demographic and Health Survey. The specific objectives were to: investigate the individual and contextual variables that may explain maternal healthcare habits; measure the individual, household and community effect on maternal healthcare attitudes in young women; assess the link between young women’s characteristics and the use of facilities for maternal healthcare; find a relationship between young women’s behavior and the community where they live; examine how the role of the local presence of healthcare facilities influences reproductive behavior, and if the specificity of services offered by healthcare facilities affects their inclination to use healthcare facilities, and measure the geographic differences that influence the propensity to seek appropriate maternal healthcare. The analysis of factors associated with maternal healthcare-seeking behavior for young women in Kenya was investigated using multilevel models. We performed three major analyses, which concerned the individual and contextual determinants influencing antenatal care (discussed in Part 6), delivery care (Part 7), and postnatal care (Part 8). Our results show that there is a significant variation in antenatal, delivery and postnatal care between communities, even if the majority of variability is explained by individual characteristics. There are differences at the women’s level on the probability of receiving antenatal care and delivering in a healthcare facility instead of at home. Moreover, community factors and availability of healthcare facilities on the territory are also crucial in influencing young women’s behavior. Therefore, policies addressed to youth’s reproductive health should also consider geographic inequalities and different types of barriers in access to healthcare facilities.