Likelihood-based estimation of the regression model with scrambled responses

A significant problem in the collection of responses to potentially sensitive questions, such as relating to illegal, immoral or embarrassing activities, is non-sampling error due to refusal to respond or false responses. Eichhorn & Hayre (1983) suggested the use of scrambled responses to reduce this form of bias. This paper considers a linear regression model in which the dependent variable is unobserved but for which the sum or product with a scrambling random variable of known distribution, is known. The performance of two likelihood-based estimators is investigated, namely of a Bayesian estimator achieved through a Markov chain Monte Carlo (MCMC) sampling scheme, and a classical maximum-likelihood estimator. These two estimators and an estimator suggested by Singh, Joarder & King (1996) are compared. Monte Carlo results show that the Bayesian estimator outperforms the classical estimators in almost all cases, and the relative performance of the Bayesian estimator improves as the responses become more scrambled.

Mechanics of raceway hysteresis in a packed bed

When a gas is introduced at high velocity through a nozzle into a packed bed, it creates a raceway in the packed bed. It has been found that the raceway size is larger when it is formed by decreasing the gas velocity from its highest value than when it is formed by increasing the gas velocity. This phenomenon is known as raceway hysteresis. A hypothesis has been oroposed to explain the hysteresis phenomenon based on a force-balance approach which includes frictional, bed-weight, and pressure forces. According to this hypothesis, the frictional force acts in different directions when the raceway is expanding and contracting. In this article, the entire packed bed has been divided into radial and Cartesian co-ordinate systems, and the forces acting on the raceway have been solved analytically for a simplified one-dimensional case. Based on the force-balance approach, a general equation has been obtained to predict the diameter of the raceway for increasing And decreasing velocities. A reasonable agreement has been found between the theoretical predictions and experimental observations. The model has also been compared with published experimental and plant data. The hysteresis mechanism in the packed beds can be described reasonably by taking into consideration the direction of frictional forces for the increasing- and decreasin-velocity cases. The effects of the particleshape factor and void fraction on the raceway hysteresis are examined.

Improved slope adjustment functions for soil erosion prediction

Numerous studies in the last 60 years have investigated the relationship between land slope and soil erosion rates. However, relatively few of these have investigated slope gradient responses: ( a) for steep slopes, (b) for specific erosion processes, and ( c) as a function of soil properties. Simulated rainfall was applied in the laboratory on 16 soils and 16 overburdens at 100 mm/h to 3 replicates of unconsolidated flume plots 3 m long by 0.8 m wide and 0.15 m deep at slopes of 20, 5, 10, 15, and 30% slope in that order. Sediment delivery at each slope was measured to determine the relationship between slope steepness and erosion rate. Data from this study were evaluated alongside data and existing slope adjustment functions from more than 55 other studies from the literature. Data and the literature strongly support a logistic slope adjustment function of the form S = A + B/[1 + exp (C - D sin theta)] where S is the slope adjustment factor and A, B, C, and D are coefficients that depend on the dominant detachment and transport processes. Average coefficient values when interill-only processes are active are A - 1.50, B 6.51, C 0.94, and D 5.30 (r(2) = 0.99). When rill erosion is also potentially active, the average slope response is greater and coefficient values are A - 1.12, B 16.05, C 2.61, and D 8.32 (r(2) = 0.93). The interill-only function predicts increases in sediment delivery rates from 5 to 30% slope that are approximately double the predictions based on existing published interill functions. The rill + interill function is similar to a previously reported value. The above relationships represent a mean slope response for all soils, yet the response of individual soils varied substantially from a 2.5-fold to a 50-fold increase over the range of slopes studied. The magnitude of the slope response was found to be inversely related ( log - log linear) to the dispersed silt and clay content of the soil, and 3 slope adjustment equations are proposed that provide a better estimate of slope response when this soil property is known. Evaluation of the slope adjustment equations proposed in this paper using independent datasets showed that the new equations can improve soil erosion predictions.

Fat as an endocrine organ: Relationship to the metabolic syndrome

Obesity and the metabolic syndrome have both reached pandemic proportions. Together they have the potential to impact on the incidence and severity of cardiovascular pathologies, with grave implications for worldwide health care systems. The metabolic syndrome is characterized by visceral obesity, insulin resistance, hypertension, chronic inflammation, and thrombotic disorders contributing to endothelial dysfunction and, subsequently, to accelerated atherosclerosis. Obesity is a key component in development of the metabolic syndrome and it is becoming increasingly clear that a central factor in this is the production by adipose cells of bioactive substances that directly influence insulin sensitivity and vascular injury. In this paper, we review advances in the understanding of biologically active molecules collectively referred to as adipokines and how dysregulated production of these factors in obese states mediates the pathogenesis of obesity associated metabolic syndrome.

