The afterglow of a gamma ray burst (GRB) can give us valuable insight into the properties of its host galaxy. To correctly interpret the spectra of the afterglow we need to have a good understanding of the foreground interstellar medium (ISM) in our own Galaxy. The common practice to correct for the foreground is to use neutral hydrogen (HI) data from the Leiden/Argentina/Bonn (LAB) survey. However, the poor spatial resolution of the single dish data may have a significant effect on the derived column densities. To investigate this, we present new high-resolution HI observations with the Australia Telescope Compact Array (ATCA) towards 4 GRBs. We combine the interferometric ATCA data with single dish data from the Galactic All Sky Survey (GASS) and derive new Galactic HI column densities towards the GRBs. We use these new foreground column densities to fit the Swift XRT X-ray spectra and calculate new intrinsic hydrogen column density values for the GRB host galaxies. We find that the new ATCA data shows higher Galactic HI column densities compared to the previous single dish data, which results in lower intrinsic column densities for the hosts. We investigate the line of sight optical depth near the GRBs and find that it may not be negligible towards one of the GRBs, which indicates that the intrinsic hydrogen column density of its host galaxy may be even lower. In addition, we compare our results to column densities derived from far-infrared data and find a reasonable agreement with the HI data.

THESEUS is a space mission concept aimed at exploiting Gamma-Ray Bursts for investigating the early Universe and at providing a substantial advancement of multi-messenger and time-domain astrophysics. These goals will be achieved through a unique combination of instruments allowing GRB and X-ray transient detection over a broad field of view (more than 1sr) with 0.5-1 arcmin localization, an energy band extending from several MeV down to 0.3 keV and high sensitivity to transient sources in the soft X-ray domain, as well as on-board prompt (few minutes) follow-up with a 0.7 m class IR telescope with both imaging and spectroscopic capabilities. THESEUS will be perfectly suited for addressing the main open issues in cosmology such as, e.g., star formation rate and metallicity evolution of the inter-stellar and intra-galactic medium up to redshift $\sim$10, signatures of Pop III stars, sources and physics of re-ionization, and the faint end of the galaxy luminosity function. In addition, it will provide unprecedented capability to monitor the X-ray variable sky, thus detecting, localizing, and identifying the electromagnetic counterparts to sources of gravitational radiation, which may be routinely detected in the late '20s / early '30s by next generation facilities like aLIGO/ aVirgo, eLISA, KAGRA, and Einstein Telescope. THESEUS will also provide powerful synergies with the next generation of multi-wavelength observatories (e.g., LSST, ELT, SKA, CTA, ATHENA).

GRB170817A, associated with the LIGO-Virgo GW170817 neutron-star merger event, lacks the short duration and hard spectrum of a Short gamma-ray burst (GRB) expected from long-standing classification models. Correctly identifying the class to which this burst belongs requires comparison with other GRBs detected by the Fermi GBM. The aim of our analysis is to classify Fermi GRBs and to test whether or not GRB170817A belongs - as suggested - to the Short GRB class. The Fermi GBM catalog provides a large database with many measured variables that can be used to explore gamma-ray burst classification. We use statistical techniques to look for clustering in a sample of 1298 gamma-ray bursts described by duration and spectral hardness. Classification of the detected bursts shows that GRB170817A most likely belongs to the Intermediate, rather than the Short GRB class. We discuss this result in light of theoretical neutron-star merger models and existing GRB classification schemes. It appears that GRB classification schemes may not yet be linked to appropriate theoretical models, and that theoretical models may not yet adequately account for known GRB class properties. We conclude that GRB170817A may not fit into a simple phenomenological classification scheme.

A significant number of the parameters of a gamma-ray burst (GRB) and its host galaxy are calculated from the afterglow. There are various methods obtaining extinction values for the necessary correction for galactic foreground. These are: galaxy counts, from HI 21 cm surveys, from spectroscopic measurements and colors of nearby Galactic stars, or using extinction maps calculated from infrared surveys towards the GRB. We demonstrate that AKARI Far-Infrared Surveyor sky surface brightness maps are useful uncovering the fine structure of the galactic foreground of GRBs. Galactic cirrus structures of a number of GRBs are calculated with a 2 arcminute resolution, and the results are compared to that of other methods.

The point sources in the Bright Source Catalogue of the AKARI Far-Infrared Surveyor (FIS) were classified based on their FIR and mid-IR fluxes and colours into young stellar object (YSO) and extragalactic source types using Quadratic Discriminant Analysis method (QDA) and Support Vector Machines (SVM). The reliability of the selection of YSO candidates is high, and the number of known YSO candidates were increased significantly, that we demonstrate in the case of the nearby open cluster IC348. Our results show that we can separate galactic and extragalactic AKARI point sources in the multidimensional space of FIR fluxes and colours with high reliability, however, differentiating among the extragalactic sub-types needs further information.

Short and long bursts were identified by the BATSE team in the early 90s. A decade ago there were some suggestions about the intermediate duration type of bursts. We are going to summarize recent analyses of the duration distributions of the Beppo-Sax and Swift data. Our conclusion is all the three satellites (CGRO, Swift, Beppo-Sax) can see the third type of the GRBs. The properties of the group members are very similar in the different data sets.

We studied the unbiased optical brightness distribution which was calculated from the survival analysis of host galaxies and its relationship with the Swift GRB data of the host galaxies observed by the Keck telescopes. Based on the sample obtained from merging the Swift GRB table and the Keck optical data we also studied the dependence of this distribution on the data of the GRBs. Finally, we compared the HGs distribution with standard galaxies distribution which is in the DEEP2 galaxies catalog.

According to the cosmological principle, Universal large-scale structure is homogeneous and isotropic. The observable Universe, however, shows complex structures even on very large scales. The recent discoveries of structures significantly exceeding the transition scale of 370 Mpc pose a challenge to the cosmological principle. We report here the discovery of the largest regular formation in the observable Universe; a ring with a diameter of 1720 Mpc, displayed by 9 gamma ray bursts (GRBs), exceeding by a factor of five the transition scale to the homogeneous and isotropic distribution. The ring has a major diameter of $43^o$ and a minor diameter of $30^o$ at a distance of 2770 Mpc in the 0.78<z<0.86 redshift range, with a probability of $2\times 10^{-6}$ of being the result of a random fluctuation in the GRB count rate. Evidence suggests that this feature is the projection of a shell onto the plane of the sky. Voids and string-like formations are common outcomes of large-scale structure. However, these structures have maximum sizes of 150 Mpc, which are an order of magnitude smaller than the observed GRB ring diameter. Evidence in support of the shell interpretation requires that temporal information of the transient GRBs be included in the analysis. This ring-shaped feature is large enough to contradict the cosmological principle. The physical mechanism responsible for causing it is unknown.

The statistical tests - done by the authors - are surveyed, which verify the null-hypothesis of the intrinsic randomness in the angular distribution of gamma-ray bursts collected at BATSE Catalog. The tests use the counts-in-cells method, an analysis of spherical harmonics, a test based on the two-point correlation function and a method based on multiscale methods. The tests suggest that the intermediate subclass of gamma-ray bursts are distributed anisotropically.

One of the most important task of the Gamma-Ray Burst field is the classification of the bursts. Many researches have proven the existence of the third kind (intermediate duration) of GRBs in the BATSE data. Recent works have analyzed BeppoSax and Swift observations and can also identify three types of GRBs in the data sets. However, the class memberships are probabilistic we have enough observed redshifts to calculate the redshift and spatial distribution of the intermediate GRBs. They are significantly farther than the short bursts and seems to be closer than the long ones.

We studied the distribution of the GRBs, observed by the Fermi satellite, in the multidimensional parameter space consisting of the duration, Fluence, Peak flux and Peak energy (if it was available). About 10% of the Fermi bursts was observed also by the Swift satellite. We did not find significant differences between the Peak flux and Peak energy of GRBs observed and not observed also by the Swift satellite. In contrast, those GRBs detected also by the Swift satellite had significantly greater Fluence and duration. We did a similar study for the GRBs detected by the Swift satellite. About 30% percent of these bursts was also measured by the Fermi satellite. We found a significant difference in the Fluence, Peak flux and Photon index but none in duration. These differences may be accounted for the different construction and observing strategy of the Fermi and Swift satellites.

Gamma-Ray Bursts (GRBs) are the most powerful cosmic explosions since the Big Bang, and thus act as signposts throughout the distant Universe. Over the last 2 decades, these ultra-luminous cosmological explosions have been transformed from a mere curiosity to essential tools for the study of high-redshift stars and galaxies, early structure formation and the evolution of chemical elements. In the future, GRBs will likely provide a powerful probe of the epoch of reionisation of the Universe, constrain the properties of the first generation of stars, and play an important role in the revolution of multi-messenger astronomy by associating neutrinos or gravitational wave (GW) signals with GRBs. Here, we describe the next steps needed to advance the GRB field, as well as the potential of GRBs for studying the Early Universe and their role in the up-coming multi-messenger revolution.

In the case of Gamma-Ray Bursts with measured redshift, we can calculate the k-correction to get the fluence and energy that were actually produced in the comoving system of the GRB. To achieve this we have to use well-fitted parameters of a GRB spectrum, available in the GCN database. The output of the calculations is the comoving isotropic energy E_iso, but this is not the endpoint: this data can be useful for estimating the {\Omega}M parameter of the Universe and for making a GRB Hubble diagram using Amati's relation.

The Fermi satellite has a particular motion during its flight which enables it to catch the gamma-ray bursts mostly well. The side-effect of this favourable feature is that the lightcurves of the GBM detectors are stressed by rapidly and extremely varying background. Before this data is processed, it needs to be separated from the background. The commonly used methods were useless for most cases of Fermi, so we developed a new technique based on the motion and orientation of the satellite. The background-free lightcurve can be used to perform statistical surveys, hence we showed the efficiency of our background-filtering method presenting a statistical analysis known from the literature.

Gamma-ray bursts are usually classified through their high-energy emission into short-duration and long-duration bursts, which presumably reflect two different types of progenitors. However, it has been shown on statistical grounds that a third, intermediate population is needed in this classification scheme, although an extensive study of the properties of this class has so far not been done. The large amount of follow-up studies generated during the Swift era allows us to have a suficient sample to attempt a study of this third population through the properties of their prompt emission and their afterglows. Our study is focused on a sample of GRBs observed by Swift during its first four years of operation. The sample contains those bursts with measured redshift since this allows us to derive intrinsic properties. Intermediate bursts are less energetic and have dimmer afterglows than long GRBs, especially when considering the X-ray light curves, which are on average one order of magnitude fainter than long bursts. There is a less significant trend in the redshift distribution that places intermediate bursts closer than long bursts. Except for this, intermediate bursts show similar properties to long bursts. In particular, they follow the Epeak vs. Eiso correlation and have, on average, positive spectral lags with a distribution similar to that of long bursts. Like long GRBs, they normally have an associated supernova, although some intermediate bursts have shown no supernova component. This study shows that intermediate bursts are different from short bursts and, in spite of sharing many properties with long bursts, there are some differences between them as well. We suggest that the physical difference between intermediate and long bursts could be that for the first the ejecta are thin shells while for the latter they are thick shells.

Recently the anisotropy of the short gamma-ray bursts detected by BATSE was announced (Vavrek et al. 2008). The impact of this discovery on cosmology is discussed. It is shown that the anisotropy found may cause the breakdown of the cosmological principle.

Earlier classification analyses found three types of gamma-ray bursts (short, long and intermediate in duration) in the BATSE sample. Recent works have shown that these three groups are also present in the RHESSI and the BeppoSAX databases. The duration distribution analysis of the bursts observed by the Swift satellite also favors the three-component model. In this paper, we extend the analysis of the Swift data with spectral information. We show, using the spectral hardness and the duration simultaneously, that the maximum likelihood method favors the three-component against the two-component model. The likelihood also shows that a fourth component is not needed.

The Gamma Ray Inverse Problem is discussed. Four methods of spectral deconvolution are presented here and applied to the BATSE's MER data type. We compare these to the Band spectra.

As GRBs trace the high-z Universe, there is a non-negligible probability of a lensing effect being imprinted on the lightcurves of the bursts. We propose to search for lensed candidates with a cross-correlation method, by looking at bursts days to years apart coming from the same part of the sky. We look for similarities and hypothesize a Singular Isothermal Sphere (SIS) model for the lens. A lensed pair would enable us to constrain the mass of the lensing object. Our search did not reveal any gravitationally lensed events.

The power spectra of 457 short BATSE bursts were analyzed, focusing on the 64 ms lightcurves' tails in the low energy bands. Using MC simulations, 22 GRBs were identified with unusually high harmonic power above 0.03 Hz. The sky distribution of these bursts shows an extraordinarily strong dipole moment with a 99.994% significance.

