Search NASA⌕ Search

SEARCH · Search NASA

Results for “sampling 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.

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

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Mechanical Performance and Mesostructure Analysis of Proton-Irradiated Fused Filament Fabrication Acrylonitrile Butadiene Styrene Material

The use of fused filament fabrication (FFF) acrylonitrile butadiene styrene (ABS) and other thermoplastics in radiation environments is beginning to be studied as possible replacements for traditionally manufactured parts and tools. Interlayer adhesion within the mesostructure has been shown in published literature to be an integral component in the strength of FFF ABS. Research of irradiated 3D printed polymers has primarily focused on the influence of gamma irradiation on the mechanical properties of FFF ABS samples, without evaluating its impact on the mesostructure of the samples. The purpose of this paper is to understand the mechanical damage caused by proton radiation on FFF ABS samples through the evaluation of the mesostructure of the ABS samples. To achieve this objective, proton radiation at 40 MeV was applied to FFF ABS samples at radiation doses up to 1.0 MGy at a dose rate of 1 MGy/hr. Following the irradiation, tensile testing was performed on the samples. Here, the fractured surfaces of the tested samples were subsequently observed using a scanning electron microscope. The mechanical test results show an increase in the ultimate tensile strength (UTS) and a decrease in the ductility of the irradiated samples. Statistical analysis on the results shows that there is a statistically significant difference in the UTS and the ductility of unirradiated and 1.0 MGy irradiated samples and between irradiated samples and 1.0 MGy samples. In addition, the percentage elongation at break has statistically significant differences in the means between irradiated samples and 1.0 MGy samples. The UTS has statistically significant differences in the means between 0.1 MGy and 1.0 MGy. The difference in the means between the unirradiated and the 1.0 MGy samples is the most significant for both the UTS and elongation at break. The scanning electron microscopy (SEM) results indicated that interlayer adhesion improved as a function of radiation dose corresponding to the increase in tensile strength. The SEM results also showed that crazing and plastic deformation were reduced; aligning with the loss in ductility observed in the tensile tests results. The proton radiation is causing these mechanical and physical changes through two mechanisms: (1) dose rate effects on ionizing radiation-induced oxidative degradation; and (2) radiation heating effects of high energy (>1 MeV) charged particles.

electron microscopy↗
DOE OSTI2025

First observation of antiproton annihilation at rest on argon in the LArIAT experiment

We report the first observation and measurement of antiproton annihilation at rest on argon track and shower multiplicities and particle identification conducted with the LArIAT experiment. Stopping antiprotons from the Fermilab Test Beam Facility’s charged particle test beam are identified using beamline instrumentation and LArIAT’s liquid argon time projection chamber (LArTPC). The charged particle multiplicity from the annihilation vertex is manually evaluated via hand scanning, yielding a mean of 3.2 ± 0.4 tracks and a standard deviation of 1.3 tracks, consistent with a semiautomated reconstruction resulting in 2.8 ± 0.4 tracks and a standard deviation of 1.2 tracks. Both methods are consistent with Monte Carlo simulations within statistical uncertainty. The shower multiplicities and particle identification for outgoing tracks are also consistent with eant4 model predictions. These results, obtained from a low-statistics sample, provide a foundation for higher-statistics studies in larger LArTPCs, which could refine modeling of intranuclear annihilation on argon and inform scenarios such as neutron-antineutron oscillations. Published by the American Physical Society 2025

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2008

Space Debris Symposium (A6.) Measurements and Space Surveillance (1.): Measurements of the Small Particle Debris Cloud from the 11 January, 2007 Chinese Anti-satellite Test

On January 11, 2007, the Chinese military conducted a test of an anti-satellite (ASAT) system, destroying their own Fengyun-1C spacecraft with an interceptor missile. The resulting hypervelocity collision created an unprecedented number of tracked debris - more than 2500 objects. These objects represent only those large enough for the US Space Surveillance Network (SSN) to track - typically objects larger than about 5-10 cm in diameter. There are expected to be even more debris objects at sizes too small to be seen and tracked by the SSN. Because of the altitude of the target satellite (865 x 845 km orbit), many of the debris are expected to have long orbital lifetimes and contribute to the orbital debris environment for decades to come. In the days and weeks following the ASAT test, NASA was able to use Lincoln Laboratory s Haystack radar on several occasions to observe portions of the ASAT debris cloud. Haystack has the capability of detecting objects down to less than one centimeter in diameter, and a large number of centimeter-sized particles corresponding to the ASAT cloud were clearly seen in the data. While Haystack cannot track these objects, the statistical sampling procedures NASA uses can give an accurate statistical picture of the characteristics of the debris from a breakup event. For years computer models based on data from ground hypervelocity collision tests (e.g., the SOCIT test) and orbital collision experiments (e.g., the P-78 and Delta-180 on-orbit collisions) have been used to predict the extent and characteristics of such hypervelocity collision debris clouds, but until now there have not been good ways to verify these models in the centimeter size regime. It is believed that unplanned collisions of objects in space similar to ASAT tests will drive the long-term future evolution of the debris environment in near-Earth space. Therefore, the Chinese ASAT test provides an excellent opportunity to test the models used to predict the future debris environment. For this study, Haystack detection events are compared to model predictions to test the model assumptions, including debris size distribution, velocity distribution, and assumptions about momentum transfer between the target and interceptor. In this paper we will present the results of these and other measurements on the size and extent of collisional breakup debris clouds.

Matney, Mark J.↗
NASA NTRSDate not supplied

A New Statistical Estimate of the Radar Coverage of the Low Earth Orbit Debris Environment

For over three decades, the NASA Orbital Debris Program Office (ODPO) has used the Goldstone Orbital Debris Radar, Haystack Ultrawideband Satellite Imaging Radar (HUSIR), and Haystack Auxiliary (HAX) radar assets to collect data on the low Earth orbit (LEO) debris environment. Each radar, with its unique beamwidth, altitude and inclination coverage, and limiting size threshold, operates in a beam park mode to statistically sample the orbital debris population in LEO. Provided that these assets are shared with other users, the orbital debris data collection is not continuous; rather, intermittent data collects are acquired and sent to the NASA ODPO. To understand the sampling process conducted by each radar over time and any related observational biases, a Statistical Confirmation of Radar Uniformity or Bias (SCRUB) code has been implemented to model the coverage of these assets for informing future operations, as well as usage of the data collected from these ground-based sensors. For this analysis, Right Ascension of the Ascending Node (RAAN) is used as a metric to measure statistical coverage. A complete survey of the LEO environment is understood to be measurements that sufficiently sample all values of RAAN for each altitude-inclination pair visible from the radar asset. Regions of incompleteness or statistical bias can help inform future observation campaigns. The SCRUB tool evaluates the coverage of the LEO environment, not by examining individual objects that may pass through a sensor’s field of view (FOV), but by determining which orbit planes pass through the FOV. Once a pointing geometry for a radar site, observation time, and range extent are configured by the user, SCRUB computes the inclinations and altitudes that are visible by the sensor. For each inclination-altitude pair, there is a distinct pair of possible RAAN values, corresponding to the ascending and descending orbit passing through that point in the sensor cone. Repeating this process for all points in the beam, and for multiple time periods during an observation window, creates a matrix of all inclination-altitude-RAAN combinations that are visible during a sensor run. This process can then be repeated for all observations within a year (for an annual survey), propagating all the RAAN values to a common epoch, typically the start of a calendar year, and combining the observations to assemble a full estimate of the RAAN coverage of the LEO environment. This RAAN coverage can then be analyzed for uniformity of sampling, within a certain inclination-altitude pair, or between larger regions of the space environment. This paper provides a general overview of the radar assets utilized by ODPO, typical analysis data products assuming circular orbits, and a discussion of the algorithms that feed between modeling and measurement operations. Following the description of the algorithm, estimates of coverage using HUSIR radar data collected from multiple years are developed and compared.

Chris L Ostrom↗
DOE OSTI2021

S0 Uncrewed Aircraft System Measurement Characterization Field Campaign Report

The lower atmosphere contributes directly to the modulation of weather and climate. Understanding the importance of this part of the atmosphere, scientists have worked to improve the representation of the atmospheric boundary layer in numerical prediction tools. Such work depends upon information from observing systems, including remote sensors, weather balloons, and research aircraft. Recent years have seen significant advances in uncrewed aircraft systems (UAS) for atmospheric research. Those efforts have provided new perspectives on atmospheric and surface conditions, particularly from smaller UAS (sUAS) platforms. This field campaign included collection of data using several sUAS, and comparing the sUAS data to those from different observing facilities at the U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) user facility’s Southern Great Plains (SGP) observatory in Oklahoma. Over the course of one week, the Black Swift Technologies and University of Nebraska teams conducted flights for the sole purpose of system evaluation and intercomparison. In addition, these teams conducted joint flight operations at the Marshall, Oklahoma Mesonet site alongside teams from the University of Colorado, University of Oklahoma, and Oklahoma State University to collect additional data for side-by-side comparison. Platforms operated by Black Swift Technologies and the University of Nebraska under this project included one fixed-wing and two rotary-wing platforms, as well as an instrumented surface vehicle. These include the Black Swift Technologies S0 UAS, the University of Nebraska M600 UAS, a University of Nebraska-operated Meteomatics Meteodrone UAS, and the University of Nebraska CoMeT (Combined Mesonet and Tracker) vehicle. Additional details on all of these platforms can be found in de Boer et al. (in prep). In total, 95 sUAS flights were conducted for a combined total of 18.6 flight hours. The S0 was operated both to the south and north of the SGP 60 m tower, while the M600 and Meteodrone were operated directly to the east of the tower. Weather conditions were generally good, with moderate winds. Flights were conducted both to follow the extra radiosondes that were launched as part of this campaign, as well as to conduct extended statistical sampling at tower instrument heights. The later flights also provided data for statistical platform intercomparison.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Foundational Principles for Large-Scale Inference: Illustrations Through Correlation Mining

When can reliable inference be drawn in the ‘‘Big Data’’ context? This article presents a framework for answering this fundamental question in the context of correlation mining, with implications for general large-scale inference. In large-scale data applications like genomics, connectomics, and eco-informatics, the data set is often variable rich but sample starved: a regime where the number n of acquired samples (statistical replicates) is far fewer than the number p of observed variables (genes, neurons, voxels, or chemical constituents). Much of recent work has focused on understanding the computational complexity of proposed methods for ‘‘Big Data.’’ Sample complexity, however, has received relatively less attention, especially in the setting when the sample size n is fixed, and the dimension p grows without bound. To address this gap, we develop a unified statistical framework that explicitly quantifies the sample complexity of various inferential tasks. Sampling regimes can be divided into several categories: 1) the classical asymptotic regime where the variable dimension is fixed and the sample size goes to infinity; 2) the mixed asymptotic regime where both variable dimension and sample size go to infinity at comparable rates; and 3) the purely high-dimensional asymptotic regime where the variable dimension goes to infinity and the sample size is fixed. Each regime has its niche but only the latter regime applies to exa-scale data dimension. We illustrate this high-dimensional framework for the problem of correlation mining, where it is the matrix of pairwise and partial correlations among the variables that are of interest. Correlation mining arises in numerous applications and subsumes the regression context as a special case. We demonstrate various regimes of correlation mining based on the unifying perspective of high-dimensional learning rates and sample complexity for different structured covariance models and different inference tasks.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

(U) A General-Purpose Code for Correlated Sampling Using Batch Statistics with MCNP6 for Fixed-Source Problems

Correlated sampling can be used to reduce the uncertainty of a difference of tallies by taking advantage of the negative covariance term in the sandwich formula. Booth first showed how correlated sampling can be applied with batch statistics using MCNP’s tally fluctuation chart (TFC) to reduce the uncertainty of a difference of tallies in fixed-source problems. Booth presented a problem in which a 1273% uncertainty in a difference was reduced to 8% by accounting for correlations. Researchers He and Su recently studied correlated sampling using the TFC in MCNP version 5. They determined that the code did not print enough digits in the TFC tally means for accurate batch statistics in some cases. After modifying the source code, they concluded that “correlated sampling can yield a standard deviation of about one magnitude smaller than that predicted by the direct, un-correlated simulation when the changes in system response are small (say about 1%), which is equivalent to saving in CPU time by a factor of 100. Such saving [sic] becomes less significant as the change in system response becomes larger.” He and Su provided the formulas needed to apply batch statistics to compute the correlated uncertainty of a difference of tallies. In this report, we follow up on their work by providing the formulas needed to apply batch statistics to compute the correlated uncertainty of a ratio of tallies and of a difference of two tallies divided by a third tally. We extend these formulas to differences and ratios of ratios. These formulas are applied to reduce the uncertainty associated with calculating a relative sensitivity. He and Su did not investigate the accuracy of their correlated sampling uncertainty estimates. We use their test problems and evaluate the accuracy of the uncertainty estimates by comparing with results obtained from random sampling, and, in simple cases, with theoretical values of the “exact” uncertainties. We find that the uncertainties obtained from batch statistics are accurate as long as at least 100 batches are used. We present a new computer code, COSUBS (COrrelated Sampling Using Batch Statistics), that reads MCNP6 TFCs and applies correlated sampling using batch statistics for the tally combinations that the user specifies. COSUBS is a very general tool that compares all TFCs for a base case and one or two perturbed cases. It computes uncertainties for ratios if given only a base case. This report is organized as follows. The equations to apply batch statistics to the difference of random tallies are reviewed in Sec. II. Section III presents the equations for applying batch statistics to a ratio of random tallies; this is useful for computing relative sensitivities using a one-sided finite difference and the relative sensitivity using the differential operator method. Section IV presents the equations for applying batch statistics to a difference of two random tallies divided by a third; this is useful for computing a relative sensitivities using a central difference. Section V presents the equations for applying batch statistics to a difference of two ratios with four random tallies. Section VI presents the equations for applying batch statistics to a one-sided finite difference estimate of the relative sensitivity of a ratio (this uses four random tallies). Section VII presents the equations for applying batch statistics to a central difference estimate of the relative sensitivity of a ratio (this uses six random tallies). Section VIII presents the equations for applying batch statistics to a sum of random tallies. Section IX discusses how to apply batch statistics using MCNP6. Section X presents COSUBS, describing its command-line options and logic. Sections XI through XVI present numerical results for various test problems. Section XVII is a summary and conclusions. Appendix A derives the theoretical Monte Carlo tally variance given certain assumptions; these variances are used to verify the batch statistics for some of the problems. Appendix B lists the MCNP6 input for the unperturbed example problem. Appendix C presents modifications made to MCNP6.3 to support this work.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1977

Optimal spatial sampling techniques for ground truth data in microwave remote sensing of soil moisture

The paper examines optimal sampling techniques for obtaining accurate spatial averages of soil moisture, at various depths and for cell sizes in the range 2.5-40 acres, with a minimum number of samples. Both simple random sampling and stratified sampling procedures are used to reach a set of recommended sample sizes for each depth and for each cell size. Major conclusions from statistical sampling test results are that (1) the number of samples required decreases with increasing depth; (2) when the total number of samples cannot be prespecified or the moisture in only one single layer is of interest, then a simple random sample procedure should be used which is based on the observed mean and SD for data from a single field; (3) when the total number of samples can be prespecified and the objective is to measure the soil moisture profile with depth, then stratified random sampling based on optimal allocation should be used; and (4) decreasing the sensor resolution cell size leads to fairly large decreases in samples sizes with stratified sampling procedures, whereas only a moderate decrease is obtained in simple random sampling procedures.

Rao, R. G. S.↗
NASA NTRS1977

Remote sensing-aided systems for snow qualification, evapotranspiration estimation, and their application in hydrologic models

The design of general remote sensing-aided methodologies was studied to provide the estimates of several important inputs to water yield forecast models. These input parameters are snow area extent, snow water content, and evapotranspiration. The study area is Feather River Watershed (780,000 hectares), Northern California. The general approach involved a stepwise sequence of identification of the required information, sample design, measurement/estimation, and evaluation of results. All the relevent and available information types needed in the estimation process are being defined. These include Landsat, meteorological satellite, and aircraft imagery, topographic and geologic data, ground truth data, and climatic data from ground stations. A cost-effective multistage sampling approach was employed in quantification of all the required parameters. The physical and statistical models for both snow quantification and evapotranspiration estimation was developed. These models use the information obtained by aerial and ground data through appropriate statistical sampling design.

Korram, S.↗
DOE OSTI2016

THE ALFALFA H α SURVEY. I. PROJECT DESCRIPTION AND THE LOCAL STAR FORMATION RATE DENSITY FROM THE FALL SAMPLE

The ALFALFA H α survey utilizes a large sample of H i-selected galaxies from the ALFALFA survey to study star formation (SF) in the local universe. ALFALFA H α contains 1555 galaxies with distances between ∼20 and ∼100 Mpc. We have obtained continuum-subtracted narrowband H α images and broadband R images for each galaxy, creating one of the largest homogeneous sets of H α images ever assembled. Our procedures were designed to minimize the uncertainties related to the calculation of the local SF rate density (SFRD). The galaxy sample we constructed is as close to volume-limited as possible, is a robust statistical sample, and spans a wide range of galaxy environments. In this paper, we discuss the properties of our Fall sample of 565 galaxies, our procedure for deriving individual galaxy SF rates, and our method for calculating the local SFRD. We present a preliminary value of log(SFRD[ M {sub ⊙} yr{sup −1} Mpc{sup −3}]) = −1.747 ± 0.018 (random) ±0.05 (systematic) based on the 565 galaxies in our Fall sub-sample. Compared to the weighted average of SFRD values around z ≈ 2, our local value indicates a drop in the global SFRD of a factor of 10.2 over that lookback time.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Parametric description of intermittent probability distribution functions in solar wind and magnetohydrodynamic turbulence

