Using individual tracking data to validate the predictions of species distribution models

The authors would like to thank the College of Life Sciences of Aberdeen University and Marine Scotland Science which funded CP's PhD project. Skate tagging experiments were undertaken as part of Scottish Government project SP004. We thank Ian Burrett for help in catching the fish and the other fishermen and anglers who returned tags. We thank José Manuel Gonzalez-Irusta for extracting and making available the environmental layers used as environmental covariates in the environmental suitability modelling procedure. We also thank Jason Matthiopoulos for insightful suggestions on habitat utilization metrics as well as Stephen C.F. Palmer, and three anonymous reviewers for useful suggestions to improve the clarity and quality of the manuscript.

Projected responses of North American grassland songbirds to climate change and habitat availability at their northern range limits in Alberta, Canada

Across North America, grassland songbirds have undergone steep population declines over recent decades, commonly attributed to agricultural intensification. Understanding the potential interactions between the impacts of climate change on the future distributions of these species and the availability of suitable vegetation for nesting can support improved risk assessments and conservation planning for this group of species. We used North American bioclimatic niche models to examine future changes in suitable breeding climate for 15 grassland songbird species at their current northern range limits along the boreal forest–prairie ecotone in Alberta, Canada. Our climate suitability projections, combined with the current distribution of native and tame pasture and cropland in Alberta, suggest that some climate-mediated range expansion of grassland songbirds in Alberta is possible. For six of the eight species projected to experience expansions of suitable climate area in Alberta, this suitable climate partly overlaps the current distribution of suitable land cover. Additionally, for more than half of the species examined, most of the area of currently suitable climate was projected to remain suitable to the end of the century, highlighting the importance of Alberta for the long-term persistence of these species. Some northern prairie-endemic species exhibited substantial projected northward shifts of both the northern and southern edges of the area of suitable climate. Baird’s Sparrow (Ammodramus bairdii) and Sprague’s Pipit (Anthus spragueii), both at-risk grassland specialists, are predicted to have limited climate stability within their current ranges, and their expansion into new areas of suitable climate may be limited by the availability of suitable land cover. Our results highlight the importance of the preservation and restoration of remaining suitable grassland habitat within areas of projected climate stability and beyond current northern range limits for the long-term persistence of many grassland songbird species in the face of climate change.

Data to Decision in a Dynamic Ocean: Robust Species Distribution Models and Spatial Decision Frameworks

Human use of the oceans is increasingly in conflict with conservation of endangered species. Methods for managing the spatial and temporal placement of industries such as military, fishing, transportation and offshore energy, have historically been post hoc; i.e. the time and place of human activity is often already determined before assessment of environmental impacts. In this dissertation, I build robust species distribution models in two case study areas, US Atlantic (Best et al. 2012) and British Columbia (Best et al. 2015), predicting presence and abundance respectively, from scientific surveys. These models are then applied to novel decision frameworks for preemptively suggesting optimal placement of human activities in space and time to minimize ecological impacts: siting for offshore wind energy development, and routing ships to minimize risk of striking whales. Both decision frameworks relate the tradeoff between conservation risk and industry profit with synchronized variable and map views as online spatial decision support systems.

For siting offshore wind energy development (OWED) in the U.S. Atlantic (chapter 4), bird density maps are combined across species with weights of OWED sensitivity to collision and displacement and 10 km2 sites are compared against OWED profitability based on average annual wind speed at 90m hub heights and distance to transmission grid. A spatial decision support system enables toggling between the map and tradeoff plot views by site. A selected site can be inspected for sensitivity to a cetaceans throughout the year, so as to capture months of the year which minimize episodic impacts of pre-operational activities such as seismic airgun surveying and pile driving.

Routing ships to avoid whale strikes (chapter 5) can be similarly viewed as a tradeoff, but is a different problem spatially. A cumulative cost surface is generated from density surface maps and conservation status of cetaceans, before applying as a resistance surface to calculate least-cost routes between start and end locations, i.e. ports and entrance locations to study areas. Varying a multiplier to the cost surface enables calculation of multiple routes with different costs to conservation of cetaceans versus cost to transportation industry, measured as distance. Similar to the siting chapter, a spatial decisions support system enables toggling between the map and tradeoff plot view of proposed routes. The user can also input arbitrary start and end locations to calculate the tradeoff on the fly.