Two classes of gamma-ray bursts were identified in the BATSE catalogs characterized by their durations. There were also some indications for the existence of a third type of gamma-ray bursts. Swift satellite detectors have different spectral sensitivity than pre-Swift ones for GRBs. Therefore in this paper we analyze the bursts' duration distribution and also the duration-hardness bivariate distribution, published in The First BAT Catalog. Similarly to the BATSE data, to explain the BAT GRBs' duration distribution three components are needed. Although, the relative frequencies of the groups are different than they were in the BATSE GRB sample, the difference in the instrument spectral sensitivities can explain this bias in a natural way. This means theoretical models may have to explain three different type of gamma-ray bursts.

Two classes of gamma-ray bursts have been identified in the BATSE catalogs characterized by durations shorter and longer than about 2 seconds. There are, however, some indications for the existence of a third type of burst. Swift satellite detectors have different spectral sensitivity than pre-Swift ones for gamma-ray bursts. Therefore it is worth to reanalyze the durations and their distribution and also the classification of GRBs. Using The First BAT Catalog the maximum likelihood estimation was used to analyzed the duration distribution of GRBs. The three log-normal fit is significantly (99.54 % probability) better than the two for the duration distribution. Monte-Carlo simulations also confirm this probability (99.2 %).

In a systematic search of the OTs at GRBs the Swift satellite determined only an upper limit of the apparent brightness in a significant fraction of cases. Combining these upper limits with the really measured OT brightness we obtained a sample well suited to survival analysis. Performing a Kaplan-Meier product limit estimation we obtained an unbiased cumulative distribution of the V visual brightness. The lg(N(V)) logarithmic cumulative distribution can be well fitted with a linear function of V in the form of lg(N(V))= 0.234 V + const. We studied the dependence of V on the gamma ray properties of the bursts. We tested the dependence on the fluence, T90 duration and peak flux. We found a dependence on the peak flux on the 99.7% significance level.

Using the discriminant analysis of the multivariate statistical analysis we compared the distribution of the physical quantities of the optically dark and bright GRBs, detected by the BAT and XRT on board of the Swift Satellite. We found that the GRBs having detected optical transients (OT) have systematically higher peak fluxes and lower HI column densities than those without OT.

We presented the results of several statistical tests of the randomness in the angular sky-distribution of gamma-ray bursts in BATSE Catalog. Thirteen different tests were presented based on Voronoi tesselation, Minimal spanning tree and Multifractal spectrum for five classes (short1, short2, intermediate, long1, long2) of gamma-ray bursts, separately. The long1 and long2 classes are distributed randomly. The intermediate subclass, in accordance with the earlier results of the authors, is distributed non-randomly. Concerning the short subclass earlier statistical tests also suggested some departure from the random distribution, but not on a high enough confidence level. The new tests presented in this article suggest also non-randomness here.

A sample of 197 long BATSE GRBs is studied statistically. In the sample 11 variables, describing for any burst the time behavior of the spectra and other quantities, are collected. The application of the factor analysis on this sample shows that five factors describe the sample satisfactorily. Both the pseudo-redshifts coming from the variability and the Amati-relation in its original form are disfavored.

In this paper we make an attempt to combine the two kinds of data from the Swift-XRT instrument (windowed timing and photon counting modes) and the from BAT. A thorough desription of the applied procedure will be given. We apply various binning techniques to the different data: Bayes blocks, exponential binning and signal-to-noise type of binning. We present a handful of lightcurves and some possible applications.

The measured redshifts of gamma-ray bursts (GRBs), which were first detected by the Swift satellite, seem to be bigger on average than the redshifts of GRBs detected by other satellites. We analyzed the redshift distribution of GRBs triggered and observed by different satellites (Swift[1], HETE2[2], BeppoSax, Ulyssses). After considering the possible biases significant difference was found at the p = 95.70% level in the redshift distributions of GRBs measured by HETE and the Swift.

We study statistically 197 long gamma-ray bursts, detected and measured in detail by the BATSE instrument of the Compton Gamma-Ray Observatory. In the sample 10 variables, describing for any burst the time behavior of the spectra and other quantities, are collected. The factor analysis method is used to find the latent random variables describing the temporal and spectral properties of GRBs. The application of this particular method to this sample indicates that five factors and the $\REpk$ spectral variable (the ratio of peak energies in the spectrum) describe the sample satisfactorily. Both the pseudo-redshifts inferred from the variability, and the Amati-relation in its original form, are disfavored.

Two classes of gamma-ray bursts have been identified in the BATSE catalogs characterized by durations shorter and longer than about 2 seconds. There are, however, some indications for the existence of a third class. Swift satellite detectors have different spectral sensitivity than pre-Swift ones for gamma-ray bursts. Therefore, it is worth to reanalyze the durations and their distribution. We analyze, the maximum likelihood estimation, the bursts duration distribution, published in The First BAT Catalog, whether it contains two, three or more groups. The three log-normal fit is significantly (99.54% probability) better than the two for the duration distribution. Monte-Carlo simulations also confirm this probability (99.2%). Similarly, in previous results we found that the fourth component is not needed. The relative frequencies of the distribution of the groups are 7% short 35% intermediate and 58% long. Although the relative frequencies of the groups are different than in the BATSE GRB sample, the difference in the instrument spectral sensitivities can explain this bias on a natural way. This means theoretical models may be needed to explain three different type of gamma-ray bursts.

Principal component analysis is a statistical method, which lowers the number of important variables in a data set. The use of this method for the bursts' spectra and afterglows is discussed in this paper. The analysis indicates that three principal components are enough among the eight ones to describe the variablity of the data. The correlation between spectral index alpha and the redshift suggests that the thermal emission component becomes more dominant at larger redshifts.

The low energy spectra of some gamma-ray bursts' show excess components beside the power-law dependence. The consequences of such a feature allows to estimate the gamma photometric redshift of the long gamma-ray bursts in the BATSE Catalog. There is good correlation between the measured optical and the estimated gamma photometric redshifts. The estimated redshift values for the long bright gamma-ray bursts are up to z=4, while for the the faint long bursts - which should be up to z=20 - the redshifts cannot be determined unambiguously with this method. The redshift distribution of all the gamma-ray bursts with known optical redshift agrees quite well with the BATSE based gamma photometric redshift distribution.

We provided separate bivariate log-normal distribution fits to the BATSE short and long burst samples using the durations and fluences. We show that these fits present an evidence for a power-law dependence between the fluence and the duration, with a statistically significant different index for the long and short groups. We argue that the effect is probably real, and the two subgroups are different physical phenomena. This may provide a potentially useful constraint for models of long and short bursts.

The origin of dark bursts - i.e. that have no observed afterglows in X-ray, optical/NIR and radio ranges - is unclear yet. Different possibilities - instrumental biases, very high redshifts, extinction in the host galaxies - are discussed and shown to be important. On the other hand, the dark bursts should not form a new subgroup of long gamma-ray bursts themselves.

Gamma-ray bursts can be divided into three groups ("short", "intermediate", "long") with respect to their durations. The third type of gamma-ray bursts - as known - has the intermediate duration. We show that the intermediate group is the softest one. An anticorrelation between the hardness and the duration is found for this subclass in contrast to the short and long groups.

We studied the statistical properties of the luminosity distance and internal extinction data of type Ia supernovae in the lists published by Tonry et al. (2003) and Barris et al. (2004). After selecting the luminosity distance in an empty Universe as a reference level we divided the sample into low $z<0.25$ and high $z \ge 0.25$ parts. We further divided these subsamples by the median of the internal extinction. Performing sign tests using the standardized residuals between the estimated logarithmic luminosity distances and those of an empty universe, on the four subsamples separately, we recognized that the residuals were distributed symmetrically in the low redshift region, independently from the internal extinction. On the contrary, the low extinction part of the data of $z \ge 0.25$ clearly showed an excess of the points with respect to an empty Universe which was not the case in the high extinction region. This diversity pointed to an interrelation between the estimated luminosity distance and internal extinction. To characterize quantitatively this interrelation we introduced a hidden variable making use of the technics of factor analysis. After subtracting that part of the residual which was explained by the hiddenmaking use of the technics of factor analysis. After subtracting that part of the residual which was explained by the hidden variable we obtained luminosity distances which were already free from interrelation with internal extinction. Fitting the corrected luminosity distances with cosmological models we concluded that the SN Ia data alone did not exclude the possibility of the $\Lambda=0$ solution.

The redshift distribution of gamma-ray bursts collected in the BATSE Catalog is compared with the star formation rate. We aim to clarify the accordance between them. We also study the case of comoving number density of bursts monotonously increasing up to redshift 6-20. A method independent of the models of the gamma-ray bursts is used. The short and the long subgroups are studied separately. The redshift distribution of the long bursts may be proportional to the star formation rate. For the short bursts this can also happen, but the proportionality is less evident. For the long bursts the monotonously increasing scenario is also less probable but still can occur. For the short bursts this alternative seems to be excluded.

Gamma-ray bursts can be divided into three groups ("short", "intermediate", "long") with respect to their durations. This classification is somewhat imprecise, since the subgroup of intermediate duration has an admixture of both short and long bursts. In this paper a physically more reasonable definition of the intermediate group is presented, using also the hardnesses of the bursts. It is shown again that the existence of the three groups is real, no further groups are needed. The intermediate group is the softest one. From this new definition it follows that 11% of all bursts belong to this group. An anticorrelation between the hardness and the duration is found for this subclass in contrast to the short and long groups. Despite this difference it is not clear yet whether this group represents a physically different phenomenon.

The CGRO/BATSE database includes many types of data such as the 16-channel continuous background or medium energy resolution burst data (CONT and MER data types). We have calculated some four hundred burst's medium energy resolution spectra and Principal Component Analysis has been applied. We found five components can describe GRBs' spectra.

The gamma-ray burst's lightcurves have been analyzed using a special wavelet transformation. The applied wavelet base is based on a typical Fast Rise-Exponential Decay (FRED) pulse. The shape of the wavelet coefficients' total distribution is determined on the observational frequency grid. Our analysis indicates that the pulses in the long bursts' high energy channel lightcurves are more FRED-like than the lower ones, independently from the actual physical time-scale.

We present an analysis indicating that there is a correlation between the fluences and the durations of gamma-ray bursts, and provide arguments that this reflects a correlation between the total emitted energies and the intrinsic durations. For the short (long) bursts the total emitted energies are roughly proportional to the first (second) power of the intrinsic duration. This difference in the energy-duration relationship is statistically significant, and may provide an interesting constraint on models aiming to explain the short and long gamma-ray bursts.

It is shown that in the duration-hardness plane the GRBs of the third intermediate subgroup are well defined. Their durations are intermediate (i.e. roughly between 2 and 10 seconds), but their hardnesses are the lowest. They are even softer than the long bursts.

The Compton Gamma Ray Observatory (CGRO) observed many types of data and one of them is the time-tagged photon events (TTE data). We use the Bayesian block analysis, using Bayesian statistics, analyses the TTE data. Our results; calculations of duration (T100), count rates (burst photon numbers in different channels) and count peaks (in 64, 16 and 4 ms). We present the duration, the peak duration and the distance between peaks distributions. Principal Component Analysis (PCA) has been also applied. The PCA shows interesting results, such as channel 4 (highest energy channel) probably is very important.

We argue that the distributions of both the intrinsic fluence and the intrinsic duration of the gamma-ray emission in gamma-ray bursts from the BATSE sample are well represented by log-normal distributions, in which the intrinsic dispersion is much larger than the cosmological time dilatation and redshift effects. We perform separate bivariate log-normal distribution fits to the BATSE short and long burst samples. The bivariate log-normal behaviour results in an ellipsoidal distribution, whose major axis determines an overall statistical relation between the fluence and the duration. We show that this fit provides evidence for a power-law dependence between the fluence and the duration, with a statistically significant different index for the long and short groups. We discuss possible biases, which might affect this result, and argue that the effect is probably real. This may provide a potentially useful constraint for models of long and short bursts.

We argue that the distributions of both the intrinsic fluence and the intrinsic duration of the gamma-ray emission in gamma-ray bursts from the BATSE sample are well represented by log-normal distributions, in which the intrinsic dispersion is much larger than the cosmological time dilatation and redshift effects. We perform separate bivariate log-normal distribution fits to the BATSE short and long burst samples. The bivariate log-normal behavior results in an ellipsoidal distribution, whose major axis determines an overall statistical relation between the fluence and the duration. We show that this fit provides evidence for a power-law dependence between the fluence and the duration, with a statistically significant different index for the long and short groups. We discuss possible biases which might affect this result, and argue that the effect is probably real. This may provide a potentially useful constraint for models of long and short bursts.

