# # - eval = FALSE ------------------------------------------------------------ # # 床上文件导入# genelist < - import.bed(~ /数据/ genelists / genes.bed) # # #进口人造石铺地面# genelist < - import.gff(~ /数据/ genelists / genes.gff) # # #出口床上修改后文件#出口。床(genelist”~ /数据/ genelists / filtered_genes.bed ") ## ---- 消息= FALSE --------------------------------------------------------- 库(BRGenomics ) ## ----------------------------------------------------------------------------- 数据(“txs_dm6_chr4”)tx4 < - txs_dm6_chr4 [c(1、10、200、300)]tx4 ## ----------------------------------------------------------------------------- tx4_pr < -促进剂(tx4上游= 50,下游= 100)tx4_pr宽度(tx4_pr ) ## ----------------------------------------------------------------------------- tx4 ## ----------------------------------------------------------------------------- genebodies (tx4开始= 300,结束= -300 ) ## ----------------------------------------------------------------------------- genebodies (tx4, -50、100、修复。结束= "开始 ") ## ----------------------------------------------------------------------------- genebodies (tx4, -100、-50、修复。结束= "开始 ") ## ----------------------------------------------------------------------------- genebodies (tx4, -1000、1000、修复。开始= "结束 ") ## ----------------------------------------------------------------------------- genebodies (tx4 0 10000,修复。开始= "结束 ") ## ----------------------------------------------------------------------------- tx < - txs_dm6_chr4[秩序(txs_dm6_chr4 gene_id美元)]#排序gene_id tx [1:10 ] ## ----------------------------------------------------------------------------- reduceByGene (tx gene_names = tx gene_id美元 ) ## ----------------------------------------------------------------------------- txs_redux < - reduceByGene (tx,gene_names = tx gene_id) txs_redux gene_id < -美元名称(txs_redux)名称(txs_redux) < -零txs_redux ## ----------------------------------------------------------------------------- tx [1:10] txs_insxt < - intersectByGene (tx gene_names = tx gene_id美元)txs_insxt(订单(名称(txs_insxt)))