Implementation of quasi-least squares With the R package qlspack

Quasi-least squares (QLS) is an alternative method for estimating the correlation parameters within the framework of generalized estimating equations (GEE) that has two main advantages over the moment estimates that are typically applied for GEE: (1) It guarantees a consistent estimate of the correlation parameter and a positive definite estimated correlation matrix, for several correlation structures; and (2) It allows for easier implementation of some correlation structures that have not yet been implemented in the framework of GEE. Furthermore, because QLS is a method in the framework of GEE, existing software can be employed within the QLS algorithm for estimation of the correlation and regression parameters. In this manuscript we describe and demonstrate the user written package qlspack that allows for implementation of QLS in R software. Our package qlspack calls up the geepack package Yan (2002) and Halekoh et al. (2006) to update the estimate of the regression parameter at the current QLS estimate of the correlation parameter; hence, geepack related functions for standard error estimation can be used after implementing qlspack

False Discovery Rate Control for High Dimensional Dependent Data with an Application to Large-Scale Genetic Association Studies

Large-scale genetic association studies are increasingly utilized for identifying novel susceptible genetic variants for complex traits, but there is little consensus on analysis methods for such data. Most commonly used methods include single SNP analysis or haplotype analysis with Bonferroni correction for multiple comparisons. Since the SNPs in typical GWAS are often in linkage disequilibrium (LD), at least locally, Bonferonni correction of multiple comparisons often leads to conservative error control and therefore lower statistical power. Motivated by an application for analysis of data from the genetic association studies, we consider the problem of false discovery rate (FDR) control under the high dimensional multivariate normal model. Using the compound decision rule framework, we develop an optimal joint oracle procedure and propose to use a marginal procedure to approximate the optimal joint optimal procedure. We show that the marginal plug-in procedure is asymptotically optimal under mild conditions. Our results indicate that the multiple testing procedure developed under the independent model is not only valid but also asymptotically optimal for the high dimensional multivariate normal data under some weak dependency. We evaluate various procedures using simulation studies and demonstrate its application to a genome-wide association study of neuroblastoma (NB). The proposed procedure identified a few more genetic variants that are potentially associated with NB than the standard p-value-based FDR controlling procedure

Genetic variants in the calcium signaling pathway genes are associated with cutaneous melanoma-specific survival

Remodeling or deregulation of the calcium signaling pathway is a relevant hallmark of cancer including cutaneous melanoma (CM). In this study, using data from a published genome-wide association study (GWAS) from The University of Texas M.D. Anderson Cancer Center, we assessed the role of 41,377 common single-nucleotide polymorphisms (SNPs) of 167 calcium signaling pathway genes in CM survival. We used another GWAS from Harvard University as the validation dataset. In the single-locus analysis, 1830 SNPs were found to be significantly associated with CM-specific survival (CMSS; P ≤ 0.050 and false-positive report probability ≤ 0.2), of which 9 SNPs were validated in the Harvard study (P ≤ 0.050). Among these, three independent SNPs (i.e. PDE1A rs6750552 T>C, ITPR1 rs6785564 A>G and RYR3 rs2596191 C>A) had a predictive role in CMSS, with a meta-analysis-derived hazards ratio of 1.52 (95% confidence interval = 1.19–1.94, P = 7.21 × 10−4), 0.49 (0.33–0.73, 3.94 × 10−4) and 0.67 (0.53–0.86, 0.0017), respectively. Patients with an increasing number of protective genotypes had remarkably improved CMSS. Additional expression quantitative trait loci analysis showed that these genotypes were also significantly associated with mRNA expression levels of the genes. Taken together, these results may help us to identify prospective biomarkers in the calcium signaling pathway for CM prognosis

A retrospective comparative study on the diagnostic efficacy and the complications: between CassiII rotational core biopsy and core needle biopsy

Accurate pathologic diagnosis and molecular classification of breast mass biopsy tissue is important for determining individualized therapy for (neo)adjuvant systemic therapies for invasive breast cancer. The CassiII rotational core biopsy system is a novel biopsy technique with a guide needle and a “stick-freeze” technology. The comprehensive assessments including the concordance rates of diagnosis and biomarker status between CassiII and core needle biopsy were evaluated in this study. Estrogen receptor (ER), progesterone receptor (PgR), human epidermal growth factor receptor 2 (HER2), and Ki67 were analyzed through immunohistochemistry. In total, 655 patients with breast cancer who underwent surgery after biopsy at Sir Run Run Shaw Hospital between January 2019 to December 2021 were evaluated. The concordance rates (CRs) of malignant surgical specimens with CassiII needle biopsy was significantly high compared with core needle biopsy. Moreover, CassiII needle biopsy had about 20% improvement in sensitivity and about 5% improvement in positive predictive value compared to Core needle biopsy. The characteristics including age and tumor size were identified the risk factors for pathological inconsistencies with core needle biopsies. However, CassiII needle biopsy was associated with tumor diameter only. The CRs of ER, PgR, HER2, and Ki67 using Cassi needle were 98.08% (kappa, 0.941; p<.001), 90.77% (kappa, 0.812; p<.001), 69.62% (kappa, 0.482; p<.001), and 86.92% (kappa, 0.552; p<.001), respectively. Post-biopsy complications with CassiII needle biopsy were also collected. The complications of CassiII needle biopsy including chest stuffiness, pain and subcutaneous ecchymosis are not rare. The underlying mechanism of subcutaneous congestion or hematoma after CassiII needle biopsy might be the larger needle diameter and the effect of temperature on coagulation function. In summary, CassiII needle biopsy is age-independent and has a better accuracy than CNB for distinguishing carcinoma in situ and invasive carcinoma

Duke EQAPOL CMV Intervention Flow Cytometry Data

This data set is collected from the External Quality Assurance Program Oversight Laboratory (EQAPOL) proficiency program by the Duke Immune Profiling Core (DIPC). This data set comes from 11 healthy volunteers who provided blood samples for flow cytometric intra-cellular cytokine staining (ICS) experiments. Blood samples from each individual was used as a negative control (“Costim”) or treated with a peptide mixture from the immunodominant cytomegalovirus (CMV) pp65 protein. Each sample contains approximately 200,000 cells for which 11 attributes, protein markers for discriminating between T cell basic, maturational and functional subsets, have been measured