Skip to content

Commit

Permalink
up
Browse files Browse the repository at this point in the history
  • Loading branch information
TuomasBorman committed Oct 2, 2023
1 parent aaf0941 commit 9ade5c3
Show file tree
Hide file tree
Showing 2 changed files with 7 additions and 8 deletions.
3 changes: 1 addition & 2 deletions DESCRIPTION
Original file line number Diff line number Diff line change
Expand Up @@ -60,8 +60,7 @@ Suggests:
patchwork,
vegan,
microbiomeDataSets,
bluster,
miaTime
bluster
Roxygen: list(markdown = TRUE)
RoxygenNote: 7.2.3
VignetteBuilder: knitr
12 changes: 6 additions & 6 deletions vignettes/miaViz.Rmd
Original file line number Diff line number Diff line change
Expand Up @@ -329,7 +329,7 @@ plotDMNFit(dmn_se, type = "laplace")

# Serial data ordination and trajectories

```{r, error=FALSE, warning=FALSE, results='hide'}
```{r, eval=FALSE}
if(!requireNamespace("devtools", quietly = TRUE)){
BiocManager::install("devtools")
}
Expand All @@ -341,7 +341,7 @@ if(!requireNamespace("miaTime", quietly = TRUE)){
Principal Coordinates Analysis using Bray-Curtis dissimilarity on the
`hitchip1006` dataset:

```{r MDS, include=FALSE}
```{r MDS, eval=FALSE}
library(miaTime)
data(hitchip1006, package = "miaTime")
tse <- hitchip1006
Expand All @@ -356,7 +356,7 @@ p

Retrieving information about all available trajectories:

```{r}
```{r, eval=FALSE}
library(dplyr)
# List subjects with two time points
Expand All @@ -369,7 +369,7 @@ table(table(tse$subject)) %>% as.data.frame() %>%

Lets look at all trajectories having two time points in the data:

```{r}
```{r, eval=FALSE}
# plot
p + geom_path(aes(x=X1, y=X2, group=subject),
arrow=arrow(length = unit(0.1, "inches")),
Expand All @@ -383,7 +383,7 @@ p + geom_path(aes(x=X1, y=X2, group=subject),

Filtering the two time point trajectories by divergence and displaying top 10%:

```{r}
```{r, eval=FALSE}
library(miaTime)
# calculating step wise divergence based on the microbial profiles
tse <- getStepwiseDivergence(tse, group = "subject", time_field = "time")
Expand Down Expand Up @@ -411,7 +411,7 @@ p + geom_path(aes(x=X1, y=X2,

Plotting an example of the trajectory with the maximum total divergence:

```{r}
```{r, eval=FALSE}
# Get subject with the maximum total divergence
selected.subject <- data.frame(reducedDim(tse), colData(tse)) %>%
group_by(subject) %>%
Expand Down

0 comments on commit 9ade5c3

Please sign in to comment.