PubMed Health⌕ Search

SEARCH · PubMed Health

Results for “bootstrap”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 829 records · Page 46Linked to original sources

Empirical study of supervised gene screening.

BACKGROUND: Microarray studies provide a way of linking variations of phenotypes with their genetic causations. Constructing predictive models using high dimensional microarray measurements usually consists of three steps: (1) unsupervised gene screening; (2) supervised gene screening; and (3) statistical model building. Supervised gene screening based on marginal gene ranking is commonly used to reduce the number of genes in the model building. Various simple statistics, such as t-statistic or signal to noise ratio, have been used to rank genes in the supervised screening. Despite of its extensive usage, statistical study of supervised gene screening remains scarce. Our study is partly motivated by the differences in gene discovery results caused by using different supervised gene screening methods. RESULTS: We investigate concordance and reproducibility of supervised gene screening based on eight commonly used marginal statistics. Concordance is assessed by the relative fractions of overlaps between top ranked genes screened using different marginal statistics. We propose a Bootstrap Reproducibility Index, which measures reproducibility of individual genes under the supervised screening. Empirical studies are based on four public microarray data. We consider the cases where the top 20%, 40% and 60% genes are screened. CONCLUSION: From a gene discovery point of view, the effect of supervised gene screening based on different marginal statistics cannot be ignored. Empirical studies show that (1) genes passed different supervised screenings may be considerably different; (2) concordance may vary, depending on the underlying data structure and percentage of selected genes; (3) evaluated with the Bootstrap Reproducibility Index, genes passed supervised screenings are only moderately reproducible; and (4) concordance cannot be improved by supervised screening based on reproducibility.

Artificial Intelligence↗

Phase analysis of circadian-related genes in two tissues.

BACKGROUND: Recent circadian clock studies using gene expression microarray in two different tissues of mouse have revealed not all circadian-related genes are synchronized in phase or peak expression times across tissues in vivo. Instead, some circadian-related genes may be delayed by 4-8 hrs in peak expression in one tissue relative to the other. These interesting biological observations prompt a statistical question regarding how to distinguish the synchronized genes from genes that are systematically lagged in phase/peak expression time across two tissues. RESULTS: We propose a set of techniques from circular statistics to analyze phase angles of circadian-related genes in two tissues. We first estimate the phases of a cycling gene separately in each tissue, which are then used to estimate the paired angular difference of the phase angles of the gene in the two tissues. These differences are modeled as a mixture of two von Mises distributions which enables us to cluster genes into two groups; one group having synchronized transcripts with the same phase in the two tissues, the other containing transcripts with a discrepancy in phase between the two tissues. For each cluster of genes we assess the association of phases across the tissue types using circular-circular regression. We also develop a bootstrap methodology based on a circular-circular regression model to evaluate the improvement in fit provided by allowing two components versus a one-component von-Mises model. CONCLUSION: We applied our proposed methodologies to the circadian-related genes common to heart and liver tissues in Storch et al. 2, and found that an estimated 80% of circadian-related transcripts common to heart and liver tissues were synchronized in phase, and the other 20% of transcripts were lagged about 8 hours in liver relative to heart. The bootstrap p-value for being one cluster is 0.063, which suggests the possibility of two clusters. Our methodologies can be extended to analyze peak expression times of circadian-related genes across more than two tissues, for example, kidney, heart, liver, and the suprachiasmatic nuclei (SCN) of the hypothalamus.

Animals↗

Accurate and robust phylogeny estimation based on profile distances: a study of the Chlorophyceae (Chlorophyta).

BACKGROUND: In phylogenetic analysis we face the problem that several subclade topologies are known or easily inferred and well supported by bootstrap analysis, but basal branching patterns cannot be unambiguously estimated by the usual methods (maximum parsimony (MP), neighbor-joining (NJ), or maximum likelihood (ML)), nor are they well supported. We represent each subclade by a sequence profile and estimate evolutionary distances between profiles to obtain a matrix of distances between subclades. RESULTS: Our estimator of profile distances generalizes the maximum likelihood estimator of sequence distances. The basal branching pattern can be estimated by any distance-based method, such as neighbor-joining. Our method (profile neighbor-joining, PNJ) then inherits the accuracy and robustness of profiles and the time efficiency of neighbor-joining. CONCLUSIONS: Phylogenetic analysis of Chlorophyceae with traditional methods (MP, NJ, ML and MrBayes) reveals seven well supported subclades, but the methods disagree on the basal branching pattern. The tree reconstructed by our method is better supported and can be confirmed by known morphological characters. Moreover the accuracy is significantly improved as shown by parametric bootstrap.

