Geogenic radon hazard was modeled based on ordered logistic regression. Our models had a classification success rate of 56%. We used k-medoids clustering to perform automatic grouping of geological units. Terrestrial gamma dose rate was modeled based on support vector machines. We took into account fault line density as predictor for geogenic radon hazard.