The impact of changes to heroin supply on blood-borne virus notifications and injecting related harms in New South Wales, Australia

Background: In early 2001 Australia experienced a sudden and unexpected disruption to heroin availability, know as the 'heroin shortage'. This 'shortage has been linked to a decrease in needle and syringe output and therefore possibly a reduction in injecting drug use. We aimed to examine changes, if any, in blood-borne viral infections and presentations for injecting related problems related to injecting drug use following the reduction heroin availability in Australia, in the context of widespread harm reduction measures. Methods: Time series analysis of State level databases on HIV, hepatitis B, hepatitis C notifications and hospital and emergency department data. Examination of changes in HIV, hepatitis B, hepatitis C notifications and hospital and emergency department admissions for injection-related problems following the onset of the heroin shortage; non-parametric curve-fitting of number of hepatitis C notifications among those aged 15 - 19 years. Results: There were no changes observed in hospital visits for injection-related problems. There was no change related to the onset heroin shortage in the number of hepatitis C notifications among persons aged 15 - 19 years, but HCV notifications have subsequently decreased in this group. No change occurred in HIV and hepatitis B notifications. Conclusion: A marked reduction in heroin supply resulted in no increase in injection-related harm at the community level. However, a delayed decrease in HCV notifications among young people may be related. These changes occurred in a setting with widespread, publicly funded harm reduction initiatives.

NMR imaging of the diffusion of water at 310 K into poly[(2-hydroxyethyl methacrylate)-co-(tetrahydrofurfuryl methacrylate)] containing vitamin B12 or aspirin

The ingress of water into copolymers of 2-hydroxyethyl methacrylate (HEMA) and tetrahydrofurfuryl methacrylate (THFMA) loaded with either one of two model drugs, ie vitamin B-12 or aspirin, was studied at 310 K using three-dimensional nuclear magnetic resonance (3D NMR) imaging. The poly(HEMA) was loaded with 5 wt% of the drugs. From the imaging profiles it was observed that incorporation of vitamin B-12 into the polymers rich in HEMA resulted in crack formation at the interface between the rubbery region and the glassy core on sorption of water, although these cracks were 'healed' behind the diffusion front. However, for the copolymers with low HEMA contents and for those containing aspirin, no evidence for similar crack formation was found. For the copolymers loaded with 5 wt% of aspirin or vitamin B-12 the values of the water diffusion coefficients, determined by curve-fitting the relative water concentration profiles from magnetic resonance imaging (MRI) measurements, were found to be smaller than those obtained from a mass uptake study. (C) 2004 Society of Chemical Industry.

Determination of Cole parameters in multiple frequency bioelectrical impedance analysis using only the measurement of impedances

Conventional bioimpedance spectrometers measure resistance and reactance over a range of frequencies and, by application of a mathematical model for an equivalent circuit (the Cole model), estimate resistance at zero and infinite frequencies. Fitting of the experimental data to the model is accomplished by iterative, nonlinear curve fitting. An alternative fitting method is described that uses only the magnitude of the measured impedances at four selected frequencies. The two methods showed excellent agreement when compared using data obtained both from measurements of equivalent circuits and of humans. These results suggest that operational equivalence to a technically complex, frequency-scanning, phase-sensitive BIS analyser could be achieved from a simple four-frequency, impedance-only analyser.

Evaluation of the indirect measures of rock brittleness and fracture toughness in rock cutting

Performance prediction models for partial face mechanical excavators, when developed in laboratory conditions, depend on relating the results of a set of rock property tests and indices to specific cutting energy (SE) for various rock types. There exist some studies in the literature aiming to correlate the geotechnical properties of intact rocks with the SE, especially for massive and widely jointed rock environments. However, those including direct and/or indirect measures of rock fracture parameters such as rock brittleness and fracture toughness, along with the other rock parameters expressing different aspects of rock behavior under drag tools (picks), are rather limited. With this study, it was aimed to investigate the relationships between the indirect measures of rock brittleness and fracture toughness and the SE depending on the results of a new and two previous linear rock cutting programmes. Relationships between the SE, rock strength parameters, and the rock index tests have also been investigated in this study. Sandstone samples taken from the different fields around Ankara, Turkey were used in the new testing programme. Detailed mineralogical analyses, petrographic studies, and rock mechanics and rock cutting tests were performed on these selected sandstone specimens. The assessment of rock cuttability was based on the SE. Three different brittleness indices (B1, B2, and B4) were calculated for sandstones samples, whereas a toughness index (T-i), being developed by Atkinson et al.(1), was employed to represent the indirect rock fracture toughness. The relationships between the SE and the large amounts of new data obtained from the mineralogical analyses, petrographic studies, rock mechanics, and linear rock cutting tests were evaluated by using bivariate correlation and curve fitting techniques, variance analysis, and Student's t-test. Rock cutting and rock property testing data that came from well-known studies of McFeat-Smith and Fowell(2) and Roxborough and Philips(3) have also been employed in statistical analyses together with the new data. Laboratory tests and subsequent analyses revealed that there were close correlations between the SE and B4 whereas no statistically significant correlation has been found between the SE and T-i. Uniaxial compressive and Brazilian tensile strengths and Shore scleroscope hardness of sandstones also exhibited strong relationships with the SE. NCB cone indenter test had the greatest influence on the SE among the other engineering properties of rocks, confirming the previous studies in rock cutting and mechanical excavation. Therefore, it was recommended to employ easy-to-use index tests of NCB cone indenter and Shore scleroscope in the estimation of laboratory SE of sandstones ranging from very low to high strengths in the absence of a rock cutting rig to measure it until the easy-to-use universal measures of the rock brittleness and especially the rock fracture toughness, being an intrinsic rock property, are developed.

Reply to Paparazzo

We would like to reply to the comments made by Paparazzo on our recent paper [1] on the “effect of curve fitting parameters on quantitative analysis of Fe0.94O and Fe2O3 using XPS”. There have been many studies on the characterisation of the properties of iron oxide surfaces. The main purpose of writing the paper was to demonstrate the extent to which the selection of input parameters for curve fitting can affect the results of the quantitative analysis, and to use this analysis to develop more consistent, reproducible and quantitative methods of analysis of these data.

Optimizing measurement of subjective amplitude of accommodation with defocus curves

PURPOSE: To determine whether letter sequences and/or lens-presentation order should be randomized when measuring defocus curves and to assess the most appropriate criterion for calculating the subjective amplitude of accommodation (AoA) from defocus curves. SETTING: Eye Clinic, School of Life & Health Sciences, Aston University, Birmingham, United Kingdom. METHODS: Defocus curves (from +3.00 diopters [D] to -3.00 D in 0.50 D steps) for 6 possible combinations of randomized or nonrandomized letter sequences and/or lens-presentation order were measured in a random order in 20 presbyopic subjects. Subjective AoA was calculated from the defocus curves by curve fitting using various published criteria, and each was correlated to subjective push-up AoA. Objective AoA was measured for comparison of blur tolerance and pupil size. RESULTS: Randomization of lens-presentation order and/or letter sequences, or lack of, did not affect the measured defocus curves (P>.05, analysis of variance). The range of defocus that maintains highest achievable visual acuity (allowing for variability of repeated measurement) was better correlated to (r = 0.84) and agreed best with ( 0.50 D) subjective push-up AoA than any other relative or absolute acuity criterion used in previous studies. CONCLUSIONS: Nonrandomized letters and lens presentation on their own did not affect subjective AoA measured by defocus curves, although their combination should be avoided. Quantification of subjective AoA from defocus curves should be standardized to the range of defocus that maintains the best achievable visual acuity.

Data methods in optometry. Part 10: non-linear regression analysis

1. The techniques associated with regression, whether linear or non-linear, are some of the most useful statistical procedures that can be applied in clinical studies in optometry. 2. In some cases, there may be no scientific model of the relationship between X and Y that can be specified in advance and the objective may be to provide a ‘curve of best fit’ for predictive purposes. In such cases, the fitting of a general polynomial type curve may be the best approach. 3. An investigator may have a specific model in mind that relates Y to X and the data may provide a test of this hypothesis. Some of these curves can be reduced to a linear regression by transformation, e.g., the exponential and negative exponential decay curves. 4. In some circumstances, e.g., the asymptotic curve or logistic growth law, a more complex process of curve fitting involving non-linear estimation will be required.

The effect of catchment characteristics on sewage settling velocity grading

A procedure has been developed which measures the settling velocity distribution of particles within a complete sewage sample. The development of the test method included observations of particle and liquid interaction using both synthetic media and sewage. Comparison studies with two other currently used settling velocity test procedures was undertaken. The method is suitable for use with either DWF or storm sewage. Information relating to the catchment characteristics of 35 No. wastewater treatment works was collected from the privatised water companies in England and Wales. 29 No. of these sites were used in an experimental programme to determine the settling velocity grading of 33 No. sewage samples. The collected data were analysed in an attempt to relate the settling velocity distribution to the characteristics of the contributing catchment. Statistical analysis of the catchment data and the measured settling velocity distributions was undertaken. A curve fitting exercise using an S-shaped curve which had the same physical characteristics as the settling velocity distributions was performed. None of these analyses found evidence that the settling velocity distribution of sewage had a significant relationship with the chosen catchment characteristics. The regression equations produced from the statistical analysis cannot be used to assist in the design of separation devices. However, a grading curve envelope was produced, the limits of which were clearly defined for the measured data set. There was no evidence of a relationship between settling velocity grading and the characteristics of the contributing catchment, particularly the catchment area. The present empirical approach to settling tank design cannot be improved upon at present by considering the variation in catchment parameters. This study has provided a basis for future research into the settling velocity measurement and should be of benefit to future workers within this field.

