The separation of the gamma-ray bursts (GRBs) into short/hard and long/soft subclasses, respectively, is well supported both theoretically and observationally. The long ones are coupled to supernovae type Ib/Ic - the short ones are connected to the merging of two neutron stars, where one or even both neutron stars can be substituted by black holes. These short GRBs - as merging binaries - can also serve as the sources of gravitation waves, and are observable as the recently detected macronovae. Since 1998 there are several statistical studies suggesting the existence of more than two subgroups. There can be a subgroup having an intermediate durations; there can be a subgroup with ultra-long durations; the long/soft subgroup itself can be divided into two subclasses with respect to the luminosity of GRBs. The authors with other collaborators provided several statistical studies in this topic. This field of the GRB-diversity is briefly surveyed in this contribution.

It is well-known that there are two types of gamma-ray bursts (GRBs): short/hard and long/soft ones, respectively. The long GRBs are coupled to supernovae, but the short ones are associated with the so called macronovae (also known as kilonovae), which can serve as the sources of gravitational waves as well. The kilonovae can arise from the merging of two neutron-stars. The neutron stars can be substituded by more massive black holes as well. Hence, the topic of gamma-ray bursts (mainly the topic of short ones) and the topic of massive binaries, are strongly connected. In this contribution, the redshifts of GRBs are studied. The surprising result - namely that the apparently fainter GRBs can be in average at smaller distances - is discussed again. In essence, the results of M\'esz\'aros et al. (2011) are studied again using newer samples of GRBs. The former result is confirmed by the newer data.

Classification of gamma-ray bursts (GRBs) into groups has been intensively studied by various statistical tests in previous years. It has been suggested that there was a distinct group of GRBs, beyond the long and short ones, with intermediate durations. However, such a group is not securely confirmed yet. Strangely, concerning the spectral hardness, the observations from the Swift and RHESSI satellites give different results. For the Swift/BAT database it is found that the intermediate-duration bursts might well be related to so-called X-ray flashes (XRFs). On the other hand, for the RHESSI dataset the intermediate-duration bursts seem to be spectrally too hard to be given by XRFs. The connection of the intermediate-duration bursts and XRFs for the BATSE database is not clear as well. The purpose of this article is to check the relation between XRFs and GRBs for the BATSE and RHESSI databases, respectively. We use an empirical definition of XRFs introduced by other authors earlier. For the RHESSI database we also use a transformation between the detected counts and the fluences based on the simulated detector response function. The purpose is to compare the hardnesses of GRBs with the definition of XRFs. There is a 1.3-4.2 % fraction of XRFs in the whole BATSE database. The vast majority of the BATSE short bursts are not XRFs because only 0.7-5.7 % of the short bursts can be given by XRFs. However, there is a large uncertainty in the fraction of XRFs among the intermediate-duration bursts. The fraction of 1-85 % of the BATSE intermediate-duration bursts can be related to XRFs. For the long bursts this fraction is between 1.0 % and 3.4 %. The uncertainties in these fractions are large, however it can be claimed that all BATSE intermediate-duration bursts cannot be given by XRFs. At least 79 % of RHESSI short bursts, at least 53 % of RHESSI intermediate-duration bursts, and at least...

The Swift satellite made a real break through with measuring simultaneously the gamma X-ray and optical data of GRBs, effectively. Although, the satellite measures the gamma, X-ray and optical properties almost in the same time a significant fractions of GRBs remain undetected in the optical domain. In a large number of cases only an upper bound is obtained. Survival analysis is a tool for studying samples where a part of the cases has only an upper (lower) limit. The obtained survival function may depend on some other variables. The Cox regression is a way to study these dependencies. We studied the dependence of the optical brightness (obtained by the UVOT) on the gamma and X-ray properties, measured by the BAT and XRT on board of the Swift satellite. We showed that the gamma peak flux has the greatest impact on the afterglow's optical brightness while the gamma photon index and the X-ray flux do not. This effect probably originates in the energetics of the jet launched from the central engine of the GRB which triggers the afterglow.

Swift satellite measurements contributed substantially to the gamma-ray burst (GRB) redshift observations through fast slewing to the source of the GRBs. Still, a large number of bursts are without redshift. We study the celestial distribution of bursts with various methods and compare them to a random catalog using Monte-Carlo simulations. We find an anisotropy in the distribution of the intermediate class of bursts and find that the short and long population are distributed isotropically.

