SEARCH · Engineering Papers

Results for “testing statistics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

NASA NTRS2005

Influence of Coronal Abundance Variations

The PI of this project was Jeff Scargle of NASA/Ames. Co-I's were Alma Connors of Eureka Scientific/Wellesley, and myself. Part of the work was subcontracted to Eureka Scientific via SAO, with Vinay Kashyap as PI. This project was originally assigned grant number NCC2-1206, and was later changed to NCC2-1350 for administrative reasons. The goal of the project was to obtain, derive, and develop statistical and data analysis tools that would be of use in the analyses of high-resolution, high-sensitivity data that are becoming available with new instruments. This is envisioned as a cross-disciplinary effort with a number of "collaborators" including some at SA0 (Aneta Siemiginowska, Peter Freeman) and at the Harvard Statistics department (David van Dyk, Rostislav Protassov, Xiao-li Meng, Epaminondas Sourlas, et al). We have developed a new tool to reliably measure the metallicities of thermal plasma. It is unfeasible to obtain high-resolution grating spectra for most stars, and one must make the best possible determination based on lower-resolution, CCD-type spectra. It has been noticed that most analyses of such spectra have resulted in measured metallicities that were significantly lower than when compared with analyses of high- resolution grating data where available (see, e.g., Brickhouse et al., 2000, ApJ 530,387). Such results have led to the proposal of the existence of so-called Metal Abundance Deficient, or "MAD" stars (e.g., Drake, J.J., 1996, Cool Stars 9, ASP Conf.Ser. 109, 203). We however find that much of these analyses may be systematically underestimating the metallicities, and using a newly developed method to correctly treat the low-counts regime at the high-energy tail of the stellar spectra (van Dyk et al. 2001, ApJ 548,224), have found that the metallicities of these stars are generally comparable to their photospheric values. The results were reported at the AAS (Sourlas, Yu, van Dyk, Kashyap, and Drake, 2000, BAAS 196, v32, #54.02), and at the conference on Statistical Challenges in Modem Astronomy (Sourlas, van Dyk, Kashyap, Drake, and Pease, 2003, SCMA 111, Eds. E.D.Feigelson, G.J.Babu, New York:Springer, p489-490). We also described the limitations of one of the most egregiously misused and misapplied statistical tests in astrophysical literature, the F-test for verifying model components (Protassov, van Dyk, Connors, Kashyap, and Siemiginowska, 2002, ApJ, 571,545). Indeed, a search through the ApJ archives turned up 170 papers in the 5 previous years that used the F-test explicitly in some form or the other, and with the vast majority of them not using it correctly! Indeed, looking at just 4 issues of the ApJ in 2001, we found 13 instances of its use, of which nine were demonstrably incorrect. Clearly, it is difficult to understate the importance of this issue. We also worked on speeding up Bayes Blocks and Sparse Bayes Blocks algorithms to make them more tractable for large searches. We also supported staistics students and postdocs in both explicit physics- model-based (spectra with tens of thousands of atomic lines) and "model-free" -- i.e. non-parametric or semi-parametric -- algorithms. Work on using more of the latter is just beginning; while using multi-scale methods for Poisson imaging has come to hition. In fact, "An Image Restoration Technique with Error Estimates", by D. Esch, A. Connors, M. Karovska, and D. van Dyk, was published by ApJ (Esch et a1.2004, ApJ, 610, 1213). The code has been delivered to M. Karovska for CXC; and is available for beta-testing upon request. The other large project we worked on was on the self-consistent modeling of logN-logs curves in the Poisson limit. logN-logs curves are a fundamental tool in the study of source populations, luminosity functions, and cosmological parameters. However, their determination is hampered by statistical effects such as the Eddington bias, incompleteness due to detection efficiency, faint source flux fluctuations, etc. We have develed a new and powerful method using the full Poisson machinery that allows us to model the logN-logs distribution of X-ray sources in a self-consistent manner. Because we properly account for all the above statistical effects, our modeling is valid over the full range of the data, and not just for strong sources, as is normally done. Using a Bayesian approach and modeling the fluxes with known functional forms such as simple or broken power-laws, and conditioning the expected photon counts on the fluxes, the background contamination, effective area, detector vignetting, and detection probability, we can delve deeply into the low counts regime and extend the usefulness of medium sensitivity surveys such as ChAMP by orders of magnitude. The built-in flexibility of the algorithm also allows a simultaneous analysis of multiple datasets. We have applied this analysis to a set a Chandra observations (Sourlas, Kashyap, Zezas, van Dyk, 2004, HEAD #8, #16.32)