Fitting genetic models to twin data with binary and ordered categorical responses: A comparison of structural equation modelling and Bayesian hierarchical models

We compare Bayesian methodology utilizing free-ware BUGS (Bayesian Inference Using Gibbs Sampling) with the traditional structural equation modelling approach based on another free-ware package, Mx. Dichotomous and ordinal (three category) twin data were simulated according to different additive genetic and common environment models for phenotypic variation. Practical issues are discussed in using Gibbs sampling as implemented by BUGS to fit subject-specific Bayesian generalized linear models, where the components of variation may be estimated directly. The simulation study (based on 2000 twin pairs) indicated that there is a consistent advantage in using the Bayesian method to detect a correct model under certain specifications of additive genetics and common environmental effects. For binary data, both methods had difficulty in detecting the correct model when the additive genetic effect was low (between 10 and 20%) or of moderate range (between 20 and 40%). Furthermore, neither method could adequately detect a correct model that included a modest common environmental effect (20%) even when the additive genetic effect was large (50%). Power was significantly improved with ordinal data for most scenarios, except for the case of low heritability under a true ACE model. We illustrate and compare both methods using data from 1239 twin pairs over the age of 50 years, who were registered with the Australian National Health and Medical Research Council Twin Registry (ATR) and presented symptoms associated with osteoarthritis occurring in joints of the hand.

Effect of adsorption deformation on thermodynamic characteristics of a fluid in slit pores at sub-critical conditions

Solvation. pressure due to adsorption of fluids in porous materials is the cause of elastic deformation of an adsorbent, which is accessible to direct experimental measurements. Such a deformation contributes to the Helmholtz free energy of the whole adsorbent-adsorbate system due to accumulation of compression or tension energy by the solid. It means that in the general case the solid has to be considered as not solely a source of the external potential field for the fluid confined in the pore volume, but also as thermodynamically nonmert component of the solid-fluid system. We present analysis of nitrogen adsorption isotherms and heat of adsorption in slit graphitic pores accounting for the adsorption deformation by means of nonlocal density functional theory. (c) 2006 Elsevier Ltd. All rights reserved.

What determines the finance-growth nexus? Empirical evidence for threshold models

This paper elaborates the notion of balanced'' financial development that is contingent on a country's general level of development. We develop an empirical framework to address this point, referring to threshold regressions and a bootstrap test for structural shift in a growth equation. We find that countries gain less from financial activity, if the latter fails to keep up with or exceeds what would follow from a balanced expansion path. These analyses contribute to the finance and growth literature in providing empirical support for the balanced'' financial development hypothesis.

Popper's Explications of Ad hocness: Circularity, Empirical Content, and Scientific Practice

Popper's explications of 'ad hoc' in relation to hypotheses and explanations turn out to be either trivial, confused or mistaken. One such explication I discuss at length is circularity; another is reduction in empirical content. I argue that non-circularity is preferable to non-ad hocness for an acceptable explanation or explanans, and I isolate some persistent errors in his analysis. Second, Popper is barking up the wrong tree in proscribing reductions in empirical content in novel hypotheses. Such reductions may constitute scientific progress. He fails to show that ad hoc hypothesis are the threat to science he claims.

A finite integral transform technique for solving the diffusion-reaction equation with Michaelis-Menten kinetics

An approximate analytical technique employing a finite integral transform is developed to solve the reaction diffusion problem with Michaelis-Menten kinetics in a solid of general shape. A simple infinite series solution for the substrate concentration is obtained as a function of the Thiele modulus, modified Sherwood number, and Michaelis constant. An iteration scheme is developed to bring the approximate solution closer to the exact solution. Comparison with the known exact solutions for slab geometry (quadrature) and numerically exact solutions for spherical geometry (orthogonal collocation) shows excellent agreement for all values of the Thiele modulus and Michaelis constant.

Refining the risk-benefit equation for thrombolysis: How to identity the low risk patient before administering thrombolytic therapy

In view of the relative risk of intracranial haemorrhage and major bleeding with thrombolytic therapy, it is important ro identify as early as possible the low risk patient who may not have a net clinical benefit from thrombolysis in the setting of acute myocardial infarction. An analysis of 5434 hospital-treated patients with myocardial infarction in the Perth MONICA study showed that age below 60 and absence of previous infarction or diabetes, shock, pulmonary oedema, cardiac arrest and Q-wave or left bundle branch block on the initial ECG identified a large group of patients with a 28 day mortality of only 1%, and one year mortality of only 2%. Identification of baseline risk in this way helps refine the risk-benefit equation for thrombolytic therapy, and may help avoid unnecessary use of thrombolysis in those unlikely to benefit.

