Skip to content

Errors in "delta_auc_in_pseudospace" function. #1

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

Open
SiyiWanggou opened this issue May 4, 2020 · 1 comment
Open

Errors in "delta_auc_in_pseudospace" function. #1

SiyiWanggou opened this issue May 4, 2020 · 1 comment

Comments

@SiyiWanggou
Copy link

Hi Team from cole-trapnell-lab,
I am using your wonderful piplines for trajectory comparison. It is really awesome and I love it so much!
Recently, I meet some errors when I use "delta_auc_in_pseudospace" function in "cds_exprs.auc" calculation. The system reports the Error like this "This tidyselect interface doesn't support predicates yet". I checked the funtion code and just wonder whether the "id" should be "gene_id" here. The original code is listed below:

cds_exprs.auc = cds_exprs %>% dplyr::select(id, Cell.Type,
Pseudotime, expectation) %>% dplyr::mutate(pseudotime_quantile = cut(Pseudotime,
breaks = seq(0, pseudotime_max, pseudotime_max/5))) %>%
unique() %>% dplyr::group_by(id, Cell.Type, pseudotime_quantile) %>%
dplyr::summarize(AUC = MESS::auc(Pseudotime, expectation,
type = "spline")) %>% dplyr::group_by(id, pseudotime_quantile) %>%
dplyr::summarize(auc_difference = (AUC[1] - AUC[2])/(AUC[1] +
AUC[2])) %>% dplyr::arrange(desc(auc_difference))

Thanks,
Siyi

@joselmc
Copy link
Collaborator

joselmc commented May 5, 2020

Hi Siyi, At that point during the function cds_exprs should contain both a gene_id and an id column due to an earlier merge. I assume the error is prior to the dplyr call, likely before the merge step, and probably due to changes or conflicts in a newer version of a package. Could you print your sessionInfo() ?

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