Scargle, Jeffrey D.↗

DOE OSTI2022

Multiscale spatial analysis of fracture arrangement and pattern reconstruction using Ripley's K-function

This work presents novel multiscale spatial data analytics using Ripley's K-function, as a measure of spatial interaction, to study one-dimensional arrangement of fractures. Fracture spatial arrangements are classified into clustered, anticlustered, or indistinguishable from random by testing statistical significance of the calculated Ripley's K-function. Characterizations of fracture arrangements are performed as a function of length scale and position. Analysis of the K-function along the study interval identifies where fracture clustering and anticlustering occur. Here, a simulation technique is also introduced here to statistically reconstruct spatial arrangements and to generate fracture realizations that are spatially similar to the fractures observed in the field. With this simulation technique, one can also fill spatial gaps in fracture measurements where data are absent, unreliable, or unused. Synthetic as well as field-measured 1D fracture datasets are used for testing and demonstration. Methods introduced in this work can be readily applied to fracture datasets observed in outcrops, borehole image logs, and cores.

58 GEOSCIENCES↗

DOE OSTI2025

PNNL-Predictive-Phenomics/ProteoMeter

ProteoMeter is a Python package that assists in the statistical analysis of global proteomics, protein post-translation modification (PTM), and limited proteolysis (LiP) data. It contains batch correction, normalization, and statistical testing methods, as well as functions that "roll up" peptide-level data to the single-site level. It has a robust user configuration system, allowing it to flexibly integrate different types of experiment designs. For basic usage, a simple configuration file provides the essential functionality. Advanced users have access to the entire statistical pipeline for fine-tuning analyses. Processed data is easily exported to many common spreadsheet and data-frame formats.

Rozum, Jordan [Pacific Northwest National Lab]↗

NASA NTRS1996

The Effects of Using Space to Teach Standard Elementary School Curriculum

This brief report and recommendation for further research brings to a formal close this effort, the original purpose of which is described in detail in The effects of using space to teach standard elementary school curriculum, Volume 1, included here as the Appendix. Volume 1 describes the project as a 3-year research program to determine the effectiveness of using space to teach. The research design is quasi experimental using standardized test data on students from Aldrin Elementary School and a District-identified 'control' school, which shall be referred to as 'School B.' Students now in fourth through sixth grades will be compared now (after one year at Aldrin) and tracked at least until the present sixth graders are through the eighth grade. Appropriate statistical tests will be applied to standardized test scores to see if Aldrin students are 'better' than School B students in areas such as: Overall academic performance; Performance in math/science; and Enrollments in math/science in middle school.

Ewell, Robert N.↗

NASA NTRS1986

The effects of expressivity and flight task on cockpit communication and resource management

The results of an investigation to develop a methodology for evaluating crew communication behavior on the flight deck and a flight simulator experiment to test the effects of crew member expressivity, as measured by the Personal Attributes Questionnarie, and flight task on crew communication and flight performance are discussed. A methodology for coding and assessing flight crew communication behavior as well as a model for predicting that behavior is advanced. Although not enough crews were found to provide valid statistical tests, the results of the study tend to indicate that crews in which the captain has high expressivity perform better than those whose captain is low in expressivity. There appears to be a strong interaction between captains and first officers along the level of command dimension of communication. The PAQ appears to identify those pilots who offer disagreements and inititate new subjects for discussion.

Jensen, R. S.↗

DOE OSTI2021

Spectral characterization of flash and high flux x-ray radiographic sources with a magnetic Compton spectrometer

In this article, we present a new analysis method applied to revitalize permanent magnet Compton spectrometers used to measure photon energy spectra in the MeV range. The inversion of the measured electron distribution to determine the original photon distribution is achieved via a method of consistent coupled radiation transport and magnetic field mapping of the input photon spectra to the measured electron distribution. The method of linear least squares was used to perform the unfolding of the electron distribution to the initial photon spectra, without any assumptions made regarding the electron distribution. We present an application of this method to data from a nominal 19.4 MeV flash radiographic source (the first axis of the Dual Axis Radiographic Hydro-Test Facility) capable of generating 500 R @ 1 m in ~60 ns and a medical therapy source (a Scanditronix M22, Microtron) capable of variable energies with nominal endpoints of 6, 10, 15, and 20 MeV and an output of ~1000–2000 R/min @ 1 m. The results provide agreement between the modeled and unfolded experimentally measured photon spectra as quantified by statistical tests, from 1.5 to 20 MeV. Experimental results are presented as well as a discussion of the novel MCNP6-based simulations and methods for reconstruction of the spectra.

47 OTHER INSTRUMENTATION↗

DOE OSTI2024

Feasibility study for test rig assessments of fish passage conditions in a Kaplan turbine

The assessment of fish passage conditions in hydroelectric turbines consists of identifying and quantifying physical magnitudes leading to increased risks of injury of fish passing through turbines in operation. Such assessments are usually carried out either with the use of computer-based methods during design or with field testing of live fish and sensors passing through prototypes. A method in between consists of test rig experimentation, which is critical for testing fish-focused design concepts and offers the opportunity for implementing the most effective design measures for improved fish survivability. However, fish-related assessments in test rigs are not sufficiently documented for industrial applications. This work presents the main findings of an experimental campaign to quantify fish-related hydraulic magnitudes in a physical model of a Kaplan turbine in a commercial test rig. Two operating conditions were tested by releasing miniaturized autonomous sensor devices (termed Sensor Fish Mini) at the turbine intake flow, passing them through the runner in motion and recovering them at the draft tube exit. During passage, time series of acceleration, absolute pressure and rotational velocity were recorded. The recordings were then interpreted to determine the magnitude and likely location of hydraulic stressors hazardous to fish. The statistical tests on the reported measurements indicated that low pressure, collisions on the runner and rotations in the draft tube were not different between the two tested operating points. On the other hand, pressure drop and collision rates on the distributor differed considerably as a function of net head. The outcomes of this investigation showed that test rig evaluations of fish-related properties with Sensor Fish Mini can contribute to an evidence-based development of turbine geometries designed for providing safer passage conditions. Future work will investigate the scaling of test rig measurements to hydraulically equivalent magnitudes in the prototype and their biological consequences.

13 HYDRO ENERGY↗

DOE OSTI2020

The look-elsewhere effect from a unified Bayesian and frequentist perspective

When searching over a large parameter space for anomalies such as events, peaks, objects, or particles, there is a large probability that spurious signals with seemingly high signi ficance will be found. This is known as the look-elsewhere effect and is prevalent throughout cosmology, (astro)particle physics, and beyond. To avoid making false claims of detection, one must account for this effect when assigning the statistical significance of an anomaly. This is typically accomplished by considering the trials factor, which is generally computed numerically via potentially expensive simulations. In this paper we develop a continuous generalization of the Bonferroni and Sidak corrections by applying the Laplace approximation to evaluate the Bayes factor, and in turn relating the trials factor to the prior-to-posterior volume ratio. Here, we use this to define a test statistic whose frequentist properties have a simple interpretation in terms of the global p-value, or statistical significance. We apply this method to various physics-based examples and show it to work well for the full range of p-values, i.e. in both the asymptotic and non-asymptotic regimes. We also show that this method naturally accounts for other model complexities such as additional degrees of freedom, generalizing Wilks' theorem. This provides a fast way to quantify statistical significance in light of the look-elsewhere effect, without resorting to expensive simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2004

A Meteorological Model's Dependence on Radiation Update Frequency

Numerical weather models are used to simulate circulations in the atmosphere including clouds and precipitation by applying a set of mathematical equations over a three-dimensional grid. The grid is composed of discrete points at which the meteorological variables are defined. As computing power continues to rise these models are being used at finer grid spacing, but they must still cover a wide range of scales. Some of the physics that must be accounted for in the model cannot be explicitly resolved, and their effects, therefore, must be estimated or "parameterized". Some of these parameterizations are computationally expensive. To alleviate the problem, they are not always updated at the time resolution of the model with the assumption being that the impact will be small. In this study, a coupled land-atmosphere model is used to assess the impact of less frequent updates of the computationally expensive radiation physics for a case on June 6, 2002, that occurred during a field experiment over the central plains known as International H20 Project (IHOP). The model was tested using both the original conditions, which were dry, and with modified conditions wherein moisture was added to the lower part of the atmosphere to produce clouds and precipitation (i.e., a wet case). For each of the conditions (i.e., dry and wet), four set of experiments were conducted wherein the model was run for a period of 24 hours and the radiation fields (including both incoming solar and outgoing longwave) were updated every 1, 3, 10, and 100 time steps. Statistical tests indicated that average quantities of surface variables for both the dry and wet cases were the same for the various update frequencies. However, spatially the results could be quite different especially in the wet case after it began to rain. The near-surface wind field was found to be different most of the time even for the dry case. In the wet case, rain intensities and average vertical profiles of heating associated with cloudy areas were found to differ for the various radiation update frequencies. The latter implies that the mean state of the model could be different as a result of not updating the radiation fields every time step and has important implications for longer term climate studies

Eastman, Joseph L.↗

DOE OSTI2020

Detectability of standing accretion shock instabilities activity in supernova neutrino signals

We introduce a novel methodology for establishing the presence of standing accretion shock instabilities (SASI) in the dynamics of a core collapse supernova from the observed neutrino event rate at water- or ice- based neutrino detectors. The methodology uses a likelihood ratio in the frequency domain as a test- statistics; it is also employed to assess the potential to estimate the frequency and the amplitude of the SASI modulations of the neutrino signal. The parameter estimation errors are consistent with the minimum possible errors as evaluated from the inverse of the Fisher information matrix, and close to the theoretical minimum for the SASI amplitude. Using results from a core-collapse simulation of a 15 solar-mass star by Kuroda et al. [Astrophys. J. 851, 62 (2017).] as a test bed for the method, we find that SASI can be identified with high confidence for a distance to the supernova of up to ~6 kpc for IceCube and up to ~3 kpc for a 0.4 Mt mass water Cherenkov detector. This methodology will aid the investigation of a future galactic supernova.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1977

Testing the accuracy of remote sensing land use maps

Some of the main aspects that need to be considered in a remote sensing sampling design are: (1) the frequency that any one land use type (on the ground) is erroneously attributed to another class by the interpreter; (2) the frequency that the wrong land use (as observed on the ground) is erroneously included in any one class by the remote sensing interpreter; (3) the proportion of all land (as determined in the field) that is mistakenly attributed by the interpreter; and (4) the determination of whether the mistakes are random (so that the overall proportions are approximately correct) or subject to a persistent bias. A sampling and statistical testing procedure is presented which allows an approximate answer to each of these aspects. The concept developed and described incorporates the probability of making incorrect interpretations at particular prescribed accuracy levels, for a certain number of errors, for a particular sample size. It is considered that this approach offers a meaningful explanation of the interpretation accuracy level of an entire remote sensing land use survey.

Vangenderen, J. L.↗

DOE OSTI2020

Statistical Modeling and Analysis of k -Layer Coverage of Two-Dimensional Materials in Inkjet Printing Processes

Two-dimensional layered materials/flakes, also known as crystalline atom-thick layer nanosheets, have recently been receiving great attention in electronics fabrication due to their unique and intriguing properties. The k-layer coverage area (i.e., the area covered by k number of overlapping layers) of the printed flake pattern significantly impacts on the properties of the printed electronics. In this work, we constructed a statistical model to describe the k-layer coverage of randomly distributed two-dimensional materials. A series of results are obtained to provide not only the expectation but also the variance of the coverage area. The boundary effects on the random flakes coverage are also studied. In addition, an approximated statistical testing approach is also developed in this work to detect abnormal coverage patterns. Furthermore, the case studies based on simulated data and real flakes images obtained from the inkjet printing process demonstrate the accuracy and effectiveness of the proposed model and analysis methods.

36 MATERIALS SCIENCE↗

DOE OSTI2018

Regression for Radioactive Source Detection

Detection of radioactive sources is an important capability that led to the deployment of networks for detecting and localizing low-level, hazardous radiation sources. It is generally expected that such networks outperform the individual detectors by intelligently fusing information from several dispersed sensors. In this paper, we develop a network detection method that uses the slope of a linear regression fit as a test statistic for detecting a point radioactive source within a field of detectors. In our regression model, we compute a least-squared linear fit between the average radiation counts at the detectors and the inverse-squared distances of known detector locations to an estimated source location. We show that the slope of this regression fit is an estimate of the source intensity and can be used as a threshold for source detection purposes. We compare performance of our proposed detection method with that of a fusion-based Sequential Probability Ratio Test (SPRT) method. For performance analyses, two datasets from the Domestic Nuclear Detection Ofces Intelligence Radiation Sensors Systems (IRSS) outdoor tests are used. Each of these tests consists of several runs of a single radioactive source moving in and out of a detector network. We present receiver operating characteristic (ROC) curves and optimal threshold values for the performance each detection method, and determine that our detection method using linear regression fit has slightly better overall performance than SPRT method.