The application of parallel processing methods to vibrational analysis

The trend in modal extraction algorithms is to use all the available frequency response functions data to obtain a global estimate of the natural frequencies, damping ratio and mode shapes. Improvements in transducer and signal processing technology allow the simultaneous measurement of many hundreds of channels of response data. The quantity of data available and the complexity of the extraction algorithms make considerable demands on the available computer power and require a powerful computer or dedicated workstation to perform satisfactorily. An alternative to waiting for faster sequential processors is to implement the algorithm in parallel, for example on a network of Transputers. Parallel architectures are a cost effective means of increasing computational power, and a larger number of response channels would simply require more processors. This thesis considers how two typical modal extraction algorithms, the Rational Fraction Polynomial method and the Ibrahim Time Domain method, may be implemented on a network of transputers. The Rational Fraction Polynomial Method is a well known and robust frequency domain 'curve fitting' algorithm. The Ibrahim Time Domain method is an efficient algorithm that 'curve fits' in the time domain. This thesis reviews the algorithms, considers the problems involved in a parallel implementation, and shows how they were implemented on a real Transputer network.

The slope of the psychometric function and non-stationarity of thresholds in spatiotemporal contrast vision

The slope of the two-interval, forced-choice psychometric function (e.g. the Weibull parameter, ß) provides valuable information about the relationship between contrast sensitivity and signal strength. However, little is known about how or whether ß varies with stimulus parameters such as spatiotemporal frequency and stimulus size and shape. A second unresolved issue concerns the best way to estimate the slope of the psychometric function. For example, if an observer is non-stationary (e.g. their threshold drifts between experimental sessions), ß will be underestimated if curve fitting is performed after collapsing the data across experimental sessions. We measured psychometric functions for 2 experienced observers for 14 different spatiotemporal configurations of pulsed or flickering grating patches and bars on each of 8 days. We found ß ˜ 3 to be fairly constant across almost all conditions, consistent with a fixed nonlinear contrast transducer and/or a constant level of intrinsic stimulus uncertainty (e.g. a square law transducer and a low level of intrinsic uncertainty). Our analysis showed that estimating a single ß from results averaged over several experimental sessions was slightly more accurate than averaging multiple estimates from several experimental sessions. However, the small levels of non-stationarity (SD ˜ 0.8 dB) meant that the difference between the estimates was, in practice, negligible.

Quantitative neuropathology:data collection and statistical analysis

The use of quantitative methods has become increasingly important in the study of neuropathology and especially in neurodegenerative disease. Disorders such as Alzheimer's disease (AD) and the frontotemporal dementias (FTD) are characterized by the formation of discrete, microscopic, pathological lesions which play an important role in pathological diagnosis. This chapter reviews the advantages and limitations of the different methods of quantifying pathological lesions in histological sections including estimates of density, frequency, coverage, and the use of semi-quantitative scores. The sampling strategies by which these quantitative measures can be obtained from histological sections, including plot or quadrat sampling, transect sampling, and point-quarter sampling, are described. In addition, data analysis methods commonly used to analysis quantitative data in neuropathology, including analysis of variance (ANOVA), polynomial curve fitting, multiple regression, classification trees, and principal components analysis (PCA), are discussed. These methods are illustrated with reference to quantitative studies of a variety of neurodegenerative disorders.

Regression analysis of ranked segment parameters for optic nerve head classification:a pilot study

