Thanks Jeff~~~ In fact I do not know how to combine and extract vectors in R.
ans<-sort(dhyper(x, m, n, k),decreasing=TRUE) rbind(ans,cumsum(ans) will show the first point that exceeds 95% threshold. The problem is: *information is lost* I can no longer identify where are the first few elements from. e.g. for 10 numbers, maybe they are from 4,5,6,7 or for 100 numbers, from 45 to 68 So to append ID's to the data for later retrieval? rbind appears to do the job but not so exactly... -- View this message in context: http://r.789695.n4.nabble.com/Retrieve-95-coverage-of-results-from-a-hypergeometric-distribution-tp4644683p4644715.html Sent from the R help mailing list archive at Nabble.com. ______________________________________________ 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.