Search NASA⌕ Search

SEARCH · Search NASA

Results for “nonparametric testing”

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 19 records

Asymptotic Efficiency of New Distribution-Free Tests of Symmetry for Generalized Skew Alternatives

The Bahadur efficiency of new nonparametric tests of symmetry recently proposed by Nikitin and Ahsanullah is calculated. In contrast to this result, where only location alternatives were discussed, in the present paper generalized skew alternatives are of interest. It is shown that the new tests are highly efficient for a large class of skew alternatives. The problem of most favorable alternatives is also studied.

97 MATHEMATICS AND COMPUTING↗

Percent-Level Test of Isotropic Expansion Using Type Ia Supernovae

We propose and implement a novel, robust, and nonparametric test of statistical isotropy of the expansion of the Universe and apply it to around 1000 type Ia supernovae from the Pantheon sample. We calculate the angular clustering of supernova magnitude residuals and compare it to the noise expected under the isotropic assumption. We also test for systematic effects and demonstrate that their effects are negligible or are already accounted for in our procedure. Here, we express our constraints as an upper limit on the rms spatial variation in the Hubble parameter at late times. For the sky smoothed with a Gaussian with FWHM = 60°, less than 1% rms spatial variation in the Hubble parameter is allowed at 99.7% confidence.

79 ASTRONOMY AND ASTROPHYSICS↗

Dietary B group vitamin intake and the bladder cancer risk: a pooled analysis of prospective cohort studies