Animals↗

Implications of molecular characters for the phylogeny of the Microbotryaceae (Basidiomycota: Urediniomycetes).

BACKGROUND: Anther smuts of the basidiomycetous genus Microbotryum on Caryophyllaceae are important model organisms for many biological disciplines. Members of Microbotryum are most commonly found parasitizing the anthers of host plants in the family Caryophyllaceae, however they can also be found on the anthers of members of the Dipsacaceae, Lamiaceae, Lentibulariaceae, and Portulacaceae. Additionally, some members of Microbotryum can be found infecting other organs of mainly Polygonaceae hosts. Based on ITS nrDNA sequences of members of almost all genera in Microbotryaceae, this study aims to resolve the phylogeny of the anther smuts and their relationship to the other members of the family of plant parasites. A multiple analysis strategy was used to correct for the effects of different equally possible ITS sequence alignments on the phylogenetic outcome, which appears to have been neglected in previous studies. RESULTS: The genera of Microbotryaceae were not clearly resolved, but alignment-independent moderate bootstrap support was achieved for a clade containing the majority of the Microbotryum species. The anther parasites appeared in two different well-supported lineages whose interrelationship remained unresolved. Whereas bootstrap support values for some clades were highly vulnerable to alignment conditions, other clades were more robustly supported. The differences in support between the different alignments were much larger than between the phylogenetic optimality criteria applied (maximum parsimony and maximum likelihood). CONCLUSION: The study confirmed, based on a larger dataset than previous work, that the anther smuts on Caryophyllaceae are monophyletic and that there exists a native North American group that diverged from the European clade before the radiation of the European species. Also a second group of anther smuts was revealed, containing parasites on Dipsacaceae, Lamiaceae, and Lentibulariaceae. At least the majority of the parasites of Asteraceae appeared as a monophylum, but delimitations of some species in this group should be reconsidered. Parasitism on Polygonaceae is likely to be the ancestral state for the Microbotryaceae on Eudicot hosts.

Base Sequence↗

Reconstructing recent human phylogenies with forensic STR loci: a statistical approach.

BACKGROUND: Forensic Short Tandem Repeat (STR) loci are effective for the purpose of individual identification, and other forensic applications. Most of these markers have high allelic variability and mutation rate because of which they have limited use in the phylogenetic reconstruction. In the present study, we have carried out a meta-analysis to explore the possibility of using only five STR loci (TPOX, FES, vWA, F13A and Tho1) to carry out phylogenetic assessment based on the allele frequency profile of 20 world population and north Indian Hindus analyzed in the present study. RESULTS: Phylogenetic analysis based on two different approaches - genetic distance and maximum likelihood along with statistical bootstrapping procedure involving 1000 replicates was carried out. The ensuing tree topologies and PC plots were further compared with those obtained in earlier phylogenetic investigations. The compiled database of 21 populations got segregated and finely resolved into three basal clusters with very high bootstrap values corresponding to three geo-ethnic groups of African, Orientals, and Caucasians. CONCLUSION: Based on this study we conclude that if appropriate and logistic statistical approaches are followed then even lesser number of forensic STR loci are powerful enough to reconstruct the recent human phylogenies despite of their relatively high mutation rates.

Forensic Anthropology↗

Resampling methods to reduce the selection bias in genetic effect estimation in genome-wide scans.