Essential to the input of these decision frameworks are distributions of the species. The two preceding chapters comprise species distribution models from two case study areas, U.S. Atlantic (chapter 2) and British Columbia (chapter 3), predicting presence and density, respectively. Although density is preferred to estimate potential biological removal, per Marine Mammal Protection Act requirements in the U.S., all the necessary parameters, especially distance and angle of observation, are less readily available across publicly mined datasets.

In the case of predicting cetacean presence in the U.S. Atlantic (chapter 2), I extracted datasets from the online OBIS-SEAMAP geo-database, and integrated scientific surveys conducted by ship (n=36) and aircraft (n=16), weighting a Generalized Additive Model by minutes surveyed within space-time grid cells to harmonize effort between the two survey platforms. For each of 16 cetacean species guilds, I predicted the probability of occurrence from static environmental variables (water depth, distance to shore, distance to continental shelf break) and time-varying conditions (monthly sea-surface temperature). To generate maps of presence vs. absence, Receiver Operator Characteristic (ROC) curves were used to define the optimal threshold that minimizes false positive and false negative error rates. I integrated model outputs, including tables (species in guilds, input surveys) and plots (fit of environmental variables, ROC curve), into an online spatial decision support system, allowing for easy navigation of models by taxon, region, season, and data provider.

For predicting cetacean density within the inner waters of British Columbia (chapter 3), I calculated density from systematic, line-transect marine mammal surveys over multiple years and seasons (summer 2004, 2005, 2008, and spring/autumn 2007) conducted by Raincoast Conservation Foundation. Abundance estimates were calculated using two different methods: Conventional Distance Sampling (CDS) and Density Surface Modelling (DSM). CDS generates a single density estimate for each stratum, whereas DSM explicitly models spatial variation and offers potential for greater precision by incorporating environmental predictors. Although DSM yields a more relevant product for the purposes of marine spatial planning, CDS has proven to be useful in cases where there are fewer observations available for seasonal and inter-annual comparison, particularly for the scarcely observed elephant seal. Abundance estimates are provided on a stratum-specific basis. Steller sea lions and harbour seals are further differentiated by ‘hauled out’ and ‘in water’. This analysis updates previous estimates (Williams & Thomas 2007) by including additional years of effort, providing greater spatial precision with the DSM method over CDS, novel reporting for spring and autumn seasons (rather than summer alone), and providing new abundance estimates for Steller sea lion and northern elephant seal. In addition to providing a baseline of marine mammal abundance and distribution, against which future changes can be compared, this information offers the opportunity to assess the risks posed to marine mammals by existing and emerging threats, such as fisheries bycatch, ship strikes, and increased oil spill and ocean noise issues associated with increases of container ship and oil tanker traffic in British Columbia’s continental shelf waters.

Starting with marine animal observations at specific coordinates and times, I combine these data with environmental data, often satellite derived, to produce seascape predictions generalizable in space and time. These habitat-based models enable prediction of encounter rates and, in the case of density surface models, abundance that can then be applied to management scenarios. Specific human activities, OWED and shipping, are then compared within a tradeoff decision support framework, enabling interchangeable map and tradeoff plot views. These products make complex processes transparent for gaming conservation, industry and stakeholders towards optimal marine spatial management, fundamental to the tenets of marine spatial planning, ecosystem-based management and dynamic ocean management.

Pileated Woodpecker Occupancy and the Occurrence and Recruitment of Key Habitat Attributes in a Managed Forest