Cordone, G.↗

NASA NTRS1984

The rotation of small asteroids

The Binzel and Mulholland (1983) sample of photoelectrically determined rotational parameters for 17 main belt asteroids with diameters less than 30 km are compared with previous observations of asteroids of that size range. Rigorous statistical tests are applied to investigate bias effects and quantify results on asteroid rotation. The samples are described and compared for rotational frequency, rotational amplitude, frequency distribution, and diameter and frequency dependence. It is concluded that the observed rotational frequency distribution can be acceptably fit by two Maxwellian distributions, which is consistent with the hypothesis that there are separate populations of slow and fast rotating asteroids. The frequency distributions of main belt asteroids less than 15 km in diameter and earth and Mars crossers do not differ significantly, but the larger mean lightcurve amplitude of these crossers is statistically significant. No significant diameter dependence on rotational frequency is found among the sample asteroids.

Binzel, R. P.↗

NASA NTRS1997

Impact Crises, Mass Extinctions, and Galactic Dynamics: A Unified Theory

A general hypothesis linking mass extinctions of life with impacts of large asteroids and comets is based on astronomical data, impact dynamics, and geological information. The waiting times of large-body impacts on the Earth, derived from the flux of Earth-crossing asteroids and comets, and the estimated size of impacts capable of causing large-scale environmental disasters predict that impacts of objects (sup 3)5 km in diameter ((sup 3)10(exp 7) Mt TNT equivalent) could be sufficient to explain the record of about 25 extinction pulses in the last 540 m.y., with the five recorded major mass extinctions related to the impacts of the largest objects of (sup 3)10 km in diameter ( (sup 3)10(exp 8) Mt events). Smaller impacts (about 10(exp 6)-10(exp 7) Mt), with significant regional and even global environmental effects, could be responsible for the lesser boundaries in the geologic record. Tests of the "kill curve" relationship for impact-induced extinctions based on new data on extinction intensities and several well-dated large impact craters suggest that major mass extinctions require large impacts, and that a step in the kill curve may exist at impacts that produce craters of -100 km diameter, with smaller impacts capable of only relatively weak extinction pulses. Single impact craters < about 60 km in diameter should not be associated with global extinction pulses detectable in the Sepkoski database (although they may explain stage and zone boundaries marked by lesser faunal turnover), but multiple impacts in that size range may produce significant stepped extinction pulses. Statistical tests of the last occurrences of species at mass-extinction boundaries are generally consistent with predictions for abrupt or stepped extinctions, and several boundaries are known to show "catastrophic" signatures of environmental disasters and biomass crash, impoverished postextinction fauna and flora dominated by stress-tolerant and opportunistic species, and gradual ecological recovery and radiation of new taxa. Isotopic and other geochemical signatures are also generally consistent with the expected after-effects of catastrophic impacts. Seven of the recognized extinction pulses are associated with concurrent (in some cases multiple) stratigraphic impact markers (e.g., layers with high Ir, shocked minerals, microtektites), and/or large, dated impact craters. Other less-well-studied crisis intervals show elevated Ir, still well below that of the K/T spike, which might be explained by low-Ir impactors, ejecta blowoff, or the sedimentary reworking and dilution of impact signatures. The best explanation for a possible periodic component of about 30 m.y. in mass extinctions and clusters of impacts is the modulation of the comet flux associated with the solar system's periodic passage through the plane of the Milky Way Galaxy. The quantitative agreement among paleontological, geological, and astronomical data suggests an important underlying unification of the processes involved.

Rampino, M.R.↗

NASA NTRSDate not supplied

Using Machine Learning to Infer Material Properties of Debris Fragments from X-ray Images in the DebriSat Project

