Statistical estimation of time-trends in two-course crop rotations

An assessment of time-trends in yield parameters is essential to the utilization of data from long-term field trials for the comparison of different crop rotations and input regimes, and the identification of sustainable production systems. The barley-vetch rotation established at Breda in northern Syria has provided the basis for estimation of the time-trends in yield data from selected treatments in a two-course crop rotation trial. The model used for the estimation accounts for the effect of rainfall, a major determinant of each annual yield value, and the first-order autocorrelation structure in the errors arising from the same plot over time. An expression for the minimum number of cycles required to detect a significant time-trend has been obtained. Results from the barley-vetch rotation under two fertilizer regimes have been discussed.     

Venezuelan Rainfall Data Analysed by Using a Bayesian Space–time Model

We consider a set of data from 80 stations in the Venezuelan state of Guárico consisting of accumulated monthly rainfall in a time span of 16 years. The problem of modelling rainfall accumulated over fixed periods of time and recorded at meteorological stations at different sites is studied by using a model based on the assumption that the data follow a truncated and transformed multivariate normal distribution. The spatial correlation is modelled by using an exponentially decreasing correlation function and an interpolating surface for the means. Missing data and dry periods are handled within a Markov chain Monte Carlo framework using latent variables. We estimate the amount of rainfall as well as the probability of a dry period by using the predictive density of the data. We considered a model based on a full second-degree polynomial over the spatial co-ordinates as well as the first two Fourier harmonics to describe the variability during the year. Predictive inferences on the data show very realistic results, capturing the typical rainfall variability in time and space for that region. Important extensions of the model are also discussed.     

Modelling spatial-temporal covariance structures in monocropping barley trials

In long-term trials, not only are individual plot errors correlated over time but there is also a consistent underlying spatial variability in field conditions. The current study sought the most appropriate covariance structure of errors correlated in three dimensions for evaluating the productivity and time-trends in the barley yield data from the monocropping system established in northern Syria. The best spatial-temporal model found reflected the contribution of autocorrelations in spatial and temporal dimensions with estimates varying with the yield variable and location. Compared with a control structure based on independent errors, this covariance structure improved the significance of the fertilizer effect and the interaction with year. Time-trends were estimated in two ways: by accounting the seasonal variable contribution in annual variability (Method 1), which is suitable for detecting significant trends in short data series; and by using the linear component of the orthogonal polynomial on time (year), which is appropriate for long series (Method 2). Method 1 strengthened time-trend detection compared with the method of Jones and Singh [J. Agri. Sci., Cambridge 135 (2000), pp. 251-259] which assumed independence of temporal errors. Most estimates of yield trends over time from fertilizer application were numerically greater than the corresponding linear trends estimated from orthogonal polynomials in time (Method 2), reflecting the effect of accounting for seasonal variables. Grain yield declined over time at the drier site in the absence of nitrogen or phosphorus application, but positive trends were observed fairly generally for straw yield and for grain yield under higher levels of fertilizer inputs. It is suggested that analyses of long-term trials on other crops and cropping systems in other agro-ecological zones could be improved by taking spatial and temporal variability into account in the data evaluation.     

Trend estimation of multivariate time series with controlled smoothness

This paper extends the univariate time series smoothing approach provided by penalized least squares to a multivariate setting, thus allowing for joint estimation of several time series trends. The theoretical results are valid for the general multivariate case, but particular emphasis is placed on the bivariate situation from an applied point of view. The proposal is based on a vector signal-plus-noise representation of the observed data that requires the first two sample moments and specifying only one smoothing constant. A measure of the amount of smoothness of an estimated trend is introduced so that an analyst can set in advance a desired percentage of smoothness to be achieved by the trend estimate. The required smoothing constant is determined by the chosen percentage of smoothness. Closed form expressions for the smoothed estimated vector and its variance-covariance matrix are derived from a straightforward application of generalized least squares, thus providing best linear unbiased estimates for the trends. A detailed algorithm applicable for estimating bivariate time series trends is also presented and justified. The theoretical results are supported by a simulation study and two real applications. One corresponds to Mexican and US macroeconomic data within the context of business cycle analysis, and the other one to environmental data pertaining to a monitored site in Scotland.     

On modelling data from degradation sample paths over time

This paper is mainly concerned with modelling data from degradation sample paths over time. It uses a general growth curve model with Box‐Cox transformation, random effects and ARMA(p, q) dependence to analyse a set of such data. A maximum likelihood estimation procedure for the proposed model is derived and future values are predicted, based on the best linear unbiased prediction. The paper compares the proposed model with a nonlinear degradation model from a prediction point of view. Forecasts of failure times with various data lengths in the sample are also compared.     

Comparison of three-dimensional profiles over time

