R 如何实现功能？

我有两个数据文件，如下所示：

gr_RNA1 <- GRanges(seqnames=RNA$chr, IRanges(start=(RNA$start)-100, end=(RNA$end)+100, names=RNA$Gene_ID))
overlaps <- findOverlaps(gr_RNA1, gr_SNP)
subsetByOver<-subsetByOverlaps(gr_RNA1, gr_SNP)
match_hit1 <- data.frame(names(gr_RNA1)[queryHits(overlaps)],names(gr_SNP)[subjectHits(overlaps)],stringsAsFactors=F)
dim(match_hit1)
[1] 17976     2

头部（RNA）

总干事（snp）

我想找到snp文件和RNA文件之间的重叠，所以我使用了Genomic R包，并完成了以下命令：

gr_RNA <- GRanges(seqnames=RNA$chr,IRanges(start=RNA$start,end=RNA$end,names=RNA$Gene_ID))
gr_SNP <- GRanges(seqnames=SNP$chr, IRanges(start=SNP$start,end=SNP$end,names=SNP$SNP_No))

overlaps <- findOverlaps(gr_RNA, gr_SNP)
subsetByOver <- subsetByOverlaps(gr_RNA, gr_SNP)
match_hit <- data.frame(names(gr_RNA)[queryHits(overlaps)],names(gr_SNP)[subjectHits(overlaps)],stringsAsFactors=F)
names(match_hit) <- c('Gene_ID','SNP')
head(match_hit)
            Gene_ID        SNP
1  ENSG00000000457.1 SNP_307301
2  ENSG00000000457.2 SNP_307307
3 ENSG00000000457.11 SNP_307365
4 ENSG00000000457.12 SNP_307387
5  ENSG00000000460.1 SNP_306845
6  ENSG00000000460.1 SNP_306846

dim(match_hit)
[1] 12287     2

---

R中的函数定义如下：

myFunction <- function(parameters) {
    #function Code
    return(result)
}

---

myFunction实际上只是把代码放在
函数（RNA=RNA，gr_SNP=gr_SNP，expand=100）{
和
}
之间，并用
expand
替换100。或者您可以将代码复制到rstudio中，并使用魔杖下的
提取功能
按钮
gr_RNA1 <- GRanges(seqnames=RNA$chr, IRanges(start=(RNA$start)-100, end=(RNA$end)+100, names=RNA$Gene_ID))
overlaps <- findOverlaps(gr_RNA1, gr_SNP)
subsetByOver<-subsetByOverlaps(gr_RNA1, gr_SNP)
match_hit1 <- data.frame(names(gr_RNA1)[queryHits(overlaps)],names(gr_SNP)[subjectHits(overlaps)],stringsAsFactors=F)
dim(match_hit1)
[1] 17976     2

myFunction <- function(parameters) {
    #function Code
    return(result)
}