Using the simulated data of Problem 2 for Genetic Analysis Workshop 14 (GAW14), we investigated the ability of three bootstrap-based resampling estimators (a shrinkage, an out-of-sample, and a weighted estimator) to reduce the selection bias for genetic effect estimation in genome-wide linkage scans. For the given marker density in the preliminary genome scans (7 cM for microsatellite and 3 cM for SNP), we found that the two sets of markers produce comparable results in terms of power to detect linkage, localization accuracy, and magnitude of test statistic at the peak location. At the locations detected in the scan, application of the three bootstrap-based estimators substantially reduced the upward selection bias in genetic effect estimation for both true and false positives. The relative effectiveness of the estimators depended on the true genetic effect size and the inherent power to detect it. The shrinkage estimator is recommended when the power to detect the disease locus is low. Otherwise, the weighted estimator is recommended.

Chromosome Mapping↗

An improved probability mapping approach to assess genome mosaicism.

BACKGROUND: Maximum likelihood and posterior probability mapping are useful visualization techniques that are used to ascertain the mosaic nature of prokaryotic genomes. However, posterior probabilities, especially when calculated for four-taxon cases, tend to overestimate the support for tree topologies. Furthermore, because of poor taxon sampling four-taxon analyses suffer from sensitivity to the long branch attraction artifact. Here we extend the probability mapping approach by improving taxon sampling of the analyzed datasets, and by using bootstrap support values, a more conservative tool to assess reliability. RESULTS: Quartets of orthologous proteins were complemented with homologs from selected reference genomes. The mapping of bootstrap support values from these extended datasets gives results similar to the original maximum likelihood and posterior probability mapping. The more conservative nature of the plotted support values allows to focus further analyses on those protein families that strongly disagree with the majority or plurality of genes present in the analyzed genomes. CONCLUSION: Posterior probability is a non-conservative measure for support, and posterior probability mapping only provides a quick estimation of phylogenetic information content of four genomes. This approach can be utilized as a pre-screen to select genes that might have been horizontally transferred. Better taxon sampling combined with subtree analyses prevents the inconsistencies associated with four-taxon analyses, but retains the power of visual representation. Nevertheless, a case-by-case inspection of individual multi-taxon phylogenies remains necessary to differentiate unrecognized paralogy and shared phylogenetic reconstruction artifacts from horizontal gene transfer events.

Bacterial Proteins↗

Performance of statistical models to predict mental health and substance abuse cost.

BACKGROUND: Providers use risk-adjustment systems to help manage healthcare costs. Typically, ordinary least squares (OLS) models on either untransformed or log-transformed cost are used. We examine the predictive ability of several statistical models, demonstrate how model choice depends on the goal for the predictive model, and examine whether building models on samples of the data affects model choice. METHODS: Our sample consisted of 525,620 Veterans Health Administration patients with mental health (MH) or substance abuse (SA) diagnoses who incurred costs during fiscal year 1999. We tested two models on a transformation of cost: a Log Normal model and a Square-root Normal model, and three generalized linear models on untransformed cost, defined by distributional assumption and link function: Normal with identity link (OLS); Gamma with log link; and Gamma with square-root link. Risk-adjusters included age, sex, and 12 MH/SA categories. To determine the best model among the entire dataset, predictive ability was evaluated using root mean square error (RMSE), mean absolute prediction error (MAPE), and predictive ratios of predicted to observed cost (PR) among deciles of predicted cost, by comparing point estimates and 95% bias-corrected bootstrap confidence intervals. To study the effect of analyzing a random sample of the population on model choice, we re-computed these statistics using random samples beginning with 5,000 patients and ending with the entire sample. RESULTS: The Square-root Normal model had the lowest estimates of the RMSE and MAPE, with bootstrap confidence intervals that were always lower than those for the other models. The Gamma with square-root link was best as measured by the PRs. The choice of best model could vary if smaller samples were used and the Gamma with square-root link model had convergence problems with small samples. CONCLUSION: Models with square-root transformation or link fit the data best. This function (whether used as transformation or as a link) seems to help deal with the high comorbidity of this population by introducing a form of interaction. The Gamma distribution helps with the long tail of the distribution. However, the Normal distribution is suitable if the correct transformation of the outcome is used.

Aged↗

Using the net benefit regression framework to construct cost-effectiveness acceptability curves: an example using data from a trial of external loop recorders versus Holter monitoring for ambulatory monitoring of "community acquired" syncope.

BACKGROUND: Cost-effectiveness acceptability curves (CEACs) describe the probability that a new treatment or intervention is cost-effective. The net benefit regression framework (NBRF) allows cost-effectiveness analysis to be done in a simple regression framework. The objective of the paper is to illustrate how net benefit regression can be used to construct a CEAC. METHODS: One hundred patients referred for ambulatory monitoring with syncope or presyncope were randomized to a one-month external loop recorder (n = 49) or 48-hour Holter monitor (n = 51). The primary endpoint was symptom-rhythm correlation during monitoring. Direct costs were calculated based on the 2003 Ontario Health Insurance Plan (OHIP) fee schedule combined with hospital case costing of labour, materials, service and overhead costs for diagnostic testing and related equipment. RESULTS: In the loop recorder group, 63.27% of patients (31/49) had symptom recurrence and successful activation, compared to 23.53% in the Holter group (12/51). The cost in US dollars for loop recording was 648.50 dollars and 212.92 dollars for Holter monitoring. The incremental cost-effectiveness ratio (ICER) of the loop recorder was 1,096 dollars per extra successful diagnosis. The probability that the loop recorder was cost-effective compared to the Holter monitor was estimated using net benefit regression and plotted on a CEAC. In a sensitivity analysis, bootstrapping was used to examine the effect of distributional assumptions. CONCLUSION: The NBRF is straightforward to use and interpret. The resulting uncertainty surrounding the regression coefficient relates to the CEAC. When the link from the regression's p-value to the probability of cost-effectiveness is tentative, bootstrapping may be used.

Bayes Theorem↗

Phylogenetic inference from homologous sequence data: minimum topological assumption, strict mutational compatibility consensus tree as the ultimate solution.

BACKGROUND: For the purposes of phylogenetic inference from molecular data sets many different methods are currently offered as alternatives for researchers in phylogenetic systematics. The vast majority of these methods are based on specific topological assumptions relating to the resultant genealogical tree. Each of these has been shown to perform effectively in special conditions and for specific data sets while yielding less reliable results in other instances. Moreover, the majority of the methods include information from homoplastic characters in spite of a universally accepted agreement in their ineffectiveness for phylogenetic inference, which may often lead to inaccuracy and inconsistency. As an alternative to such methods, a strict mutational compatibility consensus tree building method as a universally applicable and reliable method is reported. RESULTS: The analysis of a data set from a previously published experimental phylogeny demonstrates the accuracy of the strict mutational compatibility consensus tree building method and illustrates its potential for obtaining unambiguous and precise results with full resolution. CONCLUSION: The universal applicability of a simplified compatibility method in its algorithmic form for phylogenetic inference is described. Firstly, dismissal of topological assumptions creates a general potential for agreement of inferred with true phylogeny. Second, exclusion of irregular characters from analysis repeatably enables construction of consistent phylogeny. Third, a direct calculation of bootstrap proportion values for individual nodes of the resulting tree is possible rather than their empirical estimation. Finally, guidance is given for empirical assessment of the sample size necessary for full genealogical resolution and significant bootstrap proportions. REVIEWERS: This article was reviewed by Yuri I. Wolf (nominated by Eugene Koonin), Arcady Mushegian and Martijn Huynen.

Journal Article↗

Linking surveillance to action: incorporation of real-time regional data into a medical decision rule.

