Estimation of Sample Size and Power For Quantile Regression

Quantile regression (QR) was first introduced by Roger Koenker and Gilbert Bassett in 1978. It is robust to outliers which affect least squares estimator on a large scale in linear regression. Instead of modeling mean of the response, QR provides an alternative way to model the relationship between quantiles of the response and covariates. Therefore, QR can be widely used to solve problems in econometrics, environmental sciences and health sciences. Sample size is an important factor in the planning stage of experimental design and observational studies. In ordinary linear regression, sample size may be determined based on either precision analysis or power analysis with closed form formulas. There are also methods that calculate sample size based on precision analysis for QR like C.Jennen-Steinmetz and S.Wellek (2005). A method to estimate sample size for QR based on power analysis was proposed by Shao and Wang (2009). In this paper, a new method is proposed to calculate sample size based on power analysis under hypothesis test of covariate effects. Even though error distribution assumption is not necessary for QR analysis itself, researchers have to make assumptions of error distribution and covariate structure in the planning stage of a study to obtain a reasonable estimate of sample size. In this project, both parametric and nonparametric methods are provided to estimate error distribution. Since the method proposed can be implemented in R, user is able to choose either parametric distribution or nonparametric kernel density estimation for error distribution. User also needs to specify the covariate structure and effect size to carry out sample size and power calculation. The performance of the method proposed is further evaluated using numerical simulation. The results suggest that the sample sizes obtained from our method provide empirical powers that are closed to the nominal power level, for example, 80%.

An R-package for the Estimation and Testing of Multiple Covariates and Biomarker Interactions for Survival Data Based on Local Partial Likelihood

When we study the variables that a ffect survival time, we usually estimate their eff ects by the Cox regression model. In biomedical research, e ffects of the covariates are often modi ed by a biomarker variable. This leads to covariates-biomarker interactions. Here biomarker is an objective measurement of the patient characteristics at baseline. Liu et al. (2015) has built up a local partial likelihood bootstrap model to estimate and test this interaction e ffect of covariates and biomarker, but the R code developed by Liu et al. (2015) can only handle one variable and one interaction term and can not t the model with adjustment to nuisance variables. In this project, we expand the model to allow adjustment to nuisance variables, expand the R code to take any chosen interaction terms, and we set up many parameters for users to customize their research. We also build up an R package called "lplb" to integrate the complex computations into a simple interface. We conduct numerical simulation to show that the new method has excellent fi nite sample properties under both the null and alternative hypothesis. We also applied the method to analyze data from a prostate cancer clinical trial with acid phosphatase (AP) biomarker.

Time series wind power forecasting based on variant Gaussian process and TLBO

Due to the variability and stochastic nature of wind power system, accurate wind power forecasting has an important role in developing reliable and economic power system operation and control strategies. As wind variability is stochastic, Gaussian Process regression has recently been introduced to capture the randomness of wind energy. However, the disadvantages of Gaussian Process regression include its computation complexity and incapability to adapt to time varying time-series systems. A variant Gaussian Process for time series forecasting is introduced in this study to address these issues. This new method is shown to be capable of reducing computational complexity and increasing prediction accuracy. It is further proved that the forecasting result converges as the number of available data approaches innite. Further, a teaching learning based optimization (TLBO) method is used to train the model and to accelerate
the learning rate. The proposed modelling and optimization method is applied to forecast both the wind power generation of Ireland and that from a single wind farm to show the eectiveness of the proposed method.

Gas Solubility and Vapor Pressure of Electrolytes Based on Alkylcarbonates Containing Litfsi, Lifap, and LiPF6 Lithium Salt: Measurement and Prediction

Most liquid electrolytes used in commercial lithium-ion batteries are composed by alkylcarbonate mixture containing lithium salt. The decomposition of these solvents by oxidation or reduction during cycling of the cell, induce generation of gases (CO2, CH4, C2H4, CO …) increasing of pressure in the sealed cell, which causes a safety problem [1]. The prior understanding of parameters, such as structure and nature of salt, temperature pressure, concentration, salting effects and solvation parameters, which influence gas solubility and vapor pressure of electrolytes is required to formulate safer and suitable electrolytes especially at high temperature.

We present in this work the CO2, CH4, C2H4, CO solubility in different pure alkyl-carbonate solvents (PC, DMC, EMC, DEC) and their binary or ternary mixtures as well as the effect of temperature and lithium salt LiX (X = LiPF6, LiTFSI or LiFAP) structure and concentration on these properties. Furthermore, in order to understand parameters that influence the choice of the structure of the solvents and their ability to dissolve gas through the addition of a salt, we firstly analyzed experimentally the transport properties (Self diffusion coefficient (D), fluidity (h-1), and conductivity (s) and lithium transport number (tLi) using the Stock-Einstein, and extended Jones-Dole equations [2]. Furthermore, measured data for the of CO2, C2H4, CH4 and CO solubility in pure alkylcarbonates and their mixtures containing LiPF6; LiFAP; LiTFSI salt, are reported as a function of temperature and concentration in salt. Based on experimental solubility data, the Henry’s law constant of gases in these solvents and electrolytes was then deduced and compared with values predicted by using COSMO-RS methodology within COSMOthermX software. From these results, the molar thermodynamic functions of dissolution such as the standard Gibbs energy, the enthalpy, and the entropy, as well as the mixing enthalpy of the solvents and electrolytes with the gases in its hypothetical liquid state were calculated and discussed [3]. Finally, the analysis of the CO2 solubility variations with the salt addition was then evaluated by determining specific ion parameters Hi by using the Setchenov coefficients in solution. This study showed that the gas solubility is entropy driven and can been influenced by the shape, charge density, and size of the anions in lithium salt.

References

[1] S.A. Freunberger, Y. Chen, Z. Peng, J.M. Griffin, L.J. Hardwick, F. Bardé, P. Novák, P.G. Bruce, Journal of the American Chemical Society 133 (2011) 8040-8047.

[2] P. Porion, Y.R. Dougassa, C. Tessier, L. El Ouatani, J. Jacquemin, M. Anouti, Electrochimica Acta 114 (2013) 95-104.

[3] Y.R. Dougassa, C. Tessier, L. El Ouatani, M. Anouti, J. Jacquemin, The Journal of Chemical Thermodynamics 61 (2013) 32-44.

Digoxin use after diagnosis of prostate cancer and survival: a population-based cohort study

PURPOSE:
Preclinical studies have shown that digoxin exerts anticancer effects on different cancer cell lines including prostate cancer. A recent observational study has shown that digoxin use was associated with a 25% reduction in prostate cancer risk. The aim of this study was to investigate whether digoxin use after diagnosis of prostate cancer was associated with decreased prostate cancer-specific mortality.
METHODS:
A cohort of 13 134 patients with prostate cancer newly diagnosed from 1998 to 2009 was identified from English cancer registries and linked to the UK Clinical Practice Research Datalink (to provide digoxin and other prescription records) and to the Office of National Statistics mortality data (to identify 2010 prostate cancer-specific deaths). Using time-dependent Cox regression models, unadjusted and adjusted hazard ratios (HR) and 95% confidence intervals (CIs) were calculated for the association between post-diagnostic exposure to digoxin and prostate cancer-specific mortality.
RESULTS:
Overall, 701 (5%) patients with prostate cancer used digoxin after diagnosis. Digoxin use was associated with an increase in prostate cancer-specific mortality before adjustment (HR = 1.59; 95% CI 1.32-1.91), but after adjustment for confounders, the association was attenuated (adjusted HR = 1.13; 95% CI 0.93-1.37) and there was no evidence of a dose response.
CONCLUSIONS:
In this large population-based prostate cancer cohort, there was no evidence of a reduction in prostate cancer-specific mortality with digoxin use after diagnosis.

Hierarchical Demand Response for Peak Minimization Using Dantzig-Wolfe Decomposition

Demand response (DR) algorithms manipulate the energy consumption schedules of controllable loads so as to satisfy grid objectives. Implementation of DR algorithms using a centralized agent can be problematic for scalability reasons, and there are issues related to the privacy of data and robustness to communication failures. Thus, it is desirable to use a scalable decentralized algorithm for the implementation of DR. In this paper, a hierarchical DR scheme is proposed for peak minimization based on Dantzig-Wolfe decomposition (DWD). In addition, a time weighted maximization option is included in the cost function, which improves the quality of service for devices seeking to receive their desired energy sooner rather than later. This paper also demonstrates how the DWD algorithm can be implemented more efficiently through the calculation of the upper and lower cost bounds after each DWD iteration.

Attitudes and Knowledge Concerning Corneal Donation in a Population-Based Sample of Urban Chinese Adults.

PURPOSE: To better understand knowledge and attitudes concerning corneal donation among Chinese adults.
METHODS: Randomly selected residents in predetermined age strata 20 to 60+ years completed home-based questionnaires in each of 12 randomly chosen communities in Guangzhou, southern China.
RESULTS: Among 1217 selected persons, 430 (35.3%) completed the questionnaires (mean age 40.4 yrs, 57.9% female). Refusers were older (44.8 yrs, P < 0.001), but sex did not differ (52.2% female, P = 0.07). Among participants, 175 (40.7%) were willing to donate their corneas (WTD). Differences between WTD and not WTD included donation knowledge score (range, 1-12) [WTD (SD) 6.91 ± 2.21, not WTD 5.62 ± 2.43, P < 0.001]; having discussed donation (WTD 26.3%, not WTD 8.63%, P < 0.001); viewing donation as unpopular (WTD 88.0%, not WTD 96.5%, P = 0.001); and feeling donation "damages the body" (WTD 15.4%, not WTD 25.7%, P = 0.013). Associated significantly with WTD in multiple regression models were higher knowledge score [odds ratio (OR) = 1.18, 95% confidence interval (CI), 1.04-1.32, P = 0.008]; not feeling donation "damages the body" (OR = 1.91, 95% CI, 1.07-3.43, P = 0.030); and willingness to discuss donation (OR = 10.6, 95% CI, 3.35-33.9, P < 0.001). WTD did not differ by age (>60 yrs: 22/51, 43.1%; ≤60 yrs: 153/379, 40.4%, P = 0.706). Assuming all those refusing the survey would not donate, 14.4% (175/1217) were WTD for themselves, though only 7.1% (86/1217) would do so on behalf of a family member if they did not know the deceased's preference.
CONCLUSIONS: Interventions to increase knowledge and promote discussions about donation, and policies allowing widespread expression of donation preference, are needed in this setting.

New Method Based on the UNIFAC–VISCO Model for the Estimation of Ionic Liquids Viscosity Using the Experimental Data Recommended by Mathematical Gnostics

The viscosity of ionic liquids (ILs) has been modeled as a function of temperature and at atmospheric pressure using a new method based on the UNIFAC–VISCO method. This model extends the calculations previously reported by our group (see Zhao et al. J. Chem. Eng. Data 2016, 61, 2160–2169) which used 154 experimental viscosity data points of 25 ionic liquids for regression of a set of binary interaction parameters and ion Vogel–Fulcher–Tammann (VFT) parameters. Discrepancies in the experimental data of the same IL affect the quality of the correlation and thus the development of the predictive method. In this work, mathematical gnostics was used to analyze the experimental data from different sources and recommend one set of reliable data for each IL. These recommended data (totally 819 data points) for 70 ILs were correlated using this model to obtain an extended set of binary interaction parameters and ion VFT parameters, with a regression accuracy of 1.4%. In addition, 966 experimental viscosity data points for 11 binary mixtures of ILs were collected from literature to establish this model. All the binary data consist of 128 training data points used for the optimization of binary interaction parameters and 838 test data points used for the comparison of the pure evaluated values. The relative average absolute deviation (RAAD) for training and test is 2.9% and 3.9%, respectively.

Education as a Predictor of Antidepressant and Anxiolytic medication use after Bereavement: a population based record linkage study.

Purpose: Educational attainment has been shown to be positively associated with mental health and a potential buffer to stressful events. One stressful life event likely to affect everyone in their lifetime is bereavement. This paper assesses the effect of educational attainment on mental health post bereavement.
Methods: By utilising large administrative datasets, linking Census returns to death records and prescribed medication data, we analysed the bereavement exposure of 208,332 individuals aged 25-74 years. Two-level multi-level logistic regression models were constructed to determine the likelihood of antidepressant medication use (a proxy of mental ill-health) post bereavement given level of educational attainment.
Results: Individuals who are bereaved have greater antidepressant use than those who are not bereaved, with over a quarter (26.5%) of those bereaved by suicide in receipt of antidepressant medication compared to just 12.4% of those not bereaved. Within individuals bereaved by a sudden death those with a University Degree or higher qualifications are 73% less likely to be in receipt of antidepressant medication compared to those with no qualifications, after full adjustment for demographic, socio-economic and area factors (OR=0.27, 95% CI 0.09,0.75). Higher educational attainment and no qualifications have an equivalent effect for those bereaved by suicide.
Conclusions: Education may protect against poor mental health, as measured by the use of antidepressant medication, post bereavement, except in those bereaved by suicide. This is likely due to the improved cognitive, personal and psychological skills gained from time spent in education.

Statin use and survival in colorectal cancer: Results from a population-based cohort study and an updated systematic review and meta-analysis

BACKGROUND: The aim of this study was to investigate the association between statin use and survival in a population-based colorectal cancer (CRC) cohort and perform an updated meta-analysis to quantify the magnitude of any association.

METHODS: A cohort of 8391 patients with newly diagnosed Dukes' A-C CRC (2009-2012) was identified from the Scottish Cancer Registry. This cohort was linked to the Prescribing Information System and the National Records of Scotland Death Records (until January 2015) to identify 1064 colorectal cancer-specific deaths. Adjusted hazard ratios (HRs) and 95% confidence intervals (CIs) for cancer-specific mortality by statin use were calculated using time dependent Cox regression models. The systematic review included relevant studies published before January 2016. Meta-analysis techniques were used to derive combined HRs for associations between statin use and cancer-specific and overall mortality.

RESULTS: In the Scottish cohort, statin use before diagnosis (HR=0.84, 95% CI 0.75-0.94), but not after (HR=0.90, 95% CI 0.77-1.05), was associated with significantly improved cancer-specific mortality. The systematic review identified 15 relevant studies. In the meta-analysis, there was consistent (I(2)=0%,heterogeneity P=0.57) evidence of a reduction in cancer-specific mortality with statin use before diagnosis in 6 studies (n=86,622, pooled HR=0.82, 95% CI 0.79-0.86) but this association was less apparent and more heterogeneous (I(2)=67%,heterogeneity P=0.03) with statin use after diagnosis in 4 studies (n=19,152, pooled HR=0.84, 95% CI 0.68-1.04).

CONCLUSION: In a Scottish CRC cohort and updated meta-analysis there was some evidence that statin use was associated with improved survival. However, these associations were weak in magnitude and, particularly for post-diagnosis use, varied markedly between studies.

Affect-Based Effects of Simulation Games on Learning

Empirical evidence has demonstrated the benefits of using simulation games in enhancing learning especially in terms of cognitive gains. This is to be expected as the dynamism and non-linearity of simulation games are more cognitively demanding. However, the other effects of simulation games, specifically in terms of learners’ emotions, have not been given much attention and are under-investigated. This study aims to demonstrate that simulation games stimulate positive emotions from learners that help to enhance learning. The study finds that the affect-based constructs of interest, engagement and appreciation are positively correlated to learning. A stepwise multiple regression analysis shows that a model involving interest and engagement are significantly associated with learning. The emotions of learners should be considered in the development of curriculum, and the delivery of learning and teaching as positive emotions enhances learning.

MRI-based assessment of the pineal gland in a large population of children aged 0-5 years and comparison with pineoblastoma: part I, the solid gland.

INTRODUCTION: Differentiation between normal solid (non-cystic) pineal glands and pineal pathologies on brain MRI is difficult. The aim of this study was to assess the size of the solid pineal gland in children (0-5 years) and compare the findings with published pineoblastoma cases. METHODS: We retrospectively analyzed the size (width, height, planimetric area) of solid pineal glands in 184 non-retinoblastoma patients (73 female, 111 male) aged 0-5 years on MRI. The effect of age and gender on gland size was evaluated. Linear regression analysis was performed to analyze the relation between size and age. Ninety-nine percent prediction intervals around the mean were added to construct a normal size range per age, with the upper bound of the predictive interval as the parameter of interest as a cutoff for normalcy. RESULTS: There was no significant interaction of gender and age for all the three pineal gland parameters (width, height, and area). Linear regression analysis gave 99 % upper prediction bounds of 7.9, 4.8, and 25.4 mm(2), respectively, for width, height, and area. The slopes (size increase per month) of each parameter were 0.046, 0.023, and 0.202, respectively. Ninety-three percent (95 % CI 66-100 %) of asymptomatic solid pineoblastomas were larger in size than the 99 % upper bound. CONCLUSION: This study establishes norms for solid pineal gland size in non-retinoblastoma children aged 0-5 years. Knowledge of the size of the normal pineal gland is helpful for detection of pineal gland abnormalities, particularly pineoblastoma.

Scalable unstructured mesh decomposition

As the efficiency of parallel software increases it is becoming common to measure near linear speedup for many applications. For a problem size N on P processors then with software running at O(N=P ) the performance restrictions due to file i/o systems and mesh decomposition running at O(N) become increasingly apparent especially for large P . For distributed memory parallel systems an additional limit to scalability results from the finite memory size available for i/o scatter/gather operations. Simple strategies developed to address the scalability of scatter/gather operations for unstructured mesh based applications have been extended to provide scalable mesh decomposition through the development of a parallel graph partitioning code, JOSTLE [8]. The focus of this work is directed towards the development of generic strategies that can be incorporated into the Computer Aided Parallelisation Tools (CAPTools) project.

School-Based Referrals to the Juvenile Courts: Prevalence and Characteristics in Several States

Effective school discipline practices are essential to keeping schools safe and creating an optimal learning environment. However, the overreliance of exclusionary discipline often removes students from the school setting and deprives them of the opportunity to learn. Previous research has suggested that students are being introduced to the juvenile justice system through the use of school-based juvenile court referrals. In 2011, approximately 1.2 million delinquency cases were referred to the juvenile courts in the United States. Preliminary evidence suggests that an increasing number of these referrals have originated in the schools. This study investigated school-based referrals to the juvenile courts as an element of the School-to-Prison Pipeline (StPP). The likelihood of school-based juvenile court referrals and rate of dismissal of these referrals was examined in several states using data from the National Juvenile Court Data Archives. In addition, the study examined race and special education status as predictors of school-based juvenile court referrals. Descriptive statistics, logistic regression and odds ratio, were used to analyze the data, make conclusions based on the findings and recommend appropriate school discipline practices.

Nonparametric Evaluation of Quantitative Traits in Population-Based Association Studies when the Genetic Model is Unknown

Statistical association between a single nucleotide polymorphism (SNP) genotype and a quantitative trait in genome-wide association studies is usually assessed using a linear regression model, or, in the case of non-normally distributed trait values, using the Kruskal-Wallis test. While linear regression models assume an additive mode of inheritance via equi-distant genotype scores, Kruskal-Wallis test merely tests global differences in trait values associated with the three genotype groups. Both approaches thus exhibit suboptimal power when the underlying inheritance mode is dominant or recessive. Furthermore, these tests do not perform well in the common situations when only a few trait values are available in a rare genotype category (disbalance), or when the values associated with the three genotype categories exhibit unequal variance (variance heterogeneity). We propose a maximum test based on Marcus-type multiple contrast test for relative effect sizes. This test allows model-specific testing of either dominant, additive or recessive mode of inheritance, and it is robust against variance heterogeneity. We show how to obtain mode-specific simultaneous confidence intervals for the relative effect sizes to aid in interpreting the biological relevance of the results. Further, we discuss the use of a related all-pairwise comparisons contrast test with range preserving confidence intervals as an alternative to Kruskal-Wallis heterogeneity test. We applied the proposed maximum test to the Bogalusa Heart Study dataset, and gained a remarkable increase in the power to detect association, particularly for rare genotypes. Our simulation study also demonstrated that the proposed non-parametric tests control family-wise error rate in the presence of non-normality and variance heterogeneity contrary to the standard parametric approaches. We provide a publicly available R library nparcomp that can be used to estimate simultaneous confidence intervals or compatible multiplicity-adjusted p-values associated with the proposed maximum test.