In this paper, we describe an analysis for data collected on a three-dimensional spatial lattice with treatments applied at the horizontal lattice points. Spatial correlation is accounted for using a conditional autoregressive model. Observations are defined as neighbours only if they are at the same depth. This allows the corresponding variance components to vary by depth. We use the Markov chain Monte Carlo method with block updating, together with Krylov subspace methods, for efficient estimation of the model. The method is applicable to both regular and irregular horizontal lattices and hence to data collected at any set of horizontal sites for a set of depths or heights, for example, water column or soil profile data. The model for the three-dimensional data is applied to agricultural trial data for five separate days taken roughly six months apart in order to determine possible relationships over time. The purpose of the trial is to determine a form of cropping that leads to less moist soils in the root zone and beyond. We estimate moisture for each date, depth and treatment accounting for spatial correlation and determine relationships of these and other parameters over time.     

A flexible and efficient spatial interpolator for radar rainfall estimation

A key challenge in rainfall estimation is spatio-temporal variablility. Weather radars are used to estimate precipitation with high spatial and temporal resolution. Due to the inherent errors in radar estimates, spatial interpolation has been often employed to calibrate the estimates. Kriging is a simple and popular spatial interpolation method, but the method has several shortcomings. In particular, the prediction is quite unstable and often fails to be performed when sample size is small. In this paper, we proposed a flexible and efficient spatial interpolator for radar rainfall estimation, with several advantages over kriging. The method is illustrated using a real-world data set.     

Non-parametric maximum likelihood estimation of interval-censored failure time data subject to misclassification

The paper considers non-parametric maximum likelihood estimation of the failure time distribution for interval-censored data subject to misclassification. Such data can arise from two types of observation scheme; either where observations continue until the first positive test result or where tests continue regardless of the test results. In the former case, the misclassification probabilities must be known, whereas in the latter case, joint estimation of the event-time distribution and misclassification probabilities is possible. The regions for which the maximum likelihood estimate can only have support are derived. Algorithms for computing the maximum likelihood estimate are investigated and it is shown that algorithms appropriate for computing non-parametric mixing distributions perform better than an iterative convex minorant algorithm in terms of time to absolute convergence. A profile likelihood approach is proposed for joint estimation. The methods are illustrated on a data set relating to the onset of cardiac allograft vasculopathy in post-heart-transplantation patients.     

Estimation of land equivalent ratios

The estimation of the land equivalent ratios is proposed to be done by the (sum of) ratios of means of intercrop yield to sole crop yield. The bias and standard error of the estimates are obtained for large samples. Comparisons of the cropping systems have been made on the basis of these estimates and illustrated with field data.     

Use of satellite data in agricultural surveys

A major application of satellite remote sensing is the estimation of the acreage of agricultural crops. The potential for crop yield estimation using satellite remote sensing exists, but research in this area is still in its early stages. In this paper we survey the methodology for using remotely sensed data in agricultural surveys, based primarily on research conducted during the Large Area Crop Inventory Experiment (LACIE) and the follow-on program Agricultural Research and Inventory Surveys Through Aerospace Remote Sensing (AgRISTARS). The data obtained from multispectral scanner (MSS) and thematic mapper (TM) sensors onboard the Landsat series of satellites are described. Approaches for preprocessing, transferring, and modeling these data for understanding the relationship between their temporal behavior and crop growth cycles are discussed. Finally, techniques for crop identification and area and yield estimation are briefly described     

Bayesian comparative study on binary time series

In this paper, we consider the Bayesian analysis of binary time series with different priors, namely normal, Students' t, and Jeffreys prior, and compare the results with the frequentist methods through some simulation experiments and one real data on daily rainfall in inches at Mount Washington, NH. Among Bayesian methods, our results show that the Jeffreys prior perform better in most of the situations for both the simulation and the rainfall data. Furthermore, among weakly informative priors considered, Student's t prior with 7 degrees of freedom fits the data most adequately.     

Joint modeling of mixed skewed continuous and ordinal longitudinal responses: a Bayesian approach

In this paper, a joint model for analyzing multivariate mixed ordinal and continuous responses, where continuous outcomes may be skew, is presented. For modeling the discrete ordinal responses, a continuous latent variable approach is considered and for describing continuous responses, a skew-normal mixed effects model is used. A Bayesian approach using Markov Chain Monte Carlo (MCMC) is adopted for parameter estimation. Some simulation studies are performed for illustration of the proposed approach. The results of the simulation studies show that the use of the separate models or the normal distributional assumption for shared random effects and within-subject errors of continuous and ordinal variables, instead of the joint modeling under a skew-normal distribution, leads to biased parameter estimates. The approach is used for analyzing a part of the British Household Panel Survey (BHPS) data set. Annual income and life satisfaction are considered as the continuous and the ordinal longitudinal responses, respectively. The annual income variable is severely skewed, therefore, the use of the normality assumption for the continuous response does not yield acceptable results. The results of data analysis show that gender, marital status, educational levels and the amount of money spent on leisure have a significant effect on annual income, while marital status has the highest impact on life satisfaction.     

Regression analysis for discrete event history or failure time data

