Enhancer_TCGA_ATAC_1000.bed
图片.png
# -*- coding: utf-8 -*-
"""
Created on Tue Apr 16 19:38:32 2019
@author: Administrator
"""
import pandas as pd
f=open('E:\\qqdoc\\2986446306\\FileRecv\\Enhancer_TCGA_ATAC_1000.bed')
x=pd.read_table(f,header=None)
result=x.groupby(by=x[3]).size()
result.to_csv("E:\\qqdoc\\2986446306\\FileRecv\\count.csv")
Result
图片.png
网友评论