Power Analysis of Exposure Mixture Studies via Monte Carlo Simulations.

Estimating sample size and statistical power is an essential part of a good epidemiological study design. Closed-form formulas exist for simple hypothesis tests but not for advanced statistical methods designed for exposure mixture studies. Estimating power with Monte Carlo simulations is flexible and applicable to these methods. However, it is not straightforward to code a simulation for non-experienced programmers and is often hard for a researcher to manually specify multivariate associations among exposure mixtures to set up a simulation. To simplify this process, we present the R package mpower for power analysis of observational studies of environmental exposure mixtures involving recently-developed mixtures analysis methods. The components within mpower are also versatile enough to accommodate any mixtures methods that will developed in the future. The package allows users to simulate realistic exposure data and mixed-typed covariates based on public data set such as the National Health and Nutrition Examination Survey or other existing data set from prior studies. Users can generate power curves to assess the trade-offs between sample size, effect size, and power of a design. This paper presents tutorials and examples of power analysis using mpower.

crossnma: An R package to synthesize cross-design evidence and cross-format data using network meta-analysis and network meta-regression.

BACKGROUND: Although aggregate data (AD) from randomised clinical trials (RCTs) are used in the majority of network meta-analyses (NMAs), other study designs (e.g., cohort studies and other non-randomised studies, NRS) can be informative about relative treatment effects. The individual participant data (IPD) of the study, when available, are preferred to AD for adjusting for important participant characteristics and to better handle heterogeneity and inconsistency in the network. RESULTS: We developed the R package crossnma to perform cross-format (IPD and AD) and cross-design (RCT and NRS) NMA and network meta-regression (NMR). The models are implemented as Bayesian three-level hierarchical models using Just Another Gibbs Sampler (JAGS) software within the R environment. The R package crossnma includes functions to automatically create the JAGS model, reformat the data (based on user input), assess convergence and summarize the results. We demonstrate the workflow within crossnma by using a network of six trials comparing four treatments. CONCLUSIONS: The R package crossnma enables the user to perform NMA and NMR with different data types in a Bayesian framework and facilitates the inclusion of all types of evidence recognising differences in risk of bias.

synapsis: A Bioconductor Package to Automate the Analysis of Meiotic Double-Strand Break Repair and Crossover Formation.

Immunofluorescent staining is commonly used to generate images to characterize cytological phenotypes. The manual quantification of DNA double-strand breaks and their repair intermediates during meiosis using image data requires a series of subjective steps, from image selection to the counting of particular events per nucleus. Here we describe "synapsis," a bioconductor package, which includes a set of functions to automate the process of identifying meiotic nuclei and quantifying key double-strand break formation and repair events in a rapid, scalable, and reproducible workflow, and compare it to manual user quantification. The software can be extended for other applications in meiosis research, such as incorporating machine learning approaches to categorize meiotic substages.

Establishing large mammal population trends from heterogeneous count data.

Monitoring population trends is pivotal to effective wildlife conservation and management. However, wildlife managers often face many challenges when analyzing time series of census data due to heterogeneities in sampling methodology, strategy, or frequency. We present a three-step method for modeling trends from time series of count data obtained through multiple census methods (aerial or ground census and expert estimates). First, we design a heuristic for constructing credible intervals for all types of animal counts including those which come with no precision measure. Then, we define conversion factors for rendering aerial and ground counts comparable and provide values for broad classes of animals from an extant series of parallel aerial and ground censuses. Lastly, we construct a Bayesian model that takes the reconciled counts as input and estimates the relative growth rates between successive dates while accounting for their precisions. Importantly, we bound the rate of increase to account for the demographic potential of a species. We propose a flow chart for constructing credible intervals for various types of animal counts. We provide estimates of conversion factors for 5 broad classes of species. We describe the Bayesian model for calculating trends, annual rates of population increase, and the associated credible intervals. We develop a bespoke R CRAN package, popbayes, for implementing all the calculations that take the raw counts as input. It produces consistent and reliable estimates of population trends and annual rates of increase. Several examples from real populations of large African mammals illustrate the different features of our method. The approach is well-suited for analyzing population trends for heterogeneous time series and allows a principled use of all the available historical census data. The method is general and flexible and applicable to various other animal species besides African large mammals. It can readily be adapted to test predictions of various hypotheses about drivers of rates of population increase.

Person explanatory multidimensional item response theory with the instrument package in R.

