Review is devoted to estimators (introduced in astronomy by famous astronomer Ambartsumian and his followers) that are used to estimating the unknown number of flare stars and other randomly flashing objects (stars with X-ray flares, solar type stars with superflares). Some important astronomical applications of them are presented. The development of these methods in astronomy have proceeded regardless of development of analogous methods of mathematical statistics. Brief history of this development and parallels with similar statistical methods is presented.

We intensively studied the flare activity on the stellar object KIC011764567. The star was thought to be solar type, with a temperature of $T_{eff} = (5640 \pm 200)\,$K, $\log(g) = (4.3 \pm 0.3)\,$dex and a rotational period of Prot 22 d (Brown et al. 2011). High resolution spectra turn the target to an evolved object with Teff = (5300 \pm 150) K, a metalicity of $[m/H] = (-0.5 \pm 0.2)$, a surface gravity of $log(g) = (3.3 \pm 0.4)\,$dex, and a projected rotational velocity of $v sin i = (22 \pm 1)\,kms^{-1}$. Within an observing time span of 4 years we detected 150 flares in Kepler data in an energy range of $10^{36} - 10^{37}$ erg. From a dynamical Lomb-Scargle periodogram we have evidence for differential rotation as well as for stellar spot evolution and migration. Analysing the occurrence times of the flares we found hints for a periodic flare frequency cycle of $430 - 460 $d, the significance increases with an increasing threshold of the flares equivalent duration. One explanation is a very short activity cycle of the star with that period. Another possibility, also proposed by others in similar cases, is that the larger flares may be triggered by external phenomena, such as magnetically interaction with an unseen companion. Our high resolution spectra show that KIC011764567 is not a short period binary star.

The method for determination of the Supernovae (SNe) explosions frequency distribution function based on the assumption of explosions independence are offered. The method is based on assumption that the sequence of SNe explosions in an individual galaxy is a Poisson sequence. The essence of the method is in the determination of statistical moments of the frequency of the SNe explosions and subsequent determination of distribution function . The program of numerical simulation has been developed for testing the efficiency of the method. Numerical simulations show that even for a small mean number of registered SNe explosions, method allows restoring initial distribution function. The results of numerical simulations are given.