We present a catalog of the probabilistic redshift estimates and for 1366 individual Long-duration Gamma-Ray Bursts (LGRBs) detected by the Burst And Transient Source Experiment (BATSE). This result is based on a careful selection and modeling of the population distribution of 1366 BATSE LGRBs in the 5-dimensional space of redshift and the four intrinsic prompt gamma-ray emission properties: the isotropic 1024ms peak luminosity, the total isotropic emission, the spectral peak energy, as well as the intrinsic duration, while carefully taking into account the effects of sample incompleteness and the LGRB-detection mechanism of BATSE. Two fundamental plausible assumptions underlie our purely-probabilistic approach: 1. LGRBs trace, either exactly or closely, the Cosmic Star Formation Rate and 2. the joint 4-dimensional distribution of the aforementioned prompt gamma-ray emission properties is well-described by a multivariate log-normal distribution. Our modeling approach enables us to constrain the redshifts of individual BATSE LGRBs to within $0.36$ and $0.96$ average uncertainty ranges at $50\%$ and $90\%$ confidence levels, respectively. Our redshift predictions are completely at odds with the previous redshift estimates of BATSE LGRBs that were computed via the proposed phenomenological high-energy relations, specifically, the apparently-strong correlation of LGRBs' peak luminosity with the spectral peak energy, lightcurve variability, and the spectral lag. The observed discrepancies between our predictions and the previous works can be explained by the strong influence of detector threshold and sample-incompleteness in shaping these phenomenologically-proposed high-energy correlations in the literature.

We present a catalog of the redshift estimates and probability distributions for 1366 individual Long-duration Gamma-Ray Bursts (LGRBs) detected by the Burst And Transient Source Experiment (BATSE). This result is based on a careful classification and modeling of the population distribution of BATSE LGRBs in the 5-dimensional space of redshift as well as intrinsic prompt gamma-ray emission properties: peak luminosity, total isotropic emission, the spectral peak energy, and the intrinsic duration, while taking into account the detection mechanism of BATSE and sample incompleteness. The underlying assumption in our modeling approach is that LGRBs trace the Cosmic Star Formation Rate and that the joint 4-dimensional distribution of the aforementioned prompt gamma-ray emission properties follows a multivariate log-normal distribution. Our modeling approach enables us to constrain the redshifts of BATSE LGRBs to average uncertainty ranges of $0.7$ and $1.7$ at $50\%$ and $90\%$ confidence levels, respectively. Our predictions are almost entirely at odds with the previous estimates of BATSE redshifts based on the phenomenological high-energy correlations, in particular with the estimates derived from the lag-luminosity and the variability-luminosity relations. There is, however, a weak but significant correlation of strength $\sim0.26$ between our predicted redshift estimates and those derived from the hardness-brightness relations. The discrepancies between the estimates can be explained by the strong influence of sample incompleteness in shaping the phenomenologically proposed high-energy correlations in the literature. The presented catalog here can be useful for demographic studies of LGRBs and studies of individual BATSE events.

Model inadequacy and measurement uncertainty are two of the most confounding aspects of inference and prediction in quantitative sciences. The process of scientific inference (the inverse problem) and prediction (the forward problem) involve multiple steps of data analysis, hypothesis formation, model construction, parameter estimation, model validation, and finally, the prediction of the quantity of interest. This article seeks to clarify the concepts of model inadequacy and bias, measurement uncertainty, and the two traditional classes of uncertainty: aleatoric versus epistemic, as well as their relationships with each other in the process of scientific inference. Starting from basic principles of probability, we build and explain a hierarchical Bayesian framework to quantitatively deal with model inadequacy and noise in data. The methodology can be readily applied to many common inference and prediction problems in science, engineering, and statistics.

We present the results of a comprehensive study of the luminosity function, energetics, prompt gamma-ray correlations, and classification methodology of short-hard and long-soft GRBs (SGRBs and LGRBs), based on observational data in the largest catalog of GRBs available to this date: BATSE catalog of 2702 GRBs. We find that: 1. The least-biased classification method of GRBs into short and long, solely based on prompt-emission properties, appears to be the ratio of the observed spectral peak energy to the observed duration ($R=E_p/T_{90}$) with the dividing line at $R\simeq50[keV~s^{-1}]$. 2. Once data is carefully corrected for the effects of the detection threshold of gamma-ray instruments, the population distribution of SGRBs and LGRBs can be individually well described as multivariate log-normal distribution in the $4$--dimensional space of the isotropic peak gamma-ray luminosity, total isotropic gamma-ray emission, the intrinsic spectral peak energy, and the intrinsic duration. 3. Relatively large fractions of SGRBs and LGRBs with moderate-to-low spectral peak energies have been missed by BATSE detectors. 4. Relatively strong and highly significant intrinsic hardness--brightness and duration--brightness correlations likely exist in both populations of SGRBs and LGRBs, once data is corrected for selection effects. The strengths of these correlations are very similar in both populations, implying similar mechanisms at work in both GRB classes, leading to the emergence of these prompt gamma-ray correlations.

