Comparing sampling needs for variograms of soil properties computed by the method of moments and residual maximum likelihood

It has been generally accepted that the method of moments (MoM) variogram, which has been widely applied in soil science, requires about 100 sites at an appropriate interval apart to describe the variation adequately. This sample size is often larger than can be afforded for soil surveys of agricultural fields or contaminated sites. Furthermore, it might be a much larger sample size than is needed where the scale of variation is large. A possible alternative in such situations is the residual maximum likelihood (REML) variogram because fewer data appear to be required. The REML method is parametric and is considered reliable where there is trend in the data because it is based on generalized increments that filter trend out and only the covariance parameters are estimated. Previous research has suggested that fewer data are needed to compute a reliable variogram using a maximum likelihood approach such as REML, however, the results can vary according to the nature of the spatial variation. There remain issues to examine: how many fewer data can be used, how should the sampling sites be distributed over the site of interest, and how do different degrees of spatial variation affect the data requirements? The soil of four field sites of different size, physiography, parent material and soil type was sampled intensively, and MoM and REML variograms were calculated for clay content. The data were then sub-sampled to give different sample sizes and distributions of sites and the variograms were computed again. The model parameters for the sets of variograms for each site were used for cross-validation. Predictions based on REML variograms were generally more accurate than those from MoM variograms with fewer than 100 sampling sites. A sample size of around 50 sites at an appropriate distance apart, possibly determined from variograms of ancillary data, appears adequate to compute REML variograms for kriging soil properties for precision agriculture and contaminated sites. (C) 2007 Elsevier B.V. All rights reserved.

Quantifying the residual creep life of polyester mooring ropes

This paper assesses the ARELIS (Assured Residual Life Span) method for estimating residual creep life of polyester rope used in deepwater mooring lines. A statistical model has been developed to quantify the uncertainties in the method, such as the scatter in creep rupture test data and load sharing between sub-ropes. This model can be used to determine the required test load, duration and number of ARELIS tests, in order to guarantee a minimum creep life for a mooring line at its service load. Creep rupture tests have been performed to provide input for the statistical model.

Efficiency and purity control in the preparation of pure and/or aluminum-doped barium ferrites by hydrothermal methods using ferrous ions as reactants

The synthesis of hexagonal barium ferrite (BaFe12O19) was studied under hydrothermal conditions by a method in which a significant amount of ferrous chloride was introduced along side ferric chloride among the starting materials. Though all of the Fe2+ ions in the starting material were converted to Fe3+ ions in the final product, Fe2+ was confirmed to participate differently from the Fe3+ used in the conventional method in the mechanism of forming barium ferrite. Indeed the efficiency of the synthesis and the quality of the product and the lack of impurities such as Fe2O3 and BaFe2O4 were improved when Fe2+ was included. However, the amount of ferrous ions that could be included to obtain the desired product was limited with an optimum ratio of 2:8 for FeCl2/FeCl3 when only 2h of reaction time were needed. It was also found that the role of trivalent Fe3+ could be successfully replaced by Al3+. Up to 50% of their on could be replaced by Al3+ in the reactants to produce Al- doped products. It was also found that the ratio of Fe2+/M3+ could be increased in the presence of Al3+ to produce high quality barium ferrite.

Harringtonian virtue: Harrington, Machiavelli, and the method of the 'Moment'

This article presents a reinterpretation of James Harrington's writings. It takes issue with J. G. A. Pocock's reading, which treats him as importing into England a Machiavellian ‘language of political thought’. This reading is the basis of Pocock's stress on the republicanism of eighteenth-century opposition values. Harrington's writings were in fact a most implausible channel for such ideas. His outlook owed much to Stoicism. Unlike the Florentine, he admired the contemplative life; was sympathetic to commerce; and was relaxed about the threat of ‘corruption’ (a concept that he did not understand). These views can be associated with his apparent aims: the preservation of a national church with a salaried but politically impotent clergy; and the restoration of the royalist gentry to a leading role in English politics. Pocock's hypothesis is shown to be conditioned by his method; its weaknesses reflect some difficulties inherent in the notion of ‘languages of thought’.

Wavelets and the generalization of the variogram

The experimental variogram computed in the usual way by the method of moments and the Haar wavelet transform are similar in that they filter data and yield informative summaries that may be interpreted. The variogram filters out constant values; wavelets can filter variation at several spatial scales and thereby provide a richer repertoire for analysis and demand no assumptions other than that of finite variance. This paper compares the two functions, identifying that part of the Haar wavelet transform that gives it its advantages. It goes on to show that the generalized variogram of order k=1, 2, and 3 filters linear, quadratic, and cubic polynomials from the data, respectively, which correspond with more complex wavelets in Daubechies's family. The additional filter coefficients of the latter can reveal features of the data that are not evident in its usual form. Three examples in which data recorded at regular intervals on transects are analyzed illustrate the extended form of the variogram. The apparent periodicity of gilgais in Australia seems to be accentuated as filter coefficients are added, but otherwise the analysis provides no new insight. Analysis of hyerpsectral data with a strong linear trend showed that the wavelet-based variograms filtered it out. Adding filter coefficients in the analysis of the topsoil across the Jurassic scarplands of England changed the upper bound of the variogram; it then resembled the within-class variogram computed by the method of moments. To elucidate these results, we simulated several series of data to represent a random process with values fluctuating about a mean, data with long-range linear trend, data with local trend, and data with stepped transitions. The results suggest that the wavelet variogram can filter out the effects of long-range trend, but not local trend, and of transitions from one class to another, as across boundaries.

Wiggle-match dating of tree-ring sequences

Given the non-monotonic form of the radiocarbon calibration curve, the precision of single C-14 dates on the calendar timescale will always be limited. One way around this limitation is through comparison of time-series, which should exhibit the same irregular patterning as the calibration curve. This approach can be employed most directly in the case of wood samples with many years growth present (but not able to be dated by dendrochronology), where the tree-ring series of unknown date can be compared against the similarly constructed C-14 calibration curve built from known-age wood. This process of curve-fitting has come to be called "wiggle-matching." In this paper, we look at the requirements for getting good precision by this method: sequence length, sampling frequency, and measurement precision. We also look at 3 case studies: one a piece of wood which has been independently dendrochronologically dated, and two others of unknown age relating to archaeological activity at Silchester, UK (Roman) and Miletos, Anatolia (relating to the volcanic eruption at Thera).

Vibrational intensities VI: ethylene and its deutero-isotopes

Absolute intensity measurements have been made on the fundamental vibrations of ethylene and four of its deuteroisotopes. The bands were pressure broadened with nitrogen at 50 atmos, and the intensities were determined by the method of Wilson and Wells except that the observed optical density was integrated against logv rather than v. Normal coordinates have been calculated, and the intensities have been interpreted in terms of quantities (∂p/∂Si) giving the change in dipole moment with respect to each internal symmetry coordinate. Data from the different isotopic species have been used to eliminate ambiguities in the interpretation. Effective bond moments are calculated for each symmetry coordinate.

Extending Zelterman’s approach for robust estimation of population size to zero-truncated clustered data

Estimation of population size with missing zero-class is an important problem that is encountered in epidemiological assessment studies. Fitting a Poisson model to the observed data by the method of maximum likelihood and estimation of the population size based on this fit is an approach that has been widely used for this purpose. In practice, however, the Poisson assumption is seldom satisfied. Zelterman (1988) has proposed a robust estimator for unclustered data that works well in a wide class of distributions applicable for count data. In the work presented here, we extend this estimator to clustered data. The estimator requires fitting a zero-truncated homogeneous Poisson model by maximum likelihood and thereby using a Horvitz-Thompson estimator of population size. This was found to work well, when the data follow the hypothesized homogeneous Poisson model. However, when the true distribution deviates from the hypothesized model, the population size was found to be underestimated. In the search of a more robust estimator, we focused on three models that use all clusters with exactly one case, those clusters with exactly two cases and those with exactly three cases to estimate the probability of the zero-class and thereby use data collected on all the clusters in the Horvitz-Thompson estimator of population size. Loss in efficiency associated with gain in robustness was examined based on a simulation study. As a trade-off between gain in robustness and loss in efficiency, the model that uses data collected on clusters with at most three cases to estimate the probability of the zero-class was found to be preferred in general. In applications, we recommend obtaining estimates from all three models and making a choice considering the estimates from the three models, robustness and the loss in efficiency. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)

A system for studying epithelial-stromal interactions reveals distinct inductive abilities of stromal cells from benign prostatic hyperplasia and prostate cancer

The development of normal and abnormal glandular structures in the prostate is controlled at the endocrine and paracrine levels by reciprocal interactions between epithelium and stroma. To study these processes it is useful to have an efficient method of tissue acquisition for reproducible isolation of cells from defined histologies. Here we assessed the utility of a standardized system for acquisition and growth of prostatic cells from different regions of the prostate with different pathologies, and we compared the abilities of stromal cells from normal peripheral zone (PZ-S), benign prostatic hyperplasia (BPH-S), and cancer (CA-S) to induce the growth of a human prostatic epithelial cell line (BPH-1) in vivo. Using the tissue recombination method, we showed that grafting stromal cells (from any histology) alone, or BPH-1 epithelial cells alone produced no visible grafts. Recombining PZ-S with BPH-1 cells also produced no visible grafts (n = 15). Recombining BPH-S with BPH-1 cells generated small, well-organized and sharply demarcated grafts approximately 3-4 mm in diameter (n = 9), demonstrating a moderate inductive ability of BPH-S. Recombining CA-S with BPH-1 cells generated highly disorganized grafts that completely surrounded the host kidney and invaded into adjacent renal tissue, demonstrating induction of an aggressive phenotype. We conclude that acquisition of tissue from toluidine blue dye stained specimens is an efficient method to generate high quality epithelial and/or stromal cultures. Stromal cells derived by this method from areas of BPH and cancer induce epithelial cell growth in vivo which mimics the natural history of these diseases.

Comparison of new microemulsion prepared "Pt-in-Ceria" catalyst with conventional "Pt-on-Ceria" catalyst for water-gas shift reaction

New "Pt-in-CeO2" catalyst prepared by microemulsion method is shown to give higher activity for a water-gas shift reaction but with no formation of CH4, the side product from hydrogenation of carbon oxides using a hydrogen-rich reformate as compared to conventional "Pt-on-CeO2" catalysts. Detailed characterization by DRIFT analysis and temperature programmed reduction presented in this work clearly suggest the ceria coverage on Pt inhibits the metal from forming a strong CO adsorption.

Fecal water as a non-invasive biomarker in nutritional intervention: comparison of preparation methods and refinement of different endpoints

The assessment of cellular effects by the aqueous phase of human feces (fecal water, FW) is a useful biomarker approach to study cancer risks and protective activities of food. In order to refine and develop the biomarker, different protocols of preparing FW were compared. Fecal waters were prepared by 3 methods: (A) direct centrifugation; (B) extraction of feces in PBS before centrifugation; and (C) centrifugation of lyophilized and reconstituted feces. Genotoxicity was determined in colon cells using the Comet assay. Selected samples were investigated for additional parameters related to carcinogenesis. Two of 7 FWs obtained by methods A and B were similarly genotoxic. Method B, however, yielded higher volumes of FW, allowing sterile filtration for long-term culture experiments. Four of 7 samples were non-genotoxic when prepared according to all 3 methods. FW from lyophilized feces and from fresh samples were equally genotoxic. FWs modulated cytotoxicity, paracellular permeability, and invasion, independent of their genotoxicity. All 3 methods of FW preparation can be used to assess genotoxicity. The higher volumes of FWobtained by preparation method B greatly enhance the perspectives of measuring different types of biological parameters and using these to disclose activities related to cancer development.