OBJECTIVE: Broadly, to create a bidirectional communication link between public health surveillance and clinical practice. Specifically, to measure the impact of integrating public health surveillance data into an existing clinical prediction rule. We incorporate data about recent local trends in meningitis epidemiology into a prediction model differentiating aseptic from bacterial meningitis. DESIGN AND MEASUREMENTS: Retrospective analysis of a cohort of all 696 children with meningitis admitted to a large urban pediatric hospital from 1992 to 2000. We modified a published bacterial meningitis score by adding a new epidemiological context adjustor variable. We examined 540 possible rules for this adjustor, varying both the number of aseptic meningitis cases that needed to be seen, and the recent time window in which they were seen. We performed sensitivity analyses with each of 540 possibilities in order to identify the optimal rule--namely, the one that included the most cases of aseptic meningitis without missing additional cases of bacterial meningitis, as compared with the published prediction model. We used bootstrap methods to validate this new score. RESULTS: The optimal rule was found to be: "at least four cases of aseptic meningitis in the previous 10 days." The epidemiological context adjustor based on surveillance of recent cases of meningitis allowed the correct identification of an additional 47 cases (7%) of aseptic meningitis without missing any additional cases of bacterial meningitis. The epidemiological context adjustor was validated, showing significance in 84% of 1,000 bootstrap samples. CONCLUSION: Epidemiological contextual information can improve the performance of a clinical prediction rule. We provide a methodological framework for leveraging regional surveillance data to improve medical decision-making.

Child↗

Utilizing predictions of early prostate-specific antigen failure to optimize patient selection for adjuvant systemic therapy trials.

PURPOSE: Prostate-specific antigen (PSA) failure within 2 years after radical prostatectomy (RP) has been shown to be a clinically significant predictor of distant failure. This study was performed to estimate 2-year PSA failure rates on the basis of readily available clinical and pathologic factors to identify patients for whom effective adjuvant systemic therapy is needed. PATIENTS AND METHODS: A Cox regression multivariable analysis was used to determine whether the percentage of positive prostate biopsies, PSA level, and the pathologic findings at RP in 1,728 men provided clinically relevant information about PSA outcome after RP. A bootstrapping technique with 2,000 replications was used to provide 95% confidence intervals for the predicted 2-year PSA failure rates, which were determined on the basis of the independent clinical and pathologic predictors of PSA outcome. RESULTS: The independent predictors of time to PSA failure included a percentage of positive prostate biopsies of greater than 34% (P: < or =.009), PSA level greater than 10 ng/mL (P: < or =.01), seminal vesicle invasion (P: =. 02), prostatectomy Gleason score of 8 to 10 (P: =.04), and positive surgical margins (P: =.0001). Predictions of 2-year PSA failure rates and bootstrap estimates of the 95% confidence intervals were arranged in a tabular format, stratified by independent clinical and pathologic predictors of PSA outcome. CONCLUSION: Patients who are most likely to benefit from effective adjuvant systemic therapy after RP can be identified using readily available clinical and pathologic data.

Biopsy↗

Population pharmacokinetics of the novel anticancer agent E7070 during four phase I studies: model building and validation.

PURPOSE: N-(3-Chloro-7-indolyl)-1,4-benzenedisulfonamide (E7070) is a novel sulfonamide anticancer agent currently in phase II clinical development for the treatment of solid tumors. Four phase I studies have been finalized, with E7070 administered at four different treatment schedules to identify the maximum-tolerated dose and the dose-limiting toxicities. Pharmacokinetic analyses of all studies revealed E7070 to have nonlinear pharmacokinetics. A population pharmacokinetic model was designed and validated to describe the pharmacokinetics of E7070 at all four treatment schedules and to identify the possible influences of patient characteristics on the pharmacokinetic parameters. PATIENTS AND METHODS: Plasma concentration-time data of all patients (n = 143) were fitted to several pharmacokinetic models using NONMEM. Seventeen covariables were investigated for their relation with individual pharmacokinetic parameters. A bootstrap procedure was performed to check the validity of the model. RESULTS: The data were best described using a three-compartment model with nonlinear distribution to a peripheral compartment and two parallel pathways of elimination from the central compartment: a linear and a saturable pathway. Body-surface area (BSA) was significantly correlated to both the volume of distribution of the central compartment and to the maximal elimination capacity. The fits of 500 bootstrap replicates of the data set demonstrated the robustness of the developed population pharmacokinetic model. CONCLUSION: A population pharmacokinetic model has been designed and validated that accurately describes the data of four phase I studies with E7070. Furthermore, it has been demonstrated that BSA-guided dosing for E7070 is important.

Antineoplastic Agents↗

Postoperative nomogram for disease-specific survival after an R0 resection for gastric carcinoma.

