A hundred years ago (1920) in the auditorium of the Smithsonian Institution's U.S. National Museum there were two lectures under the auspices of the George Ellery Hale Lecture series, what has come to be called the 'Great Debate'. In the debate, Harlow Shapley and Heber Curtis argued over the 'Scale of the Universe'. Curtis argued that the Universe is composed of many galaxies like our own and they are relatively small. Shapley argued that the Universe was composed of only one big Galaxy. In Shapley's model, our Sun was far from the center of this great island Universe.

Studying the GRBs' gamma-ray spectra may reveal some physical information of Gamma-ray Bursts. The Fermi satellite observed more than two thousand GRBs. The FERMIGBRST catalog contains GRB parameters (peak energy, spectral indices, intensity) estimated for both the total emission (fluence), and the emission during the interval of the peak flux. We found a relationship with linear discriminant analysis between the spectral categories and the model independent physical data. We compared the Swift and Fermi spectral types. We found a connection between the Fermi fluence spectra and the Swift spectra but the result of the peak flux spectra can be disputable. We found that those GRBs which were observed by both Swift and Fermi can similarly discriminate as the complete Fermi sample. We concluded that the common observation probably did not mean any trace of selection effects in the spectral behavior of GRBs.

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.

The Fermi GBM catalog provides a large database with many measured variables that can be used to explore and verify gamma-ray burst classification results. We have used Principal Component Analysis and statistical clustering techniques to look for clustering in a sample of 801 gamma-ray bursts described by sixteen classification variables. The analysis recovers what appears to be the Short class and two long-duration classes that differ from one another via peak flux, with negligible variations in fluence, duration and spectral hardness. Neither class has properties entirely consistent with the Intermediate GRB class. Spectral hardness has been a critical Intermediate class property. Rather than providing spectral hardness, Fermi GBM provides a range of fitting variables for four different spectral models; it is not intuitive how these variables can be used to support or disprove previous GRB classification results.

In the first part of this paper we summarize the discovery and the early study of the gamma-ray bursts. The second part studies the spatial and sky distribution of the bursters. The brightness distribution is also studied. In the conclusion we discuss that these are supporting the cosmological origin of the gamma-ray bursts.

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.

In 2011 the Nobel Prize in Physics was awarded for the 1998 discovery of the nonzero cosmological constant. This discovery is very important and surely worth to receive the Nobel Prize. However, years earlier several papers had been published (Paal, Horvath, & Lukacs 1992; Holba et al. 1992, Holba et al. 1994) about a very similar discovery from observational data.

We analyze pulse properties of Short gamma-ray bursts (GRBs) from a new catalog containing 434 pulses from 387 BATSE Time-Tagged Event (TTE) GRBs. Short GRB pulses exhibit correlated properties of duration, fluence, hardness, and amplitude, and they evolve hard-to-soft while undergoing similar triple-peaked light curves similar to those found in Long/Intermediate bursts. We classify pulse light curves using their temporal complexities, demonstrating that Short GRB pulses exhibit a range of complexities from smooth to highly variable. Most of the bright, hard, chaotic emission seen in complex pulses seems to represent a separate highly-variable emission component. Unlike Long/Intermediate bursts, as many as 90\% of Short GRBs are single-pulsed. However, emission in Short multi-pulsed bursts is coupled such that the first pulse's duration is a predictor of both the interpulse separation and subsequent pulse durations. These results strongly support the idea that external shocks produce the prompt emission seen in Short GRBs. The similarities between the triple-peaked structures and spectral evolution of Long, Short, and Intermediate GRBs then suggests that external shocks are responsible for the prompt emission observed in all GRB classes. In addition to these findings, we identify a new type of gamma-ray transient in which peak amplitudes occur at the end of the burst rather than at earlier times. Some of these "Crescendo" bursts are preceded by rapid-fire "Staccato" pulses, whereas the remaining are preceded by a variable episode that could be unresolved staccato pulses.

