Simultaneous estimation of hunting pressure, harvest and hunter success rates using WinBUGS

The use of hierarchical Bayesian spatial models in the analysis of ecological data is increasingly prevalent. The implementation of these models has been heretofore limited to specifically written software that required extensive programming knowledge to create. The advent of WinBUGS provides access to Bayesian hierarchical models for those without the programming expertise to create their own models and allows for the more rapid implementation of new models and data analysis. This facility is demonstrated here using data collected by the Missouri Department of Conservation for the Missouri Turkey Hunting Survey of 1996. Three models are considered, the first uses the collected data to estimate the success rate for individual hunters at the county level and incorporates a conditional autoregressive (CAR) spatial effect. The second model builds upon the first by simultaneously estimating the success rate and harvest at the county level, while the third estimates the success rate and hunting pressure at the county level. These models are discussed in detail as well as their implementation in WinBUGS and the issues arising therein. Future areas of application for WinBUGS and the latest developments in WinBUGS are discussed as well.

Temporal Characterization of Ossification of the Crania in Australian Subadults: New standards for Age Estimation using Computed Tomography

After attending this presentation, attendees will gain awareness of the ontogeny of cranial maturation, specifically: (1) the fusion timings of primary ossification centers in the basicranium; and (2) the temporal pattern of closure of the anterior fontanelle, to develop new population-specific age standards for medicolegal death investigation of Australian subadults. This presentation will impact the forensic science community by demonstrating the potential of a contemporary forensic subadult Computed Tomography (CT) database of cranial scans and population data, to recalibrate existing standards for age estimation and quantify growth and development of Australian children. This research welcomes a study design applicable to all countries faced with paucity in skeletal repositories. Accurate assessment of age-at-death of skeletal remains represents a key element in forensic anthropology methodology. In Australian casework, age standards derived from American reference samples are applied in light of scarcity in documented Australian skeletal collections. Currently practitioners rely on antiquated standards, such as the Scheuer and Black1 compilation for age estimation, despite implications of secular trends and population variation. Skeletal maturation standards are population specific and should not be extrapolated from one population to another, while secular changes in skeletal dimensions and accelerated maturation underscore the importance of establishing modern standards to estimate age in modern subadults. Despite CT imaging becoming the gold standard for skeletal analysis in Australia, practitioners caution the application of forensic age standards derived from macroscopic inspection to a CT medium, suggesting a need for revised methodologies. Multi-slice CT scans of subadult crania and cervical vertebrae 1 and 2 were acquired from 350 Australian individuals (males: n=193, females: n=157) aged birth to 12 years. The CT database, projected at 920 individuals upon completion (January 2014), comprises thin-slice DICOM data (resolution: 0.5/0.3mm) of patients scanned since 2010 at major Brisbane Childrens Hospitals. DICOM datasets were subject to manual segmentation, followed by the construction of multi-planar and volume rendering cranial models, for subsequent scoring. The union of primary ossification centers of the occipital bone were scored as open, partially closed or completely closed; while the fontanelles, and vertebrae were scored in accordance with two stages. Transition analysis was applied to elucidate age at transition between union states for each center, and robust age parameters established using Bayesian statistics. In comparison to reported literature, closure of the fontanelles and contiguous sutures in Australian infants occur earlier than reported, with the anterior fontanelle transitioning from open to closed at 16.7±1.1 months. The metopic suture is closed prior to 10 weeks post-partum and completely obliterated by 6 months of age, independent of sex. Utilizing reverse engineering capabilities, an alternate method for infant age estimation based on quantification of fontanelle area and non-linear regression with variance component modeling will be presented. Closure models indicate that the greatest rate of change in anterior fontanelle area occurs prior to 5 months of age. This study complements the work of Scheuer and Black1, providing more specific age intervals for union and temporal maturity of each primary ossification center of the occipital bone. For example, dominant fusion of the sutura intra-occipitalis posterior occurs before 9 months of age, followed by persistence of a hyaline cartilage tongue posterior to the foramen magnum until 2.5 years; with obliteration at 2.9±0.1 years. Recalibrated age parameters for the atlas and axis are presented, with the anterior arch of the atlas appearing at 2.9 months in females and 6.3 months in males; while dentoneural, dentocentral and neurocentral junctions of the axis transitioned from non-union to union at 2.1±0.1 years in females and 3.7±0.1 years in males. These results are an exemplar of significant sexual dimorphism in maturation (p<0.05), with girls exhibiting union earlier than boys, justifying the need for segregated sex standards for age estimation. Studies such as this are imperative for providing updated standards for Australian forensic and pediatric practice and provide an insight into skeletal development of this population. During this presentation, the utility of novel regression models for age estimation of infants will be discussed, with emphasis on three-dimensional modeling capabilities of complex structures such as fontanelles, for the development of new age estimation methods.

