Variance estimation in ranked set sampling using a concomitant variable

We propose a nonparametric variance estimator when ranked set sampling (RSS) and judgment post stratification (JPS) are applied by measuring a concomitant variable. Our proposed estimator is obtained by conditioning on observed concomitant values and using nonparametric kernel regression.

Survey estimation and marginal effects in Stata

Stata is a general purpose software package that has become popular among various disciplines such as epidemiology, economics, or social sciences. Users like Stata for its scientific approach, its robustness and reliability, and the ease with which its functionality can be extended by user written programs. In this talk I will first give a brief overview of the functionality of Stata and then discuss two specific features: survey estimation and predictive margins/marginal effects. Most surveys are based on complex samples that contain multiple sampling stages, are stratified or clustered, and feature unequal selection probabilities. Standard estimators can produce misleading results in such samples unless the peculiarities of the sampling plan are taken into account. Stata offers survey statistics for complex samples for a wide variety of estimators and supports several variance estimation procedures such as linearization, jackknife, and balanced repeated replication (see Kreuter and Valliant, 2007, Stata Journal 7: 1-21). In the talk I will illustrate these features using applied examples and I will also show how user written commands can be adapted to support complex samples. Complex can also be the models we fit to our data, making it difficult to interpret them, especially in case of nonlinear or non-additive models (Mood, 2010, European Sociological Review 26: 67-82). Stata provides a number of highly useful commands to make results of such models accessible by computing and displaying predictive margins and marginal effects. In my talk I will discuss these commands provide various examples demonstrating their use.

3G network traffic sources measurement and analysis

In this paper, we show statistical analyses of several types of traffic sources in a 3G network, namely voice, video and data sources. For each traffic source type, measurements were collected in order to, on the one hand, gain better understanding of the statistical characteristics of the sources and, on the other hand, enable forecasting traffic behaviour in the network. The latter can be used to estimate service times and quality of service parameters. The probability density function, mean, variance, mean square deviation, skewness and kurtosis of the interarrival times are estimated by Wolfram Mathematica and Crystal Ball statistical tools. Based on evaluation of packet interarrival times, we show how the gamma distribution can be used in network simulations and in evaluation of available capacity in opportunistic systems. As a result, from our analyses, shape and scale parameters of gamma distribution are generated. Data can be applied also in dynamic network configuration in order to avoid potential network congestions or overflows. Copyright © 2013 John Wiley & Sons, Ltd.

LORENZ: Stata module to estimate and display Lorenz curves and concentration curves

Lorenz estimates Lorenz and concentration curves from individual-level data and, optionally, displays the results in a graph. Relative as well as generalized, absolute, unnormalized, or custom-normalized Lorenz or concentration curves are supported, and tools for computing contrasts between different subpopulations or outcome variables are provided. Variance estimation for complex samples is fully supported.

Estimating Lorenz and concentration curves in Stata

Lorenz and concentration curves are widely used tools in inequality research. In this paper I present a new Stata command called -lorenz- that estimates Lorenz and concentration curves from individual-level data and, optionally, displays the results in a graph. The -lorenz- command supports relative as well as generalized, absolute, unnormalized, or custom-normalized Lorenz or concentration curves, and provides tools for computing contrasts between different subpopulations or outcome variables. Variance estimation for complex samples is fully supported.

DEVCON: Stata module to apply the deviation contrast transform to estimation results

devcon transforms the coefficients of 0/1 dummy variables so that they reflect deviations from the "grand mean" rather than deviations from the reference category (the transformed coefficients are equivalent to those obtained by the so called "effects coding") and adds the coefficient for the reference category. The variance-covariance matrix of the estimates is transformed accordingly. The transformed estimated can be used with post estimation procedures. In particular, devcon can be used to solve the identification problem for dummy variable effects in the so-called Blinder-Oaxaca decomposition (see the oaxaca package).

Estimation of pelvic tilt on anteroposterior X-rays--a comparison of six parameters

OBJECTIVE: To compare six different parameters described in literature for estimation of pelvic tilt on an anteroposterior pelvic radiograph and to create a simple nomogram for tilt correction of prosthetic cup version in total hip arthroplasty. DESIGN: Simultaneous anteroposterior and lateral pelvic radiographs are taken routinely in our institution and were analyzed prospectively. The different parameters (including three distances and three ratios) were measured and compared to the actual pelvic tilt on the lateral radiograph using simple linear regression analysis. PATIENTS: One hundred and four consecutive patients (41 men, 63 women with a mean age of 31.7 years, SD 9.2 years, range 15.7-59.1 years) were studied. RESULTS: The strongest correlation between pelvic tilt and one of the six parameters for both men and women was the distance between the upper border of the symphysis and the sacrococcygeal joint. The correlation coefficient was 0.68 for men (P<0.001) and 0.61 for women (P<0.001). Based on this linear correlation, a nomogram was created that enables fast, tilt-corrected cup version measurements in clinical routine use. CONCLUSION: This simple method for correcting variations in pelvic tilt on plain radiographs can potentially improve the radiologist's ability to diagnose and interpret malformations of the acetabulum (particularly acetabular retroversion and excessive acetabular overcoverage) and post-operative orientation of the prosthetic acetabulum.

Deadspace estimation from CO2 versus molar mass measurements in infants

BACKGROUND: Estimation of respiratory deadspace is often based on the CO2 expirogram, however presence of the CO2 sensor increases equipment deadspace, which in turn influences breathing pattern and calculation of lung volume. In addition, it is necessary to correct for the delay between the sensor and flow signals. We propose a new method for estimation of effective deadspace using the molar mass (MM) signal from an ultrasonic flowmeter device, which does not require delay correction. We hypothesize that this estimation is correlated with that calculated from the CO2 signal using the Fowler method. METHODS: Breath-by-breath CO2, MM and flow measurements were made in a group of 77 term-born healthy infants. Fowler deadspace (Vd,Fowler) was calculated after correcting for the flow-dependent delay in the CO2 signal. Deadspace estimated from the MM signal (Vd,MM) was defined as the volume passing through the flowhead between start of expiration and the 10% rise point in MM. RESULTS: Correlation (r = 0.456, P < 0.0001) was found between Vd,MM and Vd,Fowler averaged over all measurements, with a mean difference of -1.4% (95% CI -4.1 to 1.3%). Vd,MM ranged from 6.6 to 11.4 ml between subjects, while Vd,Fowler ranged from 5.9 to 12.0 ml. Mean intra-measurement CV over 5-10 breaths was 7.8 +/- 5.6% for Vd,MM and 7.8 +/- 3.7% for Vd,Fowler. Mean intra-subject CV was 6.0 +/- 4.5% for Vd,MM and 8.3 +/- 5.9% for Vd,Fowler. Correcting for the CO2 signal delay resulted in a 12% difference (P = 0.022) in Vd,Fowler. Vd,MM could be obtained more frequently than Vd,Fowler in infants with CLD, with a high variability. CONCLUSIONS: Use of the MM signal provides a feasible estimate of Fowler deadspace without introducing additional equipment deadspace. The simple calculation without need for delay correction makes individual adjustment for deadspace in FRC measurements possible. This is especially important given the relative large range of deadspace seen in this homogeneous group of infants.

Children's performance estimation in mathematics and science tests over a school year: A pilot study

The metacognitve ability to accurately estimate ones performance in a test, is assumed to be of central importance for initializing task-oriented effort. In addition activating adequate problem-solving strategies, and engaging in efficient error detection and correction. Although school children's' ability to estimate their own performance has been widely investigated, this was mostly done under highly-controlled, experimental set-ups including only one single test occasion. Method: The aim of this study was to investigate this metacognitive ability in the context of real achievement tests in mathematics. Developed and applied by a teacher of a 5th grade class over the course of a school year these tests allowed the exploration of the variability of performance estimation accuracy as a function of test difficulty. Results: Mean performance estimations were generally close to actual performance with somewhat less variability compared to test performance. When grouping the children into three achievement levels, results revealed higher accuracy of performance estimations in the high achievers compared to the low and average achievers. In order to explore the generalization of these findings, analyses were also conducted for the same children's tests in their science classes revealing a very similar pattern of results compared to the domain of mathematics. Discussion and Conclusion: By and large, the present study, in a natural environment, confirmed previous laboratory findings but also offered additional insights into the generalisation and the test dependency of students' performances estimations.

