tag:blogger.com,1999:blog-6232819486261696035.post5795034677568715970..comments2023-09-25T09:01:44.323-05:00Comments on Getting Genetics Done: Prune GWAS data in RStephen Turnerhttp://www.blogger.com/profile/06656711316726116187noreply@blogger.comBlogger1125tag:blogger.com,1999:blog-6232819486261696035.post-82509708965950390462011-03-29T15:39:43.953-05:002011-03-29T15:39:43.953-05:00Probably faster to pick (one chromosome) like
...Probably faster to pick (one chromosome) like<br /><br /> pickSNPs1chr <- function(pos, dist=100000)<br /> {<br /> res <- logical(length(pos)); i <- 1L; d <- dist<br /> while ((i <- which.min((pos-d)^2)) != length(pos)) {<br /> if (pos[i] < d) i <- i + 1<br /> res[i] <- TRUE<br /> d <- pos[i] + dist<br /> }<br /> which(res)<br /> }Anonymousnoreply@blogger.com