Online estimation of Allan variance parameters

A new online method is presented for estimation of the angular randomwalk and rate randomwalk coefficients of inertial measurement unit gyros and accelerometers. In the online method, a state-space model is proposed, and recursive parameter estimators are proposed for quantities previously measured from offline data techniques such as the Allan variance method. The Allan variance method has large offline computational effort and data storage requirements. The technique proposed here requires no data storage and computational effort of approximately 100 calculations per data sample.

On-line estimation of Allan variance parameters

A new online method is presented for estimation of the angular random walk and rate random walk coefficients of IMU (inertial measurement unit) gyros and accelerometers. The online method proposes a state space model and proposes parameter estimators for quantities previously measured from off-line data techniques such as the Allan variance graph. Allan variance graphs have large off-line computational effort and data storage requirements. The technique proposed here requires no data storage and computational effort of O(100) calculations per data sample.

Short-data recursive HMM parameter estimation for rapid vision-based aircraft heading estimation

Rapid recursive estimation of hidden Markov Model (HMM) parameters is important in applications that place an emphasis on the early availability of reasonable estimates (e.g. for change detection) rather than the provision of longer-term asymptotic properties (such as convergence, convergence rate, and consistency). In the context of vision- based aircraft (image-plane) heading estimation, this paper suggests and evaluates the short-data estimation properties of 3 recursive HMM parameter estimation techniques (a recursive maximum likelihood estimator, an online EM HMM estimator, and a relative entropy based estimator). On both simulated and real data, our studies illustrate the feasibility of rapid recursive heading estimation, but also demonstrate the need for careful step-size design of HMM recursive estimation techniques when these techniques are intended for use in applications where short-data behaviour is paramount.

Is retinal shape different in Asians and Caucasians? Estimation from peripheral refraction and peripheral eye length methods

Purpose:Race appears to be associated with myopiogenesis, with East Asians showing high myopia prevalence. Considering structural variations in the eye, it is possible that retinal shapes are different between races. The purpose of this study was to quantify and compare retinal shapes between racial groups using peripheral refraction (PR) and peripheral eye lengths (PEL). Methods:A Shin-Nippon SRW5000 autorefractor and a Haag-Streit Lenstar LS900 biometer measured PR and PEL, respectively, along horizontal (H) and vertical (V) fields out to ±35° in 5° steps in 29 Caucasian (CA), 16 South Asian (SA) and 23 East Asian (EA) young adults (spherical equivalent range +0.75D to –5.00D in all groups). Retinal vertex curvature Rv and asphericity Q were determined from two methods: a) PR (Dunne): The Gullstrand-Emsley eye was modified according to participant’s intraocular lengths and anterior cornea curvature. Ray-tracing was performed at each angle through the stop, altering cornea asphericity until peripheral astigmatism matched experimental measurements. Retinal curvature and hence retinal co-ordinate intersection with the chief ray were altered until sagittal refraction matched its measurement. b) PEL: Ray-tracing was performed at each angle through the anterior corneal centre of curvature of the Gullstrand-Emsley eye. Ignoring lens refraction, retinal co-ordinates relative to the fovea were determined from PEL and trigonometry. From sets of retinal co-ordinates, conic retinal shapes were fitted in terms of Rv and Q. Repeated-measures ANOVA were conducted on Rv and Q, and post hoc t-tests with Bonferroni correction were used to compare races. Results:In all racial groups both methods showed greater Rv for the horizontal than for the vertical meridian and greater Rv for myopes than emmetropes. Rv was greater in EA than in CA (P=0.02), with Rv for SA being intermediate and not significantly different from CA and EA. The PEL method provided larger Rv than the PR method: PEL: EA vs CA 87±13 vs 83±11 m-1 (H), 79±13 vs 72±14 m-1 (V); PR: EA vs CA 79±10 vs 67±10 m-1 (H), 71±17 vs 66±12 m-1 (V). Q did not vary significantly with race. Conclusions:Estimates of Rv, but not of Q, varied significantly with race. The greater Rv found in EA than in CA and the comparatively high prevalence rate of myopia in many Asian countries may be related.

Dewatering of microalgal culture for biodiesel production : exploring polymer flocculation and tangential flow filtration

