-
Notifications
You must be signed in to change notification settings - Fork 22
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Merging qname from BAM file to readId.x #435
Comments
Hi there, Could you paste the head() of the two tables you are trying to merge so that I can get a better idea of the problem. Kind Regards, |
Hi, Thank you for your response! The tables that I am trying to merge are the unique_read_info which contains information from the bam file and the merged_data. The merged_data contains metadata(se.multiSample)$readToTranscriptMaps[[1]] which I had merged with the fullLengthCounts. Here are the head() of the two tables. The qname and the readID.x are not matching up when I run se.multiSample, but when I run each sample individually and not together, it does match up. If it is helpful, here is the code I ran for the se.multiSample: |
Hi, Thanks for sharing these. It looks like you are running bambu correctly as I do not see any issues there. I just have a few more questions. Is unique_read_info from 1 bam file (ie. RNA_cell_naive1.bam)? And does that match with the first file name in this vector When you try merge unique_read_info and merged data what is the output? Is it an error or perhaps an empty table? In the examples you show, could you find if How did you merge the full length counts and the read to transcript map, they do not have any columns that directly key into each other. There are ways around this but it is a bit messy as reads can have multiple equal matches? What is your final goal output that you want to produce/question you want to ask? Maybe I can suggest an alternative way of producing it, if I know how. Kind Regards, |
Hi! After looking more thoroughly into the code, I was able to resolve the merging issue! Essentially, I merged the metadata with the rowData of the multisample. This provided me with the columns needed to merge with the full length counts. Because of the merging of different data frames, it did get complicated and I had to be more careful with the process. Thank you so much for your support! |
Hi!
I ran bambu analysis using se.multiSample (I have 3 replicates) and was having trouble merging the readIDs from the bam files to the output of se.multiSample. I have been trying to merge by the "readId.x" from the metadata(se.multiSample)$readToTranscriptMaps[[1]] with the "qname" from the bam file. Any suggestions or thoughts would be really helpful, thank you!!
The text was updated successfully, but these errors were encountered: