Estimating the mean and standard deviation of concentration distributions using censored samples

Environmental data sets of pollutant concentrations in air, water, and soil frequently include unquantified sample values reported only as being below the analytical method detection limit. These values, referred to as censored values, should be considered in the estimation of distribution parameters as each represents some value of pollutant concentration between zero and the detection limit. Most of the currently accepted methods for estimating the population parameters of environmental data sets containing censored values rely upon the assumption of an underlying normal (or transformed normal) distribution. This assumption can result in unacceptable levels of error in parameter estimation due to the unbounded left tail of the normal distribution. With the beta distribution, which is bounded by the same range of a distribution of concentrations, $\rm\lbrack0\le x\le1\rbrack,$ parameter estimation errors resulting from improper distribution bounds are avoided. This work developed a method that uses the beta distribution to estimate population parameters from censored environmental data sets and evaluated its performance in comparison to currently accepted methods that rely upon an underlying normal (or transformed normal) distribution. Data sets were generated assuming typical values encountered in environmental pollutant evaluation for mean, standard deviation, and number of variates. For each set of model values, data sets were generated assuming that the data was distributed either normally, lognormally, or according to a beta distribution. For varying levels of censoring, two established methods of parameter estimation, regression on normal ordered statistics, and regression on lognormal ordered statistics, were used to estimate the known mean and standard deviation of each data set. The method developed for this study, employing a beta distribution assumption, was also used to estimate parameters and the relative accuracy of all three methods were compared. For data sets of all three distribution types, and for censoring levels up to 50%, the performance of the new method equaled, if not exceeded, the performance of the two established methods. Because of its robustness in parameter estimation regardless of distribution type or censoring level, the method employing the beta distribution should be considered for full development in estimating parameters for censored environmental data sets. ^

Hyalotekite, (Ba,Pb,K)(4)(Ca,Y)(2)Si-8(B,Be)(2)(Si,B)(2)O28F, a Tectosilicate Related to Scapolite: New Structure Refinement, Phase Transitions and a Short-Range Ordered 3B Superstructure

Hyalotekite, a framework silicate of composition (Ba,Pb,K)(4)(Ca,Y)(2)Si-8(B,Be)(2) (Si,B)(2)O28F, is found in relatively high-temperature(greater than or equal to 500 degrees C) Mn skarns at Langban, Sweden, and peralkaline pegmatites at Dara-i-Pioz, Tajikistan. A new paragenesis at Dara-i-Pioz is pegmatite consisting of the Ba borosilicates leucosphenite and tienshanite, as well as caesium kupletskite, aegirine, pyrochlore, microcline and quartz. Hyalotekite has been partially replaced by barylite and danburite. This hyalotekite contains 1.29-1.78 wt.% Y2O3, equivalent to 0.172-0.238 Y pfu or 8-11% Y on the Ca site; its Pb/(Pb+Ba) ratio ranges 0.36-0.44. Electron microprobe F contents of Langban and Dara-i-Pioz hyalotekite range 1.04-1.45 wt.%, consistent with full occupancy of the F site. A new refinement of the structure factor data used in the original structural determination of a Langban hyalotekite resulted in a structural formula, (Pb1.96Ba1.86K0.18)Ca-2(B1.76Be0.24)(Si1.56B0.44)Si8O28F, consistent with chemical data and all cations with positive-definite thermal parameters, although with a slight excess of positive charge (+57.14 as opposed to the ideal +57.00). An unusual feature of the hyalotekite framework is that 4 of 28 oxygens are non-bridging; by merging these 4 oxygens into two, the framework topology of scapolite is obtained. The triclinic symmetry of hyalotekite observed at room temperature is obtained from a hypothetical tetragonal parent structure via a sequence of displacive phase transitions. Some of these transitions are associated with cation ordering, either Pb-Ba ordering in the large cation sites, or B-Be and Si-B ordering on tetrahedral sites. Others are largely displacive but affect the coordination of the large cations (Pb, Ba, K, Ca). High-resolution electron microscopy suggests that the undulatory extinction characteristic of hyalotekite is due to a fine mosaic microstructure. This suggests that at least one of these transitions occurs in nature during cooling, and that it is first order with a large volume change. A diffuse superstructure observed by electron diffraction implies the existence of a further stage of short-range cation ordering which probably involves both (Pb,K)-Ba and (BeSi,BB)-BSi.

Regression of the inflammatory microenvironment of the peritoneal cavity in women with endometriosis by GnRHa treatment.

OBJECTIVE To investigate the effect of gonadotropin-releasing hormone analogues (GnRHa) on the peritoneal fluid microenvironment in women with endometriosis. STUDY DESIGN Peritoneal fluid was collected from 85 women with severe endometriosis (rAFS stage III and IV) during laparoscopic surgery during the proliferative phase. Prior to surgery clinical data were collected. The concentrations of specific markers for endometriosis in the peritoneal fluid were determined using an ELISA and a comparison between peritoneal fluid markers in women using GnRHa and no hormonal treatment was performed using a non-parametric Mann-Whitney U test. RESULTS The study included peritoneal fluid from 39 patients who had been administered GnRHa (Zoladex(®)) in the three months prior to surgery and 46 from women with no hormonal treatment in this period. Concentrations of IL-8, PAPP-A, glycodelin-A and midkine were significantly reduced in the GnRHa treatment group compared to women receiving no hormonal treatment. RANTES, MCP-1, ENA-78, TNF-α, OPG, IP-10 and defensin showed no significant change between the two groups. CONCLUSIONS GnRHa mediate a significant regression in the inflammatory nature of the peritoneal microenvironment in women with endometriosis.

Gynaecologists' awareness of bone healthcare in Switzerland

PRINCIPLES To assess gynaecologists' awareness of bone healthcare in women and the prevalence of application of national recommendations on bone healthcare in Switzerland. METHODS During the annual meeting of the Swiss Society of Gynaecology and Obstetrics 2012, the Swiss Association against Osteoporosis (SVGO) performed standardised interviews with conference participants (n = 210). Questions addressed responsibility for bone healthcare, and whether diagnostic procedures, initiation of bone-specific treatment and follow-up were performed in accordance with SVGO recommendations. RESULTS The majority of respondents were aged 30-50 years (60%), female (70%) and working as board-certified gynaecologists (69%). Ninety-three percent of respondents considered care for bone health as part of the gynaecologist's expertise. As diagnostic procedures, 44% recommended performing bone densitometry (DXA) only, 34% ordered additional laboratory testing. Seventy-two percent of respondents initiated a bone-specific treatment. Predictors for not performing diagnostic procedures and not initiating a bone-specific treatment were physician's age below 30, being a trainee gynaecologist, and working at a university clinic. Particularly, young trainee gynaecologists working at a university clinic were especially unlikely to initiate a bone-specific treatment (regression coefficient = -2.68; odds ratio [OR] 0.069, 95% confidence interval [CI] 0.01-0.61; p = 0.16). Follow-ups were performed by 77% of respondents, but were less likely to be by female physicians (OR 0.27, 95% CI 0.09-0.84; p = 0.024). CONCLUSIONS Although the majority of board-certified gynaecologists follow national recommendations on bone healthcare, current medical training in obstetrics and gynaecology does not sufficiently cover the subject of women's health. However, since this is a small study our findings may not reflect the findings in the total population.

Independent measurement of biogenic silica in sediments by FTIR spectroscopy and PLS regression