Gamma-ray bursts are the most luminous events in the Universe. Going beyond the short-long classification scheme we work in the context of three burst populations with the third group of intermediate duration and softest spectrum. We are looking for physical properties which discriminate the intermediate duration bursts from the other two classes. We use maximum likelihood fits to establish group memberships in the duration-hardness plane. To confirm these results we also use k-means and hierarchical clustering. We use Monte-Carlo simulations to test the significance of the existence of the intermediate group and we find it with 99.8% probability. The intermediate duration population has a significantly lower peak-flux (with 99.94% significance). Also, long bursts with measured redshift have higher peak-fluxes (with 98.6% significance) than long bursts without measured redshifts. As the third group is the softest, we argue that we have {related} them with X-ray flashes among the gamma-ray bursts. We give a new, probabilistic definition for this class of events.

After the discovery of the anisotropy in the sky-distribution of intermediate gamma-ray bursts recently also the distribution of the short gamma-ray bursts is proven to be anisotropic. The impact of these behaviors on the validity of the cosmological principle is shortly discussed.

A sample of 427 gamma-ray bursts (GRBs), measured by the RHESSI satellite, is studied statistically with respect to duration and hardness ratio. Standard statistical tests are used, such as $\chi^2$, F-test and the maximum likelihood ratio test, in order to compare the number of GRB groups in the RHESSI database with that of the BATSE database. Previous studies based on the BATSE Catalog claim the existence of an intermediate GRB group, besides the long and short groups. Using only the GRB duration $T_{90}$ as information and $\chi^2$ or F-test, we have not found any statistically significant intermediate group in the RHESSI data. However, maximum likelihood ratio test reveals a significant intermediate group. Also using the 2-dimensional hardness / $T_{90}$ plane, the maximum likelihood analysis reveals a significant intermediate group. Contrary to the BATSE database, the intermediate group in the RHESSI data-set is harder than the long one. The existence of an intermediate group follows not only from the BATSE data-set, but also from the RHESSI one.

We studied the complete randomness of the angular distribution of gamma-ray bursts (GRBs) detected by BATSE. Since GRBs seem to be a mixture of objects of different physical nature we divided the BATSE sample into 5 subsamples (short1, short2, intermediate, long1, long2) based on their durations and peak fluxes and studied the angular distributions separately. We used three methods, Voronoi tesselation, minimal spanning tree and multifractal spectra to search for non-randomness in the subsamples. To investigate the eventual non-randomness in the subsamples we defined 13 test-variables (9 from the Voronoi tesselation, 3 from the minimal spanning tree and one from the multifractal spectrum). Assuming that the point patterns obtained from the BATSE subsamples are fully random we made Monte Carlo simulations taking into account the BATSE's sky-exposure function. The MC simulations enabled us to test the null hypothesis i.e. that the angular distributions are fully random. We tested the randomness by binomial test and introducing squared Euclidean distances in the parameter space of the test-variables. We concluded that the short1, short2 groups deviate significantly (99.90%, 99.98%) from the fully randomness in the distribution of the squared Euclidean distances but it is not the case at the long samples. At the intermediate group the squared Euclidean distances also give significant deviation (98.51%).

Until 6 October 2005 sixteen redshifts have been measured of long gamma-ray bursts discovered by the Swift satellite. Further 45 redshifts have been measured of the long gamma-ray bursts discovered by other satellites. Here we perform five statistical tests comparing the redshift distributions of these two samples assuming - as the null hypothesis - identical distribution for the two samples. Three tests (Student's $t$-test, Mann-Whitney test, Kolmogorov-Smirnov test) reject the null hypothesis on the significance levels between 97.19 and 98.55%. Two different comparisons of the medians show extreme $(99.78-99.99994)$% significance levels of rejection. This means that the redshifts of the Swift sample and the redshifts of the non-Swift sample are distributed differently - in the Swift sample the redshifts are on average larger. This statistical result suggests that the long GRBs should on average be at the higher redshifts of the Swift sample.

It is known that the soft tail of the gamma-ray bursts' spectra show excesses from the exact power-law dependence. In this article we show that this departure can be detected in the peak flux ratios of different BATSE DISCSC energy channels. This effect allows to estimate the redshift of the bright long gamma-ray bursts in the BATSE Catalog. A verification of these redshifts is obtained for the 8 GRB which have both BATSE DISCSC data and measured optical spectroscopic redshifts. There is good correlation between the measured and esti redshifts, and the average error is $\Delta z \approx 0.33$. The method is similar to the photometric redshift estimation of galaxies in the optical range, hence it can be called as "gamma photometric redshift estimation". The estimated redshifts for the long bright gamma-ray bursts are up to $z \simeq 4$. For the the faint long bursts - which should be up to $z \simeq 20$ - the redshifts cannot be determined unambiguously with this method.