help for fine mappting

Hi

Did you test my suggestions? If not, why not? If yes, in what respect they did not work?

sel <- map$POS >= ref$POS_start[1] & map$POS < ref$POS_end[1]
result1 <- sum( map$post_prob[sel] )
result2 <- length( map$post_prob[sel] )
result3 <- min( map$p[sel] )

should give you desired values. It is up to you how do you want to organise them, as from your examples I do not have faintest idea what you want to do.

And keep your responds to r help list, I cc?d it.

Cheers
Petr

From: greg holly [mailto:mak.hholly at gmail.com]
Sent: Thursday, June 16, 2016 3:06 PM
To: PIKAL Petr <petr.pikal at precheza.cz>
Subject: Re: [R] help for fine mappting

Hi PIKAL;

Thanks so much your writing. I am sorry if I could not explain precisely. All information in ref file are exist in map file. So they are in common. Ref file has about 560 and map file has 27 million rows.That is CHR column common in both and all value given ref$POS_start & ref$POS_end columns  are exist in map$POS.

Thanks in advance,

Greg

Hi