Purpose: To determine whether curve-fitting analysis of the ranked segment distributions of topographic optic nerve head (ONH) parameters, derived using the Heidelberg Retina Tomograph (HRT), provide a more effective statistical descriptor to differentiate the normal from the glaucomatous ONH. Methods: The sample comprised of 22 normal control subjects (mean age 66.9 years; S.D. 7.8) and 22 glaucoma patients (mean age 72.1 years; S.D. 6.9) confirmed by reproducible visual field defects on the Humphrey Field Analyser. Three 10°-images of the ONH were obtained using the HRT. The mean topography image was determined and the HRT software was used to calculate the rim volume, rim area to disc area ratio, normalised rim area to disc area ratio and retinal nerve fibre cross-sectional area for each patient at 10°-sectoral intervals. The values were ranked in descending order, and each ranked-segment curve of ordered values was fitted using the least squares method. Results: There was no difference in disc area between the groups. The group mean cup-disc area ratio was significantly lower in the normal group (0.204 ± 0.16) compared with the glaucoma group (0.533 ± 0.083) (p < 0.001). The visual field indices, mean deviation and corrected pattern S.D., were significantly greater (p < 0.001) in the glaucoma group (-9.09 dB ± 3.3 and 7.91 ± 3.4, respectively) compared with the normal group (-0.15 dB ± 0.9 and 0.95 dB ± 0.8, respectively). Univariate linear regression provided the best overall fit to the ranked segment data. The equation parameters of the regression line manually applied to the normalised rim area-disc area and the rim area-disc area ratio data, correctly classified 100% of normal subjects and glaucoma patients. In this study sample, the regression analysis of ranked segment parameters method was more effective than conventional ranked segment analysis, in which glaucoma patients were misclassified in approximately 50% of cases. Further investigation in larger samples will enable the calculation of confidence intervals for normality. These reference standards will then need to be investigated for an independent sample to fully validate the technique. Conclusions: Using a curve-fitting approach to fit ranked segment curves retains information relating to the topographic nature of neural loss. Such methodology appears to overcome some of the deficiencies of conventional ranked segment analysis, and subject to validation in larger scale studies, may potentially be of clinical utility for detecting and monitoring glaucomatous damage. © 2007 The College of Optometrists.

Short fatigue crack growth behaviour:data analysis effects

Consideration of the influence of test technique and data analysis method is important for data comparison and design purposes. The paper highlights the effects of replication interval, crack growth rate averaging and curve-fitting procedures on crack growth rate results for a Ni-base alloy. It is shown that an upper bound crack growth rate line is not appropriate for use in fatigue design, and that the derivative of a quadratic fit to the a vs N data looks promising. However, this type of averaging, or curve fitting, is not useful in developing an understanding of microstructure/crack tip interactions. For this purpose, simple replica-to-replica growth rate calculations are preferable. © 1988.

Cortical degeneration in frontotemporal lobar degeneration with TDP-43 proteinopathy caused by progranulin gene mutation

Familial frontotemporal lobar degeneration with transactive response (TAR) DNA-binding protein of 43 kDa (TDP-43) proteinopathy (FTLD-TDP) is most commonly caused by progranulin (GRN) gene mutation. To characterize cortical degeneration in these cases, changes in density of the pathology across the cortical laminae of the frontal and temporal lobe were studied in seven cases of FTLD-TDP with GRN mutation using quantitative analysis and polynomial curve fitting. In 50% of gyri studied, neuronal cytoplasmic inclusions (NCI) exhibited a peak of density in the upper cortical laminae. Most frequently, neuronal intranuclear inclusions (NII) and dystrophic neurites (DN) exhibited a density peak in lower and upper laminae, respectively, glial inclusions (GI) being distributed in low densities across all laminae. Abnormally enlarged neurons (EN) were distributed either in the lower laminae or were more uniformly distributed across the cortex. The distribution of all neurons present varied between cases and regions, but most commonly exhibited a bimodal distribution, density peaks occurring in upper and lower laminae. Vacuolation primarily affected the superficial laminae and density of glial cell nuclei increased with distance across the cortex from pia mater to white matter. The densities of the NCI, GI, NII, and DN were not spatially correlated. The laminar distribution of the pathology in GRN mutation cases was similar to previously reported sporadic cases of FTLD-TDP. Hence, pathological changes initiated by GRN mutation, and by other causes in sporadic cases, appear to follow a parallel course resulting in very similar patterns of cortical degeneration in FTLD-TDP.

Determining the change of Brillouin frequency shift by using the similarity matching method

We have proposed a similarity matching method (SMM) to obtain the change of Brillouin frequency shift (BFS), in which the change of BFS can be determined from the frequency difference between detecting spectrum and selected reference spectrum by comparing their similarity. We have also compared three similarity measures in the simulation, which has shown that the correlation coefficient is more accurate to determine the change of BFS. Compared with the other methods of determining the change of BFS, the SMM is more suitable for complex Brillouin spectrum profiles. More precise result and much faster processing speed have been verified in our simulation and experiments. The experimental results have shown that the measurement uncertainty of the BFS has been improved to 0.72 MHz by using the SMM, which is almost one-third of that by using the curve fitting method, and the speed of deriving the BFS change by the SMM is 120 times faster than that by the curve fitting method.