File tree 2 files changed +4
-4
lines changed
2 files changed +4
-4
lines changed Original file line number Diff line number Diff line change @@ -206,7 +206,7 @@ def to_model(bam_path, output):
206
206
ins_r , del_r = modeller .indel_matrix_to_choices (indel_matrix_r , read_length )
207
207
208
208
logger .info ("Calculating insert size distribution" )
209
- # insert_size = int(np.mean(insert_size_dist))
209
+ # insert_size = int(np.mean(insert_size_dist))
210
210
hist_insert_size = modeller .insert_size (template_length_dist , read_length )
211
211
212
212
write_to_file (
Original file line number Diff line number Diff line change @@ -14,8 +14,8 @@ def insert_size(template_length_dist, read_length):
14
14
distributin. Uses 1D kernel density estimation.
15
15
16
16
Args:
17
- insert_size_distribution (list): list of insert sizes from aligned
18
- read pairs
17
+ template_length_dist (list): List of template lengths from bam file.
18
+ read_length (int): The length of the read.
19
19
20
20
Returns:
21
21
1darray: a cumulative density function
@@ -28,7 +28,7 @@ def insert_size(template_length_dist, read_length):
28
28
max_mask = tld < 2000
29
29
tld = tld [max_mask ]
30
30
31
- isd = tld - (2 * read_length ) # convert to insert size
31
+ isd = tld - (2 * read_length ) # convert to insert size
32
32
33
33
kde = stats .gaussian_kde (isd , bw_method = 0.2 / np .std (isd , ddof = 1 ))
34
34
x_grid = np .linspace (min (isd ), max (isd ), 2000 )
You can’t perform that action at this time.
0 commit comments