We present the new R package instrument to perform Bayesian estimation of person explanatory multidimensional item response theory. The package implements an exploratory multidimensional item response theory model and a higher-order multidimensional item response theory model, a type of confirmatory multidimensional item response theory. Explanation of person parameters is accomplished by fixed and random effect linear regression models. Estimation is carried out using Hamiltonian Monte Carlo in Stan. In this article, we provide a detailed description of the models; we use the instrument package to demonstrate fitting explanatory item response models with fixed and random effects (i.e., mixed modeling) of person parameters in R; and, we perform a simulation study to evaluate the performance of our implementation of the models.

CleanUpRNAseq: An R/Bioconductor Package for Detecting and Correcting DNA Contamination in RNA-Seq Data.

RNA sequencing (RNA-seq) has become a standard method for profiling gene expression, yet genomic DNA (gDNA) contamination carried over to the sequencing library poses a significant challenge to data integrity. Detecting and correcting this contamination is vital for accurate downstream analyses. Particularly, when RNA samples are scarce and invaluable, it becomes essential not only to identify but also to correct gDNA contamination to maximize the data's utility. However, existing tools capable of correcting gDNA contamination are limited and lack thorough evaluation. To fill the gap, we developed CleanUpRNAseq, which offers a comprehensive set of functionalities for identifying and correcting gDNA-contaminated RNA-seq data. Our package offers three correction methods for unstranded RNA-seq data and a dedicated approach for stranded data. Through rigorous validation on published RNA-seq datasets with known levels of gDNA contamination and real-world RNA-seq data, we demonstrate CleanUpRNAseq's efficacy in detecting and correcting detrimental levels of gDNA contamination across diverse library protocols. CleanUpRNAseq thus serves as a valuable tool for post-alignment quality assessment of RNA-seq data and should be integrated into routine workflows for RNA-seq data analysis. Its incorporation into OneStopRNAseq should significantly bolster the accuracy of gene expression quantification and differential expression analysis of RNA-seq data.

Research Trends and Visualization of Cerebrospinal Fluid Dynamics (2013-2023).

OBJECTIVE: This study aims to analyze CSF dynamics using VOSviewer, CiteSpace, and the Bibliometrix R-package software to identify research hotspots and future directions. METHODS: Search by Web of Science Core Collection Database for related literature on CSF dynamics from 2013 to 2023. Bibliometric and visual analysis of data on number of citations, number of publications, most productive countries and institutions, important authors and journals, time of publication, popular topics and keywords were performed by CiteSpace and VOSviewer. RESULTS AND CONCLUSION: In the field of cerebrospinal fluid (CSF) dynamics, there is an observable upward trend in annual publications. The United States, Japan, and Germany rank as the top three countries in terms of publication output. Copenhagen University, Idaho University, and Zurich University are the leading institutions in terms of research publications.The most frequently published literature over the past decade is "Coupled electrophysiological, hemodynamic, and cerebrospinal fluid oscillations in human sleep." The most prolific authors in this field are Martin and Bryn. Active authors have formed several stable research teams, indicating collaboration, which could benefit from further enhancement across different teams.Key research keywords such as cerebrospinal fluid, hydrocephalus, dynamics, brain, flow, cerebrospinal-fluid, pressure, CSF flow, MRI highlight the focal areas in CSF dynamics research. These keywords represent both the current focus and forefront of research in this field.

Dietaryindex: A User-Friendly and Versatile R Package for Standardizing Dietary Pattern Analysis in Epidemiological and Clinical Studies.