The paper deals with discrete-time regression models to analyze multistate—multiepisode models for event history data or failure time data collected in follow-up studies, retrospective studies, or longitudinal panels. The models are applicable if the events are not dated exactly but only a time interval is recorded. The models include individual specific parameters to account for unobserved heterogeneity. The explantory variables may be time-varying and random with distributions depending on the observed history of the process. Different estimation procedures are considered: Estimation of structural as well as individual specific parameters by maximization of a joint likelihood function, estimation of the structural parameters by maximization of a conditional likelihood function conditioning on a set of sufficient statistics for the individual specific parameters, and estimation of the structural parameters by maximization of a marginal likelihood function assuming that the individual specific parameters follow a distribution. The advantages and limitations of the different approaches are discussed.     

Quantile smoothing in financial time series

Various parametric models have been designed to analyze volatility in time series of financial market data. For maximum likelihood estimation these parametric methods require the assumption of a known conditional distribution. In this paper we examine the conditional distribution of daily DAX returns with the help of nonparametric methods. We use kernel estimators for conditional quantiles resulting from a kernel estimation of conditional distributions. This work was financially supported by the Deutsche Forschungsgemeinschaft     

Space–time calibration of radar rainfall data

Motivated by a specific problem concerning the relationship between radar reflectance and rainfall intensity, the paper develops a space–time model for use in environmental monitoring applications. The model is cast as a high dimensional multivariate state space time series model, in which the cross-covariance structure is derived from the spatial context of the component series, in such a way that its interpretation is essentially independent of the particular set of spatial locations at which the data are recorded. We develop algorithms for estimating the parameters of the model by maximum likelihood, and for making spatial predictions of the radar calibration parameters by using realtime computations. We apply the model to data from a weather radar station in Lancashire, England, and demonstrate through empirical validation the predictive performance of the model.     

Forecasting spring flow time series

Summary The paper deals with a statistical analysis, carried out to define the underlying reason of some of the damage observed in many buildings of a southern Italian town. Engineering considerations, substantiated by specific measurements, attributed them to the lowering of the groundwater table in the area below the building locations. Due to two coinciding events which occurred in the preceding years, i.e. a persistent drought and the start up of a system of wells, it was not possible to define the cause of the former phenomenon. As shutting down the wells could generate additional problems, an accurate picture of the whole situation was necessary, before taking any action. By taking advantage of some fragmentary data belonging to the flow of a spring located in the area and on the basis of the knowledge of the rainfall data recorded in the Italian hydrographic service directory, two models have been developed which reproduce the spring flow time series in relation to the rainfall recorded in the surrounding area. By comparing the spring flow predictions with the actual data it has been possible to highlight the main role played by the wells.     

A note on the estimation methods of crop production at the block level

This paper presents a method of estimation of crop-production statistics at smaller geographical levels like a community development block (generally referred to as a block) to make area-specific plans for agricultural development programmes in India. Using available district-level data on crop yield from crop-cutting experiments and data on auxiliary variables from various administrative sources, a suitable regression model is fitted. The fitted model is then used to predict the crop production at the block level. Some scaled estimators are also developed using predicted estimates. An empirical study is also carried out to judge the merits of the proposed estimators.     

The Bayesian choice of crop variety and fertilizer dose

Recent contributions to the theory of optimizing fertilizer doses in agricultural crop production have introduced Bayesian ideas to incorporate information on crop yield from several environments and on soil nutrients from a soil test, but they have not used a fully Bayesian formulation. We present such a formulation and demonstrate how the resulting Bayes decision procedure can be evaluated in practice by using Markov chain Monte Carlo methods. The approach incorporates expert knowledge of the crop and of regional and local soil conditions and allows a choice of crop variety as well as of fertilizer level. Alternative dose–response functions are expressed in terms of a common interpretable set of parameters to facilitate model comparisons and the specification of prior distributions. The approach is illustrated with a set of yield data from spring barley nitrogen–response trials and is found to be robust to changes in the dose–response function and the prior distribution for indigenous soil nitrogen.     

Additive mixed effect model for recurrent gap time data

Gap times between recurrent events are often of primary interest in medical and observational studies. The additive hazards model, focusing on risk differences rather than risk ratios, has been widely used in practice. However, the marginal additive hazards model does not take the dependence among gap times into account. In this paper, we propose an additive mixed effect model to analyze gap time data, and the proposed model includes a subject-specific random effect to account for the dependence among the gap times. Estimating equation approaches are developed for parameter estimation, and the asymptotic properties of the resulting estimators are established. In addition, some graphical and numerical procedures are presented for model checking. The finite sample behavior of the proposed methods is evaluated through simulation studies, and an application to a data set from a clinic study on chronic granulomatous disease is provided.     

Maximum likelihood estimation of linear continuous time long memory processes with discrete time data

Summary.  We develop a new class of time continuous autoregressive fractionally integrated moving average (CARFIMA) models which are useful for modelling regularly spaced and irregu-larly spaced discrete time long memory data. We derive the autocovariance function of a stationary CARFIMA model and study maximum likelihood estimation of a regression model with CARFIMA errors, based on discrete time data and via the innovations algorithm. It is shown that the maximum likelihood estimator is asymptotically normal, and its finite sample properties are studied through simulation. The efficacy of the approach proposed is demonstrated with a data set from an environmental study.