We present an independent calibration model for the determination of biogenic silica (BSi) in sediments, developed from analysis of synthetic sediment mixtures and application of Fourier transform infrared spectroscopy (FTIRS) and partial least squares regression (PLSR) modeling. In contrast to current FTIRS applications for quantifying BSi, this new calibration is independent from conventional wet-chemical techniques and their associated measurement uncertainties. This approach also removes the need for developing internal calibrations between the two methods for individual sediments records. For the independent calibration, we produced six series of different synthetic sediment mixtures using two purified diatom extracts, with one extract mixed with quartz sand, calcite, 60/40 quartz/calcite and two different natural sediments, and a second extract mixed with one of the natural sediments. A total of 306 samples—51 samples per series—yielded BSi contents ranging from 0 to 100 %. The resulting PLSR calibration model between the FTIR spectral information and the defined BSi concentration of the synthetic sediment mixtures exhibits a strong cross-validated correlation ( R2cv = 0.97) and a low root-mean square error of cross-validation (RMSECV = 4.7 %). Application of the independent calibration to natural lacustrine and marine sediments yields robust BSi reconstructions. At present, the synthetic mixtures do not include the variation in organic matter that occurs in natural samples, which may explain the somewhat lower prediction accuracy of the calibration model for organic-rich samples.

Assessment of tumor regression of esophageal adenocarcinomas after neoadjuvant chemotherapy: comparison of 2 commonly used scoring approaches.

Histopathologic determination of tumor regression provides important prognostic information for locally advanced gastroesophageal carcinomas after neoadjuvant treatment. Regression grading systems mostly refer to the amount of therapy-induced fibrosis in relation to residual tumor or the estimated percentage of residual tumor in relation to the former tumor site. Although these methods are generally accepted, currently there is no common standard for reporting tumor regression in gastroesophageal cancers. We compared the application of these 2 major principles for assessment of tumor regression: hematoxylin and eosin-stained slides from 89 resection specimens of esophageal adenocarcinomas following neoadjuvant chemotherapy were independently reviewed by 3 pathologists from different institutions. Tumor regression was determined by the 5-tiered Mandard system (fibrosis/tumor relation) and the 4-tiered Becker system (residual tumor in %). Interobserver agreement for the Becker system showed better weighted κ values compared with the Mandard system (0.78 vs. 0.62). Evaluation of the whole embedded tumor site showed improved results (Becker: 0.83; Mandard: 0.73) as compared with only 1 representative slide (Becker: 0.68; Mandard: 0.71). Modification into simplified 3-tiered systems showed comparable interobserver agreement but better prognostic stratification for both systems (log rank Becker: P=0.015; Mandard P=0.03), with independent prognostic impact for overall survival (modified Becker: P=0.011, hazard ratio=3.07; modified Mandard: P=0.023, hazard ratio=2.72). In conclusion, both systems provide substantial to excellent interobserver agreement for estimation of tumor regression after neoadjuvant chemotherapy in esophageal adenocarcinomas. A simple 3-tiered system with the estimation of residual tumor in % (complete regression/1% to 50% residual tumor/>50% residual tumor) maintains the highest reproducibility and prognostic value.

Impact of histologic subtype on cancer-specific survival in patients with renal cell carcinoma and tumor thrombus

BACKGROUND Although different prognostic factors for patients with renal cell carcinoma (RCC) and vena cava tumor thrombus (TT) have been studied, the prognostic value of histologic subtype in these patients remains unclear. OBJECTIVE We analyzed the impact of histologic subtype on cancer-specific survival (CSS). DESIGN, SETTINGS, AND PARTICIPANTS We retrospectively analyzed the records of 1774 patients with RCC and TT who underwent radical nephrectomy and tumor thrombectomy from 1971 to 2012 at 22 US and European centers. OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS Multivariable ordered logistic and Cox regression models were used to quantify the impact of tumor histology on CSS. RESULTS AND LIMITATIONS Overall 5-yr CSS was 53.4% (confidence interval [CI], 50.5-56.2) in the entire group. TT level (according to the Mayo classification of macroscopic venous invasion in RCC) was I in 38.5% of patients, II in 30.6%, III in 17.3%, and IV in 13.5%. Histologic subtypes were clear cell renal cell carcinoma (cRCC) in 89.9% of patients, papillary renal cell carcinoma (pRCC) in 8.5%, and chromophobe RCC in 1.6%. In univariable analysis, pRCC was associated with a significantly worse CSS (p<0.001) compared with cRCC. In multivariable analysis, the presence of pRCC was independently associated with CSS (hazard ratio: 1.62; CI, 1.01-2.61; p<0.05). Higher TT level, positive lymph node status, distant metastasis, and fat invasion were also independently associated with CSS. CONCLUSIONS In our multi-institutional series, we found that patients with pRCC and vena cava TT who underwent radical nephrectomy and tumor thrombectomy had significantly worse cancer-specific outcomes when compared with patients with other histologic subtypes of RCC. We confirmed that higher TT level and fat invasion were independently associated with reduced CSS.