The DebriSat project is a collaboration effort with the NASA Orbital Debris Program Office, the U.S. Space Force Space Systems Command Center, The Aerospace Corporation, and the University of Florida. To date, over 200,000 fragments from this ground-based, hypervelocity impact experiment have been collected, and processing is underway to determine their physical characteristics, such as material, shape, color, characteristic length, and average cross-sectional area. The x-ray process is primarily used to identify the location of the fragments and estimated size for extraction, so that these physical characteristics can be assessed. This paper proposes a machine learning-based approach to characterize materials from x-ray images of debris fragments embedded in soft-catch foam used in the DebriSat project. The novel methodology discussed in this paper will highlight the use of x-ray imagery data to characterize these fragments without extraction or a human-in-the-loop. Both supervised and unsupervised machine learning techniques are utilized with this approach to infer the physical parameters of the fragments embedded in the soft-catch foam panels used in the impact experiment based on x-ray images of the foam panels. Additionally, 3D reconstructions of the extracted fragments are created with images taken from two different angles using the structure from motion (SfM) method. The characteristic lengths and shape from the 3D reconstruction, alongside the physical characteristics of the debris, are used in the inference of the material type. To develop and test the approach, a dataset of x-ray images of debris fragments of varying sizes and materials is collected. Supervised learning methods such as convolutional neural networks (CNNs), support vector machines (SVM), decision trees, and random forest classifiers are used due to the high-dimensional feature spaces of the debris and nonlinear decision boundaries for material categorization. Given the limited pre-labeled data of embedded debris materials smaller than 10 mm, unsupervised machine learning techniques such as clustering algorithms and autoencoders are used, in addition to supervised learning methods. The clustering algorithms group similar fragments together based on their physical properties, and autoencoders reduce the dimensionality of the x ray images and extract relevant features. The performance of the proposed approach's is analyzed using a range of statistical methods, including confusion matrices, receiver operating characteristic curves, and precision-recall curves. The results are compared with those obtained using a baseline approach that relies on manual identification and classification of debris fragments. To evaluate the effectiveness of different machine learning methods, statistical tests such as t-tests, ANOVA, and cross-validation are performed, comparing the performance of CNNs, SVMs, clustering algorithms, and autoencoders. Additional analysis needs to be conducted to identify any sources of bias or variability that may affect the results, such as variations in imaging conditions or fragmentation patterns. Other topics explored are limitations, refinements, and the potential use of semi-supervised learning techniques, such as self-training to label unlabeled datasets and co-training using x-ray images taken from two different angles as two different models.

Saik Anam Siam↗

NASA NTRS2023

Using Machine Learning to Infer Material Properties of Debris Fragments from X-ray Images in the DebriSat Project

Saik Anam Siam↗

DOE OSTI2019

Smoking, Alcohol, and Biliary Tract Cancer Risk: A Pooling Project of 26 Prospective Studies

Abstract Background Tobacco and alcohol are well-established risk factors for numerous cancers, yet their relationship to biliary tract cancers remains unclear. Methods We pooled data from 26 prospective studies to evaluate associations of cigarette smoking and alcohol consumption with biliary tract cancer risk. Study-specific hazard ratios (HRs) and 95% confidence intervals (CIs) for associations with smoking and alcohol consumption were calculated. Random-effects meta-analysis produced summary estimates. All statistical tests were two-sided. Results Over a period of 38 369 156 person-years of follow-up, 1391 gallbladder, 758 intrahepatic bile duct, 1208 extrahepatic bile duct, and 623 ampulla of Vater cancer cases were identified. Ever, former, and current smoking were associated with increased extrahepatic bile duct and ampulla of Vater cancers risk (eg, current vs never smokers HR = 1.69, 95% CI = 1.34 to 2.13 and 2.22, 95% CI = 1.69 to 2.92, respectively), with dose-response effects for smoking pack-years, duration, and intensity (all Ptrend < .01). Current smoking and smoking intensity were also associated with intrahepatic bile duct cancer (eg, >40 cigarettes per day vs never smokers HR = 2.15, 95 % CI = 1.15 to 4.00; Ptrend = .001). No convincing association was observed between smoking and gallbladder cancer. Alcohol consumption was only associated with intrahepatic bile duct cancer, with increased risk for individuals consuming five or more vs zero drinks per day (HR = 2.35, 95%CI = 1.46 to 3.78; Ptrend = .04). There was evidence of statistical heterogeneity among several cancer sites, particularly between gallbladder cancer and the other biliary tract cancers. Conclusions Smoking appears to increase the risk of developing all biliary tract cancers except gallbladder cancer. Alcohol may increase the risk of intrahepatic bile duct cancer. Findings highlight etiologic heterogeneity across the biliary tract.

McGee, Emma E.↗