-
Notifications
You must be signed in to change notification settings - Fork 5
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
Odd result in pseudotime #73
Comments
Hi, Is the code for the trajectory calculation exactly as per here? Jack |
I've just re-ran this code on a fresh installation of dyno and everything looked fine
|
@jackbibby1 I ran the code same like the example. One difference between the example data set and mine is that I cannot split the data by anything. It is from a bunch of diseased samples. I know that there is a trajectory as I have checked it with monocle. In your example code when you run the plot command it says "root cell or milestone not provided, trying first outgoing milestone_id Could you let me know how to procees? Thanks! |
If you've ran monocle already then you can split the cells across that trajectory and use it as an input to SCPA -- it doesn't have to be slingshot. Otherwise, this is more likely a question for the Slingshot authors and asking what they would advise |
|
This is likely not the best or most elegant way to do it, but it should work. You can split the trajectory into x number of parts across all values, then use these to annotate your Seurat object or whatever you're using e.g.
... then use the trajectory groups in SCPA |
@jackbibby1 ...I have been trying to dun your pseudotime codes in my own data. The data entered as per the instruction, but when I run the plot command this is the figure that comes
plot_dimred(model_n4,
"pseudotime",
pseudotime = calculate_pseudotime(model_n4),
hex_cells = F,
plot_trajectory = T,
size_cells = 1, alpha_cells = 0.8) +
theme(aspect.ratio = 1)
root cell or milestone not provided, trying first outgoing milestone_id
Using '3' as root
Warning message:
In calculate_pseudotime(model_n4) :
Trajectory is not rooted. Add a root to the trajectory using dynwrap::add_root(). This will result in an error in future releases.
When I added a root cell
root_cell <- "celltype1"
model_n4 <- add_root(model_n4, root_cell)
Error in add_root(model_n4, root_cell) : Invalid root_cell_id
I notice that the "dataset_n4" has cell id as barcodes AAACCAACACAGCCAT-1.
Could you help me fix this? Thanks :)
The text was updated successfully, but these errors were encountered: