Spatially explicit structural equation modeling

Structural equation modeling (SEM) is a powerful statistical approach for the testing of networks of direct and indirect theoretical causal relationships in complex data sets with intercorrelated dependent and independent variables. SEM is commonly applied in ecology, but the spatial information commonly found in ecological data remains difficult to model in a SEM framework. Here we propose a simple method for spatially explicit SEM (SE-SEM) based on the analysis of variance/covariance matrices calculated across a range of lag distances. This method provides readily interpretable plots of the change in path coefficients across scale and can be implemented using any standard SEM software package. We demonstrate the application of this method using three studies examining the relationships between environmental factors, plant community structure, nitrogen fixation, and plant competition. By design, these data sets had a spatial component, but were previously analyzed using standard SEM models. Using these data sets, we demonstrate the application of SE-SEM to regularly spaced, irregularly spaced, and ad hoc spatial sampling designs and discuss the increased inferential capability of this approach compared with standard SEM. We provide an R package, sesem, to easily implement spatial structural equation modeling.

Modeling meat yield based on measurements of body traits in genetically improved giant freshwater prawn (GFP) Macrobrachium rosenbergii

A single-generation dataset consisting of 1,730 records from a selection program for high growth rate in giant freshwater prawn (GFP, Macrobrachium rosenbergii) was used to derive prediction equations for meat weight and meat yield. Models were based on body traits [body weight, total length and abdominal width (AW)] and carcass measurements (tail weight and exoskeleton-off weight). Lengths and width were adjusted for the systematic effects of selection line, male morphotypes and female reproductive status, and for the covariables of age at slaughter within sex and body weight. Body and meat weights adjusted for the same effects (except body weight) were used to calculate meat yield (expressed as percentage of tail weight/body weight and exoskeleton-off weight/body weight). The edible meat weight and yield in this GFP population ranged from 12 to 15 g and 37 to 45 %, respectively. The simple (Pearson) correlation coefficients between body traits (body weight, total length and AW) and meat weight were moderate to very high and positive (0.75–0.94), but the correlations between body traits and meat yield were negative (−0.47 to −0.74). There were strong linear positive relationships between measurements of body traits and meat weight, whereas relationships of body traits with meat yield were moderate and negative. Step-wise multiple regression analysis showed that the best model to predict meat weight included all body traits, with a coefficient of determination (R 2) of 0.99 and a correlation between observed and predicted values of meat weight of 0.99. The corresponding figures for meat yield were 0.91 and 0.95, respectively. Body weight or length was the best predictor of meat weight, explaining 91–94 % of observed variance when it was fitted alone in the model. By contrast, tail width explained a lower proportion (69–82 %) of total variance in the single trait models. It is concluded that in practical breeding programs, improvement of meat weight can be easily made through indirect selection for body trait combinations. The improvement of meat yield, albeit being more difficult, is possible by genetic means, with 91 % of the variation in the trait explained by the body and carcass traits examined in this study.

Development and validation of the Emergency Department Assessment of Chest pain Score and 2h accelerated diagnostic protocol

Objective Risk scores and accelerated diagnostic protocols can identify chest pain patients with low risk of major adverse cardiac event who could be discharged early from the ED, saving time and costs. We aimed to derive and validate a chest pain score and accelerated diagnostic protocol (ADP) that could safely increase the proportion of patients suitable for early discharge. Methods Logistic regression identified statistical predictors for major adverse cardiac events in a derivation cohort. Statistical coefficients were converted to whole numbers to create a score. Clinician feedback was used to improve the clinical plausibility and the usability of the final score (Emergency Department Assessment of Chest pain Score [EDACS]). EDACS was combined with electrocardiogram results and troponin results at 0 and 2 h to develop an ADP (EDACS-ADP). The score and EDACS-ADP were validated and tested for reproducibility in separate cohorts of patients. Results In the derivation (n = 1974) and validation (n = 608) cohorts, the EDACS-ADP classified 42.2% (sensitivity 99.0%, specificity 49.9%) and 51.3% (sensitivity 100.0%, specificity 59.0%) as low risk of major adverse cardiac events, respectively. The intra-class correlation coefficient for categorisation of patients as low risk was 0.87. Conclusion The EDACS-ADP identified approximately half of the patients presenting to the ED with possible cardiac chest pain as having low risk of short-term major adverse cardiac events, with high sensitivity. This is a significant improvement on similar, previously reported protocols. The EDACS-ADP is reproducible and has the potential to make considerable cost reductions to health systems.

Exploring the role of energy, trade and financial development in explaining economic growth in South Africa: A revisit

South Africa is an emerging and industrializing economy which is experiencing remarkable progress. We contend that amidst the developments in the economy, the role of energy, trade openness and financial development are critical. In this article, we revisit the pivotal role of these factors. We use the ARDL bounds [72], the Bayer and Hanck [11] cointegration techniques, and an extended Cobb–Douglas framework, to examine the long-run association with output per worker over the sample period 1971–2011. The results support long-run association between output per worker, capital per worker and the shift parameters. The short-run elasticity coefficients are as follows: energy (0.24), trade (0.07), financial development (−0.03). In the long-run, the elasticity coefficients are: trade openness (0.05), energy (0.29), and financial development (−0.04). In both the short-run and the long-run, we note the post-2000 period has a marginal positive effect on the economy. The Toda and Yamamoto [91] Granger causality results show that a unidirectional causality from capital stock and energy consumption to output; and from capital stock to trade openness; a bidirectional causality between trade openness and output; and absence (neutrality) of any causality between financial development and output thus indicating that these two variables evolve independent of each other.

Patient reported adherence to hypertension treatment: A revalidation study

Background Adherence to hypertension management in patients with hypertension is known to influence their blood pressure control. It is important to measure patients’ adherence behaviours to assist with designing appropriate interventions to improve blood pressure control. Aims The purposes of this study were to use confirmatory factor analysis to revalidate the Therapeutic Adherence Scale for Hypertensive Patients (TASHP), and to calculate the cut-off score for classifying adherence behaviours into two groups: satisfactory and low adherence behaviours. Methods Systematic random sampling was used to recruit patients with hypertension in China. Demographic characteristics, the TASHP and blood pressure were collected. The psychometric tests of the TASHP included: construct validity, criteria-related validity, internal reliability, and split-half reliability. The area under the receiver operating characteristics curve and Youden index were used to identify the cut-off score of the TASHP for blood pressure control. Results This study involved 366 patients. Confirmatory factor analysis supported the four-component structure of the TASHP proposed in the original scale development study. The TASHP has a satisfactory internal reliability (Cronbach’s α > 0.7) and a satisfactory split-half reliability (Spearman–Brown coefficients > 0.7). The patients with overall scores of the TASHP ⩾ 109 points were considered to have satisfactory adherence behaviours. Conclusion The TASHP is a validated and reliable instrument to measure the adherence to hypertension management in Chinese patients with hypertension. The cut-off score of 109 points can be considered as an effective measure to classify the level of adherence into satisfactory and low adherence behaviours.

Determination of imidacloprid in paddy water and soil by liquid chromatography electrospray ionization-tandem mass spectrometry

A method for the determination of imidacloprid in paddy water and soil was developed using liquid chromatography electrospray ionization-tandem mass spectrometry (LC/ESI-MS/MS). Separation of imidacloprid was carried out on a Shimadzu C18 column (150 mm × 4.6 mm, 4.6 μm) with an acetonitrile?water (50 : 50, v/v) mobile phase containing 0.1% of acetic acid. The flow rate was 0.3 mL/min in isocratic mode. The product ion at 209 m/z was selected for quantification in multiple-reaction monitoring scan mode. Imidacloprid residues in soil were extracted by a solid-liquid extraction method with acetonitrile. Water samples were filtered and directly injected for analysis without extraction. Detection limits of 0.5 μg/kg and 0.3 μg/L were achieved for soil and water samples, respectively. The method had recoveries of 90 ± 2% (n = 4) for soil samples and 100 ± 2% (n = 4) for water samples. A linear relationship was observed throughout the investigated range of concentrations (1-200 μg/L), with the correlation coefficients ranging from 0.999 to 1.000. © Pleiades Publishing, Ltd., 2010.

Bearing capacity of cold-formed unlipped channels with restrained flanges under EOF and IOF load cases

Bearing failure is a form of localized failure that occurs when thin-walled cold-formed steel sections are subjected to concentrated loads or support reactions. To determine the bearing capacity of cold-formed channel sections, a unified design equation with different bearing coefficients is given in the current North American specification AISI S100 and the Australian/New Zealand standard AS/NZS 4600. However, coefficients are not available for unlipped channel sections that are normally fastened to supports through their flanges. Eurocode 3 Part 1.3 includes bearing capacity equations for different load cases, but does not distinguish between fastened and unfastened support conditions. Therefore, an experimental study was conducted to determine the bearing capacities of these sections as used in floor systems. Twenty-eight web bearing tests on unlipped channel sections with restrained flanges were conducted under End One Flange (EOF) and Interior One Flange (IOF) load cases. Using the results from this study, a new equation was proposed within the AISI S100 and AS/NZS 4600 guidelines to determine the bearing capacities of cold-formed unlipped channels with flanges fastened to supports. A new design rule was also proposed based on the direct strength method.

Acoustic classification of Australian anurans using syllable features

Acoustic classification of anurans (frogs) has received increasing attention for its promising application in biological and environment studies. In this study, a novel feature extraction method for frog call classification is presented based on the analysis of spectrograms. The frog calls are first automatically segmented into syllables. Then, spectral peak tracks are extracted to separate desired signal (frog calls) from background noise. The spectral peak tracks are used to extract various syllable features, including: syllable duration, dominant frequency, oscillation rate, frequency modulation, and energy modulation. Finally, a k-nearest neighbor classifier is used for classifying frog calls based on the results of principal component analysis. The experiment results show that syllable features can achieve an average classification accuracy of 90.5% which outperforms Mel-frequency cepstral coefficients features (79.0%).

Acoustic feature extraction using perceptual wavelet packet decomposition for frog call classification

Frog protection has become increasingly essential due to the rapid decline of its biodiversity. Therefore, it is valuable to develop new methods for studying this biodiversity. In this paper, a novel feature extraction method is proposed based on perceptual wavelet packet decomposition for classifying frog calls in noisy environments. Pre-processing and syllable segmentation are first applied to the frog call. Then, a spectral peak track is extracted from each syllable if possible. Track duration, dominant frequency and oscillation rate are directly extracted from the track. With k-means clustering algorithm, the calculated dominant frequency of all frog species is clustered into k parts, which produce a frequency scale for wavelet packet decomposition. Based on the adaptive frequency scale, wavelet packet decomposition is applied to the frog calls. Using the wavelet packet decomposition coefficients, a new feature set named perceptual wavelet packet decomposition sub-band cepstral coefficients is extracted. Finally, a k-nearest neighbour (k-NN) classifier is used for the classification. The experiment results show that the proposed features can achieve an average classification accuracy of 97.45% which outperforms syllable features (86.87%) and Mel-frequency cepstral coefficients (MFCCs) feature (90.80%).

Application of image processing techniques for frog call classification

Frogs have received increasing attention due to their effectiveness for indicating the environment change. Therefore, it is important to monitor and assess frogs. With the development of sensor techniques, large volumes of audio data (including frog calls) have been collected and need to be analysed. After transforming the audio data into its spectrogram representation using short-time Fourier transform, the visual inspection of this representation motivates us to use image processing techniques for analysing audio data. Applying acoustic event detection (AED) method to spectrograms, acoustic events are firstly detected from which ridges are extracted. Three feature sets, Mel-frequency cepstral coefficients (MFCCs), AED feature set and ridge feature set, are then used for frog call classification with a support vector machine classifier. Fifteen frog species widely spread in Queensland, Australia, are selected to evaluate the proposed method. The experimental results show that ridge feature set can achieve an average classification accuracy of 74.73% which outperforms the MFCCs (38.99%) and AED feature set (67.78%).

The Post-Illumination Pupil Response (PIPR)

Purpose The post-illumination pupil response (PIPR) has been quantified using four metrics, but the spectral sensitivity of only one is known; here we determine the other three. To optimize the human PIPR measurement, we determine the protocol producing the largest PIPR, the duration of the PIPR, and the metric(s) with the lowest coefficient of variation. Methods The consensual pupil light reflex (PLR) was measured with a Maxwellian view pupillometer. - Experiment 1: Spectral sensitivity of four PIPR metrics [plateau, 6 s, area under curve (AUC) early and late recovery] was determined from a criterion PIPR to a 1s pulse and fitted with Vitamin A1 nomogram (λmax = 482nm). - Experiment 2: The PLR was measured as a function of three stimulus durations (1s, 10s, 30s), five irradiances spanning low to high melanopsin excitation levels (retinal irradiance: 9.8 to 14.8 log quanta.cm-2.s-1), and two wavelengths, one with high (465nm) and one with low (637nm) melanopsin excitation. Intra and inter-individual coefficients of variation (CV) were calculated. Results The melanopsin (opn4) photopigment nomogram adequately describes the spectral sensitivity of all four PIPR metrics. The PIPR amplitude was largest with 1s short wavelength pulses (≥ 12.8 log quanta.cm-2.s-1). The plateau and 6s PIPR showed the least intra and inter-individual CV (≤ 0.2). The maximum duration of the sustained PIPR was 83.0±48.0s (mean±SD) for 1s pulses and 180.1±106.2s for 30s pulses (465nm; 14.8 log quanta.cm-2.s-1). Conclusions All current PIPR metrics provide a direct measure of the intrinsic melanopsin photoresponse. To measure progressive changes in melanopsin function in disease, we recommend that the PIPR be measured using short duration pulses (e.g., ≤ 1s) with high melanopsin excitation and analyzed with plateau and/or 6s metrics. Our PIPR duration data provide a baseline for the selection of inter-stimulus intervals between consecutive pupil testing sequences.

The post-illumination pupil response of melanopsin expressing retinal ganglion cells

Purpose The post-illumination pupil response (PIPR) has been quantified in the literature by four metrics. The spectral sensitivity of only one metric is known and this study quantifies the other three. To optimize the measurement of the PIPR in humans, we also determine the stimulus protocol producing the largest PIPR, the duration of the PIPR, and the metric(s) with the lowest coefficient of variation. Methods The consensual pupil light reflex (PLR) was measured with a Maxwellian view pupillometer (35.6° diameter stimulus). - Experiment 1: Spectral sensitivity of four PIPR metrics [plateau, 6 s, area under curve (AUC) early and late recovery] was determined from a criterion PIPR (n = 2 participants) to a 1 s pulse at five wavelengths (409-592nm) and fitted with Vitamin A nomogram (ƛmax = 482 nm). - Experiment 2: The PLR was measured in five healthy participants [29 to 42 years (mean = 32.6 years)] as a function of three stimulus durations (1 s, 10 s, 30 s), five irradiances spanning low to high melanopsin excitation levels (retinal irradiance: 9.8 to 14.8 log quanta.cm-2.s-1), and two wavelengths, one with high (465 nm) and one with low (637 nm) melanopsin excitation. Intra and inter-individual coefficients of variation (CV) were calculated. Results The melanopsin (opn4) photopigment nomogram adequately described the spectral sensitivity derived from all four PIPR metrics. The largest PIPR amplitude was observed with 1 s short wavelength pulses (retinal irradiance ≥ 12.8 log quanta.cm-2.s-1). Of the 4 PIPR metrics, the plateau and 6 s PIPR showed the least intra and inter-individual CV (≤ 0.2). The maximum duration of the sustained PIPR was 83.4 ± 48.0 s (mean ± SD) for 1 s pulses and 180.1 ± 106.2 s for 30 s pulses (465 nm; 14.8 log quanta.cm-2.s-1). Conclusions All current PIPR metrics provide a direct measure of intrinsic melanopsin retinal ganglion cell function. To measure progressive changes in melanopsin function in disease, we recommend that the intrinsic melanopsin response should be measured using a 1 s pulse with high melanopsin excitation and the PIPR should be analyzed with the plateau and/or 6 s metrics. That the PIPR can have a sustained constriction for as long as 3 minutes, our PIPR duration data provide a baseline for the selection of inter-stimulus intervals between consecutive pupil testing sequences.

Pulse-echo ultrasound transit time spectroscopy: A comparison of experimental measurements and simulation prediction