BACKGROUND: Few standardized and open-source tools exist for calculating multiple dietary pattern indexes from dietary intake data in epidemiological and clinical studies. OBJECTIVE: The primary aim is to develop and validate dietaryindex, a user-friendly and versatile R package that standardizes the calculation of dietary indexes. METHODS: Dietaryindex can calculate multiple dietary indexes of high interest in research, including Healthy Eating Index (HEI) - 2020, Alternative Healthy Eating Index 2010, Dietary Approaches to Stop Hypertension Index, Alternate Mediterranean Diet Score, Dietary Inflammatory Index, American Cancer Society 2020 dietary index, and Planetary Health Diet Index from the EAT-Lancet Commission. The package includes generic dietary index calculation functions that accept any dietary assessment with preprocessed serving sizes of food groups and nutrients, as defined by the research group that developed each index. For ease of use and to eliminate any need for data preprocessing, dietaryindex also offers one-step functions that directly reference common datasets and tools, including the National Health and Nutrition Examination Survey (NHANES), Automated Self-Administered 24-Hour Dietary Assessment Tool (ASA24), Diet History Questionnaire III (DHQ3), and Block Food Frequency Questionnaire. At least two independent researchers validated the serving size definitions and scoring algorithms of dietaryindex. RESULTS: In our validation process, dietaryindex demonstrated full accuracy (100%) in all generic functions with two-decimal rounding precision in comparison to hand-calculated results. Similarly, using NHANES 2017-2018 data and ASA24 and DHQ3 example data, the HEI2015 outputs from dietaryindex aligned (99.95% - 100%) with results using the SAS codes from the National Cancer Institute. CONCLUSIONS: Dietaryindex is a user-friendly, versatile, and validated informatics tool for standardized dietary index calculations. We have open-sourced all the validation files and codes with detailed tutorials on GitHub (https://github.com/jamesjiadazhan/dietaryindex).

SurvdigitizeR: an algorithm for automated survival curve digitization.

BACKGROUND: Decision analytic models and meta-analyses often rely on survival probabilities that are digitized from published Kaplan-Meier (KM) curves. However, manually extracting these probabilities from KM curves is time-consuming, expensive, and error-prone. We developed an efficient and accurate algorithm that automates extraction of survival probabilities from KM curves. METHODS: The automated digitization algorithm processes images from a JPG or PNG format, converts them in their hue, saturation, and lightness scale and uses optical character recognition to detect axis location and labels. It also uses a k-medoids clustering algorithm to separate multiple overlapping curves on the same figure. To validate performance, we generated survival plots form random time-to-event data from a sample size of 25, 50, 150, and 250, 1000 individuals split into 1,2, or 3 treatment arms. We assumed an exponential distribution and applied random censoring. We compared automated digitization and manual digitization performed by well-trained researchers. We calculated the root mean squared error (RMSE) at 100-time points for both methods. The algorithm's performance was also evaluated by Bland-Altman analysis for the agreement between automated and manual digitization on a real-world set of published KM curves. RESULTS: The automated digitizer accurately identified survival probabilities over time in the simulated KM curves. The average RMSE for automated digitization was 0.012, while manual digitization had an average RMSE of 0.014. Its performance was negatively correlated with the number of curves in a figure and the presence of censoring markers. In real-world scenarios, automated digitization and manual digitization showed very close agreement. CONCLUSIONS: The algorithm streamlines the digitization process and requires minimal user input. It effectively digitized KM curves in simulated and real-world scenarios, demonstrating accuracy comparable to conventional manual digitization. The algorithm has been developed as an open-source R package and as a Shiny application and is available on GitHub: https://github.com/Pechli-Lab/SurvdigitizeR and https://pechlilab.shinyapps.io/SurvdigitizeR/ .

FnR: R package for computing inbreeding and numerator relationship coefficients.

BACKGROUND: Inbreeding and relationship coefficients are essential for conservation and breeding programs. Whether dealing with a small conserved population or a large commercial population, monitoring the inbreeding rate and designing mating plans that minimize the inbreeding rate and maximize the effective population size is important. Free, open-source, and efficient software may greatly contribute to conservation and breeding programs and help students and researchers. Efficient methods exist for calculating inbreeding coefficients. Therefore, an efficient way of calculating the numerator relationship coefficients is via the inbreeding coefficients. i.e., the relationship coefficient between parents is twice the inbreeding coefficient of their progeny. A dummy progeny is introduced where no progeny exists for a pair of individuals. Calculating inbreeding coefficients is very fast, and finding whether a pair of individuals has a progeny and picking one from multiple progenies is computationally more demanding. Therefore, the R package introduces a dummy progeny for any pair of individuals whose relationship coefficient is of interest, whether they have a progeny or not. RESULTS: Runtime and peak memory usage were benchmarked for calculating relationship coefficients between two sets of 250 and 800 animals (200,000 dummy progenies) from a pedigree of 2,721,252 animals. The program performed efficiently (200,000 relationship coefficients, which involved calculating 2,721,252 + 200,000 inbreeding coefficients) within 3:45 (mm:ss). Providing the inbreeding coefficients (for real animals), the runtime was reduced to 1:08. Furthermore, providing the diagonal elements of D in A = TDT ' (d), the runtime was reduced to 54s. All the analyses were performed on a machine with a total memory size of 1 GB. CONCLUSIONS: The R package FnR is free and open-source software with implications in conservation and breeding programs. It proved to be time and memory efficient for large populations and many dummy progenies. Calculation of inbreeding coefficients can be resumed for new animals in the pedigree. Thus, saving the latest inbreeding coefficient estimates is recommended. Calculation of d coefficients (from scratch) was very fast, and there was limited value in storing those for future use.

Modeling the Depth of Cellular Dormancy from RNA-Sequencing Data.

High-throughput transcriptome RNA sequencing is a powerful tool for understanding dynamic biological processes. Here, we present a computational framework, implemented in an R package QDSWorkflow, to characterize heterogeneous cellular dormancy depth using RNA-sequencing data from bulk samples and single cells.

Mime: A flexible machine-learning framework to construct and visualize models for clinical characteristics prediction and feature selection.

The widespread use of high-throughput sequencing technologies has revolutionized the understanding of biology and cancer heterogeneity. Recently, several machine-learning models based on transcriptional data have been developed to accurately predict patients' outcome and clinical response. However, an open-source R package covering state-of-the-art machine-learning algorithms for user-friendly access has yet to be developed. Thus, we proposed a flexible computational framework to construct a machine learning-based integration model with elegant performance (Mime). Mime streamlines the process of developing predictive models with high accuracy, leveraging complex datasets to identify critical genes associated with prognosis. An in silico combined model based on de novo PIEZO1-associated signatures constructed by Mime demonstrated high accuracy in predicting the outcomes of patients compared with other published models. Furthermore, the PIEZO1-associated signatures could also precisely infer immunotherapy response by applying different algorithms in Mime. Finally, SDC1 selected from the PIEZO1-associated signatures demonstrated high potential as a glioma target. Taken together, our package provides a user-friendly solution for constructing machine learning-based integration models and will be greatly expanded to provide valuable insights into current fields. The Mime package is available on GitHub (https://github.com/l-magnificence/Mime).

Interactive medical and safety monitoring in clinical trials with clinDataReview: a validated and open-source reporting tool.

Continuous medical and safety monitoring of subject data during a clinical trial is a critical part of evaluating the safety of trial participants and as such is governed by protocol procedures and regulatory guidelines to meet the trial's intended objectives. We present an open-source validated graphical tool (clinDataReview R package) which provides access to the trial data with drill-down to individual patient profiles. The tool incorporates functionalities that facilitate detection of error and data inconsistencies requiring follow-up. It supports regular medical monitoring and oversight as well as safety monitoring committees with interactive tables and listings alongside graphical visualizations of the primary safety data in reports. An implementation example is given where the tool is used to deliver validated outputs following FDA/EMA guidelines. As such, this tool enables a more efficient, interactive, and reproducible review of safety data collected during an ongoing clinical trial.

GasanalyzeR: advancing reproducible research using a new R package for photosynthesis data workflows.

The analysis of photosynthetic traits has become an integral part of plant (eco-)physiology. Many of these characteristics are not directly measured, but calculated from combinations of several, more direct, measurements. The calculations of such derived variables are based on underlying physical models and may use additional constants or assumed values. Commercially available gas-exchange instruments typically report such derived variables, but the available implementations use different definitions and assumptions. Moreover, no software is currently available to allow a fully scripted and reproducible workflow that includes importing data, pre-processing and recalculating derived quantities. The R package gasanalyzer aims to address these issues by providing methods to import data from different instruments, by translating photosynthetic variables to a standardized nomenclature, and by optionally recalculating derived quantities using standardized equations. In addition, the package facilitates performing sensitivity analyses on variables or assumptions used in the calculations to allow researchers to better assess the robustness of the results. The use of the package and how to perform sensitivity analyses are demonstrated using three different examples.

allofus: an R package to facilitate use of the All of Us Researcher Workbench.

OBJECTIVES: Despite easy-to-use tools like the Cohort Builder, using All of Us Research Program data for complex research questions requires a relatively high level of technical expertise. We aimed to increase research and training capacity and reduce barriers to entry for the All of Us community through an R package, allofus. In this article, we describe functions that address common challenges we encountered while working with All of Us Research Program data, and we demonstrate this functionality with an example of creating a cohort of All of Us participants by synthesizing electronic health record and survey data with time dependencies. TARGET AUDIENCE: All of Us Research Program data are widely available to health researchers. The allofus R package is aimed at a wide range of researchers who wish to conduct complex analyses using best practices for reproducibility and transparency, and who have a range of experience using R. Because the All of Us data are transformed into the Observational Medical Outcomes Partnership Common Data Model (OMOP CDM), researchers familiar with existing OMOP CDM tools or who wish to conduct network studies in conjunction with other OMOP CDM data will also find value in the package. SCOPE: We developed an initial set of functions that solve problems we experienced across survey and electronic health record data in our own research and in mentoring student projects. The package will continue to grow and develop with the All of Us Research Program. The allofus R package can help build community research capacity by increasing access to the All of Us Research Program data, the efficiency of its use, and the rigor and reproducibility of the resulting research.

PCAS: An Integrated Tool for Multi-Dimensional Cancer Research Utilizing Clinical Proteomic Tumor Analysis Consortium Data.

Proteomics offers a robust method for quantifying proteins and elucidating their roles in cellular functions, surpassing the insights provided by transcriptomics. The Clinical Proteomic Tumor Analysis Consortium database, enriched with comprehensive cancer proteomics data including phosphorylation and ubiquitination profiles, alongside transcriptomics data from the Genomic Data Commons, allow for integrative molecular studies of cancer. The ProteoCancer Analysis Suite (PCAS), our newly developed R package and Shinyapp, leverages these resources to facilitate in-depth analyses of proteomics, phosphoproteomics, and transcriptomics, enhancing our understanding of the tumor microenvironment through features like immune infiltration and drug sensitivity analysis. This tool aids in identifying critical signaling pathways and therapeutic targets, particularly through its detailed phosphoproteomic analysis. To demonstrate the functionality of the PCAS, we conducted an analysis of GAPDH across multiple cancer types, revealing a significant upregulation of protein levels, which is consistent with its important biological and clinical significance in tumors, as indicated in our prior research. Further experiments were used to validate the findings performed using the tool. In conclusion, the PCAS is a powerful and valuable tool for conducting comprehensive proteomic analyses, significantly enhancing our ability to uncover oncogenic mechanisms and identify potential therapeutic targets in cancer research.

Alternative splicing coupled to nonsense-mediated decay coordinates downregulation of non-neuronal genes in developing mouse neurons.

BACKGROUND: The functional coupling between alternative pre-mRNA splicing (AS) and the mRNA quality control mechanism called nonsense-mediated decay (NMD) can modulate transcript abundance. Previous studies have identified several examples of such a regulation in developing neurons. However, the systems-level effects of AS-NMD in this context are poorly understood. RESULTS: We developed an R package, factR2, which offers a comprehensive suite of AS-NMD analysis functions. Using this tool, we conducted a longitudinal analysis of gene expression in pluripotent stem cells undergoing induced neuronal differentiation. Our analysis uncovers hundreds of AS-NMD events with significant potential to regulate gene expression. Notably, this regulation is significantly overrepresented in specific functional groups of developmentally downregulated genes. Particularly strong association with gene downregulation is detected for alternative cassette exons stimulating NMD upon their inclusion into mature mRNA. By combining bioinformatic analyses with CRISPR/Cas9 genome editing and other experimental approaches we show that NMD-stimulating cassette exons regulated by the RNA-binding protein PTBP1 dampen the expression of their genes in developing neurons. We also provided evidence that the inclusion of NMD-stimulating cassette exons into mature mRNAs is temporally coordinated with NMD-independent gene repression mechanisms. CONCLUSIONS: Our study provides an accessible workflow for the discovery and prioritization of AS-NMD targets. It further argues that the AS-NMD pathway plays a widespread role in developing neurons by facilitating the downregulation of functionally related non-neuronal genes.

A decade of insight: bibliometric analysis of gut microbiota's role in osteoporosis (2014-2024).

Purpose: Osteoporosis represents a profound challenge to public health, underscoring the critical need to dissect its complex etiology and identify viable targets for intervention. Within this context, the gut microbiota has emerged as a focal point of research due to its profound influence on bone metabolism. Despite this growing interest, the literature has yet to see a bibliometric study addressing the gut microbiota's contribution to both the development and management of osteoporosis. This study aims to fill this gap through an exhaustive bibliometric analysis. Our objective is to uncover current research hotspots, delineate key themes, and identify future research trends. In doing so, we hope to provide direction for future studies and the development of innovative treatment methods. Methods: Relevant publications in this field were retrieved from the Web of Science Core Collection database. We used VOSviewer, CiteSpace, an online analysis platform and the R package "Bibliometrix" for bibliometric analysis. Results: A total of 529 publications (including 351 articles and 178 reviews) from 61 countries, 881 institutions, were included in this study. China leads in publication volume and boast the highest cumulative citation. Shanghai Jiao Tong University and Southern Medical University are the leading research institutions in this field. Nutrients contributed the largest number of articles, and J Bone Miner Res is the most co-cited journal. Of the 3,166 scholars who participated in the study, Ohlsson C had the largest number of articles. Li YJ is the most co-cited author. "Probiotics" and "inflammation" are the keywords in the research. Conclusion: This is the first bibliometric analysis of gut microbiota in osteoporosis. We explored current research status in recent years and identified frontiers and hot spots in this research field. We investigate the impact of gut microbiome dysregulation and its associated inflammation on OP progression, a topic that has garnered international research interest in recent years. Additionally, our study delves into the potential of fecal microbiota transplantation or specific dietary interventions as promising avenues for future research, which can provide reference for the researchers who focus on this research filed.

Global trends and hotspots in pain associated with bipolar disorder in the last 20 years: a bibliometric analysis.

Purpose: The prevalence of comorbid pain and Bipolar Disorder in clinical practice continues to be high, with an increasing number of related publications. However, no study has used bibliometric methods to analyze the research progress and knowledge structure in this field. Our research is dedicated to systematically exploring the global trends and focal points in scientific research on pain comorbidity with bipolar disorder from 2003 to 2023, with the goal of contributing to the field. Methods: Relevant publications in this field were retrieved from the Web of Science core collection database (WOSSCC). And we used VOSviewer, CiteSpace, and the R package "Bibliometrix" for bibliometric analysis. Results: A total of 485 publications (including 360 articles and 125 reviews) from 66 countries, 1019 institutions, were included in this study. Univ Toront and Kings Coll London are the leading research institutions in this field. J Affect Disorders contributed the largest number of articles, and is the most co-cited journal. Of the 2,537 scholars who participated in the study, Stubbs B, Vancampfort D, and Abdin E had the largest number of articles. Stubbs B is the most co-cited author. "chronic pain," "neuropathic pain," "psychological pain" are the keywords in the research. Conclusion: This is the first bibliometric analysis of pain-related bipolar disorder. There is growing interest in the area of pain and comorbid bipolar disorder. Focusing on different types of pain in bipolar disorder and emphasizing pain management in bipolar disorder are research hotspots and future trends. The study of pain related bipolar disorder still has significant potential for development, and we look forward to more high-quality research in the future.

CureAuxSP: An R package for estimating mixture cure models with auxiliary survival probabilities.

BACKGROUND AND OBJECTIVE: There is a rising interest in exploiting aggregate information from external medical studies to enhance the statistical analysis of a modestly sized internal dataset. Currently available software packages for analyzing survival data with a cure fraction ignore the potentially available auxiliary information. This paper aims at filling this gap by developing a new R package CureAuxSP that can include subgroup survival probabilities extracted outside into an interested internal survival dataset. METHODS: The newly developed R package CureAuxSP provides an efficient approach for information synthesis under the mixture cure models, including Cox proportional hazards mixture cure model and the accelerated failure time mixture cure model as special cases. It focuses on synthesizing subgroup survival probabilities at multiple time points and the underlying method development lies in the control variate technique. Evaluation of homogeneity assumption based on a test statistic can be automatically carried out by our package and if heterogeneity does exist, the original outputs can be further refined adaptively. RESULTS: The R package CureAuxSP provides a main function SMC.AxuSP() that helps us adaptively incorporate external subgroup survival probabilities into the analysis of an internal survival data. We also provide another function Print.SMC.AuxSP() for printing the results with a better presentation. Detailed usages are described, and implementations are illustrated with numerical examples, including a simulated dataset with a well-designed data generating process and a real breast cancer dataset. Substantial efficiency gain can be observed by our results. CONCLUSIONS: Our R package CureAuxSP can make the wide applications of utilizing auxiliary information possible. It is anticipated that the performance of mixture cure models can be improved for the survival data with a cure fraction, especially for those with small sample sizes.