Non-invasive and non-occlusive blood pressure estimation via a chest sensor

The clinical demand for a device to monitor Blood Pressure (BP) in ambulatory scenarios with minimal use of inflation cuffs is increasing. Based on the so-called Pulse Wave Velocity (PWV) principle, this paper introduces and evaluates a novel concept of BP monitor that can be fully integrated within a chest sensor. After a preliminary calibration, the sensor provides non-occlusive beat-by-beat estimations of Mean Arterial Pressure (MAP) by measuring the Pulse Transit Time (PTT) of arterial pressure pulses travelling from the ascending aorta towards the subcutaneous vasculature of the chest. In a cohort of 15 healthy male subjects, a total of 462 simultaneous readings consisting of reference MAP and chest PTT were acquired. Each subject was recorded at three different days: D, D+3 and D+14. Overall, the implemented protocol induced MAP values to range from 80 ± 6 mmHg in baseline, to 107 ± 9 mmHg during isometric handgrip maneuvers. Agreement between reference and chest-sensor MAP values was tested by using intraclass correlation coefficient (ICC = 0.78) and Bland-Altman analysis (mean error = 0.7 mmHg, standard deviation = 5.1 mmHg). The cumulative percentage of MAP values provided by the chest sensor falling within a range of ±5 mmHg compared to reference MAP readings was of 70%, within ±10 mmHg was of 91%, and within ±15mmHg was of 98%. These results point at the fact that the chest sensor complies with the British Hypertension Society (BHS) requirements of Grade A BP monitors, when applied to MAP readings. Grade A performance was maintained even two weeks after having performed the initial subject-dependent calibration. In conclusion, this paper introduces a sensor and a calibration strategy to perform MAP measurements at the chest. The encouraging performance of the presented technique paves the way towards an ambulatory-compliant, continuous and non-occlusive BP monitoring system.

Simple Estimation of Incident HIV Infection Rates in Notification Cohorts Based on Window Periods of Algorithms for Evaluation of Line-Immunoassay Result Patterns

Background Tests for recent infections (TRIs) are important for HIV surveillance. We have shown that a patient's antibody pattern in a confirmatory line immunoassay (Inno-Lia) also yields information on time since infection. We have published algorithms which, with a certain sensitivity and specificity, distinguish between incident (< = 12 months) and older infection. In order to use these algorithms like other TRIs, i.e., based on their windows, we now determined their window periods. Methods We classified Inno-Lia results of 527 treatment-naïve patients with HIV-1 infection < = 12 months according to incidence by 25 algorithms. The time after which all infections were ruled older, i.e. the algorithm's window, was determined by linear regression of the proportion ruled incident in dependence of time since infection. Window-based incident infection rates (IIR) were determined utilizing the relationship ‘Prevalence = Incidence x Duration’ in four annual cohorts of HIV-1 notifications. Results were compared to performance-based IIR also derived from Inno-Lia results, but utilizing the relationship ‘incident = true incident + false incident’ and also to the IIR derived from the BED incidence assay. Results Window periods varied between 45.8 and 130.1 days and correlated well with the algorithms' diagnostic sensitivity (R2 = 0.962; P<0.0001). Among the 25 algorithms, the mean window-based IIR among the 748 notifications of 2005/06 was 0.457 compared to 0.453 obtained for performance-based IIR with a model not correcting for selection bias. Evaluation of BED results using a window of 153 days yielded an IIR of 0.669. Window-based IIR and performance-based IIR increased by 22.4% and respectively 30.6% in 2008, while 2009 and 2010 showed a return to baseline for both methods. Conclusions IIR estimations by window- and performance-based evaluations of Inno-Lia algorithm results were similar and can be used together to assess IIR changes between annual HIV notification cohorts.

