Estimating equations with nonignorably missing response data

Troxel, Lipsitz, and Brennan (1997, Biometrics 53, 857-869) considered parameter estimation from survey data with nonignorable nonresponse and proposed weighted estimating equations to remove the biases in the complete-case analysis that ignores missing observations. This paper suggests two alternative modifications for unbiased estimation of regression parameters when a binary outcome is potentially observed at successive time points. The weighting approach of Robins, Rotnitzky, and Zhao (1995, Journal of the American Statistical Association 90, 106-121) is also modified to obtain unbiased estimating functions. The suggested estimating functions are unbiased only when the missingness probability is correctly specified, and misspecification of the missingness model will result in biases in the estimates. Simulation studies are carried out to assess the performance of different methods when the covariate is binary or normal. For the simulation models used, the relative efficiency of the two new methods to the weighting methods is about 3.0 for the slope parameter and about 2.0 for the intercept parameter when the covariate is continuous and the missingness probability is correctly specified. All methods produce substantial biases in the estimates when the missingness model is misspecified or underspecified. Analysis of data from a medical survey illustrates the use and possible differences of these estimating functions.

Short-term responses of leaf growth rate to water deficit scale up to whole-plant and crop levels: An integrated modelling approach in maize.

Physiological and genetic studies of leaf growth often focus on short-term responses, leaving a gap to whole-plant models that predict biomass accumulation, transpiration and yield at crop scale. To bridge this gap, we developed a model that combines an existing model of leaf 6 expansion in response to short-term environmental variations with a model coordinating the development of all leaves of a plant. The latter was based on: (1) rates of leaf initiation, appearance and end of elongation measured in field experiments; and (2) the hypothesis of an independence of the growth between leaves. The resulting whole-plant leaf model was integrated into the generic crop model APSIM which provided dynamic feedback of environmental conditions to the leaf model and allowed simulation of crop growth at canopy level. The model was tested in 12 field situations with contrasting temperature, evaporative demand and soil water status. In observed and simulated data, high evaporative demand reduced leaf area at the whole-plant level, and short water deficits affected only leaves developing during the stress, either visible or still hidden in the whorl. The model adequately simulated whole-plant profiles of leaf area with a single set of parameters that applied to the same hybrid in all experiments. It was also suitable to predict biomass accumulation and yield of a similar hybrid grown in different conditions. This model extends to field conditions existing knowledge of the environmental controls of leaf elongation, and can be used to simulate how their genetic controls flow through to yield.

Response of minimally processed carambola to chemical treatments and low-oxygen atmospheres

Fresh-cut carambola (Averrhoa carambola L.) has limited marketability due to cut-surface browning. The effect of chemical treatments (ascorbic acid, citric acid and Ca-EDTA), controlled atmosphere (0.4-20.3% O2) and the association of these processes was investigated. Post-cutting dip and low-oxygen atmospheres did not prevent discoloration or improve sensory and physicochemical parameters. However, ascorbic acid (0.5% and 1%) dips reduced polyphenol oxidase (PPO) activity during storage at 4.5 °C, with 1% ascorbic acid inducing the lowest activity. Although cut-surface browning of 'Maha' slices was not relevant, carambola slices treated with 1% ascorbic acid in association with 0.4% oxygen did not present significant browning or loss of visual quality for up to 12 days, 3 days longer than low oxygen alone (0.4% O2), thus, their quality can be significantly improved by combining both treatments.

A fuzzy dynamic flood routing model for natural channels

