Révision b20637ed src/R/nucleominer.R
| b/src/R/nucleominer.R | ||
|---|---|---|
| 767 | 767 |
} |
| 768 | 768 |
non_inter_fuzzy = rec_substract_region(region1, region2) |
| 769 | 769 |
if (is.null(non_inter_fuzzy)) {return(non_inter_fuzzy)}
|
| 770 |
return(collapse_regions(non_inter_fuzzy)) |
|
| 770 |
tmp_ulist = unlist(non_inter_fuzzy) |
|
| 771 |
tmp_names = names(tmp_ulist)[1:4] |
|
| 772 |
non_inter_fuzzy = data.frame(matrix(tmp_ulist, ncol=4, byrow=TRUE), stringsAsFactors=FALSE) |
|
| 773 |
names(non_inter_fuzzy) = tmp_names |
|
| 774 |
non_inter_fuzzy$chr = as.character(non_inter_fuzzy$chr) |
|
| 775 |
non_inter_fuzzy$chr = as.numeric(non_inter_fuzzy$chr) |
|
| 776 |
non_inter_fuzzy$lower_bound = as.numeric(non_inter_fuzzy$lower_bound) |
|
| 777 |
non_inter_fuzzy$upper_bound = as.numeric(non_inter_fuzzy$upper_bound) |
|
| 778 |
non_inter_fuzzy = non_inter_fuzzy[order(non_inter_fuzzy$lower_bound),] |
|
| 779 |
return(non_inter_fuzzy) |
|
| 771 | 780 |
} |
| 772 | 781 |
|
| 773 | 782 |
union_regions = function(# Aggregate regions that intersect themnselves. |
| ... | ... | |
| 841 | 850 |
### Utils to reformat an "apply manipulated" list of regions |
| 842 | 851 |
regions ##< a list of regions |
| 843 | 852 |
) {
|
| 844 |
regions = do.call("rbind", regions)
|
|
| 853 |
regions = do.call(rbind, regions)
|
|
| 845 | 854 |
regions$chr = as.character(regions$chr) |
| 846 | 855 |
regions$chr = as.numeric(regions$chr) |
| 847 | 856 |
regions$lower_bound = as.numeric(regions$lower_bound) |
Formats disponibles : Unified diff