In the article a test is developed, which allows to test the null-hypothesis of the intrinsic randomness in the angular distribution of gamma-ray bursts collected at the Current BATSE Catalog. The method is a modified version of the well-known counts-in-cells test, and fully eliminates the non-uniform sky-exposure function of BATSE instrument. Applying this method to the case of all gamma-ray bursts no intrinsic non-randomness is found. The test also did not find intrinsic non-randomnesses for the short and long gamma-ray bursts, respectively. On the other hand, using the method to the new intermediate subclass of gamma-ray bursts, the null-hypothesis of the intrinsic randomness for 181 intermediate gamma-ray bursts is rejected on the 96.4% confidence level. Taking 92 dimmer bursts from this subclass itself, we obtain the surprising result: This "dim" subclass of the intermediate subclass has an intrinsic non-randomness on the 99.3% confidence level. On the other hand, the 89 "bright" bursts show no intrinsic non-randomness.

The counts-in-cells and the two-point angular correlation function method are used to test the randomnesses in the angular distributions of both the all gamma-ray bursts collected at BATSE Catalog, and also their three subclasses ("short", "intermediate", "long"). The methods elimate the non-zero sky-exposure function of BATSE instrument. Both tests suggest intrinsic non-randomnesses for the intermediate subclass; for the remaining three cases only the correlation function method. The confidence levels are between 95% and 99.9%. Separating the GRBs into two parts ("dim half" and "bright half", respectively) we obtain the result that the "dim" half shows a non-randomness on the 99.3% confidence level from the counts-in-cells test.

The anisotropy of the sky distribution of 2025 gamma-ray bursts (GRBs) collected in Current BATSE catalog is confirmed. It is shown that the quadrupole term being proportional to $\sim \sin 2b \sin l$ is non-zero with a probability 99.9%. The occurrence of this anisotropy term is then supported by the binomial test even with the probability 99.97%. It is also argued that this anisotropy cannot be caused exclusively by instrumental effects due to the non-uniform sky exposure of BATSE instrument; there should exist also some intrinsic anisotropy in the angular distribution of GRBs. Separating GRBs into short and long subclasses, it is shown that the 251 short ones are distributed anisotropically, but the 681 long ones seem to be distributed still isotropically. The 2-sample Kolmogorov-Smirnov test shows that they are distributed differently with a 98.7% probability. The character of anisotropy suggests that the cosmological origin of short GRBs further holds, and there is no evidence for their Galactical origin. The work in essence contains the key ideas and results of a recently published paper (\cite{balazs}), to which the new result following from the 2-sample Kolmogorov-Smirnov test is added, too.

The isotropy of gamma-ray bursts collected in Current BATSE catalog is studied. It is shown that the quadrupole term being proportional to $\sim \sin 2b \sin l$ is non-zero with a probability 99.9%. The occurrence of this anisotropy term is then confirmed by the binomial test even with the probability 99.97 %. Hence, the sky distribution of all known gamma-ray bursts is anisotropic. It is also argued that this anisotropy cannot be caused exclusively by instrumental effects due to the nonuniform sky exposure of BATSE instrument. Separating the GRBs into short and long subclasses, it is shown that the short ones are distributed anisotropically, but the long ones seem to be distributed still isotropically. The character of anisotropy suggests that the cosmological origin of short GRBs further holds, and there is no evidence for their Galactical origin.

We have carried out a principal component analysis for 625 gamma-ray bursts in the BATSE 3B catalog for which non-zero values exist for the nine measured variables. This shows that only two out of the three basic quantities of duration, peak flux and fluence are independent, even if this relation is strongly affected by instrumental effects, and these two account for 91.6% of the total information content. The next most important variable is the fluence in the fourth energy channel (at energies above 320 keV). This has a larger variance and is less correlated with the fluences in the remaining three channels than the latter correlate among themselves. Thus a separate consideration of the fourth channel, and increased attention on the related hardness ratio $H43$ appears useful for future studies. The analysis gives the weights for the individual measurements needed to define a single duration, peak flux and fluence. It also shows that, in logarithmic variables, the hardness ratio $H32$ is significantly correlated with peak flux, while $H43$ is significantly anticorrelated with peak flux. The principal component analysis provides a potentially useful tool for estimating the improvement in information content to be achieved by considering alternative variables or performing various corrections on available measurements