Background: Conventional biodiesel production relies on trans-esterification of lipids extracted from vegetable crops. However, the use of valuable vegetable food stocks as raw material for biodiesel production makes it an unfeasibly expensive process. Used cooking oil is a finite resource and requires extra downstream processing, which affects the amount of biodiesel that can be produced and the economics of the process. Lipids extracted from microalgae are considered an alternative raw material for biodiesel production. This is primarily due to the fast growth rate of these species in a simple aquaculture environment. However, the dilute nature of microalgae culture puts a huge economic burden on the dewatering process especially on an industrial scale. This current study explores the performance and economic viability of chemical flocculation and tangential flow filtration (TFF) for the dewatering of Tetraselmis suecicamicroalgae culture. Results: Results show that TFF concentrates the microalgae feedstock up to 148 times by consuming 2.06 kWh m-3 of energy while flocculation consumes 14.81 kWhm-3 to concentrate the microalgae up to 357 times. Economic evaluation demonstrates that even though TFF has higher initial capital investment than polymer flocculation, the payback period for TFF at the upper extreme ofmicroalgae revenue is ∼1.5 years while that of flocculation is ∼3 years. Conclusion: These results illustrate that improved dewatering levels can be achieved more economically by employing TFF. The performances of these two techniques are also compared with other dewatering techniques.

Engineering cellulose nanofibre suspensions to control filtration resistance and sheet permeability

This study examines and quantifies the effect of adding polyelectrolytes to cellulose nanofibre suspensions on the gel point of cellulose nanofibre suspensions, which is the lowest solids concentration at which the suspension forms a continuous network. The lower the gel point, the faster the drainage time to produce a sheet and the higher the porosity of the final sheet formed. Two new techniques were designed to measure the dynamic compressibility and the drainability of nanocellulose–polyelectrolyte suspensions. We developed a master curve which showed that the independent variable controlling the behaviour of nanocellulose suspensions and its composite is the structure of the flocculated suspension which is best quantified as the gel point. This was independent of the type of polyelectrolyte used. At an addition level of 2 mg/g of nanofibre, a reduction in gel point over 50 % was achieved using either a high molecular weight (13 MDa) linear cationic polyacrylamide (CPAM, 40 % charge), a dendrimer polyethylenimine of high molecular weight of 750,000 Da (HPEI) or even a low molecular weight of 2000 Da (LPEI). There was no significant difference in the minimum gel point achieved, despite the difference in polyelectrolyte morphology and molecular weight. In this paper, we show that the gel point controls the flow through the fibre suspension, even when comparing fibre suspensions with solids content above the gel point. A lower gel point makes it easier for water to drain through the fibre network,reducing the pressure required to achieve a given dewatering rate and reducing the filtering time required to form a wet laid sheet. We further show that the lower gel point partially controls the structure of the wet laid sheet after it is dried. Halving the gel point increased the air permeability of the dry sheet by 37, 46 and 25 %, when using CPAM, HPEI and LPEI, respectively. The resistance to liquid flow was reduced by 74 and 90 %, when using CPAM and LPEI. Analysing the paper formed shows that sheet forming process and final sheet properties can be engineered and controlled by adding polyelectrolytes to the nanofibre suspension.

Reliable estimation via simulation

Let a and s denote the inter arrival times and service times in a GI/GI/1 queue. Let a (n), s (n) be the r.v.s, with distributions as the estimated distributions of a and s from iid samples of a and s of sizes n. Let w be a r.v. with the stationary distribution lr of the waiting times of the queue with input (a, s). We consider the problem of estimating E [w~], tx > 0 and 7r via simulations when (a (n), s (n)) are used as input. Conditions for the accuracy of the asymptotic estimate, continuity of the asymptotic variance and uniformity in the rate of convergence to the estimate are obtained. We also obtain rates of convergence for sample moments, the empirical process and the quantile process for the regenerative processes. Robust estimates are also obtained when an outlier contaminated sample of a and s is provided. In the process we obtain consistency, continuity and asymptotic normality of M-estimators for stationary sequences. Some robustness results for Markov processes are included.

Sediment concentration prediction and statistical evaluation for annual load estimation