Efficient estimation of the total number of acini in adult rat lung

Pulmonary airways are subdivided into conducting and gas-exchanging airways. An acinus is defined as the small tree of gas-exchanging airways, which is fed by the most distal purely conducting airway. Until now a dissector of five consecutive sections or airway casts were used to count acini. We developed a faster method to estimate the number of acini in young adult rats. Right middle lung lobes were critical point dried or paraffin embedded after heavy metal staining and imaged by X-ray micro-CT or synchrotron radiation-based X-rays tomographic microscopy. The entrances of the acini were counted in three-dimensional (3D) stacks of images by scrolling through them and using morphological criteria (airway wall thickness and appearance of alveoli). Segmentation stopper were placed at the acinar entrances for 3D visualizations of the conducting airways. We observed that acinar airways start at various generations and that one transitional bronchiole may serve more than one acinus. A mean of 5612 (±547) acini per lung and a mean airspace volume of 0.907 (±0.108) μL per acinus were estimated. In 60-day-old rats neither the number of acini nor the mean acinar volume did correlate with the body weight or the lung volume.

Image registration for triggered and non-triggered DTI of the human kidney: Reduced variability of diffusion parameter estimation

BACKGROUND: To investigate if non-rigid image-registration reduces motion artifacts in triggered and non-triggered diffusion tensor imaging (DTI) of native kidneys. A secondary aim was to determine, if improvements through registration allow for omitting respiratory-triggering. METHODS: Twenty volunteers underwent coronal DTI of the kidneys with nine b-values (10-700 s/mm2 ) at 3 Tesla. Image-registration was performed using a multimodal nonrigid registration algorithm. Data processing yielded the apparent diffusion coefficient (ADC), the contribution of perfusion (FP ), and the fractional anisotropy (FA). For comparison of the data stability, the root mean square error (RMSE) of the fitting and the standard deviations within the regions of interest (SDROI ) were evaluated. RESULTS: RMSEs decreased significantly after registration for triggered and also for non-triggered scans (P < 0.05). SDROI for ADC, FA, and FP were significantly lower after registration in both medulla and cortex of triggered scans (P < 0.01). Similarly the SDROI of FA and FP decreased significantly in non-triggered scans after registration (P < 0.05). RMSEs were significantly lower in triggered than in non-triggered scans, both with and without registration (P < 0.05). CONCLUSION: Respiratory motion correction by registration of individual echo-planar images leads to clearly reduced signal variations in renal DTI for both triggered and particularly non-triggered scans. Secondarily, the results suggest that respiratory-triggering still seems advantageous.J. Magn. Reson. Imaging 2014. (c) 2014 Wiley Periodicals, Inc.

Localization and Segmentation of 3D Intervertebral Discs in MR Images by Data Driven Estimation

This paper addresses the problem of fully-automatic localization and segmentation of 3D intervertebral discs (IVDs) from MR images. Our method contains two steps, where we first localize the center of each IVD, and then segment IVDs by classifying image pixels around each disc center as foreground (disc) or background. The disc localization is done by estimating the image displacements from a set of randomly sampled 3D image patches to the disc center. The image displacements are estimated by jointly optimizing the training and test displacement values in a data-driven way, where we take into consideration both the training data and the geometric constraint on the test image. After the disc centers are localized, we segment the discs by classifying image pixels around disc centers as background or foreground. The classification is done in a similar data-driven approach as we used for localization, but in this segmentation case we are aiming to estimate the foreground/background probability of each pixel instead of the image displacements. In addition, an extra neighborhood smooth constraint is introduced to enforce the local smoothness of the label field. Our method is validated on 3D T2-weighted turbo spin echo MR images of 35 patients from two different studies. Experiments show that compared to state of the art, our method achieves better or comparable results. Specifically, we achieve for localization a mean error of 1.6-2.0 mm, and for segmentation a mean Dice metric of 85%-88% and a mean surface distance of 1.3-1.4 mm.