Earthquakes are among the most destructive natural disasters, which cause damage on a huge scale not only to life but to property also. Radon is one of the precursory phenomena that exist in connection to the occurrence of earthquakes and may have potential in forecasting these hazards. The data set in this study contains the observations of radon from August 01, 2014, to January 31, 2015 collected in Sobra city, northern Pakistan. Weibull, Gamma, Log-Normal, Log-Logistic, and Pareto probability distribution were fitted over the radon on its original scale and a log scale. Log-logistic best fits the radon on both scales. Log-logistic based Tukey control charts suggest a central limit of 2465.459 Bq/m3 and upper control limit of 12704.080 Bq/m3 on the original radon data, while they suggest a central limit of 7.759, upper control limit of 11.080 and lower control limit of 4.524 on the log scale. The Tukey control charts reveal several anomalies that were compared with earthquake occurrence. There are five earthquakes that occurred during the same period as the radon monitoring program, having magnitude ranges between 4.9 5.5, with the ratio between strain radius and distance to the epicenter greater than or equal to 1. Among several radon anomalies exceeding the upper control limit, the five are identified as earthquake precursors. The results of this study demonstrate that for earthquakes, seismic events show a correlation with increasing concentrations of radon gas before their occurrence.