We consider the development of statistical models for prediction of constituent concentration of riverine pollutants, which is a key step in load estimation from frequent flow rate data and less frequently collected concentration data. We consider how to capture the impacts of past flow patterns via the average discounted flow (ADF) which discounts the past flux based on the time lapsed - more recent fluxes are given more weight. However, the effectiveness of ADF depends critically on the choice of the discount factor which reflects the unknown environmental cumulating process of the concentration compounds. We propose to choose the discount factor by maximizing the adjusted R-2 values or the Nash-Sutcliffe model efficiency coefficient. The R2 values are also adjusted to take account of the number of parameters in the model fit. The resulting optimal discount factor can be interpreted as a measure of constituent exhaustion rate during flood events. To evaluate the performance of the proposed regression estimators, we examine two different sampling scenarios by resampling fortnightly and opportunistically from two real daily datasets, which come from two United States Geological Survey (USGS) gaging stations located in Des Plaines River and Illinois River basin. The generalized rating-curve approach produces biased estimates of the total sediment loads by -30% to 83%, whereas the new approaches produce relatively much lower biases, ranging from -24% to 35%. This substantial improvement in the estimates of the total load is due to the fact that predictability of concentration is greatly improved by the additional predictors.

Comparative analysis of traffic state estimation ─ Cumulative counts-based and trajectory-based methods

The Macroscopic Fundamental Diagram (MFD) relates space-mean density and flow. Since the MFD represents the area-wide network traffic performance, studies on perimeter control strategies and network-wide traffic state estimation utilising the MFD concept have been reported. Most previous works have utilised data from fixed sensors, such as inductive loops, to estimate the MFD, which can cause biased estimation in urban networks due to queue spillovers at intersections. To overcome the limitation, recent literature reports the use of trajectory data obtained from probe vehicles. However, these studies have been conducted using simulated datasets; limited works have discussed the limitations of real datasets and their impact on the variable estimation. This study compares two methods for estimating traffic state variables of signalised arterial sections: a method based on cumulative vehicle counts (CUPRITE), and one based on vehicles’ trajectory from taxi Global Positioning System (GPS) log. The comparisons reveal some characteristics of taxi trajectory data available in Brisbane, Australia. The current trajectory data have limitations in quantity (i.e., the penetration rate), due to which the traffic state variables tend to be underestimated. Nevertheless, the trajectory-based method successfully captures the features of traffic states, which suggests that the trajectories from taxis can be a good estimator for the network-wide traffic states.

Site classification and estimation of surface level seismic hazard using geophysical data and probabilistic approach

This paper presents the site classification of Bangalore Mahanagar Palike (BMP) area using geophysical data and the evaluation of spectral acceleration at ground level using probabilistic approach. Site classification has been carried out using experimental data from the shallow geophysical method of Multichannel Analysis of Surface wave (MASW). One-dimensional (1-D) MASW survey has been carried out at 58 locations and respective velocity profiles are obtained. The average shear wave velocity for 30 m depth (Vs(30)) has been calculated and is used for the site classification of the BMP area as per NEHRP (National Earthquake Hazards Reduction Program). Based on the Vs(30) values major part of the BMP area can be classified as ``site class D'', and ``site class C'. A smaller portion of the study area, in and around Lalbagh Park, is classified as ``site class B''. Further, probabilistic seismic hazard analysis has been carried out to map the seismic hazard in terms spectral acceleration (S-a) at rock and the ground level considering the site classes and six seismogenic sources identified. The mean annual rate of exceedance and cumulative probability hazard curve for S. have been generated. The quantified hazard values in terms of spectral acceleration for short period and long period are mapped for rock, site class C and D with 10% probability of exceedance in 50 years on a grid size of 0.5 km. In addition to this, the Uniform Hazard Response Spectrum (UHRS) at surface level has been developed for the 5% damping and 10% probability of exceedance in 50 years for rock, site class C and D These spectral acceleration and uniform hazard spectrums can be used to assess the design force for important structures and also to develop the design spectrum.

Estimation of red blood cell lifespan from alveolar carbon monoxide measurements

