|
| ||
|
|
||
Actual Observations (actual + observation)
Selected AbstractsApplication of robust procedures for estimation of breeding values in multiple-trait random regression test-day modelJOURNAL OF ANIMAL BREEDING AND GENETICS, Issue 1 2007J. Jamrozik Summary Robust procedures for estimation of breeding values were applied to multiple-trait random regression test-day (TD) model to reduce the influence of outliers on inferences. Robust estimation methods consisted of correcting selected observations (defined as outliers) in the process of solving mixed-model equations in such a way that ,new' observations gave residuals (actual observation minus predicted) within k residual standard deviations for a given day in milk in 305-day lactation. Data were 980 503 TD records on 63 346 Canadian Jersey cows. Milk, fat, protein and somatic cell score in the first three lactations were analysed jointly in the model that included fixed herd-TD effect and regressions within region,age,season of calving, and regressions with random coefficients for animal genetic and permanent environmental effects. All regressions were orthogonal polynomials of order 4. Robust procedures for k = 1.5, 2.0, 2.5, 2.75 and 3.0 were contrasted with the regular best linear unbiased prediction (BLUP) method in terms of numbers and distributions of outliers, and estimated breeding values (EBV) of animals. Distributions of outliers were similar across traits and lactations. Early days in milk (from 5 to 15) were associated with larger frequency of outliers compared with the remaining part of lactation. Several, computationally simple, robust methods (for k > 2.0) reduced the influence of outlier observations in the model and improved the overall model performance. Differences in rankings of animals from robust evaluations were small compared with the regular BLUP method. No clear associations between changes in EBV (rankings) of top animals from different methods and the occurrence of outliers were detected. [source] Improving the analysis of movement data from marked individuals through explicit estimation of observer heterogeneityJOURNAL OF AVIAN BIOLOGY, Issue 1 2010Fränzi Korner-Nievergelt Ring re-encounter data, in particular ring recoveries, have made a large contribution to our understanding of bird movements. However, almost every study based on ring re-encounter data has struggled with the bias caused by unequal observer distribution. Re-encounter probabilities are strongly heterogeneous in space and over time. If this heterogeneity can be measured or at least controlled for, the enormous number of ring re-encounter data collected can be used effectively to answer many questions. Here, we review four different approaches to account for heterogeneity in observer distribution in spatial analyses of ring re-encounter data. The first approach is to measure re-encounter probability directly. We suggest that variation in ring re-encounter probability could be estimated by combining data whose re-encounter probabilities are close to one (radio or satellite telemetry) with data whose re-encounter probabilities are low (ring re-encounter data). The second approach is to measure the spatial variation in re-encounter probabilities using environmental covariates. It should be possible to identify powerful predictors for ring re-encounter probabilities. A third approach consists of the comparison of the actual observations with all possible observations using randomization techniques. We encourage combining such randomisations with ring re-encounter models that we discuss as a fourth approach. Ring re-encounter models are based on the comparison of groups with equal re-encounter probabilities. Together these four approaches could improve our understanding of bird movements considerably. We discuss their advantages and limitations and give directions for future research. [source] Undersampled radial MRI with multiple coils.MAGNETIC RESONANCE IN MEDICINE, Issue 6 2007Iterative image reconstruction using a total variation constraint Abstract The reconstruction of artifact-free images from radially encoded MRI acquisitions poses a difficult task for undersampled data sets, that is for a much lower number of spokes in k-space than data samples per spoke. Here, we developed an iterative reconstruction method for undersampled radial MRI which (i) is based on a nonlinear optimization, (ii) allows for the incorporation of prior knowledge with use of penalty functions, and (iii) deals with data from multiple coils. The procedure arises as a two-step mechanism which first estimates the coil profiles and then renders a final image that complies with the actual observations. Prior knowledge is introduced by penalizing edges in coil profiles and by a total variation constraint for the final image. The latter condition leads to an effective suppression of undersampling (streaking) artifacts and further adds a certain degree of denoising. Apart from simulations, experimental results for a radial spin-echo MRI sequence are presented for phantoms and human brain in vivo at 2.9 T using 24, 48, and 96 spokes with 256 data samples. In comparison to conventional reconstructions (regridding) the proposed method yielded visually improved image quality in all cases. Magn Reson Med 57:1086,1098, 2007. © 2007 Wiley-Liss, Inc. [source] Mesoscale simulations of atmospheric flow and tracer transport in Phoenix, ArizonaMETEOROLOGICAL APPLICATIONS, Issue 3 2006Ge Wang Abstract Large urban centres located within confining rugged or complex terrain can frequently experience episodes of high concentrations of lower atmospheric pollution. Metropolitan Phoenix, Arizona (United States), is a good example, as the general population is occasionally subjected to high levels of lower atmospheric ozone, carbon monoxide and suspended particulate matter. As a result of dramatic but continuous increase in population, the accompanying environmental stresses and the local atmospheric circulation that dominates the background flow, an accurate simulation of the mesoscale pollutant transport across Phoenix and similar urban areas is becoming increasingly important. This is particularly the case in an airshed, such as that of Phoenix, where the local atmospheric circulation is complicated by the complex terrain of the area. Within the study presented here, a three-dimensional time-dependent mesoscale meteorological model (HOTMAC) is employed for simulation of lower-atmospheric flow in Phoenix, for both winter and summer case-study periods in 1998. The specific purpose of the work is to test the model's ability to replicate the atmospheric flow based on the actual observations of the lower-atmospheric wind profile and known physical principles. While a reasonable general agreement is found between the model-produced flow and the observed one, the simulation of near-surface wind direction produces a much less accurate representation of actual conditions, as does the simulation of wind speed over 1,000 metres above the surface. Using the wind and turbulence output from the mesoscale model, likely particle plume trajectories are simulated for the case-study periods using a puff dispersion model (RAPTAD). Overall, the results provide encouragement for the efforts towards accurately simulating the mesoscale transport of lower-atmospheric pollutants in environments of complex terrain. Copyright © 2006 John Wiley & Sons, Ltd. [source] Statistical interpretation of NWP products in IndiaMETEOROLOGICAL APPLICATIONS, Issue 1 2002Parvinder Maini Although numerical weather prediction (NWP) models provide an objective forecast, poor representation of local topography and other features in these models, necessitates statistical interpretation (SI) of NWP products in terms of local weather. The Perfect Prognostic Method (PPM) is one of the techniques for accomplishing this. At the National Center for Medium Range Weather Forecasting, PPM models for precipitation (quantitative, probability, yes/no) and maximum/minimum temperatures are developed for monsoon season by using analyses from the European Centre for Medium-Range Weather Forecasts. The SI forecast is then obtained by using these PPM models and output from the operational NWP model at the Center. Direct model output (DMO) obtained from the NWP model and the SI forecast are verified against the actual observations. The present study shows the verification scores obtained during the 1997 monsoon season for 10 locations in India. The results show that the SI forecast has good skill and is an improvement over DMO. Copyright © 2002 Royal Meteorological Society. [source] A review on the use of the adjoint method in four-dimensional atmospheric-chemistry data assimilationTHE QUARTERLY JOURNAL OF THE ROYAL METEOROLOGICAL SOCIETY, Issue 576 2001K.-Y. Wang Abstract In this paper we review a theoretical formulation of the adjoint method to be used in four-dimensional (4D) chemistry data assimilation. The goal of the chemistry data assimilation is to combine an atmospheric-chemistry model and actual observations to produce the best estimate of the chemistry of the atmosphere. The observational dataset collected during the past decades is an unprecedented expansion of our knowledge of the atmosphere. The exploitation of these data is the best way to advance our understanding of atmospheric chemistry, and to develop chemistry models for chemistry-climate prediction. The assimilation focuses on estimating the state of the chemistry in a chemically and dynamically consistent manner (if the model allows online interactions between chemistry and dynamics). In so doing, we can: produce simultaneous and chemically consistent estimates of all species (including model parameters), observed and unobserved; fill in data voids; test the photochemical theories used in the chemistry models. In this paper, the Hilbert space is first formulated from the geometric structure of the Banach space, followed by the development of the adjoint operator in Hilbert space. The principle of the adjoint method is described, followed by two examples which show the relationship of the gradient of the cost function with respect to the output vector and the gradient of the cost function with respect to the input vector. Applications to chemistry data assimilation are presented for both continuous and discrete cases. The 4D data variational adjoint method is then tested in the assimilation of stratospheric chemistry using a simple catalytic ozone-destruction mechanism, and the test results indicate that the performance of the assimilation method is good. [source] Seismic reflection imaging of active offshore faults in the Gulf of Corinth: their seismotectonic significanceBASIN RESEARCH, Issue 4 2002A. Stefatos ABSTRACT High resolution seismic reflection surveys over one of the most active and rapidly extending regions in the world, the Gulf of Corinth, have revealed that the gulf is a complex asymmetric graben whose geometry varies significantly along its length. A detailed map of the offshore faults in the gulf shows that a major fault system of nine distinct faults limits the basin to the south. The northern Gulf appears to be undergoing regional subsidence and is affected by an antithetic major fault system consisting of eight faults. All these major faults have been active during the Quaternary. Uplifted coastlines along their footwalls, growth fault patterns and thickening of sediment strata toward the fault planes indicate that some of these offshore faults on both sides of the graben are active up to present. Our data ground-truth recent models and provides actual observations of the distribution of variable deformation rates in the Gulf of Corinth. Furthermore they suggest that the offshore faults should be taken into consideration in explaining the high extension rates and the uplift scenarios of the northern Peloponnesos coast. The observed coastal uplift appears to be the result of the cumulative effect of deformation accommodated by more than one fault and therefore, average uplift rates deduced from raised fossil shorelines, should be treated with caution when used to infer individual fault slip rates. Seismic reflection profiling is a vital tool in assessing seismic hazard and basin-formation in areas of active extension. [source] | |||||||||||||