A fuzzy dynamic flood routing model (FDFRM) for natural channels is presented, wherein the flood wave can be approximated to a monoclinal wave. This study is based on modification of an earlier published work by the same authors, where the nature of the wave was of gravity type. Momentum equation of the dynamic wave model is replaced by a fuzzy rule based model, while retaining the continuity equation in its complete form. Hence, the FDFRM gets rid of the assumptions associated with the momentum equation. Also, it overcomes the necessity of calculating friction slope (S-f) in flood routing and hence the associated uncertainties are eliminated. The fuzzy rule based model is developed on an equation for wave velocity, which is obtained in terms of discontinuities in the gradient of flow parameters. The channel reach is divided into a number of approximately uniform sub-reaches. Training set required for development of the fuzzy rule based model for each sub-reach is obtained from discharge-area relationship at its mean section. For highly heterogeneous sub-reaches, optimized fuzzy rule based models are obtained by means of a neuro-fuzzy algorithm. For demonstration, the FDFRM is applied to flood routing problems in a fictitious channel with single uniform reach, in a fictitious channel with two uniform sub-reaches and also in a natural channel with a number of approximately uniform sub-reaches. It is observed that in cases of the fictitious channels, the FDFRM outputs match well with those of an implicit numerical model (INM), which solves the dynamic wave equations using an implicit numerical scheme. For the natural channel, the FDFRM Outputs are comparable to those of the HEC-RAS model.

Long-term response of corn to limited irrigation and crop rotations

Dwindling water supplies for irrigation are prompting alternative management choices by irrigators. Limited irrigation, where less water is applied than full crop demand, may be a viable approach. Application of limited irrigation to corn was examined in this research. Corn was grown in crop rotations with dryland, limited irrigation, or full irrigation management from 1985 to 1999. Crop rotations included corn following corn (continuous corn), corn following wheat, followed by soybean (wheat-corn-soybean), and corn following soybean (corn-soybean). Full irrigation was managed to meet crop evapotranspiration requirements (ETc). Limited irrigation was managed with a seasonal target of no more than 150 mm applied. Precipitation patterns influenced the outcomes of measured parameters. Dryland yields had the most variation, while fully irrigated yields varied the least. Limited irrigation yields were 80% to 90%> of fully irrigated yields, but the limited irrigation plots received about half the applied water. Grain yields were significantly different among irrigation treatments. Yields were not significantly different among rotation treatments for all years and water treatments. For soil water parameters, more statistical differences were detected among the water management treatments than among the crop rotation treatments. Economic projections of these management practices showed that full irrigation produced the most income if water was available. Limited irrigation increased income significantly from dryland management.

PMU measurement based dynamic load modeling using SVC devices in online environment

Online dynamic load modeling has become possible with the availability of Static Voltage Compensator (SVC) and Phasor Measurement Unit (PMU) devices. The power of the load response to the small random bounded voltage fluctuations caused from SVC can be measured by PMU for modelling purposes. The aim of this paper is to illustrate the capability of identifying an aggregated load model from high voltage substation level in the online environment. The induction motor is used as the main test subject since it contributes the majority of the dynamic loads. A test system representing simple electromechanical generator model serving dynamic loads through the transmission network is used to verify the proposed method. Also, dynamic load with multiple induction motors are modeled to achieve a better realistic load representation.

Probabilistic Seismic Design of Pile Foundations in Non Liquefiable Soil by Response Spectrum Approach

The behavior of pile foundations in non liquefiable soil under seismic loading is considerably influenced by the variability in the soil and seismic design parameters. Hence, probabilistic models for the assessment of seismic pile design are necessary. Deformation of pile foundation in non liquefiable soil is dominated by inertial force from superstructure. The present study considers a pseudo-static approach based on code specified design response spectra. The response of the pile is determined by equivalent cantilever approach. The soil medium is modeled as a one-dimensional random field along the depth. The variability associated with undrained shear strength, design response spectrum ordinate, and superstructure mass is taken into consideration. Monte Carlo simulation technique is adopted to determine the probability of failure and reliability indices based on pile failure modes, namely exceedance of lateral displacement limit and moment capacity. A reliability-based design approach for the free head pile under seismic force is suggested that enables a rational choice of pile design parameters.

Optimal dynamic inversion-based semi-active control of benchmark bridge using MR dampers

A new two-stage state feedback control design approach has been developed to monitor the voltage supplied to magnetorheological (MR) dampers for semi-active vibration control of the benchmark highway bridge. The first stage contains a primary controller, which provides the force required to obtain a desired closed-loop response of the system. In the second stage, an optimal dynamic inversion (ODI) approach has been developed to obtain the amount of voltage to be supplied to each of the MR dampers such that it provides the required force prescribed by the primary controller. ODI is formulated by optimization with dynamic inversion, such that an optimal voltage is supplied to each damper in a set. The proposed control design has been simulated for both phase-I and phase-II study of the recently developed benchmark highway bridge problem. The efficiency of the proposed controller is analyzed in terms of the performance indices defined in the benchmark problem definition. Simulation results demonstrate that the proposed approach generally reduces peak response quantities over those obtained from the sample semi-active controller, although some response quantities have been seen to be increasing. Overall, the proposed control approach is quite competitive as compared with the sample semi-active control approach.