Statistical studies of gamma-ray burst (GRB) spectra may result in important information on the physics of GRBs. The Fermi GBM catalog contains GRB parameters (peak energy, spectral indices, intensity) estimated fitting the gamma-ray SED of the total emission (fluence, flnc), and during the time of the peak flux pflx. Using contingency tables we studied the relationship of the models best fitting pflx and flnc time intervals. Our analysis revealed an ordering of the spectra into a power law - Comptonized - smoothly broken power law - Band series. This result was further supported by a correspondence analysis (CA) of the pflx and flnc spectra categorical variables. We performed a linear discriminant analysis (LDA) to find a relationship between categorical (spectral) and model independent physical data. LDA resulted in highly significant physical differences among the spectral types, that is more pronounced in the case of the pflx spectra, than for the flnc spectra. We interpreted this difference as caused by the temporal variation of the spectrum during the outburst. This spectral variability is confirmed by the differences in the low energy spectral index and peak energy, between the pflx and flnc spectra. We found that the synchrotron radiation is significant in GBM spectra. The mean low energy spectral index is close to the canonical value of {\alpha} = -2/3 during the peak flux. However, $\alpha$ is ~ -0.9 for the spectra of the fluences. We interpret this difference as showing that the effect of cooling is important only for the fluence spectra.

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.

Decades ago two classes of gamma-ray bursts were identified and delineated as having durations shorter and longer than about 2 s. Subsequently indications also supported the existence of a third class. Using maximum likelihood estimation we analyze the duration distribution of 888 Swift BAT bursts observed before October 2015. Fitting three log-normal functions to the duration distribution of the bursts provides a better fit than two log-normal distributions, with 99.9999% significance. Similarly to earlier results, we found that a fourth component is not needed. The relative frequencies of the distribution of the groups are 8% for short, 35% for intermediate and 57% for long bursts which correspond to our previous results. We analyse the redshift distribution for the 269 GRBs of the 888 GRBs with known redshift. We find no evidence for the previously suggested difference between the long and intermediate GRBs' redshift distribution. The observed redshift distribution of the 20 short GRBs differs with high significance from the distributions of the other groups.

Several large structures, including the Sloan Great Wall, the Huge Large Quasar Group, and a large gamma-ray burst cluster referred to as the Hercules-Corona Borealis Great Wall, appear to exceed the maximum structural size predicted by Universal inflationary models. The existence of very large structures such as these might necessitate cosmological model modifications. Gamma-ray bursts are the most luminous sources found in nature. They are associated with the stellar endpoints of massive stars and are found in and near distant galaxies. Since they are viable indicators of the dense part of the Universe containing normal matter, the spatial distribution of gamma-ray bursts can serve as tracers of Universal large-scale structure. An increased sample size of gamma-ray bursts with known redshift provides us with the opportunity to validate or invalidate the existence of the Hercules-Corona Borealis Great Wall. Nearest-neighbour tests are used to search the larger sample for evidence of clustering and a bootstrap point-radius method is used to estimate the angular cluster size. The potential influence of angular sampling biasing is studied to determine the viability of the results. The larger gamma-ray burst database further supports the existence of a statistically significant gamma-ray burst cluster at 1.6 < z < 2.1 with an estimated angular size of 2000-3000 Mpc. Although small number statistics limit our angular resolution and do not rule out the existence of adjacent and/or line-of-sight smaller structures, these structures must still clump together in order for us to see the large gamma-ray burst cluster detected here. This cluster provides support for the existence of very large-scale universal heterogeneities.

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.

Swift's remarkable ability to quickly localize gamma-ray bursts has led to the accumulation of a sizable burst sample for which both angular locations and redshifts are measured. This sample has become large enough that it can potentially be used to probe angular anisotropies indicative of large-scale universal structure. In a previous work, a large clustering of gamma-ray bursts at redshift z about 2 was reported in the general direction of the constellations of Hercules and Corona Borealis. Since that report, a 42 per cent increase in the number of z about 2 gamma-ray bursts has been observed, warranting an updated analysis. Surprisingly, the cluster is more pronounced now than it was when it was first reported.

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 Fermi GBM Catalog has been recently published. Previous classification analyses of the BATSE, RHESSI, BeppoSAX, and Swift databases found three types of gamma-ray bursts. Now we analyzed the GBM catalog to classify the GRBs. PCA and Multiclustering analysis revealed three groups. Validation of these groups, in terms of the observed variables, shows that one of the groups coincides with the short GRBs. The other two groups split the long class into a bright and dim part, as defined by the peak flux. Additional analysis is needed to determine whether this splitting is only a mathematical byproduct of the analysis or has some real physical meaning.