Thesis (Master's)--University of Washington, 2016-08

As rodovias e as aves de rapina nocturnas relação entre a densidade, o habitat e a mortalidade

As estradas e tráfego inerente surgem como a criação antrópica mais conspícua e penetrante na paisagem natural, sendo considerados os principais agentes causadores de fragmentação e destruição de habitats, assim como representam um obstáculo físico sem precedentes, limitando as relações directas entre os indivíduos, por diminuição da frequência de dispersão e aumento da mortalidade por atropelamento, impedindo o fluxo natural de genes e suscitando o aumento de fenómenos de inbreeding e perda de heterozigotia. Todos os impactes deletérios associados às rodovias são claramente perceptíveis em vertebrados, onde as aves de rapinas nocturnas não são excepção. Uma vez que estas rapaces beneficiam das suas bermas e orlas, como locais de poiso, nidificação ou como corredores de dispersão através da paisagem, são frequentemente vítimas de mortalidade por atropelamento em estradas, sendo esta problemática considerada actualmente uma das mais recentes e importantes formas de mortalidade não natural em rapinas nocturnas e vinculada como um dos maiores problemas de conservação que afecta este grupo. Não obstante, esse mútuo efeito de atracção/repulsa das estradas a estas rapaces, as rodovias criam uma barreira específica que limita a dinâmica, comportamento e densidade populacional das espécies residentes, reconhecendo-se que o isolamento daí resultante, pode comprometer a viabilidade populacional a longo prazo, podendo mesmo conduzir a altos riscos de extinção das populações locais devido a efeitos estocásticos. Mediante esta problemática, este trabalho debruçou-se sobre um único objectivo principal: a avaliação do impacte das rodovias e do tráfego, na densidade das aves de rapina nocturnas. Este estudo foi efectuado na região Alentejana, abrangendo uma área de cariz tipicamente mediterrânico, delimitada pelas localidades de Montemor-o-Novo, Arraiolos e Évora, sendo seccionada por 143 quilómetros de estradas, divididas em autoestrada, rodovias com elevada e reduzida densidade de tráfego. A monitorização das rapinas nocturnas foi conduzida em dois anos amostrais (2005 e 2007), tendo sido focalizada sobretudo em duas espécies de Strigiformes, a Coruja­ do-mato Strix aluco e o Mocho-galego Athene noctua, recorrendo ao uso de playbacks com reprodução de vocalizações de indivíduos conspecíficos. Foram usadas 32 variáveis explicativas integradas em três grandes grupos: variáveis de estrada, métricas da paisagem, uso do solo, tendo sido analiticamente testadas, recorrendo à aplicação de Modelos Lineares Generalizados. Os principais resultados obtidos demonstram que as variáveis de estrada, aliadas à densidade de tráfego e ruído inerente à sua circulação, são provavelmente, responsáveis por um comportamento de repulsa das espécies de aves de rapina nocturnas em estudo, apresentando estas densidades mais elevadas longe de áreas antropicamente perturbadas e, portanto, de menor qualidade que se encontram adjacentes às rodovias. Todavia a presença de habitat favorável a estas rapaces é provavelmente o descritor com maior poder estatístico no que concerne à sua distribuição e densidade, sendo os montados densos e a presença de zonas agrícolas de sequeiro, positivamente correlacionadas, respectivamente com a densidade de casais reprodutores de Coruja-do-mato e Mocho-galego. Mediante a observação dos resultados será veemente a aplicação de medidas de mitigação específicas, que fundamentalmente considerem o afastamento dos efectivos populacionais longe das estradas e tráfego, conservando e assegurando as características estruturais, requisitos e qualidade dos habitats, de modo a incrementar e garantir a viabilidade e densidade das populações, fidelizando a territorialidade e permanência destas aves nestas áreas. Adicionalmente deverá investir-se na conectividade das manchas de paisagem fragmentada pelas rodovias, criando opções de conservação estratégicas, em zonas ecologicamente mais sensíveis, que não somente minimizem o efeito de repulsa reconhecido nestas aves, mas também os níveis de mortalidade por atropelamento, tornando os ecossistemas mais funcionais para a sobrevivência destes rapaces. ABSTRACT; Roads and traffic are the most conspicuous and pervasive human creation, being the great responsible for fragmentation and habitats destruction, reducing animal movement through landscape, which implies decrease of gene flow and loss of variability that can fragmented populations, thus reducing their sizes and densities. All deleterious impacts associated with roads are clearly visible in vertebrates, where owls aren't exception, being frequent victims of road mortality, since they can use roadside habitats and edges for hunting, nesting or dispersal corridors through the landscape, being nowadays one of the most recent and important causes of non­natural mortality in owls and has been recognized as one of the largest conservation problems affecting this group. However, the attractive and avoidance effect of roads and his edges on owls creates a barrier effect that limits dynamics, behaviour and breeding density of resident species, recognizing that possible isolation, could compromise populations survival, make them more vulnerable to high risks of local extinction due to stochastic effects. Despite that, several authors suggest that owls use roads to hunt, as marginal habitats, or for navigation corridors through the landscape ln this context, the major aims of this study was to verify if there are negative effects of roads on the density of owls, considering traffic as an influencing factor. This study was conducted in Alentejo, covering a typical Mediterranean area, including three main localities Montemor-o-Novo, Évora and Arraiolos, being sectioned by 143 km of roads, including highway, roads with high traffic density, and the remaining with low traffic density. The owl census was conducted in two sampling years (2005 and 2007) and was focused mainly on Little owl Athene noctua and Tawny Owl Strix aluco species, using the playback technique, with conspecific calls. We used 32 explanatory variables, mainly included in three groups: road variables, landscape metrics and land use, having been analytically tested, with application of Generalised Linear Models. The main results show that noisy roads with high traffic density are probably the most responsible for the avoidance behaviour of owls, under the study area, showing density depression near high anthropogenic disturbed areas adjacent to roads. However, the presence of habitat quality to these birds is probably the descriptor with greater statistical power, considering its distribution and density, with the dense oak woodland and croplands and arable lands, positively correlated, respectively with Tawny owl and Little owl density. ln consequence, the great conservation effort should be done in order to keep breeding populations away from roads and traffic, ensuring the structural features, requirements and quality of its habitats in order to enhance and ensure the viability and density of owl's populations in these areas. ln addition, it is important to invest in connectivity between roadside fragmented patches, creating strategic conservation options, in sensitive areas, which minimize the avoidance effect, recognized in owls, but also road-kill levels, making ecosystems more functional to survival of these top predators.

Can we use niche models of indicator species to predict the distribution of endangered communities?

European-wide conservation policies are based on the identification of priority habitats. However, research on conservation biogeography often relies on the results and projections of species distribution models to assess species' vulnerability to global change. We assess whether the distribution and structure of threatened communities can be predicted by the suitability of the environmental conditions for their indicator species. We present some preliminary results elucidating if using species distribution models of indicator species at a regional scale is a valid approach to predict these endangered communities. Dune plant assemblages, affected by severe conditions, are excellent models for studying possible interactions among their integrating species and the environment. We use data from an extensive survey of xerophytic inland sand dune scrub communities from Portugal, one of the most threatened habitat types of Europe. We identify indicator shrub species of different types of communities, model their geographical response to the environment, and evaluate whether the output of these niche models are able to predict the distribution of each type of community in a different region.

Environmental variables, habitat discontinuity and life history shaping the genetic structure of Pomatoschistus marmoratus

Coastal lagoons are semi-isolated ecosystems exposed to wide fluctuations of environmental conditions and showing habitat fragmentation. These features may play an important role in separating species into different populations, even at small spatial scales. In this study, we evaluate the concordance between mitochondrial (previous published data) and nuclear data analyzing the genetic variability of Pomatoschistus marmoratus in five localities, inside and outside the Mar Menor coastal lagoon (SE Spain) using eight microsatellites. High genetic diversity and similar levels of allele richness were observed across all loci and localities, although significant genic and genotypic differentiation was found between populations inside and outside the lagoon. In contrast to the FST values obtained from previous mitochondrial DNA analyses (control region), the microsatellite data exhibited significant differentiation among samples inside the Mar Menor and between lagoonal and marine samples. This pattern was corroborated using Cavalli-Sforza genetic distances. The habitat fragmentation inside the coastal lagoon and among lagoon and marine localities could be acting as a barrier to gene flow and contributing to the observed genetic structure. Our results from generalized additive models point a significant link between extreme lagoonal environmental conditions (mainly maximum salinity) and P. marmoratus genetic composition. Thereby, these environmental features could be also acting on genetic structure of coastal lagoon populations of P. marmoratus favoring their genetic divergence. The mating strategy of P. marmoratus could be also influencing our results obtained from mitochondrial and nuclear DNA. Therefore, a special consideration must be done in the selection of the DNA markers depending on the reproductive strategy of the species.

Influência da gestão do habitat na distribuição e abundância do coelho-bravo no Sítio Natura 2000 Monchique

O coelho-bravo, devido à sua importância ecológica e económica, tem sido alvo de diversos planos de gestão e vários esforços têm sido empreendidos no sentido de contrariar o decréscimo das suas populações. Este estudo foi realizado em três zonas de caça do Sítio Monchique e o principal objectivo é determinar se as medidas de gestão implementadas influenciam a distribuição e abundância da espécie na área de estudo. A abundância relativa foi interpolada com o método "Inverso do Peso da Distância" {IDW), e as relações entre presença de coelho e os descritores ambientais foram analisadas através de Modelos Lineares Generalizados (GLM). Os resultados da modelação estatística mostraram que as medidas de melhoria de habitat parecem ter sido determinantes para um aumento da área de distribuição do coelho-bravo nos locais intervencionados. São propostas novas medidas de gestão, cujo objectivo será promover a continuação do aumento da ocorrência e abundância da espécie neste local. /ABSTRACT: The wild rabbit, due to its ecological and economical role, has been the target of several management plans and considerable efforts have been made to enhance its populations. This study was held in three game estates located inside Monchique Natura 2000. Site and aims to determine if the habitat management actions implemented in the study area influence rabbit distribution and abundance. The relative abundance was interpolated to all study area with lnverse Distance Weight method {IDW} and the relationships between rabbit presence and the environmental descriptors were evaluated with Generalized Linear Models (GLM). The results of the statistical modelling showed that the management actions seem to have contributed significantly to an enhancement on the rabbit occurrence in the studied game estates. Several new management actions are proposed with the aim to continue to increase rabbit occurrence and abundance in this site.

Evaluating Karner Blue Butterfly (Lycaeides Melissa Samuelis Nabokov) Habitat Selection in the State of Wisconsin, USA

The federally endangered Karner blue butterfly (Lycaeides melissa samuelis Nabokov) persists in rare oak/pine grassland communities spanning across the Great Lakes region, relying on host plant wild blue lupine (Lupinus perennis). Conservation efforts since 1992 have led to the development of several programs that restore and monitor habitat. This study aims to evaluate Karner blue habitat selection in the state of Wisconsin and develop high-resolution tools for use in conservation efforts. Spatial predictive models developed during this study accurately predicted potential habitat across state properties based on soils and canopy cover, and identified ~51-100% of Karner blue occurrences based on lupine and shrub/tree cover, and focal nectar plant abundance. When evaluated relative to American bison (Bison bison), Karner blues and lupine were more likely to occur in areas of low disturbance, but aggregated where bison were recently present in areas of moderate/high disturbance. Lupine C:N ratio increased relative to cover of shrubs/trees and focal nectar plant abundance and decreased relative to cover of groundlitter. Karner blue density increased with lupine C:N ratio, decreased with nitrogen content, and was not related to phenolic levels. We strongly suggest that areas of different soil textures must be managed differently and that maintenance techniques should generate a mix of shrubs/tree cover (10-45%), groundlitter cover (~10-40%), >5% cover of lupine, and establish an abundance of focal nectar plants. This study provides unique tools for use in conservation and should aid in focusing management efforts and recovery of this species.

Combining multispecies home range and distribution models aids assessment of MPA effectiveness

Marine protected areas (MPAs) are today's most important tools for the spatial management and conservation of marine species. Yet, the true protection that they provide to individual fish is unknown, leading to uncertainty associated with MPA effectiveness. In this study, conducted in a recently established coastal MPA in Portugal, we combined the results of individual home range estimation and population distribution models for 3 species of commercial importance and contrasting life histories to infer (1) the size of suitable areas where they would be fully protected and (2) the vulnerability to fishing mortality of each species. Results show that the relationship between MPA size and effective protection is strongly modulated by both the species' home range and the distribution of suitable habitat inside and outside the MPA. This approach provides a better insight into the true potential of MPAs in effectively protecting marine species, since it can reveal the size and location of the areas where protection is most effective and a clear, quantitative estimation of the vulnerability to fishing throughout an entire MPA.

The mesophotic zone of the Mediterranean Sea: spatial definition, biodiversity, and predictive models

The mesophotic zone is frequently defined as ranging between 30-40 and 150 m depth. However, these borders are necessarily imprecise due to variations in the penetration of light along the water column related to local factors. Moreover, density of data on mesophotic ecosystems vary along geographical distance, with temperate latitudes largely less explored than tropical situations. This is the case of the Mediterranean Sea, where information on mesophotic ecosystems is largely lower with respect to tropical situations. The lack of a clear definition of the borders of the mesophotic zone may represent a problem when information must be transferred to the policy that requires a coherent spatial definition to plan proper management and conservation measures. The present thesis aims at providing information on the spatial definition of the mesophotic zone in the Mediterranean Sea, its biodiversity and distribution of its ecosystems. The first chapter analyzes information on mesophotic ecosystems in the Mediterranean Sea to identify gaps in the literature and map the mesophotic zone in the Mediterranean Sea using light penetration estimated from satellite data. In the second chapter, different visual techniques to study mesophotic ecosystems are compared to identify the best analytical method to estimate diversity and habitat extension. In the third chapter, a set of Remotely Operated vehicles (ROV) surveys performed on mesophotic assemblages in the Mediterranean Sea are analyzed to describe their taxonomic and functional diversity and environmental factors influencing their structure. A Habitat Suitability Model is run in the fourth chapter to map the distribution of areas suitable for the presence of deep-water oyster reefs in the Adriatic-Ionian area. The fifth chapter explores the mesophotic zone in the northern Gulf of Mexico providing its spatial and vertical extension of the mesophotic zone and information on the diversity associated with mesophotic ecosystems.

Characterization Of Microsatellite Markers Developed From Prosopis Rubriflora And Prosopis Ruscifolia (leguminosae - Mimosoideae), Legume Species That Are Used As Models For Genetic Diversity Studies In Chaquenian Areas Under Anthropization In South America.

Prosopis rubriflora and Prosopis ruscifolia are important species in the Chaquenian regions of Brazil. Because of the restriction and frequency of their physiognomy, they are excellent models for conservation genetics studies. The use of microsatellite markers (Simple Sequence Repeats, SSRs) has become increasingly important in recent years and has proven to be a powerful tool for both ecological and molecular studies. In this study, we present the development and characterization of 10 new markers for P. rubriflora and 13 new markers for P. ruscifolia. The genotyping was performed using 40 P. rubriflora samples and 48 P. ruscifolia samples from the Chaquenian remnants in Brazil. The polymorphism information content (PIC) of the P. rubriflora markers ranged from 0.073 to 0.791, and no null alleles or deviation from Hardy-Weinberg equilibrium (HW) were detected. The PIC values for the P. ruscifolia markers ranged from 0.289 to 0.883, but a departure from HW and null alleles were detected for certain loci; however, this departure may have resulted from anthropic activities, such as the presence of livestock, which is very common in the remnant areas. In this study, we describe novel SSR polymorphic markers that may be helpful in future genetic studies of P. rubriflora and P. ruscifolia.

Experiences From The Brazilian Atlantic Forest: Ecological Findings And Conservation Initiatives.

The Brazilian Atlantic Forest hosts one of the world's most diverse and threatened tropical forest biota. In many ways, its history of degradation describes the fate experienced by tropical forests around the world. After five centuries of human expansion, most Atlantic Forest landscapes are archipelagos of small forest fragments surrounded by open-habitat matrices. This 'natural laboratory' has contributed to a better understanding of the evolutionary history and ecology of tropical forests and to determining the extent to which this irreplaceable biota is susceptible to major human disturbances. We share some of the major findings with respect to the responses of tropical forests to human disturbances across multiple biological levels and spatial scales and discuss some of the conservation initiatives adopted in the past decade. First, we provide a short description of the Atlantic Forest biota and its historical degradation. Secondly, we offer conceptual models describing major shifts experienced by tree assemblages at local scales and discuss landscape ecological processes that can help to maintain this biota at larger scales. We also examine potential plant responses to climate change. Finally, we propose a research agenda to improve the conservation value of human-modified landscapes and safeguard the biological heritage of tropical forests.

Censored Linear Regression Models For Irregularly Observed Longitudinal Data Using The Multivariate-t Distribution.

In acquired immunodeficiency syndrome (AIDS) studies it is quite common to observe viral load measurements collected irregularly over time. Moreover, these measurements can be subjected to some upper and/or lower detection limits depending on the quantification assays. A complication arises when these continuous repeated measures have a heavy-tailed behavior. For such data structures, we propose a robust structure for a censored linear model based on the multivariate Student's t-distribution. To compensate for the autocorrelation existing among irregularly observed measures, a damped exponential correlation structure is employed. An efficient expectation maximization type algorithm is developed for computing the maximum likelihood estimates, obtaining as a by-product the standard errors of the fixed effects and the log-likelihood function. The proposed algorithm uses closed-form expressions at the E-step that rely on formulas for the mean and variance of a truncated multivariate Student's t-distribution. The methodology is illustrated through an application to an Human Immunodeficiency Virus-AIDS (HIV-AIDS) study and several simulation studies.

Augmented Mixed Models For Clustered Proportion Data.

Often in biomedical research, we deal with continuous (clustered) proportion responses ranging between zero and one quantifying the disease status of the cluster units. Interestingly, the study population might also consist of relatively disease-free as well as highly diseased subjects, contributing to proportion values in the interval [0, 1]. Regression on a variety of parametric densities with support lying in (0, 1), such as beta regression, can assess important covariate effects. However, they are deemed inappropriate due to the presence of zeros and/or ones. To evade this, we introduce a class of general proportion density, and further augment the probabilities of zero and one to this general proportion density, controlling for the clustering. Our approach is Bayesian and presents a computationally convenient framework amenable to available freeware. Bayesian case-deletion influence diagnostics based on q-divergence measures are automatic from the Markov chain Monte Carlo output. The methodology is illustrated using both simulation studies and application to a real dataset from a clinical periodontology study.