The transient response of certain third-ordernon-linear systems

In this paper a method of solving certain third-order non-linear systems by using themethod of ultraspherical polynomial approximation is proposed. By using the method of variation of parameters the third-order equation is reduced to three partial differential equations. Instead of being averaged over a cycle, the non-linear functions are expanded in ultraspherical polynomials and with only the constant term retained, the equations are solved. The results of the procedure are compared with the numerical solutions obtained on a digital computer. A degenerate third-order system is also considered and results obtained for the above system are compared with numerical results obtained on the digital computer. There is good agreement between the results obtained by the proposed method and the numerical solution obtained on digital computer.

Quantifying the response of cotton production in eastern Australia to climate change

Abstract The paper evaluates the effect of future climate change (as per the CSIRO Mk3.5 A1FI future climate projection) on cotton yield in Southern Queensland and Northern NSW, eastern Australia by using of the biophysical simulation model APSIM (Agricultural Production Systems sIMulator). The simulations of cotton production show that changes in the influential meteorological parameters caused by climate change would lead to decreased future cotton yields without the effect of CO2 fertilisation. By 2050 the yields would decrease by 17 %. Including the effects of CO2 fertilisation ameliorates the effect of decreased water availability and yields increase by 5.9 % by 2030, but then decrease by 3.6 % in 2050. Importantly, it was necessary to increase irrigation amounts by almost 50 % to maintain adequate soil moisture levels. The effect of CO2 was found to have an important positive impact of the yield in spite of deleterious climate change. This implies that the physiological response of plants to climate change needs to be thoroughly understood to avoid making erroneous projections of yield and potentially stifling investment or increasing risk.

Growth analysis and biomass partitioning of Cyperus iria in response to rice planting density and nitrogen rate

Cyperus iria is a weed of rice with widespread occurrence throughout the world. Because of concerns about excessive and injudicious use of herbicides, cultural weed management approaches that are safe and economical are needed. Developing such approaches will require a better understanding of weed biology and ecology, as well as of weed response to increases in crop density and nutrition. Knowledge of the effects of nitrogen (N) fertilizer on crop-weed competitive interactions could also help in the development of integrated weed management strategies. The present study was conducted in a screenhouse to determine the effects of rice planting density (0, 5, 10, and 20 plants pot−1) and N rate (0, 50, 100, and 150 kg ha−1) on the growth of C. iria. Tiller number per plant decreased by 73–88%, leaf number by 85–94%, leaf area by 85–98%, leaf biomass by 92–99%, and inflorescence biomass by 96–99% when weed plants were grown at 20 rice plants pot−1 (i.e., 400 plants m−2) compared with weed plants grown alone. All of these parameters increased when N rates were increased. On average, weed biomass increased by 118–389% and rice biomass by 121–275% with application of 50–150 kg N ha−1, compared to control. Addition of N favored weed biomass production relative to rice biomass. Increased N rates reduced the root-to-shoot weight ratio of C. iria. Rice interference reduced weed growth and biomass and completely suppressed C. iria when no N was applied at high planting densities (i.e., 20 plants pot−1). The weed showed phenotypic plasticity in response to N application, and the addition of N increased the competitive ability of the weed over rice at densities of 5 and 10 rice plants pot−1 compared with 20 plants pot−1. The results of the present study suggest that high rice density (i.e., 400 plants m−2) can help suppress C. iria growth even at high N rates (150 kg ha−1).

Dynamic Instabilities in Slender Space Launch Vehicles under Propulsive Thrust and Aerodynamic Forces