Prognostic value of histopathological regression in 850 neoadjuvantly treated oesophagogastric adenocarcinomas.

BACKGROUND Recently, histopathological tumour regression, prevalence of signet ring cells, and localisation were reported as prognostic factors in neoadjuvantly treated oesophagogastric (junctional and gastric) cancer. This exploratory retrospective study analyses independent prognostic factors within a large patient cohort after preoperative chemotherapy including clinical and histopathological factors. METHODS In all, 850 patients presenting with oesophagogastric cancer staged cT3/4 Nany cM0/x were treated with neoadjuvant chemotherapy followed by resection in two academic centres. Patient data were documented in a prospective database and retrospectively analysed. RESULTS Of all factors prognostic on univariate analysis, only clinical response, complications, ypTNM stage, and R category were independently prognostic (P<0.01) on multivariate analysis. Tumour localisation and signet ring cells were independently prognostic only when investigator-dependent clinical response evaluation was excluded from the multivariate model. Histopathological tumour regression correlates with tumour grading, Laurén classification, clinical response, ypT, ypN, and R categories but was not identified as an independent prognostic factor. Within R0-resected patients only surgical complications and ypTNM stage were independent prognostic factors. CONCLUSIONS Only established prognostic factors like ypTNM stage, R category, and complications were identified as independent prognostic factors in resected patients after neoadjuvant chemotherapy. In contrast, histopathological tumour regression was not found as an independent prognostic marker.

Plotting regression coefficients and other estimates

Graphical display of regression results has become increasingly popular in presentations and in scientific literature because graphs are often much easier to read than tables. Such plots can be produced in Stata by the marginsplot command (see [R] marginsplot). However, while marginsplot is versatile and flexible, it has two major limitations: it can only process results left behind by margins (see [R] margins), and it can handle only one set of results at a time. In this article, I introduce a new command called coefplot that overcomes these limitations. It plots results from any estimation command and combines results from several models into one graph. The default behavior of coefplot is to plot markers for coefficients and horizontal spikes for confidence intervals. However, coefplot can also produce other types of graphs. I illustrate the capabilities of coefplot by using a series of examples.

Sheaf representations of MV-algebras and lattice-ordered abelian groups via duality

We study representations of MV-algebras -- equivalently, unital lattice-ordered abelian groups -- through the lens of Stone-Priestley duality, using canonical extensions as an essential tool. Specifically, the theory of canonical extensions implies that the (Stone-Priestley) dual spaces of MV-algebras carry the structure of topological partial commutative ordered semigroups. We use this structure to obtain two different decompositions of such spaces, one indexed over the prime MV-spectrum, the other over the maximal MV-spectrum. These decompositions yield sheaf representations of MV-algebras, using a new and purely duality-theoretic result that relates certain sheaf representations of distributive lattices to decompositions of their dual spaces. Importantly, the proofs of the MV-algebraic representation theorems that we obtain in this way are distinguished from the existing work on this topic by the following features: (1) we use only basic algebraic facts about MV-algebras; (2) we show that the two aforementioned sheaf representations are special cases of a common result, with potential for generalizations; and (3) we show that these results are strongly related to the structure of the Stone-Priestley duals of MV-algebras. In addition, using our analysis of these decompositions, we prove that MV-algebras with isomorphic underlying lattices have homeomorphic maximal MV-spectra. This result is an MV-algebraic generalization of a classical theorem by Kaplansky stating that two compact Hausdorff spaces are homeomorphic if, and only if, the lattices of continuous [0, 1]-valued functions on the spaces are isomorphic.