Skip to main content

Return the range of a dataframe not within a range of another dataframe

I have the following main dataframe:

first.seqnames  first.start first.end   first.width first.strand    second.seqnames second.start    second.end  second.width    second.strand
0   chr1    346212  346975  7   *   chr1    10882   10888   7   *
1   chr1    3135476 3136100 2   *   chr1    10890   10891   2   *
2   chr1    11473   11484   12  *   chr1    10893   10904   12  *
3   chr1    5388140 5388730 2   *   chr1    11096   11097   2   *
4   chr1    346213  346984  68  *   chr1    11202   11269   68  *

I want to return the rows of above dataframe that don't exist within the range of the following dataframe:

first.seqnames  first.start first.end   3   4   5
3503    chr1    346213  346984  .   0   .
3504    chr1    3135466 3136202 .   0   .
3505    chr1    3190760 3191377 .   0   .
3506    chr1    3354604 3355258 .   0   .
3507    chr1    5388136 5388749 .   0   .

Here, the first dataframe's 'first.start' and 'first.end' should not exist within the range(346213, 346984),.......

I've tried the following code which creates memory and time complexity. Even the result isn't accurate. Here, some of the df1 ranges exist exactly between the df2 ranges and some overlaps. In case of overlaps, the ranges can be ignored.

def range_subset(range1, range2):
    """Whether range1 is a subset of range2."""
    if not range1:
        return True  # empty range is subset of anything
    if not range2:
        return False  # non-empty range can't be subset of empty range
    if len(range1) > 1 and range1.step % range2.step:
        return False  # must have a single value or integer multiple step
    return range1.start in range2 and range1[-1] in range2

for a,b in zip(df1['first.start'], df1['first.end']):
    for i,j in zip(df2['first.start'], df2['first.end']):
        if(range_subset(range(a, b), range(i, j)) == True):
            print(a,b)

Output:

first.seqnames  first.start first.end   first.width first.strand    second.seqnames second.start    second.end  second.width    second.strand
0   chr1    346212  346975  7   *   chr1    10882   10888   7   *
2   chr1    11473   11484   12  *   chr1    10893   10904   12  *


source https://stackoverflow.com/questions/71227508/return-the-range-of-a-dataframe-not-within-a-range-of-another-dataframe

Comments

Popular posts from this blog

ValueError: X has 10 features, but LinearRegression is expecting 1 features as input

So, I am trying to predict the model but its throwing error like it has 10 features but it expacts only 1. So I am confused can anyone help me with it? more importantly its not working for me when my friend runs it. It works perfectly fine dose anyone know the reason about it? cv = KFold(n_splits = 10) all_loss = [] for i in range(9): # 1st for loop over polynomial orders poly_order = i X_train = make_polynomial(x, poly_order) loss_at_order = [] # initiate a set to collect loss for CV for train_index, test_index in cv.split(X_train): print('TRAIN:', train_index, 'TEST:', test_index) X_train_cv, X_test_cv = X_train[train_index], X_test[test_index] t_train_cv, t_test_cv = t[train_index], t[test_index] reg.fit(X_train_cv, t_train_cv) loss_at_order.append(np.mean((t_test_cv - reg.predict(X_test_cv))**2)) # collect loss at fold all_loss.append(np.mean(loss_at_order)) # collect loss at order plt.plot(np.log(al...

Sorting large arrays of big numeric stings

I was solving bigSorting() problem from hackerrank: Consider an array of numeric strings where each string is a positive number with anywhere from to digits. Sort the array's elements in non-decreasing, or ascending order of their integer values and return the sorted array. I know it works as follows: def bigSorting(unsorted): return sorted(unsorted, key=int) But I didnt guess this approach earlier. Initially I tried below: def bigSorting(unsorted): int_unsorted = [int(i) for i in unsorted] int_sorted = sorted(int_unsorted) return [str(i) for i in int_sorted] However, for some of the test cases, it was showing time limit exceeded. Why is it so? PS: I dont know exactly what those test cases were as hacker rank does not reveal all test cases. source https://stackoverflow.com/questions/73007397/sorting-large-arrays-of-big-numeric-stings

How to load Javascript with imported modules?

I am trying to import modules from tensorflowjs, and below is my code. test.html <!DOCTYPE html> <html lang="en"> <head> <meta charset="UTF-8"> <title>Document</title </head> <body> <script src="https://cdn.jsdelivr.net/npm/@tensorflow/tfjs@2.0.0/dist/tf.min.js"></script> <script type="module" src="./test.js"></script> </body> </html> test.js import * as tf from "./node_modules/@tensorflow/tfjs"; import {loadGraphModel} from "./node_modules/@tensorflow/tfjs-converter"; const MODEL_URL = './model.json'; const model = await loadGraphModel(MODEL_URL); const cat = document.getElementById('cat'); model.execute(tf.browser.fromPixels(cat)); Besides, I run the server using python -m http.server in my command prompt(Windows 10), and this is the error prompt in the console log of my browser: Failed to loa...