ABSTRACT In this work, we find empirical evidence that the scale-dependent statistical properties of solar wind and magnetohydrodynamic (MHD) turbulence can be described in terms of a family of parametric probability distribution functions (PDFs) known as Normal Inverse Gaussian (NIG). Understanding these PDFs is one of the most important goals in turbulence theory, as they are inherently connected to the intermittent properties of solar wind turbulence. We investigate the properties of PDFs of Elsasser increments based on a large statistical sample from solar wind observations and high-resolution numerical simulations of MHD turbulence. In order to measure the PDFs and their corresponding properties, three experiments are presented: fast and slow solar wind for experimental data and a simulation of reduced MHD (RMHD) turbulence. Conditional statistics on a 23-yr-long sample of WIND data near 1 au and high-resolution pseudo-spectral simulation of steadily driven RMHD turbulence on a $2048^3$ mesh are used to construct scale-dependent PDFs. The empirical PDFs are fitted to NIG distributions, which depend on four free parameters. Our analysis shows that NIG distributions accurately capture the evolution of the PDFs, with scale-dependent parameters, from large scales characterized by a Gaussian distribution, turning to exponential tails within the inertial range and stretched exponentials at dissipative scales. We also show that empirically-measured NIG parameters exhibit well-defined scaling properties that are similar across the three empirical data sets, which may be indicative of universal behaviour.

Astronomy & Astrophysics↗
NASA NTRS2002

Evolution of Galaxies and the Star Formation Rate in the Infrared

A central goal of extragalactic observational astronomy is to understand how normal galaxies evolve with redshift, and particularly when galaxies formed their stars. While optical and rest-frame UV (ultraviolet) observations have begun to address these issues, the interpretation of such data is particularly challenging because of the sensitivity to dust obscureness (at optical and UV wavelengths). The absorbed light is re-radiated at IR (infrared) wavelengths, hence the optimal indicators of the star formation rate (SFR) is at a rest-frame of (lambda) (is approximately equal to) 60 microns. The SIRTF (space infrared telescope facility) mission will revolutionize the study of the global evolution of the SFR by providing mass-selected, complete samples of galaxies and fares estimators of the SFR. This research program is to study the SFR using statistical samples of galaxies in the local universe, at intermediate redshifts, and set the stage for continuing studies up to z = 5.

Pahre, Michael, A.↗
NASA NTRS2021

New giant planet beyond the snow line for an extended MOA exoplanet microlens sample

Characterizing a planet detected by microlensing is hard if the planetary signal is weak or the lens-source relative trajectory is far from caustics. However, statistical analyses of planet demography must include those planets to accurately determine occurrence rates. As part of a systematic modelling effort in the context of a >10-yr retrospective analysis of MOA’s survey observations to build an extended MOA statistical sample, we analyse the light curve of the planetary microlensing event MOA-2014-BLG-472. This event provides weak constraints on the physical parameters of the lens, as a result of a planetary anomaly occurring at low magnification in the light curve. We use a Bayesian analysis to estimate the properties of the planet, based on a refined Galactic model and the assumption that all Milky Way’s stars have an equal planet-hosting probability. We find that a lens consisting of a 1.9(+2.2,−1.2)M(J) giant planet orbiting a 0.31(+0.36,−0.19)Mꙩ host at a projected separation of 0.75±0.24au is consistent with the observations and is most likely, based on the Galactic priors. The lens most probably lies in the Galactic bulge, at 7.2(+0.6,−1.7)kpc from Earth. The accurate measurement of the measured planet-to-host star mass ratio will be included in the next statistical analysis of cold planet demography detected by microlensing.

Clément Ranc↗
NASA NTRS2015

Optimization of Microphone Locations for Acoustic Liner Impedance Eduction

Two impedance eduction methods are explored for use with data acquired in the NASA Langley Grazing Flow Impedance Tube. The first is an indirect method based on the convected Helmholtz equation, and the second is a direct method based on the Kumaresan and Tufts algorithm. Synthesized no-flow data, with random jitter to represent measurement error, are used to evaluate a number of possible microphone locations. Statistical approaches are used to evaluate the suitability of each set of microphone locations. Given the computational resources required, small sample statistics are employed for the indirect method. Since the direct method is much less computationally intensive, a Monte Carlo approach is employed to gather its statistics. A comparison of results achieved with full and reduced sets of microphone locations is used to determine which sets of microphone locations are acceptable. For the indirect method, each array that includes microphones in all three regions (upstream and downstream hard wall sections, and liner test section) provides acceptable results, even when as few as eight microphones are employed. The best arrays employ microphones well away from the leading and trailing edges of the liner. The direct method is constrained to use microphones opposite the liner. Although a number of arrays are acceptable, the optimum set employs 14 microphones positioned well away from the leading and trailing edges of the liner. The selected sets of microphone locations are also evaluated with data measured for ceramic tubular and perforate-over-honeycomb liners at three flow conditions (Mach 0.0, 0.3, and 0.5). They compare favorably with results attained using all 53 microphone locations. Although different optimum microphone locations are selected for the two impedance eduction methods, there is significant overlap. Thus, the union of these two microphone arrays is preferred, as it supports usage of both methods. This array contains 3 microphones in the upstream hard wall section, 14 microphones opposite the liner, and 3 microphones in the downstream hard wall section.

Jones, M. G.↗
NASA NTRS2013

Orbital Debris: Past, Present, and Future

In the early days of spaceflight, the gBig Sky h theory was the near universally accepted paradigm for dealing with collisions of orbiting objects. This theory was also used during the early years of the aviation industry. Just as it did in aviation, the gBig Sky h theory breaks down as more and more objects accumulate in the environment. Fortunately, by the late 1970 fs some visionaries in NASA and the US Department of Defense (DoD) realized that trends in the orbital environment would inevitably lead to increased risks to operational spacecraft from collisions with other orbiting objects. The NASA Orbital Debris Program was established at and has been conducted at Johnson Space Center since 1979. At the start of 1979, fewer than 5000 objects were being tracked by the US Space Surveillance Network and very few attempts had been made to sample the environment for smaller sizes. Today, the number of tracked objects has quadrupled. Ground ]based and in situ measurements have statistically sampled the LEO environment over most sizes and mitigation guidelines and requirements are common among most space faring nations. NASA has been a leader, not only in defining the debris environment, but in promoting awareness of the issues in the US and internationally, and in providing leadership in developing policies to address the issue. This paper will discuss in broad terms the evolution of the NASA debris program from its beginnings to its present broad range of debris related research. The paper will discuss in some detail current research topics and will attempt to predict future research trends.

Stansbery, Gene↗
DOE OSTI2017

Statistical interpretation of the Amerithrax “morph” assay results

In the Amerithrax investigation PCR–based “morph assays” were used to link the anthrax letters with the RMR–1029 flask at USAMRIID. Quantitative data reported for several of these assays are not consistent with Poisson sampling statistics, but instead exhibit “Taylor's Law” behavior where the variance greatly exceeds the mean. A plausible statistical model for this behavior can explain the large number of observed negative and “inconclusive” findings, and implies a high likelihood that a repository sample could contain a “morph” mutant at concentrations well above the nominal detection limit but nonetheless give a negative or inconclusive test result. A Bayesian framework relates the assay results to the probability that a sample actually contains all four morph mutants, even though it tested negative for at least one. The analysis implies that the observed false negative rate actually does not significantly weaken the conclusion that the morph assays correctly excluded all but the stocks derived from RMR–1029 as possible sources of the letter powders, at least when the test results were unambiguous. In conclusion, these findings expand upon and resolve some of the issues cited in recent reviews, and indicate the importance of developing a rigorous statistical framework for interpreting “morph” assay data.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1998

Implementation and Testing of Turbulence Models for the F18-HARV Simulation

This report presents three methods of implementing the Dryden power spectral density model for atmospheric turbulence. Included are the equations which define the three methods and computer source code written in Advanced Continuous Simulation Language to implement the equations. Time-history plots and sample statistics of simulated turbulence results from executing the code in a test program are also presented. Power spectral densities were computed for sample sequences of turbulence and are plotted for comparison with the Dryden spectra. The three model implementations were installed in a nonlinear six-degree-of-freedom simulation of the High Alpha Research Vehicle airplane. Aircraft simulation responses to turbulence generated with the three implementations are presented as plots.

Yeager, Jessie C.↗