Greetings to all -- I have a set of ~1300 radon measurements from homes across British Columbia in Canada. We have split these measurements into five groups according the the tectonic belt in which they fall, with N ranging from ~160 to ~600 measurements per group. The values are roughly log-normally distributed (1) overall and (2) within each group.
For those areas where values are low I would like to be able to estimate the probability of finding values greater than 200 Bq/m3 (the new safety standard) even though the maximum measured value is 174 Bq/m3. I know this is a question of extreme value distributions, and the reading I have done so far suggests that there are multiple approaches (maximum likelihood, probability-weighted moments,...). I am wondering if someone could point me towards R packages/functions appropriate for these data? I assume that it would evir and I have read through the documentation, but I would appreciate some guidance if someone with more experience is willing to offer it. Many thanks, Sarah ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
