getting average score of a GRanges over another GRanges object. summarizeScores(loc.gr,score.gr,score.col)

summarizeScoresOnRegions.R

.libPaths("/work2/gschub/altuna/RlibsDev")
require(data.table)

summarizeScores<-function(loc.gr,score.gr,score.col){

  ov <- as.matrix(findOverlaps(loc.gr,score.gr))
  
  dt=data.table(id=ov[,1],score=values(score.gr)[ov[,2],which(names(values(score.gr))==score.col)])
  dt=dt[,list(av=mean(score)),by=id]
  
  res=loc.gr[dt$id,]
  mcols(res)=DataFrame(mcols(res),score=dt$av)
  res
}