A mechanics based linear analysis of the problem of dynamic instabilities in slender space launch vehicles is undertaken. The flexible body dynamics of the moving vehicle is studied in an inertial frame of reference, including velocity induced curvature effects, which have not been considered so far in the published literature. Coupling among the rigid-body modes, the longitudinal vibrational modes and the transverse vibrational modes due to asymmetric lifting-body cross-section are considered. The model also incorporates the effects of aerodynamic forces and the propulsive thrust of the vehicle. The effects of the coupling between the combustion process (mass variation, developed thrust etc.) and the variables involved in the flexible body dynamics (displacements and velocities) are clearly brought out. The model is one-dimensional, and it can be employed to idealised slender vehicles with complex shapes. Computer simulations are carried out using a standard eigenvalue problem within h-p finite element modelling framework. Stability regimes for a vehicle subjected to propulsive thrust are validated by comparing the results from published literature. Numerical simulations are carried out for a representative vehicle to determine the instability regimes with vehicle speed and propulsive thrust as the parameters. The phenomena of static instability (divergence) and dynamic instability (flutter) are observed. The results at low Mach number match closely with the results obtained from previous models published in the literature.

Brillouin Scattering by Phonons Renormalized by Two-Magnon Excitations in Ferromagnetic Dielectrics

Brillouin scattering by one-phonon-two-magnon interacting excitations in ferromagnetic dielectrics is discussed. The basic light scattering mechanism is taken to be the modulation of the density-dependent optical dielectric polarizability of the medium by the dynamic strain field generated by the longitudinal acoustic (LA) phonons. The renormalization effects arising from the scattering of phonons by the two-magnon creation-annihilation processes have, however, been taken into account. Via these interactions, the Brillouin components corresponding to the two-magnon excitations are reflected indirectly in the spectrum of the phonon scattered light as line-broadening of the otherwise relatively sharp Brillouin doublet. The present mechanism is shown to be dominant in a clean saturated ferromagnetic dielectric with large magneto-strictive coupling constant, and with the magnetic ions in the orbitally quenched states. Following the linear response theory, an expression has been derived for the spectral density of the scattered light as a function of temperature, scattering angle, and the strength of the externally applied magnetic field. Some estimates are given for the line-width and line-shift of the Brillouin components for certain typical choice of parameters involved. The results are discussed in relation to some available calculations on the ultrasonic attenuation in ferromagnetic insulators at low temperatures.

Study of third-order non-linear systems- Variation of parameters approach

In this paper, the transient response of a third-order non-linear system is obtained by first reducing the given third-order equation to three first-order equations by applying the method of variation of parameters. On the assumption that the variations of amplitude and phase are small, the functions are expanded in ultraspherical polynomials. The expansion is restricted to the constant term. The resulting equations are solved to obtain the response of the given third-order system. A numerical example is considered to illustrate the method. The results show that the agreement between the approximate and digital solution is good thus vindicating the approximation.

Automatic Path Planning and Control Design for Autonomous Landing of UAVs using Dynamic Inversion

In this paper a nonlinear control has been designed using the dynamic inversion approach for automatic landing of unmanned aerial vehicles (UAVs), along with associated path planning. This is a difficult problem because of light weight of UAVs and strong coupling between longitudinal and lateral modes. The landing maneuver of the UAV is divided into approach, glideslope and flare. In the approach UAV aligns with the centerline of the runway by heading angle correction. In glideslope and flare the UAV follows straight line and exponential curves respectively in the pitch plane with no lateral deviations. The glideslope and flare path are scheduled as a function of approach distance from runway. The trajectory parameters are calculated such that the sink rate at touchdown remains within specified bounds. It is also ensured that the transition from the glideslope to flare path is smooth by ensuring C-1 continuity at the transition. In the outer loop, the roll rate command is generated by assuring a coordinated turn in the alignment segment and by assuring zero bank angle in the glideslope and flare segments. The pitch rate command is generated from the error in altitude to control the deviations from the landing trajectory. The yaw rate command is generated from the required heading correction. In the inner loop, the aileron, elevator and rudder deflections are computed together to track the required body rate commands. Moreover, it is also ensured that the forward velocity of the UAV at the touch down remains close to a desired value by manipulating the thrust of the vehicle. A nonlinear six-DOF model, which has been developed from extensive wind-tunnel testing, is used both for control design as well as to validate it.