A model is presented here that is capable of simultaneously describing the luminosity function and the underlying joint population distribution of the prompt spectral and temporal parameters of Gamma-Ray Bursts (GRBs), subject to the detection threshold of gamma-ray instruments--in particular, BATSE and Fermi. Based on 2130 GRB prompt emission data in the BATSE catalog, I show that the population properties of the two classes of GRBs -- Long and Short durations -- bear striking similarities in the 4-dimensional space of prompt parameters: peak luminosity (Liso), total isotropic emission (Eiso), time-integrated spectral peak energy (Ep,z) and the prompt duration (T90,z). The two well-known Amati (Eiso-Ep,z) and Yonetoku (Liso-Ep,z) relations are shown to be highly affected by selection effects, undermining the legitimacy of their frequent uses in the studies of Dark Energy's equation of state and different cosmological models. In particular, I show that the slope of the Amati relation is likely in the range 0.2-0.3, corresponding to a Pearson's correlation strength of rho=0.58\pm0.04. This predicted slope is significantly (>16sigma) less than the value currently perceived by the GRB community (~0.56). I argue that similar Eiso-Ep,z and Liso-Ep,z relations with approximately same strength and significance should exist in the population of Short GRBs. Also predicted by the model is the strong positive correlation of the prompt emission duration (e.g., T90,z) with Eiso and Liso in both classes of Short and Long GRBs.

It is proposed that the luminosity function, the rest-frame spectral correlations and distributions of cosmological Long-duration (Type-II) Gamma-Ray Bursts (LGRBs) may be very well described as multivariate log-normal distribution. This result is based on careful selection, analysis and modeling of LGRBs' temporal and spectral variables in the largest catalog of Gamma-Ray Bursts available to date: 2130 BATSE GRBs, while taking into account the detection threshold and possible selection effects. Constraints on the joint rest-frame distribution of the isotropic peak luminosity (Liso), total isotropic emission (Eiso), the time-integrated spectral peak energy (Epkz) and duration (T90z) of LGRBs are derived. The presented analysis provides evidence for a relatively large fraction of LGRBs that have been missed by BATSE detector with Eiso extending down to ~ 10^49 [erg] and observed spectral peak energies (Epk) as low as ~ 5 [keV]. LGRBs with rest-frame duration T90z < 1 [s] or observer-frame duration T90 < 2 [s] appear to be rare events (<0.1% chance of occurrence). The model predicts a fairly strong but highly significant correlation (rho=0.58 \pm 0.04) between Eiso & Epkz of LGRBs. Also predicted are strong correlations of Liso & Eiso with T90z and moderate correlation between Liso & Epkz. The strength and significance of the correlations found, encourage the search for underlying mechanisms, though undermine their capabilities as probes of Dark Energy's equation of state at high redshifts. The presented analysis favors -- but does not necessitate -- a cosmic rate for BATSE LGRBs tracing metallicity evolution consistent with a cutoff ~ 0.2-0.5, assuming no luminosity-redshift evolution.

GRB satellites are relatively inefficient detectors of dim hard bursts because they trigger on photon counts, which are number-biased against hard photons. Therefore, for example, given two bursts of identical peak luminosity near the detection threshold, a dim soft burst will be preferentially detected over a dim hard burst. This detector bias can create or skew an apparent correlation where increasingly hard GRBs appear increasingly bright. Although such correlations may be obfuscated by a middle step where GRBs need to be bright enough to have their actual redshifts determined, it is found that the bias is generally pervasive. This result is derived here through simulations convolving a wide variety of possible GRB brightnesses and spectra with the BATSE Large Area Detectors (LAD) detection thresholds. The presented analyses indicate that the rest-frame $\nu F_{\nu}$ spectrum peak energy of long-duration GRBs, $\epi$, is not a good cosmological standard candle without significant corrections for selection effects. Therefore, the appearance of $\epi$ in seeming correlations such as the Amati ($E_{iso}-\epi$), Ghirlanda ($E_{\gamma}-\epi$), and $L_{iso}-\epi$ relations is statistically real but strongly influenced by so far uncalibrated GRB detector thresholds.

Simple hardness ratios are found to be a good estimator for the spectral peak energy in Gamma-Ray Bursts (GRBs). Specifically, a high correlation strength is found between the $\nu F_{\nu}$ peak in the spectrum of BATSE GRBs, $\epo$, and the hardness of GRBs, $\hr$, as defined by the fluences in channels 3 and 4, divided by the combined fluences in channels 1 and 2 of the BATSE Large Area Detectors. The correlation is independent of the type of the burst, whether Long-duration GRB (LGRB) or Short-duration (SGRB) and remains almost linear over the wide range of the BATSE energy window (20-2000 KeV). Based on Bayes theorem and Markov Chain Monte Carlo techniques, we also present multivariate analyses of the observational data while accounting for data truncation and sample-incompleteness. Prediction intervals for the proposed \hrep ~relation are derived. Results and further simulations are used to compute $\epo$ estimates for nearly the entire BATSE catalog: 2130 GRBs. These results may be useful for investigating the cosmological utility of the spectral peak in GRBs intrinsic luminosity estimates.