PURPOSE: Few published studies have addressed individual patient risk after R0 resection for gastric cancer. We developed and internally validated a nomogram that combines these factors to predict the probability of 5-year gastric cancer-specific survival on the basis of 1,039 patients treated at a single institution. METHODS: Nomogram predictor variables included age, sex, primary site (distal one-third, middle one-third, gastroesophageal junction, and proximal one-third), Lauren histotype (diffuse, intestinal, mixed), number of positive lymph nodes resected, number of negative lymph nodes resected, and depth of invasion. Death as a result of gastric cancer was the predicted end point. The concordance index was used as an accuracy measure, with bootstrapping to correct for optimistic bias. Calibration plots were constructed. RESULTS: Gastric cancer-specific survival at 5 years was 50%. A nomogram was constructed on the basis of a Cox regression model. The bootstrap-corrected concordance index was 0.80. When compared with the predictive ability of American Joint Committee on Cancer stage, the nomogram discrimination was superior (P <.001). Nomogram calibration appeared to be excellent. CONCLUSION: A nomogram was developed to predict 5-year disease-specific survival after R0 resection for gastric cancer. This tool should be useful for patient counseling, follow-up scheduling, and clinical trial eligibility determination.

Adult↗

The addition of interleukin-6 soluble receptor and transforming growth factor beta1 improves a preoperative nomogram for predicting biochemical progression in patients with clinically localized prostate cancer.

PURPOSE: Several preoperative prostate cancer nomograms have been developed that predict risk of progression using pretreatment prostate-specific antigen (PSA) level, clinical stage, and biopsy Gleason grade. We describe the development and performance of a new nomogram. The nomogram adds new markers to the standard clinical predictors that reflect the biologic behavior of prostate cancer: pretreatment plasma levels of interleukin-6 soluble receptor (IL6SR) and transforming growth factor beta1 (TGF-beta1). PATIENTS AND METHODS: Between November 7, 1994 and December 22, 1997, 714 patients with stage cT1c to cT3a prostate cancer and no prior therapy were treated with radical prostatectomy at the Methodist Hospital, Houston TX. Plasma levels of IL6SR and TGF-beta1 were measured in banked preoperative plasma. With these data, a nomogram was developed to predict the probability of PSA progression within 5 years of surgery. The nomogram was validated with bootstrapping to assess its discrimination and calibration performance. RESULTS: In the multivariable Cox model, PSA (P =.004), IL6SR (P <.001), TGF-beta1 (P <.001), primary Gleason grade (P <.002), and secondary Gleason grade (P =.029) were associated with PSA progression, whereas clinical stage (P =.696) was not. The nomogram seemed to be well calibrated and had a bootstrap-corrected area under the receiver operating characteristic curve (ie, concordance index) of 0.83. For comparison, a nomogram that omitted IL6SR and TGF-beta1 achieved a concordance index of only 0.75. CONCLUSION: We found that pretreatment plasma levels of IL6SR and TGF-beta1 improved the ability to predict biochemical progression by a prognostically substantial margin. A nomogram including the pretreatment levels of these molecular markers, along with standard clinical markers, has been developed and internally validated.

Adult↗

Health-related quality of life parameters as prognostic factors in a nonmetastatic breast cancer population: an international multicenter study.

PURPOSE: The purpose of this research was to evaluate whether baseline health-related quality of life (HRQOL) parameters are prognostic factors for survival in locally advanced breast cancer patients. Although the literature highlights the important role of HRQOL parameters in predicting survival in advanced metastatic disease, little evidence exists for earlier stages. PATIENTS AND METHODS: The overall sample consisted of 448 patients randomly assigned to receive cyclophosphamide, epirubicin, and fluorouracil versus epirubicin, cyclophosphamide, and granulocyte colony-stimulating factor. Patients were enrolled in 12 countries. HRQOL baseline scores were assessed using the European Organization for Research and Treatment of Cancer Quality of Life Questionnaire C30. The Cox proportional hazards regression model was used for both univariate and multivariate analyses of survival. In addition, a bootstrap resampling technique was used to assess the stability of the outcomes. Bootstrap results were then applied for model averaging purposes as a means to account for the observed model selection uncertainty. RESULTS: The final multivariate model retained inflammatory breast cancer (T4d) as the only factor predicting overall survival (OS) with a hazard ratio of 1.375 (95% CI, 1.027 to 1.840; P =.03). The presence of inflammatory breast cancer lowers the median survival time from 6.6 to 4.2 years (36% reduction). None of the preselected HRQOL variables were prognostic for OS or disease-free survival, in either the univariate or multivariate analysis. CONCLUSION: Our findings suggest that baseline HRQOL parameters have no prognostic value in a nonmetastatic breast cancer population.

