在墙内搞生物信息学偶尔会遇到一些莫名其妙的问题,比如突然间有天就没办法使用下面这个函数:
?...怪不得总是有些人问到使用它的各种失败,各种报错,因为我大部分时间都是在墙外所以根本就没办法重复出求助者的错误。...简单看了加源代码,发现是download_KEGG('hsa') 这个函数的问题,这个函数每次使用都会去下载这个数据,如果某些地方某些电脑无法访问这个KEGG的官方网站,那么这个包就没办法使用了。...它描述了由有限个物件中抽出n个物件,成功抽出指定种类的物件的个数(不归还)
gened=intersect(names(GeneID2Path),gene_up) ## all genes been..." "))
return(c(thisP,p,OddsRatio,exp_count,k,n,M,N))
}))
colnames(kegg_r)=c('path_id','p','OddsRatio