Truncation errors in finite difference models for solute transport equation with first-order reaction

The truncation errors associated with finite difference solutions of the advection-dispersion equation with first-order reaction are formulated from a Taylor analysis. The error expressions are based on a general form of the corresponding difference equation and a temporally and spatially weighted parametric approach is used for differentiating among the various finite difference schemes. The numerical truncation errors are defined using Peclet and Courant numbers and a new Sink/Source dimensionless number. It is shown that all of the finite difference schemes suffer from truncation errors. Tn particular it is shown that the Crank-Nicolson approximation scheme does not have second order accuracy for this case. The effects of these truncation errors on the solution of an advection-dispersion equation with a first order reaction term are demonstrated by comparison with an analytical solution. The results show that these errors are not negligible and that correcting the finite difference scheme for them results in a more accurate solution. (C) 1999 Elsevier Science B.V. All rights reserved.

Empirical catchment-wide rainfall erosivity models for two rivers in the humid tropics of Australia

The concept of rainfall erosivity is extended to the estimation of catchment sediment yield and its variation over time. Five different formulations of rainfall erosivity indices, using annual, monthly and daily rainfall data, are proposed and tested on two catchments in the humid tropics of Australia. Rainfall erosivity indices, using simple power functions of annual and daily rainfall amounts, were found to be adequate in describing the interannual and seasonal variation of catchment sediment yield. The parameter values of these rainfall erosivity indices for catchment sediment yield are broadly similar to those for rainfall erosivity models in relation to the R-factor in the Universal Soil Loss Equation.

Muscle fibre types and size distribution in sub-antarctic notothenioid fishes

The presumptive tonic muscles fibres of Cottoperca gobio, Champsocephalus esox, Harpagifer bispinis, Eleginops maclovinus, Patagonothen tessellata, P. cornucola and Paranotothenia magellanica stained weakly or were unstained for glycogen, lipid, succinic dehydrogenase (SDHase) and myosin ATPase (mATPase) activity. Slow, intermediate and fast twitch muscle fibres, distinguished on the basis of the pH stability of their mATPases, showed intense, moderate and low staining activity for SDHase, respectively. Slow fibres were the major component of the pectoral fin adductor profundis muscle. The proportion of different muscle fibre types varied from the proximal to distal end of the muscle, but showed relatively little variation between species. The myotomes contained a lateral superficial strip of red muscle composed of presumptive tonic, slow twitch and intermediate fibres, thickening to a major wedge at the horizontal septum. All species also had characteristic secondary dorsal and ventral wedges of red muscle. The relative abundance and localization of muscle fibre types in the red muscle varied between species and with body size in the protandric hermaphrodite E. maclovinus. The frequency distribution of diameters for fast twitch muscle fibres, the major component of deep white muscle, was determined in fish of a range of body sizes. The absence of fibres <20 mu m diameter was used as a criterion for the cessation of muscle fibre recruitment. Fibre recruitment had stopped in P, tessellata of 13.8 cm L-T and E, maclovinus of 32.8 cm L-T, equivalent to 49 and 36.5% of their recorded maximum sizes respectively. As a result in 20-cm P. tessellata, the maximum fibre diameter was 300 mu m and 36% of fibres were in excess of 200 mu m The unusually large maximum fibre diameter, the general arrangement of the red muscle layer and the extreme pH lability of the mATPase of fast twitch fibres are all common characters of the sub-Antarctic and Antarctic Notothenioids, including Cottoperca gobio, the suggested sister group to the Notothenidae. (C) 2000 The Fisheries Society of the British Isles.

Comments on "Exact solution for horizontal redistribution by general similarity"

No abstract

General properties of cosmological models with an isotropic singularity

Much of the published work regarding the Isotropic Singularity is performed under the assumption that the matter source for the cosmological model is a barotropic perfect fluid, or even a perfect fluid with a gamma-law equation of state. There are, however, some general properties of cosmological models which admit an Isotropic Singularity, irrespective of the matter source. In particular, we show that the Isotropic Singularity is a point-like singularity and that vacuum space-times cannot admit an Isotropic Singularity. The relationships between the Isotropic Singularity, and the energy conditions, and the Hubble parameter is explored. A review of work by the authors, regarding the Isotropic Singularity, is presented.