Gamma-ray burst sources are distributed with a high level of isotropy, which is compatible with either a cosmological origin or an extended Galactic halo origin. The brightness distribution is another indicator used to characterize the spatial distribution in distance. In this paper the author discusses detailed fits of the BATSE gamma-ray burst peak-flux distributions with Friedmann models taking into account possible density evolution and standard candle luminosity functions. A chi-square analysis is used to estimate the goodness of the fits and the author derives the significance level of limits on the density evolution and luminosity function parameters. Cosmological models provide a good fit over a range of parameter space which is physically reasonable

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.

Context. Research over the past three decades has revolutionized cosmology while supporting the standard cosmological model. However, the cosmological principle of Universal homogeneity and isotropy has always been in question, since structures as large as the survey size have always been found each time the survey size has increased. Until 2013, the largest known structure in our Universe was the Sloan Great Wall, which is more than 400 Mpc long located approximately one billion light years away. Aims. Gamma-ray bursts are the most energetic explosions in the Universe. As they are associated with the stellar endpoints of massive stars and are found in and near distant galaxies, they are viable indicators of the dense part of the Universe containing normal matter. The spatial distribution of gamma-ray bursts can thus help expose the large scale structure of the Universe. Methods. As of July 2012, 283 GRB redshifts have been measured. Subdividing this sample into nine radial parts, each containing 31 GRBs, indicates that the GRB sample having 1.6 < z < 2.1 differs significantly from the others in that 14 of the 31 GRBs are concentrated in roughly 1/8 of the sky. A two-dimensional Kolmogorov-Smirnov test, a nearest-neighbour test, and a Bootstrap Point-Radius Method explore the significance of this clustering. Results. All tests used indicate that there is a statistically significant clustering of the GRB sample at 1.6 < z < 2.1. Furthermore, this angular excess cannot be entirely attributed to known selection biases, making its existence due to chance unlikely. Conclusions. This huge structure lies ten times farther away than the Sloan Great Wall, at a distance of approximately ten billion light years. The size of the structure defined by these GRBs is about 2000-3000 Mpc, or more than six times the size of the Sloan Great Wall.

Research over the past three decades has revolutionized the field of cosmology while supporting the standard cosmological model. However, the cosmological principle of Universal homogeneity and isotropy has always been in question, since structures as large as the survey size have always been found as the survey size has increased. Until now, the largest known structure in our Universe is the Sloan Great Wall (SGW), which is more than 400 Mpc long and located approximately one billion light-years away. Here we report the discovery of a structure at least six times larger than the Sloan Great Wall that is suggested by the distribution of gamma-ray bursts (GRBs). Gamma-ray bursts are the most energetic explosions in the Universe. They are associated with the stellar endpoints of massive stars and are found in and near distant galaxies. Therefore, they are very good indicators of the dense part of the Universe containing normal matter. As of July 2012, 283 GRB redshifts have been measured. If one subdivides this GRB sample into nine radial parts and compares the sky distributions of these subsamples (each containing 31 GRBs), one can observe that the fourth subsample (1.6 < z < 2.1) differs significantly from the others in that many of the GRBs are concentrated in the same angular area of the sky. Using the two-dimensional Kolmogorov-Smirnov test, the significance of this observation is found to be less than 0.05 per cent. Fourteen out of the 31 Gamma-Ray Bursts in this redshift band are concentrated in approximately 1/8 of the sky. The binomial probability to find such a deviation is p=0.0000055. This huge structure lies ten times farther away than the Sloan Great Wall, at a distance of approximately ten billion light-years. The size of the structure defined by these GRBs is about 2000-3000 Mpc, or more than six times the size of the largest known object (SGW) in the Universe.

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.

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.

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.

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.

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.

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.