Adult↗

A random sampling approach for robust estimation of tissue-to-plasma ratio from extremely sparse data.

his study was performed to develop a new nonparametric approach for the estimation of robust tissue-to-plasma ratio from extremely sparsely sampled paired data (ie, one sample each from plasma and tissue per subject). Tissue-to-plasma ratio was estimated from paired/unpaired experimental data using independent time points approach, area under the curve (AUC) values calculated with the naïve data averaging approach, and AUC values calculated using sampling based approaches (eg, the pseudoprofile-based bootstrap [PpbB] approach and the random sampling approach [our proposed approach]). The random sampling approach involves the use of a 2-phase algorithm. The convergence of the sampling/resampling approaches was investigated, as well as the robustness of the estimates produced by different approaches. To evaluate the latter, new data sets were generated by introducing outlier(s) into the real data set. One to 2 concentration values were inflated by 10% to 40% from their original values to produce the outliers. Tissue-to-plasma ratios computed using the independent time points approach varied between 0 and 50 across time points. The ratio obtained from AUC values acquired using the naive data averaging approach was not associated with any measure of uncertainty or variability. Calculating the ratio without regard to pairing yielded poorer estimates. The random sampling and pseudoprofile-based bootstrap approaches yielded tissue-to-plasma ratios with uncertainty and variability. However, the random sampling approach, because of the 2-phase nature of its algorithm, yielded more robust estimates and required fewer replications. Therefore, a 2-phase random sampling approach is proposed for the robust estimation of tissue-to-plasma ratio from extremely sparsely sampled data.

Algorithms↗

Using plasma topotecan pharmacokinetics to estimate topotecan exposure in cerebrospinal fluid of children with medulloblastoma.

The purpose of this study was to estimate ventricular cerebrospinal fluid (vCSF) topotecan lactone (TPT) exposures in pediatric medulloblastoma patients from plasma concentration-time data by using a pharmacokinetic (PK) model. We studied children with high-risk medulloblastoma who received pharmacokinetically guided TPT (target plasma area under the concentration-time curve [AUC], 120-160 ng/ml-h) and obtained serial vCSF samples to assess TPT exposure. Population pharmacokinetic parameters were determined by using linear mixed-effects modeling via the two-stage approach. We simulated TPT vCSF exposure duration at plasma TPT AUC values of 120 to 200 ng/ml-h and determined percentages of studies meeting or exceeding the vCSF exposure duration threshold (EDT) of 1 ng/ml for 8 h. We then used bootstrap methods to estimate variability in vCSF EDT. Eighteen PK studies were conducted in six patients (median age, 4.8 years). In these patients, seven of nine studies attaining target plasma TPT AUC achieved the vCSF EDT. Given a plasma TPT AUC of 120 ng/ml-h, the median percentage of results meeting or exceeding EDT was 78% (95% CI, 61%-100%). One patient (four studies) with tumor blockage of CSF flow, which can alter CSF pharmacokinetics, was removed, and the bootstrap analysis was repeated. In this subset, a median 93% (95% CI, 79%-100%) of studies achieved vCSF EDT. Increasing plasma TPT AUC values resulted in increased ability to achieve vCSF EDT. We demonstrated that a plasma PK model could estimate vCSF TPT concentrations. Further, our results indicate that the TPT vCSF EDT can be achieved in more than 80% of studies targeted to a plasma TPT AUC of 120 ng/ml-h.

Adolescent↗