@@ -54,19 +54,19 @@ def interpolate_map(map_gz, pgen):
54
54
Reading in map file (format has columns [chrom, SNP, cM-pos, bp])
55
55
"""
56
56
if (map_gz [:- 3 ] == ".gz" ) :
57
- maps = pd .read_table (map_gz , header = None , compression = 'gzip' , delim_whitespace = True )
57
+ maps = pd .read_table (map_gz , header = None , compression = 'gzip' , sep = " \\ s+" )
58
58
else :
59
- maps = pd .read_table (map_gz , header = None , delim_whitespace = True )
59
+ maps = pd .read_table (map_gz , header = None , sep = " \\ s+" )
60
60
cm_pos_map = maps [2 ].values .astype (np .float64 )
61
61
phys_pos_map = maps [3 ].values .astype (np .float64 )
62
62
pvar = pgen .replace ("pgen" , "pvar" )
63
63
bim = pgen .replace ("pgen" , "bim" )
64
64
65
65
physical_positions = None
66
66
if os .path .isfile (bim ):
67
- physical_positions = np .array (pd .read_table (bim , delim_whitespace = True , header = None , comment = '#' )[3 ]).astype (np .float64 )
67
+ physical_positions = np .array (pd .read_table (bim , sep = " \\ s+" , header = None , comment = '#' )[3 ]).astype (np .float64 )
68
68
elif os .path .isfile (pvar ):
69
- physical_positions = np .array (pd .read_table (pvar , delim_whitespace = True , header = None , comment = '#' )[1 ]).astype (np .float64 )
69
+ physical_positions = np .array (pd .read_table (pvar , sep = " \\ s+" , header = None , comment = '#' )[1 ]).astype (np .float64 )
70
70
else :
71
71
raise RuntimeError (f"Can't find { bim } or { pvar } " )
72
72
@@ -87,10 +87,10 @@ def get_map_from_bim(pgen, rho):
87
87
cm_out = None
88
88
physical_positions = None
89
89
if os .path .isfile (bim ):
90
- physical_positions = np .array (pd .read_table (bim , delim_whitespace = True , header = None , comment = '#' )[3 ]).astype (int )
90
+ physical_positions = np .array (pd .read_table (bim , sep = " \\ s+" , header = None , comment = '#' )[3 ]).astype (int )
91
91
cm_out = rho * 100 * physical_positions
92
92
elif os .path .isfile (pvar ):
93
- physical_positions = np .array (pd .read_table (pvar , delim_whitespace = True , header = None , comment = '#' )[1 ]).astype (int )
93
+ physical_positions = np .array (pd .read_table (pvar , sep = " \\ s+" , header = None , comment = '#' )[1 ]).astype (int )
94
94
cm_out = rho * 100 * physical_positions
95
95
else :
96
96
raise RuntimeError (f"Can't find { bim } or { pvar } " )
@@ -101,5 +101,5 @@ def get_map_from_bim(pgen, rho):
101
101
return cm_out , physical_positions
102
102
103
103
def parse_demography (demography ):
104
- d = pd .read_table (demography , delim_whitespace = True , header = None )
104
+ d = pd .read_table (demography , sep = " \\ s+" , header = None )
105
105
return list (d [0 ]), list (d [1 ])
0 commit comments