Abstract Purpose Diet may play an essential role in the aetiology of bladder cancer (BC). The B group complex vitamins involve diverse biological functions that could be influential in cancer prevention. The aim of the present study was to investigate the association between various components of the B group vitamin complex and BC risk. Methods Dietary data were pooled from four cohort studies. Food item intake was converted to daily intakes of B group vitamins and pooled multivariate hazard ratios (HRs), with corresponding 95% confidence intervals (CIs), were obtained using Cox-regression models. Dose–response relationships were examined using a nonparametric test for trend. Results In total, 2915 BC cases and 530,012 non-cases were included in the analyses. The present study showed an increased BC risk for moderate intake of vitamin B1 (HR B1 : 1.13, 95% CI: 1.00–1.20). In men, moderate intake of the vitamins B1, B2, energy-related vitamins and high intake of vitamin B1 were associated with an increased BC risk (HR (95% CI): 1.13 (1.02–1.26), 1.14 (1.02–1.26), 1.13 (1.02–1.26; 1.13 (1.02–1.26), respectively). In women, high intake of all vitamins and vitamin combinations, except for the entire complex, showed an inverse association (HR (95% CI): 0.80 (0.67–0.97), 0.83 (0.70–1.00); 0.77 (0.63–0.93), 0.73 (0.61–0.88), 0.82 (0.68–0.99), 0.79 (0.66–0.95), 0.80 (0.66–0.96), 0.74 (0.62–0.89), 0.76 (0.63–0.92), respectively). Dose–response analyses showed an increased BC risk for higher intake of vitamin B1 and B12. Conclusion Our findings highlight the importance of future research on the food sources of B group vitamins in the context of the overall and sex-stratified diet.

60 APPLIED LIFE SCIENCES↗

Influence of tumor location on the intensity-modulated radiation therapy plan of helical tomotherapy

Given the design of the Helical TomoTherapy device, the patient's central axis is routinely aligned with the machine's rotational axis to prevent the patient's body from colliding with the machine walls. However, for treatment of tumors located away from the patient's central axis, this position may not be optimal as the adequate radiation dose may not reach the affected site. Our study aimed to investigate the influence of tumor location on dose quality and delivery efficiency of tomotherapy plans. A phantom and 15 patients were selected for this study. Two plans, A and B, were implemented for each case. In plan A, the patient's central axis was aligned with the machine's rotational axis, whereas in plan B, the center of the planning target volume (PTV) was aligned with the machine's rotational axis. Both plans were optimized with the same planning parameters, and the dose quality of the plans was evaluated using dosimetrics. The delivery efficiency was determined from delivery time and monitor units (MUs). A paired t-test or nonparametric Wilcoxon signed-rank test was performed for statistical comparison. In the phantom study, the median delivery times were 358 and 336 seconds for plans A and B, respectively, and this difference was significant (p = 0.005). In the patient study, the median delivery times were 348 and 317 seconds for plans A and B, respectively, and this difference was also significant (p = 0.001). The dose qualities of both plans for each patient were nearly identical. No significant differences were found in the conformal index, heterogeneity index, and mean dose delivered to normal tissue between the plans. Both phantom and patient studies showed that for normal-sized patients, the delivery time reduced as the distance between the PTV and the patient's central axis increased when the PTV center was aligned with the machine axis. In conclusion, aligning the PTV center with the machine's rotational axis by shifting the patient during tomotherapy reduces the delivery time without compromising the dose quality of intensity-modulated radiation therapy.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Model-based approach to investigate equipment-induced error in pressure-waveform derived hemodynamic measurements

Advanced hemodynamic monitoring systems have provided less invasive methods for estimating pressure-derived measurements such as pressure-derived cardiac output (CO) measurements. These devices apply algorithms to arterial pressure waveforms recorded via pressure recording components that transmit the pressure signal to a pressure monitor. While standards have been developed for pressure monitoring equipment, it is unclear how the equipment-induced error can affect secondary measurements from pressure waveforms. We propose an approach for modelling different components of a pressure monitoring system and use this model-based approach to investigate the effect of different pressure recording configurations on pressure-derived hemodynamic measurements. The proposed model-based approach is a three-step process. (1) Modelling the response of pressure recording components using bench tests; (2) verifying the identified models through nonparametric equivalence tests; and (3) assessing the effects of pressure recording components on pressure-derived measurements. To delineate the application of this approach, we performed a series of model-based analyses to quantify the combined effect of a wide range of tubing configurations with various damping ratios and natural frequencies and monitors with different bandwidths on pressure waveforms and CO measurements by six pulse contour algorithms. Model-based results show the error in pressure-derived CO measurements because of tubing configurations with different natural frequencies and damping ratios. Tubing configurations with low natural frequencies (<23 Hz) altered characteristics of pressure waveforms in a way that affected the CO measurement, some by as much as 20%. Our method can serve as a tool to quantify the performance of pressure recording systems with different dynamic properties. In conclusion, this approach can be applied to investigate the effects of physiologic signal recording configurations on various pressure-derived hemodynamic measurements.

42 ENGINEERING↗

Anomalies in the topology of the temperature fluctuations in the cosmic microwave background: An analysis of the NPIPE and FFP10 data releases

We present a topological analysis of the temperature fluctuation maps from the Planck 2020 Data Release 4 NPIPE dataset and the Planck 2018 Data Release 3 FFP10 dataset. We performed a multiscale analysis in terms of the homology characteristics of the maps, invoking relative homology to account for the analysis in the presence of masks. We performed our analysis for a range of smoothing scales spanning sub- and super-horizon scales corresponding to a full width at half maximum (FWHM) of 5',10',20',40',80',160',320', and 640', and employed simulations based on the standard model for comparison, which assumes the initial fluctuation field to be an isotropic and homogeneous Gaussian random field. Examining the behavior of topological components, represented by the 0D homology group, we find the observations to be approximately 2σ or less deviant from the simulations for all resolutions and scales for the NPIPE dataset. For the FFP10 dataset, we detect a 2.96σ deviation between the observations and simulations at N = 128, FWHM = 80'. For the topological loops, represented by the first homology group, the simulations and observations are consistent within 2σ for most resolutions and scales for both the datasets. However, for the NPIPE dataset, we observe a high deviation between the observation and simulations in the number of loops at FWHM = 320', but at a low dimensionless threshold ν = –2.5. Under a Gaussian assumption, this would amount to a deviation of ~4σ. However, the distribution in this bin is manifestly non-Gaussian and does not obey Poisson statistics either. In the absence of a true theoretical understanding, we simply note that the significance is higher than what may be resolved by 600 simulations, yielding an empirical p-value of at most 0.0016. Specifically in this case, our tests indicate that the numbers arise from a statistically stable regime, despite being based on small numbers. For the FFP10 dataset, the differences are not as strong as for the NPIPE dataset, indicating a 2.77σ deviation at this resolution and threshold. The Euler characteristic, which is the alternating sum of the ranks of relative homology groups, reflects the deviations in the components and loops. To assess the significance of combined levels for a given scale, we employed the empirical and theoretical versions of the χ 2 test as well as the nonparametric Tukey depth test. Although all statistics exhibit a stable distribution, we favor the empirical version of the χ 2 test in the final interpretation, as it indicates the most conservative differences. For the NPIPE dataset, we find that the components and loops differ at more than 95%, but agree within the 99% confidence level with respect to the base model at N = 32, FWHM = 320'. The Euler characteristic at this resolution displays a per mil deviation. In contrast, the FFP10 dataset shows that the observations are consistent with the base model within the 95% confidence level, at this and smaller scales. This is consistent with the observations of the Planck analysis pipeline via Minkowski functionals. For the largest smoothing scale, N = 16, FWHM = 640', both datasets exhibit an anomalous behavior of the loops, where FFP10 data exhibit a deviation that is larger by an order of magnitude than that of the NPIPE dataset. In contrast, the values for the topological components and the Euler characteristic agree between observations and model to within a confidence level of 99%. However, for the largest scales, the statistics are based on low numbers and may have to be regarded with caution. Even though both datasets exhibit mild to significant discrepancies, they also exhibit contrasting behaviors at various instances. Therefore, we do not find it feasible to convincingly accept or reject the null hypothesis. Disregarding the large-scale anomalies that persist at similar scales in WMAP and Planck, observations of the cosmic microwave background are largely consistent with the standard cosmological model within 2σ.

79 ASTRONOMY AND ASTROPHYSICS↗

Probing parity violation with the four-point correlation function of BOSS galaxies

Parity-violating physics in the early universe can leave detectable traces in late-time observables. While vector- and tensor-type parity violation can be observed in the B -modes of the cosmic microwave background, scalar-type signatures are visible only in the four-point correlation function (4PCF) and beyond. This work presents a blind test for parity violation in the 4PCF of the baryon oscillation spectroscopic survey (BOSS) CMASS sample, considering galaxy separations in the range [20; 160]h -1 Mpc The parity-odd 4PCF contains no contributions from standard Λ CDM physics and can be efficiently measured using recently developed estimators. Data are analyzed using both a nonparametric rank test (comparing the BOSS 4PCFs to those of realistic simulations) and a compressed χ 2 analysis, with the former avoiding the assumption of a Gaussian likelihood. These find similar results, with the rank test giving a detection probability of 99.6% (2.9 σ). This provides significant evidence for parity violation, from either cosmological sources or systematics. We perform a number of systematic tests: although these do not reveal any observational artifacts, we cannot exclude the possibility that our detection is caused by the simulations not faithfully representing the statistical properties of the BOSS data. Our measurements can be used to constrain physical models of parity violation. As an example, we consider a coupling between the inflaton and a U (1) gauge field and place bounds on the latter’s energy density, which are several orders of magnitude stronger than those previously reported. Upcoming probes such as DESI and Euclid will reveal whether our detection of parity violation is due to new physics, and strengthen the bounds on a variety of models.

79 ASTRONOMY AND ASTROPHYSICS↗

A Nonparametric Method for the Inference of Halo Occupation Distributions

The galaxy–halo connection traces processes by which galaxies form and evolve. The halo occupation distribution (HOD) describes the relationship between galaxies and their host dark matter haloes. Measurements of the galaxy two-point correlation function (2PCF) allow us to extract information about the HODs of observed galaxy samples. Several parametric HOD models have been proposed in the literature, but the choice of parameterization restricts the space of possible HODs. To resolve this issue, we introduce a nonparametric HOD fitting method in which we train an emulator to learn the mappings among the galaxy 2PCF, physical properties used to select galaxy samples, and the HOD, all obtained from simulated past light cones constructed with the Santa Cruz semianalytic model. Implementing this emulator within a likelihood analysis framework, we derive constraints on the HOD of a galaxy sample when provided with a measurement of its 2PCF. Using the emulator to accelerate likelihood evaluations, we test the nonparametric HOD approach on a set of 2PCFs for mock galaxy samples drawn from the TNG100-1 simulation and selected above threshold values of stellar mass and star formation rate. Our framework is able to recover TNG100-1 HODs within 0.2 dex. We use the TNG100-1 mocks to tune the reported uncertainties to estimate those expected in the analysis of observations. Comparing to parametric HOD modelling routines applied to the same mock galaxy samples, our approach consistently infers the HOD with comparable or greater precision and accuracy.

Kennedy, Jacob [Rutgers Univ., Piscataway, NJ (Uni↗

Diffusion and perfusion MR parameters to assess preoperative short-course radiotherapy response in locally advanced rectal cancer: a comparative explorative study among Standardized Index of Shape by DCE-MRI, intravoxel incoherent motion- and diffusion kurtosis imaging-derived parameters

Purpose: To assess preoperative short-course radiotherapy (SCR) tumor response in locally advanced rectal cancer (LARC) by means of Standardized Index of Shape (SIS) by dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI), apparent diffusion coefficient (ADC), intravoxel incoherent motion (IVIM) and diffusion kurtosis imaging (DKI) parameters derived from diffusion-weighted MRI (DW-MRI). Materials and methods: Thirty-four patients with LARC who underwent MRI scans before and after SCR followed by delayed surgery, retrospectively, were enrolled. SIS, ADC, IVIM parameters [tissue diffusion (D{sub t}), pseudo-diffusion (D{sub p}), perfusion fraction (f{sub p})] and DKI parameters [mean diffusivity (MD), mean of diffusional kurtosis (MK)] were calculated for each patient. IVIM parameters were estimated using two methods, namely conventional biexponential fitting (CBFM) and variable projection (VARPRO). After surgery, the pathological TNM and tumor regression grade (TRG) were estimated. For each parameter, percentage changes between before and after SCR were evaluated. Furthermore, an artificial neural network was trained for outcome prediction. Nonparametric sample tests and receiver operating characteristic curve (ROC) analysis were performed. Results: Fifteen patients were classified as responders (TRG ≤ 2) and 19 as not responders (TRG > 3). Seven patients had TRG 1 (pathological complete response, pCR). Mean and standard deviation values of pre-treatment CBFM D{sub p} and mean value of VARPRO D{sub p} pre-treatment showed statistically significant differences to predict pCR. (p value at Mann–Whitney test was 0.05, 0.03 and 0.008, respectively.) Exclusively SIS percentage change showed significant differences between responder and non-responder patients after SCR (p value << 0.001) and to assess pCR after SCR (p value << 0.001). The best results to predict pCR were obtained by VARPRO Fp mean value pre-treatment with area under ROC of 0.84, a sensitivity of 96.4%, a specificity of 71.4%, a positive predictive value (PPV) of 92.9%, a negative predictive value (NPV) of 83.3% and an accuracy of 91.2%. The best results to assess after treatment complete pathological response were obtained by SIS with an area under ROC of 0.89, a sensitivity of 85.7%, a specificity of 92.6%, a PPV of 75.0%, a NPV of 96.1% and an accuracy of 91.2%. Moreover, the best results to differentiate after treatment responders vs. non-responders were obtained by SIS with an area under ROC of 0.94, a sensitivity of 93.3%, a specificity of 84.2%, a PPV of 82.4%, a NPV of 94.1% and an accuracy of 88.2%. Promising initial results were obtained using a decision tree tested with all ADC, IVIM and DKI extracted parameter: we reached high accuracy to assess pathological complete response after SCR in LARC (an accuracy of 85.3% to assess pathological complete response after SCR using VARPRO D{sub p} mean value post-treatment, ADC standard deviation value pre-treatment, MD standard deviation value post-treatment). Conclusion: SIS is a hopeful DCE-MRI angiogenic biomarker to assess preoperative treatment response after SCR with delayed surgery. Furthermore, an important prognostic role was obtained by VARPRO F{sub p} mean value pre-treatment and by a decision tree composed by diffusion parameters derived by DWI and DKI to assess pathological complete response.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Object-Based Evaluation of Dynamical and Statistical Downscaled Precipitation Products over CONUS

High-resolution precipitation data, generated through dynamical downscaling (DD) or statistical downscaling (SD) of global climate model output, provide critical information for regional climate assessment and adaptation planning. Most downscaling development and validation have focused on accurate gridscale precipitation construction and ignored the spatial structure of precipitation across model grids and at the event scale. However, many applications, e.g., hydrologic modeling and the analysis using the downscaled precipitation, require a reasonable representation of the spatial structure of precipitation within watersheds. Therefore, a set of standard metrics to evaluate the representation of the spatial structure of individual storms across diverse downscaled precipitation products is desired. To address this need, we conducted an object-based evaluation of precipitation in decades-long DD and SD products over the contiguous United States (CONUS). Specifically, we evaluate their ability to reproduce various features of precipitation objects in the observations: total volume, precipitation area, peak intensity, and spatial structure. Multiple metrics (bias, Perkins score, and nonparametric statistical tests) are used to quantify model performance. Our evaluation reveals notable variations in performance among individual products across different climate zones and seasons, as well as between extreme and nonextreme events. In general, most DD products exhibit balanced performance across the four precipitation object features, while SD products vary more significantly in their performance across products. Based on this comprehensive evaluation, we provide guidance on choosing downscaled products for specific regions, seasons, and precipitation object features. These findings and recommendations can inform precipitation-relevant modeling and analysis over CONUS, guide future downscaling technique developments, and provide actionable information for climate impact assessment and adaptation.

Downscaling↗

Object-Based Evaluation of Dynamical and Statistical Downscaled Precipitation Products over CONUS

High-resolution precipitation data, generated through dynamical downscaling (DD) or statistical downscaling (SD) of global climate model output, provide critical information for regional climate assessment and adaptation planning. Most downscaling development and validation have focused on accurate gridscale precipitation construction and ignored the spatial structure of precipitation across model grids and at the event scale. However, many applications, e.g., hydrologic modeling and the analysis using the downscaled precipitation, require a reasonable representation of the spatial structure of precipitation within watersheds. Therefore, a set of standard metrics to evaluate the representation of the spatial structure of individual storms across diverse downscaled precipitation products is desired. To address this need, we conducted an object-based evaluation of precipitation in decades-long DD and SD products over the contiguous United States (CONUS). Specifically, we evaluate their ability to reproduce various features of precipitation objects in the observations: total volume, precipitation area, peak intensity, and spatial structure. Multiple metrics (bias, Perkins score, and nonparametric statistical tests) are used to quantify model performance. Our evaluation reveals notable variations in performance among individual products across different climate zones and seasons, as well as between extreme and nonextreme events. In general, most DD products exhibit balanced performance across the four precipitation object features, while SD products vary more significantly in their performance across products. Based on this comprehensive evaluation, we provide guidance on choosing downscaled products for specific regions, seasons, and precipitation object features. These findings and recommendations can inform precipitation-relevant modeling and analysis over CONUS, guide future downscaling technique developments, and provide actionable information for climate impact assessment and adaptation.

Environmental sciences↗

Evaluating various composite sampling modes for detecting pathogenic SARS-CoV-2 virus in raw sewage

Inadequate sampling approaches to wastewater analyses can introduce biases, leading to inaccurate results such as false negatives and significant over- or underestimation of average daily viral concentrations, due to the sporadic nature of viral input. To address this challenge, we conducted a field trial within the University of Tennessee residence halls, employing different composite sampling modes that encompassed different time intervals (1 h, 2 h, 4 h, 6 h, and 24 h) across various time windows (morning, afternoon, evening, and late-night). Our primary objective was to identify the optimal approach for generating representative composite samples of SARS-CoV-2 from raw wastewater. Utilizing reverse transcription-quantitative polymerase chain reaction, we quantified the levels of SARS-CoV-2 RNA and pepper mild mottle virus (PMMoV) RNA in raw sewage. Our findings consistently demonstrated that PMMoV RNA, an indicator virus of human fecal contamination in water environment, exhibited higher abundance and lower variability compared to pathogenic SARS-CoV-2 RNA. Significantly, both SARS-CoV-2 and PMMoV RNA exhibited greater variability in 1 h individual composite samples throughout the entire sampling period, contrasting with the stability observed in other time-based composite samples. Through a comprehensive analysis of various composite sampling modes using the Quade Nonparametric ANCOVA test with date, PMMoV concentration and site as covariates, we concluded that employing a composite sampler during a focused 6 h morning window for pathogenic SARS-CoV-2 RNA is a pragmatic and cost-effective strategy for achieving representative composite samples within a single day in wastewater-based epidemiology applications. This method has the potential to significantly enhance the accuracy and reliability of data collected at the community level, thereby contributing to more informed public health decision-making during a pandemic.

sampling timing↗

Observed Spatiotemporal Changes in the Mechanisms of Extreme Water Available for Runoff in the Western United States

We used a quality-controlled Snowpack Telemetry (SNOTEL) data set (1979–2017) combined with the nonparametric regional Kendall test (RKT) to examine changes in annual maximum water available for runoff at the land surface under four hydrometeorological conditions (snowmelt only, rain-on-snow, all melt with or without rainfall, and all melt plus rainfall on snow-free ground) over the mountainous regions of the western U.S. Our RKT analyses indicated significant declining trends at point (SNOTEL) scale under all four conditions. The annual maximum of all melt plus rainfall decreased primarily in the southwestern U.S., while the frequency of rain-on-snow events increased significantly in the northwestern U.S. The rate of annual maximum snowmelt only decreased in 10 of 11 ecological regions across the western U.S. These results suggest that under a warming climate, the four types of annual maximum water available for runoff have been substantially decreasing at local scale across the western U.S.

54 ENVIRONMENTAL SCIENCES↗

Spatially Dependent Multiple Testing Under Model Misspecification, With Application to Detection of Anthropogenic Influence on Extreme Climate Events

The Weather Risk Attribution Forecast (WRAF) is a forecasting tool that uses output from global climate models to make simultaneous attribution statements about whether and how greenhouse gas emissions have contributed to extreme weather across the globe. However, in conducting a large number of simultaneous hypothesis tests, the WRAF is prone to identifying false “discoveries.” A common technique for addressing this multiple testing problem is to adjust the procedure in a way that controls the proportion of true null hypotheses that are incorrectly rejected, or the false discovery rate (FDR). Unfortunately, generic FDR procedures suffer from low power when the hypotheses are dependent, and techniques designed to account for dependence are sensitive to misspecification of the underlying statistical model. In this paper, we develop a Bayesian decision-theoretical approach for dependent multiple testing and a nonparametric hierarchical statistical model that flexibly controls false discovery and is robust to model misspecification. We illustrate the robustness of our procedure to model error with a simulation study, using a framework that accounts for generic spatial dependence and allows the practitioner to flexibly specify the decision criteria. Finally, we apply our procedure to several seasonal forecasts and discuss implementation for the WRAF workflow. Lastly, supplementary materials for this article, including a standardized description of the materials available for reproducing the work, are available as an online supplement.

54 ENVIRONMENTAL SCIENCES↗

Strong Influence of Irrigation on Water Budget and Land Surface Temperature in Indian Subcontinental River Basins

Irrigation modulates the terrestrial water budget, especially in intensively irrigated regions such as Indian sub-continent river basins, where 35% of agricultural land has been irrigated during 2010. Here, we investigate the impact of irrigation on water budget components and land surface temperature (LST), as well as seasonal and annual trends of irrigation water demand in response to climate variability over the Indian sub-continental region. We use the Variable Infiltration Capacity (VIC) model coupled with an irrigation scheme (VIC-IRR) to represent the irrigation activities. The VIC-IRR model was calibrated and evaluated against in-situ and remote sensing observations including streamflow from India-Water Resource Information System (WRIS), Terrestrial Water Storage Change (TWSC) from Gravity Recovery and Climate Experiment (GRACE), LST from Advanced Along-Track Scanning Radiometer (AATSR), and observed gravimetric soil moisture data from India Meteorology Department (IMD). We conducted two simulations (i.e., with and without irrigation) over 18 Indian sub-continental basins at 0.25° spatial resolution from 1951 to 2012. The comparison between the two simulations indicates that as a result of irrigation, evapotranspiration (ET) increases in all the basins. We find that mean annual ET of Indus and Ganges basins increases by 47% and 12%, respectively due to irrigation. After irrigation, basin averaged total runoff (TR) was increased in the most irrigated basins (i.e., Sabarmati, Ganges, and Indus) by 4.62%, 11.87% and 7.06% in the monsoon, pre-monsoon, and post-monsoon seasons, respectively. The land surface temperature was cooled down in range of 0.2-0.65 °C in the Indian sub-continental basins especially in Ganges and Indus basins. Compare to other basins; irrigation water demand was high in Indus and Ganges basins during the post-monsoon season. Using the nonparametric Mann–Kendall trend test we find that in the pre-monsoon and monsoon seasons over the study period, irrigation demand decreased by 8.25 mm/year in 13 basins. Whereas for the post-monsoon season, irrigation water demand increased in 12 out of the 18 basins. Changes in irrigation demands are largely driven by changes in soil moisture due to decline in the summer monsoon precipitation.

54 ENVIRONMENTAL SCIENCES↗

Quantifying Model Form Uncertainty of Epidemic Forecasting Models from Incidence Data

Objective We present a mathematical framework for non-parametric estimation of the force of infection, together with statistical upper and lower confidence bands. The resulting estimates allow to assess how well simpler models, such as SEIR, fit the observed time series of incidence data.IntroductionUncertainty Quantification (UQ), the ability to quantify the impact of sample-to-sample variations and model misspecification on predictions and forecasts, is a critical aspect of disease surveillance. While quantifying the impact of stochastic uncertainty in the data is well understood, quantifying the impact of model misspecification is significantly harder. For the latter, one needs a "universal model" to which more restrictive parametric models are compared too.MethodsThis talk presents a useful modeling framework for time series of incidence data from contagious diseases that enables one to identify and quantify the impact of model form uncertainty. Specifically, we propose to focus on estimating the timedependent force of infection. The latter is a universal parameters for all contagious disease model. Using a machine learning technique for estimating monotone functions, i.e., isotonic regression and its variants, one can estimate the force ofinfection without addtional assumptions. We note that most contagious disease model do satisfy this monotonicity assumption, due to a combination of factors: depletion of susceptibles, implementation of mitigation strategies, behavior change, etc. Comparing the resulting "non-parametric" estimate with parametric estimates, obtained by fitting an SEIR for example, can reveal model deficiencies and help quantify model form uncertainties.Finally, we discuss how ideas from "strict bound theory" can be used to develop upper and lower uncertainty bands for force of infection that acknowledge the intrinsic stochasticity in the data.ResultsWe demonstrate the application of the methodology to weekly Influenza Like Illness (ILI) incidence data from France andcompare the results to fitted SIR and SEIR models. This comparison can be seen as a nonparametric goodness of fit test, providing one with tools to do simple model selection.ConclusionsWe present a novel and flexible model to statistically describe the force of infection as a function of time. Comparing the fit to incidence data of that model with the fit of simpler parametric models enables the quantification of model form uncertainty and associated model selection.

Hengartner, Nicholas↗

Normative Ranges for, and Interrater Reliability of, Rotational Vestibular and Balance Tests in U.S. Military Service Members and Veterans

Purpose: The objectives of this study were to (a) describe normative ranges—expressed as reference intervals (RIs)—for vestibular and balance function tests in a cohort of Service Members and Veterans (SMVs) and (b) to describe the interrater reliability of these tests. Method: As part of the Defense and Veterans Brain Injury Center (DVBIC)/Traumatic Brain Injury Center of Excellence 15-year Longitudinal Traumatic Brain Injury (TBI) Study, participants completed the following: vestibulo-ocular reflex suppression, visual-vestibular enhancement, subjective visual vertical, subjective visual horizontal, sinusoidal harmonic acceleration, the computerized rotational head impulse test (crHIT), and the sensory organization test. RIs were calculated using nonparametric methods and interrater reliability was assessed using intraclass correlation coefficients between three audiologists who independently reviewed and cleaned the data. Results: Reference populations for each outcome measure comprised 40 to 72 individuals, 19 to 61 years of age, who served either as noninjured controls (NIC) or injured controls (IC) in the 15-year study; none had a history of TBI or blast exposure. A subset of 15 SMVs from the NIC, IC, and TBI groups were included in the interrater reliability calculations. RIs are reported for 27 outcome measures from the seven rotational vestibular and balance tests. Interrater reliability was considered excellent for all tests except the crHIT, which was found to have good interrater reliability. Conclusion: This study provides clinicians and scientists with important information regarding normative ranges and interrater reliability for rotational vestibular and balance tests in SMVs.

Audiology & Speech-Language Pathology↗