Question: CustomproDB JunctionType function
0
gravatar for georges.bedran
3.3 years ago by
georges.bedran0 wrote:

hello,

i'm currently working on a proteogenomics project, and i am using customProDB to generate a custom DB, it's a great R package and i find it really helpful.

During the process and while trying to generate my novel junctions peptides fasta, I believe that I ran into a little bug.

at line 131 in the function JunctionType.R "if(length(index_NA) >= 0) matchid_new <- matchid[-index_NA, ]" the if condition will have a true value either if (length(index_NA) is equal to or greater than 0, in the case when (length(index_NA)  is equal to 0, R is returning a null  value for matchid_new, thus returning a wrong list of jun_type.

for example, in my case it was giving me a list full of junctions that connect 2 known exons (97000 junctions that connect 2 known exons).

and after changing the code to :

if(length(index_NA) > 0) matchid_new <- matchid[-index_NA, ]

else {matchid_new <- matchid}

i got a list of 95000 known junctions and 2000 new alternative splicing junctions.

i am using R version 3.3.0 and customProDB version 1.10.0

thank you.

customprodb • 457 views
ADD COMMENTlink modified 3.3 years ago by xiaojing.wang50 • written 3.3 years ago by georges.bedran0
Answer: CustomproDB JunctionType function
0
gravatar for xiaojing.wang
3.3 years ago by
United States
xiaojing.wang50 wrote:
Thanks for pointing this out. It must be a typo. I’ll update the package soon. Xiaojing From: georges.bedran [bioc] [mailto:noreply@bioconductor.org] Sent: Thursday, June 30, 2016 5:51 AM To: Wang, Xiaojing <xiaojing.wang@vanderbilt.edu> Subject: [bioc] CustomproDB JunctionType function Activity on a post you are following on support.bioconductor.org<https: support.bioconductor.org=""> User georges.bedran<https: support.bioconductor.org="" u="" 11010=""/> wrote Question: CustomproDB JunctionType function<https: support.bioconductor.org="" p="" 84473=""/>: hello, i'm currently working on a proteogenomics project, and i am using customProDB to generate a custom DB, it's a great R package and i find it really helpful. During the process and while trying to generate my novel junctions peptides fasta, I believe that I ran into a little bug. at line 131 in the function JunctionType.R "if(length(index_NA) >= 0) matchid_new <- matchid[-index_NA, ]" the if condition will have a true value either if (length(index_NA) is equal to or greater than 0, in the case when (length(index_NA) is equal to 0, R is returning a null value for matchid_new, thus returning a wrong list of jun_type. for example, in my case it was giving me a list full of junctions that connect 2 known exons (97000 junctions that connect 2 known exons). and after changing the code to : if(length(index_NA) > 0) matchid_new <- matchid[-index_NA, ] else {matchid_new <- matchid} i got a list of 95000 known junctions and 2000 new alternative splicing junctions. i am using R version 3.3.0 and customProDB version 1.10.0 thank you. ________________________________ Post tags: customprodb You may reply via email or visit CustomproDB JunctionType function
ADD COMMENTlink written 3.3 years ago by xiaojing.wang50
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 284 users visited in the last hour