Considering ultrasound propagation through complex composite media as an array of parallel sonic rays, a comparison of computer simulated prediction with experimental data has previously been reported for transmission mode (where one transducer serves as transmitter, the other as receiver) in a series of ten acrylic step-wedge samples, immersed in water, exhibiting varying degrees of transit time inhomogeneity. In this study, the same samples were used but in pulse-echo mode, where the same ultrasound transducer served as both transmitter and receiver, detecting both ‘primary’ (internal sample interface) and ‘secondary’ (external sample interface) echoes. A transit time spectrum (TTS) was derived, describing the proportion of sonic rays with a particular transit time. A computer simulation was performed to predict the transit time and amplitude of various echoes created, and compared with experimental data. Applying an amplitude-tolerance analysis, 91.7±3.7% of the simulated data was within ±1 standard deviation (STD) of the experimentally measured amplitude-time data. Correlation of predicted and experimental transit time spectra provided coefficients of determination (R2) ranging from 100.0% to 96.8% for the various samples tested. The results acquired from this study provide good evidence for the concept of parallel sonic rays. Further, deconvolution of experimental input and output signals has been shown to provide an effective method to identify echoes otherwise lost due to phase cancellation. Potential applications of pulse-echo ultrasound transit time spectroscopy (PE-UTTS) include improvement of ultrasound image fidelity by improving spatial resolution and reducing phase interference artefacts.

Wavelet-based multiresolution analysis of Wivenhoe Dam water temperatures

Water temperature measurements from Wivenhoe Dam offer a unique opportunity for studying fluctuations of temperatures in a subtropical dam as a function of time and depth. Cursory examination of the data indicate a complicated structure across both time and depth. We propose simplifying the task of describing these data by breaking the time series at each depth into physically meaningful components that individually capture daily, subannual, and annual (DSA) variations. Precise definitions for each component are formulated in terms of a wavelet-based multiresolution analysis. The DSA components are approximately pairwise uncorrelated within a given depth and between different depths. They also satisfy an additive property in that their sum is exactly equal to the original time series. Each component is based upon a set of coefficients that decomposes the sample variance of each time series exactly across time and that can be used to study both time-varying variances of water temperature at each depth and time-varying correlations between temperatures at different depths. Each DSA component is amenable for studying a certain aspect of the relationship between the series at different depths. The daily component in general is weakly correlated between depths, including those that are adjacent to one another. The subannual component quantifies seasonal effects and in particular isolates phenomena associated with the thermocline, thus simplifying its study across time. The annual component can be used for a trend analysis. The descriptive analysis provided by the DSA decomposition is a useful precursor to a more formal statistical analysis.

Population structure, mortality and growth of Pinna nobilis Linnaeus, 1758 (Mollusca, Bivalvia) at different depths in Moraira bay (Alicante, Western Mediterranean)

An investigation to characterize the causes of Pinna nobilis population structure in Moraira bay (Western Mediterranean) was developed. Individuals of two areas of the same Posidonia meadow, located at different depths (A1, -13 and A2, -6 m), were inventoried, tagged, their positions accurately recorded and monitored from July 1997 to July 2002. On each area, different aspects of population demography were studied (i.e. spatial distribution, size structure, displacement evidences, mortality, growth and shell orientation). A comparison between both groups of individuals was carried out, finding important differences between them. In A1, the individuals were more aggregated and mean and maximum size were higher (A1, 10.3 and A2, 6 individuals/100 m(2); A1, x = 47.2 +/- 9.9; A2, x = 29.8 +/- 7.4 cm, P < 0.001, respectively). In A2, growth rate and mortality were higher, the latter concentrated on the largest individuals, in contrast to A1, where the smallest individuals had the higher mortality rate [A1, L = 56.03(1 - e(-0.17t)); A2, L = 37.59(1 - e(-0.40t)), P < 0.001; mean annual mortality A1: 32 dead individuals out of 135, 23.7% and A2: 16 dead individuals out of 36, 44.4%, and total mortality coefficients (z), z(A1(-30)) = 0.28, z(A1(31-45)) = 0.05, z(A1(46-)) = 0.08; z(A2(-30)) = 0.15, z(A2(31-45)) = 0.25]. A common shell orientation N-S, coincident with the maximum shore exposure, was observed in A2. Spatial distribution in both areas showed not enough evidence to discard a random distribution of the individuals, despite the greater aggregation on the deeper area (A1) (A1, chi(2) = 0.41, df = 3, P > 0.5, A2, chi(2)= 0.98, df = 2 and 0.3 < P < 0.5). The obtained results have demonstrated that the depth-related size segregation usually shown by P. nobilis is mainly caused by differences in mortality and growth among individuals located at different depths, rather than by the active displacement of individuals previously reported in the literature. Furthermore, dwarf individuals are observed in shallower levels and as a consequence, the relationship between size and age are not comparable even among groups of individuals inhabiting the same meadow at different depths. The final causes of the differences on mortality and growth are also discussed.