search

nghiavtr / FuSeq

GNU General Public License v3.0
32 stars 12 forks source link

erro in FuSeq.R #11

Open reykeryao opened 4 years ago

reykeryao commented 4 years ago

I have two dataset, one runs fine, but the second one gave this error message

Checking possible paralogs... Post processing mapped reads (MR)...'select()' returned 1:many mapping between keys and columns Error in if (sum(res$crt.tx3.gap == 1) == 0) { : missing value where TRUE/FALSE needed Calls: postProcessMappedRead -> refineJunctionBreak In addition: There were 50 or more warnings (use warnings() to see the first 50) Execution halted

unique379r commented 4 years ago

Hi I also got such am error once, i was using k31 for 50nt long reads then i changed it to k21 and ran it again, was fixed. Please check if this is your case too. I am not sure, else author will reply.

reykeryao commented 4 years ago

Thanks for the quick replay, I used k21 for GRCh38 references, I will try using k31 to see if it happens.

Thanks Jun

On Oct 3, 2019, at 2:45 PM, unique379r notifications@github.com wrote:

Hi I also got such a error once, i was using k31 for 50nt long reads then i changed it to k21 and ran it again, was fixed. Please check if this is your case too. I am not sure, else author will reply.

— You are receiving this because you authored the thread. Reply to this email directly, view it on GitHub https://github.com/nghiavtr/FuSeq/issues/11?email_source=notifications&email_token=ABXFCL4MX2K34F5PCTVX3W3QMZDURA5CNFSM4I5HAZB2YY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGOEAJLYPA#issuecomment-538098748, or mute the thread https://github.com/notifications/unsubscribe-auth/ABXFCLYO4XOFJ6SX5KREX6TQMZDURANCNFSM4I5HAZBQ.

reykeryao commented 4 years ago

Don't know why, but after switching k-mer size back to 31 solved the problem. Thanks