Skip to content

Commit

Permalink
Browse files Browse the repository at this point in the history
  • Loading branch information
mbcole committed Jun 20, 2016
2 parents fb508fa + de6a40e commit 1b1dfdd
Show file tree
Hide file tree
Showing 2 changed files with 256 additions and 3 deletions.
7 changes: 4 additions & 3 deletions vignettes/slingshot.Rmd
Original file line number Diff line number Diff line change
Expand Up @@ -151,15 +151,16 @@ We can then pick out the top genes for each lineage and visualize their expressi
topgenes1 <- names(sort(gam.pval[[1]], decreasing = FALSE))[1:100]
heatdata1 <- vargenes[rownames(vargenes) %in% topgenes1, order(crv[[1]]$pseudotime, na.last = NA)]
heatclus1 <- clus[order(crv[[1]]$pseudotime, na.last = NA)]
plotHeatmap(heatdata1, clusterSamples=FALSE, sampleData = data.frame(Cluster = heatclus1))
ce1 <- clusterExperiment(heatdata1, heatclus1, transformation=identity)
plotHeatmap(ce1, clusterSamplesData="orderSamplesValue")
topgenes2 <- names(sort(gam.pval[[2]], decreasing = FALSE))[1:100]
heatdata2 <- vargenes[rownames(vargenes) %in% topgenes2, order(crv[[2]]$pseudotime, na.last = NA)]
heatclus2 <- clus[order(crv[[2]]$pseudotime, na.last = NA)]
plotHeatmap(heatdata2, clusterSamples=FALSE, sampleData = data.frame(Cluster = heatclus2))
ce2 <- clusterExperiment(heatdata2, heatclus2, transformation=identity)
plotHeatmap(ce2, clusterSamplesData="orderSamplesValue")
```


# Session Info

```{r session}
Expand Down
Loading

0 comments on commit 1b1dfdd

Please sign in to comment.