GRB 080319B is one of the brightest and most extensively sampled bursts. It has good coverage at many wavelengths. Here we present the optical observations of the Palomar 60 inch telescope, which spans a long time interval after the burst. We augment the optical dataset with freely available Swift BAT and XRT observations reduced by us. We also compare our conclusions with the published parameters from the rich literature about this burst.

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.

Aims: We present multiwavelength observations of one of the most distant gamma-ray bursts detected so far, GRB 080913. Based on these observations, we consider whether it could be classified as a short-duration GRB and discuss the implications for the progenitor nature and energy extraction mechanisms. Methods: Multiwavelength X-ray, near IR and millimetre observations were made between 20.7 hours and 16.8 days after the event. Results: Whereas a very faint afterglow was seen at the 3.5m CAHA telescope in the nIR, the X-ray afterglow was clearly detected in both Swift and XMM-Newton observations. An upper limit is reported in the mm range. We have modeled the data assuming a collimated $\theta_0$ $\gtrsim$ 3$^\circ$ blast wave with an energy injection at 0.5 days carrying $5\sim 10^{52}$ erg or approximately 12 times the initial energy of the blast wave. We find that GRB 080913 shares many of the gamma-ray diagnostics with the more recent burst GRB 090423 for being classified as short had they ocurred at low redshift. If the progenitor were a compact binary merger, it is likely composed by a NS and BH. The Blandford-Znajek (BZ) mechanism is the preferred one to extract energy from the central, maximally-rotating BH. Both the magnetic field close to the event horizon (B) and the BH mass ($M_{bh}$) are restricted within a relatively narrow range, such that $(B / 3\times 10^{16} \rm{G}) (M_{bh} / 7 M_\odot) \sim 1$. Similar constraints on the central BH hold for collapsar-like progenitor systems if the BZ-mechanism works for the system at hand.

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.

The BeppoSAX Catalog has been very recently published. In this paper we analyze - using the Maximum Likelihood (ML) method - the duration distribution of the 1003 GRBs listed in the catalog with duration. The ML method can identify the long and the intermediate duration groups. The short population of the bursts is identified only at a 96% significance level. MC simulation has been also applied and gives a similar significance level; 95%. However, the existence of the short bursts is not a scientific question after the Compton Gamma-Ray Observatory's observation. Our minor result is this well-known fact that in the BeppoSAX data the short bursts are underrepresented, mainly caused by the different triggering system. Our major result is the identification of the intermediate group in the BeppoSAX data. Therefore, four different satellites (CGRO, Swift, RHESSI and BeppoSAX) observed the intermediate type Gamma-Ray Burst.

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.

Two classes of Gamma Ray Bursts have been identified so far, characterized by T90 durations shorter and longer than approximately 2 seconds. It was shown that the BATSE 3B data allow a good fit with three Gaussian distributions in log T90. In the same Volume in ApJ. another paper suggested that the third class of GRBs is may exist. Using the full BATSE catalog here we present the maximum likelihood estimation, which gives us 0.5% probability to having only two subclasses. The MC simulation confirms this probability.

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.

Two classes of gamma-ray bursts have been identified so far, characterized by durations shorter and longer than approximately 2 seconds. In 1998 two independent papers indicated the existence of the third class of the bursts roughly duration between 2 and 10 seconds. In this paper, using the full BATSE Catalog, the maximum likelihood estimation is presented, which gives a 0.5% probability to having only two subclasses. The Monte-Carlo simulation confirms this probability, too.

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.

Two classes of Gamma Ray Bursts have been identified so far, characterized by T_90 durations shorter and longer than approximately 2 seconds. We show here that the BATSE 3B data allow a good fit with three Gaussian distributions in log T_90. The chi^2 statistic indicates a 40 % probability for two Gaussians, whereas the three-Gaussian fit probability is 98 %. Using another statistical method, it is argued that the probability that the third class is a random fluctuation is less than 0.02 %.

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

We discuss detailed fits of the BATSE and PVO gamma-ray burst peak-flux distributions with Friedman models taking into account possible density evolution and standard candle or power law luminosity functions. A chi-square analysis is used to estimate the goodness of the fits and we derive the significance level of limits on the density evolution and luminosity function parameters. Cosmological models provide a good fit over a range of parameter space which is physically reasonable.