Measurement of alveolar carbon monoxide (CO) presents a facile technique to estimate the lifespan, L, of red blood cells (RBCs) in vivo. Several recent studies employ this technique and calculate L (in days) using the expression, L = 13.8 (Hb)/P-CO(end), where (Hb) is the concentration (in g/dL) of hemoglobin in blood, and P-CO(end) is the endogenous production of CO (in ppm). Implicit in this calculation is the assumption that the fraction, f, of endogenous CO production due to RBC turnover is a constant equal to 0.7, which yields the expected RBC lifespan, L approximate to 120 days, in normal controls. In anemic patients, however, enhanced RBC turnover may increase f substantially above 0.7. The above expression then overestimates L. Here, we deriv an alternative tive expression, L = 3390[Hb]/322P(CO (end)-110, that accounts explicitly for the dependence of f on the rate of RBC turnover and thereby provides more accurate estimates of L without requiring additional measurements. Using the latter expression, we recalculate L from recent measurements on hepatitis C virus infected patients undergoing treatment with ribavirin. We find that our estimates of L in these patients (39 +/- 13 days) are significantly lower than current estimates (46 +/- 14 days), indicating that ribavirin affects RBC survival more severely than expected from current studies. Our expression for L is simple to employ in a clinical setting and would render the broadly applicable technique of alveolar CO measurement for the estimation of RBC lifespan more accurate.

Estimation of Cavitation Erosion with Incubation Periods and Material Properties

A detailed study of the normalized correlations between the incubation period tc and the properties of various materials tested in a rotating disk device indicates that, at very high intensities, the strength properties influence the duration of tc. The analysis of extensive data from other laboratories for cavitation and liquid impingement erosion also indicates that, while both energy and strength properties influence the duration of tc, the latter ones predominate for a majority of cases. A fatigue-type failure occurs during tc. For estimating the time required to pierce a metallic specimen in a rotating device a relationship tp = 160 tc0.44 is proposed. A detailed study of normalized correlations between erosion resistance (inverse of erosion rate) and tc values of different materials tested in the rotating disk shows that correlations are good. Analysis of data from eight other investigators clearly points out the validity and the usefulness of this type of prediction.

Molecular basis of the kidney filtration barrier : Role of the nephrin protein complex

The kidney filtration barrier consists of fenestrated endothelial cell layer, glomerular basement membrane and slit diaphragm (SD), the specialized junction between glomerular viscelar epithelial cells (podocytes). Podocyte injury is associated with the development of proteinuria, and if not reversed the injury will lead to permanent deterioration of the glomerular filter. The early events are characterized by disruption of the integrity of the SD, but the molecular pathways involved are not fully understood. Congenital nephrotic syndrome of the Finnish type (CNF) is caused by mutations in NPHS1, the gene encoding the SD protein nephrin. Lack of nephrin results in loss of the SD and massive proteinuria beginning before birth. Furthermore, nephrin expression is decreased in acquired human kidney diseases including diabetic nephropathy. This highlights the importance of nephrin and consequently SD in regulating the kidney filtration function. However, the precise molecular mechanism of how nephrin is involved in the formation of the SD is unknown. This thesis work aimed at clarifying the role of nephrin and its interaction partners in the formation of the SD. The purpose was to identify novel proteins that associate with nephrin in order to define the essential molecular complex required for the establishment of the SD. The aim was also to decipher the role of novel nephrin interacting proteins in podocytes. Nephrin binds to nephrin-like proteins Neph1 and Neph2, and to adherens junction protein P-cadherin. These interactions have been suggested to play a role in the formation of the SD. In this thesis work, we identified densin as a novel interaction partner for nephrin. Densin was localized to the SD and it was shown to bind to adherens junction protein beta-catenin. Furthermore, densin was shown to behave in a similar fashion as adherens junction proteins in cell-cell contacts. These results indicate that densin may play a role in cell adhesion and, therefore, may contribute to the formation of the SD together with nephrin and adherens junction proteins. Nephrin was also shown to bind to Neph3, which has been previously localized to the SD. Neph3 and Neph1 were shown to induce cell adhesion alone, whereas nephrin needed to trans-interact with Neph1 or Neph3 from the opposite cell surface in order to make cell-cell contacts. This was associated with the decreased tyrosine phosphorylation of nephrin. These data extend the current knowledge of the molecular composition of the nephrin protein complex at the SD and also provide novel insights of how the SD may be formed. This thesis work also showed that densin was up-regulated in the podocytes of CNF patients. Neph3 was up-regulated in nephrin deficient mouse kidneys, which share similar podocyte alterations and lack of the SD as observed in CNF patients podocytes. These data suggest that densin and Neph3 may have a role in the formation of morphological alterations in podocytes detected in CNF patients. Furthermore, this thesis work showed that deletion of beta-catenin specifically from adult mouse podocytes protected the mice from the development of adriamycin-induced podocyte injury and proteinuria compared to wild-type mice. These results show that beta-catenin play a role in the adriamycin induced podocyte injury. Podocyte injury is a hallmark in many kidney diseases and the changes observed in the podocytes of CNF patient share characteristics with injured podocytes observed in chronic kidney diseases. Therefore, the results obtained in this thesis work suggest that densin, Neph3 and beta-catenin participate in the molecular pathways which result in morphological alterations commonly detected in injured podocytes in kidney diseases.