Structures of Pd(CN)2and Pt(CN)2: intrinsically nanocrystalline materials?

Analysis and modeling of X-ray and neutron Bragg and total diffraction data show that the compounds referred to in the literature as “Pd(CN)2”and“Pt(CN)2” are nanocrystalline materials containing of small sheets of vertex-sharing square-planar M(CN)4 units, layered in a disordered manner with an intersheet separation of 3.44 A at 300 K. The small size of the crystallites means that the sheets’ edges form a significant fraction of each material. The Pd(CN)2 nanocrystallites studied using total neutron diffraction are terminated by water and the Pt(CN)2 nanocrystallites by ammonia, in place of half of the terminal cyanide groups, thus maintaining charge neutrality. The neutron samples contain sheets of approximate dimensions 30 A x 30 A. For sheets of the size we describe, our structural models predict compositions of Pd(CN)2-xH2O and Pt(CN)2-yNH3 (x = y = 0.29). These values are in good agreement with those obtained from total neutron diffraction and thermal analysis, and are also supported by infrared and Raman spectroscopy measurements. It is also possible to prepare related compounds Pd(CN)2-pNH3 and Pt(CN)2-qH2O, in which the terminating groups are exchanged. Additional samples showing sheet sizes in the range 10 A x 10 A (y = 0.67) to 80 A x 80 A (p = q = 0.12), as determined by X-ray diffraction, have been prepared. The related mixed-metal phase, Pd1/2Pt1/2(CN)2-qH2O(q = 0.50), is also nanocrystalline (sheet size 15 A x 15 A). In all cases, the interiors of the sheets are isostructural with those found in Ni(CN)2. Removal of the final traces of water or ammonia by heating results in decomposition of the compounds to Pd and Pt metal, or in the case of the mixed-metal cyanide, the alloy, Pd1/2Pt1/2, making it impossible to prepare the simple cyanides, Pd(CN)2, Pt(CN)2 or Pd1/2Pt1/2(CN)2, by this method.

Synthesis of broadband anti-reflection coatings for spaceflight infrared optics

A synthesis method is outlined for the design of broadband anti-reflection coatings for use in spaceborne infrared optics. The Golden Section optimisation routine is used to make a search, using designated non-absorptive dielectric thin film combinations, for the coating design which fulfils the required spectral requirements using the least number of layers and different materials. Three examples are given of coatings designed by this method : (I) 1µm to 12µm anti-reflection coating on Zinc Sulphide using Zinc Sulphide and Yttrium Fluoride thin film materials. (ii) 2µm to 14µm anti-reflection coating on Germanium using Germanium and Ytterbium Fluoride thin film materials. (iii) 6µm to 17µm anti-reflection coating on Germanium using Lead Telluride, Zinc Selenide and Barium Fluoride. The measured spectral performance of the manufactured 6µm to 17µm coating on Germanium is given. This is the anti-reflection coating for the germanium optics in the NASA Cassini Orbiter CIRS instrument.

Calibration of closed loop controllers for setting impedances in force-reflecting systems

This paper discusses a new method of impedance control that has been successfully implemented on the master robot of a teleoperation system. The method involves calibrating the robot to quantify the effect of adjustable controller parameters on the impedances along its different axes. The empirical equations relating end-effector impedance to the controller's feedback gains are obtained by performing system identification tests along individual axes of the robot. With these equations, online control of end-effector stiffness and damping is possible without having to monitor joint torques or solving complex algorithms. Hard contact conditions and compliant interfaces have been effectively demonstrated on a telemanipulation test-bed using appropriate combinations of stiffness and damping settings obtained by this method.