Skip to content
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

Questions and advice for GDCRNATools #6

Open
huangwb8 opened this issue Dec 4, 2018 · 2 comments
Open

Questions and advice for GDCRNATools #6

huangwb8 opened this issue Dec 4, 2018 · 2 comments

Comments

@huangwb8
Copy link

huangwb8 commented Dec 4, 2018

When I use it to build ceRNA network,I notice the internal dataset "pcTarget" and "lncTarget" of this package.

QUESTION:
According to your published article,you use 3 data sources including StarBase v2.0 (Li et al.,2014), miRcode (Jeggari et al., 2012), and spongeScan (Furio´ -Tarı´et al., 2016).So I was confused when I found only 137 lncRNAs in lncTarget dataset.Whether you use intersect or union of them?
Because starBase had been updated to v3.0 with 531 lncRNAs with miRNA interaction.I wonder whether to use my new data.Counld you give some suggestion?

ADVICE:
I found an error would appear if some miRNAs in the matrix of users didn't exist in the "pcTarget"/ "lncTarget" or self-defined datasets of users.And I think the next code for mirCorTestFun function would avoid this problems.Your suggestion?

mirCorTestFun <- function(lncDa,pcDa, mir,mir.expr) { if(mir %in% rownames(mir.expr)){ ## miRNA in the matrix(your raw codes) mirDa <- unlist(mir.expr[mir,]) corlm <- cor.test(lncDa, mirDa, alternative='less') corpm <- cor.test(pcDa, mirDa, alternative='less') reglm <- corlm$estimate regpm <- corpm$estimate } else { ## miRNA not exist,so correlation could not be calculated. reglm <- NA regpm <- NA } return (c(reglm, regpm)) ## lnc then pc }

Look forward to your return.Thank you!

@huangwb8
Copy link
Author

huangwb8 commented Dec 4, 2018

Moreover,for the result from gdcCEAnalysis function,could you give a suggestion about proper cut-off of parameters for a real ceRNA network?

@isha-ashi
Copy link

Hi,
I am trying to use gdcCEAnalysis function and it gives two main errors:
1.) Error in colnames<-(*tmp*, value = c("lncRNAs", "Genes", "Counts", :
attempt to set 'colnames' on an object with less than two dimensions

2.)Error in cor.test.default(lncDa, mirDa, alternative = "less") :
not enough finite observations.

I believe huangwb8 has suggested some changes regarding same. However, I am new to R and it would really help if you can help me resolve the issue.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants