Browse code

Efficiency and user experience improvement & BUG fixes for TSCAN workflow

Yichen Wang authored on 16/06/2022 18:54:41
Showing 1 changed files
... ...
@@ -2,65 +2,62 @@
2 2
 % Please edit documentation in R/runTSCAN.R
3 3
 \name{plotTSCANPseudotimeGenes}
4 4
 \alias{plotTSCANPseudotimeGenes}
5
-\title{Run plotTSCANPseudotimeGenes function to plot genes with significant 
6
-changes}
5
+\title{Plot expression changes of top features along a TSCAN pseudotime path}
7 6
 \usage{
8 7
 plotTSCANPseudotimeGenes(
9 8
   inSCE,
10 9
   pathIndex,
11 10
   direction = c("increasing", "decreasing"),
12
-  n = 10
11
+  topN = 10,
12
+  useAssay = NULL,
13
+  featureDisplay = metadata(inSCE)$featureDisplay
13 14
 )
14 15
 }
15 16
 \arguments{
16 17
 \item{inSCE}{Input \linkS4class{SingleCellExperiment} object.}
17 18
 
18
-\item{pathIndex}{Path index for which the pseudotime values should be used. 
19
-PathIndex corresponds to the terminal node of specific path from the root 
20
-node to the terminal node.}
19
+\item{pathIndex}{Path index for which the pseudotime values should be used.
20
+Should have being used in \code{\link{runTSCANDEG}}.}
21 21
 
22
-\item{direction}{Which direction to use. Choices are increasing or 
23
-decreasing.}
22
+\item{direction}{Should we show features with expression increasing or
23
+decreeasing along the increase in TSCAN pseudotime? Choices are
24
+\code{"increasing"} or \code{"decreasing"}.}
24 25
 
25
-\item{n}{An integer. Only to plot this number of top genes that are 
26
-increasing/decreasing in expression with increasing pseudotime along 
26
+\item{topN}{An integer. Only to plot this number of top genes that are
27
+increasing/decreasing in expression with increasing pseudotime along
27 28
 the path in the MST. Default 10}
29
+
30
+\item{useAssay}{A single character to specify a feature expression matrix in
31
+\code{assays} slot. The expression of top features from here will be
32
+visualized. Default \code{NULL} use the one used for
33
+\code{\link{runTSCANDEG}}.}
34
+
35
+\item{featureDisplay}{Specify the feature ID type to display. Users can set
36
+default value with \code{\link{setSCTKDisplayRow}}. \code{NULL} or
37
+\code{"rownames"} specifies the rownames of \code{inSCE}. Other character
38
+values indicates \code{rowData} variable.}
28 39
 }
29 40
 \value{
30
-A plot with the top genes that increase/decrease in expression with 
31
-increasing pseudotime along the path in the MST
41
+A \code{.ggplot} object with the facets of the top genes. Expression
42
+on y-axis, pseudotime on x-axis.
32 43
 }
33 44
 \description{
34
-A wrapper function which visualizes outputs from the 
45
+A wrapper function which visualizes outputs from the
35 46
 \code{\link{runTSCANDEG}} function. Plots the genes that increase or decrease
36 47
 in expression with increasing pseudotime along the path in the MST.
48
+\code{\link{runTSCANDEG}} has to be run in advance with using the same
49
+\code{pathIndex} of interest.
37 50
 }
38 51
 \examples{
39
-data("scExample", package = "singleCellTK")
40
-sce <- subsetSCECols(sce, colData = "type != 'EmptyDroplet'")
41
-rowData(sce)$Symbol <- rowData(sce)$feature_name
42
-rownames(sce) <- rowData(sce)$Symbol
43
-sce <- runNormalization(sce, 
44
-                        normalizationMethod = "LogNormalize", 
45
-                        useAssay = "counts", 
46
-                        outAssayName = "logcounts")
47
-sce <- runDimReduce(inSCE = sce, 
48
-                    method = "scaterPCA", 
49
-                    useAssay = "logcounts", 
50
-                    reducedDimName = "PCA")
51
-sce <- runDimReduce(inSCE = sce, 
52
-                    method = "rTSNE", 
53
-                    useReducedDim = "PCA", 
54
-                    reducedDimName = "TSNE")
55
-sce <- runTSCAN (inSCE = sce, 
56
-                 useReducedDim = "PCA", 
57
-                 seed = NULL)
58
-terminalNodes <- listTSCANTerminalNodes(sce, analysisName = "Pseudotime")
59
-sce <- runTSCANDEG(inSCE = sce, 
60
-                   pathIndex = terminalNodes[1])
61
-plotTSCANPseudotimeGenes(inSCE = sce, 
62
-                         pathIndex = terminalNodes[1], 
63
-                         direction = "increasing")
52
+data("mouseBrainSubsetSCE", package = "singleCellTK")
53
+mouseBrainSubsetSCE <- runTSCAN(inSCE = mouseBrainSubsetSCE,
54
+                                useReducedDim = "PCA_logcounts")
55
+terminalNodes <- listTSCANTerminalNodes(mouseBrainSubsetSCE)
56
+mouseBrainSubsetSCE <- runTSCANDEG(inSCE = mouseBrainSubsetSCE,
57
+                                   pathIndex = terminalNodes[1])
58
+plotTSCANPseudotimeGenes(mouseBrainSubsetSCE,
59
+                         pathIndex = terminalNodes[1],
60
+                         useAssay = "logcounts")
64 61
 }
65 62
 \author{
66 63
 Nida Pervaiz
Browse code

TSCAN example case bug fixed

nida pervaiz authored on 26/05/2022 14:22:30
Showing 1 changed files
... ...
@@ -15,13 +15,16 @@ plotTSCANPseudotimeGenes(
15 15
 \arguments{
16 16
 \item{inSCE}{Input \linkS4class{SingleCellExperiment} object.}
17 17
 
18
-\item{pathIndex}{Path index for which the pseudotime values should be used. PathIndex 
19
-corresponds to the terminal node of specific path from the root node to the terminal node.}
18
+\item{pathIndex}{Path index for which the pseudotime values should be used. 
19
+PathIndex corresponds to the terminal node of specific path from the root 
20
+node to the terminal node.}
20 21
 
21
-\item{direction}{Which direction to use. Choices are increasing or decreasing.}
22
+\item{direction}{Which direction to use. Choices are increasing or 
23
+decreasing.}
22 24
 
23
-\item{n}{An integer. Only to plot this number of top genes that are increasing/decreasing 
24
-in expression with increasing pseudotime along the path in the MST. Default 10}
25
+\item{n}{An integer. Only to plot this number of top genes that are 
26
+increasing/decreasing in expression with increasing pseudotime along 
27
+the path in the MST. Default 10}
25 28
 }
26 29
 \value{
27 30
 A plot with the top genes that increase/decrease in expression with 
... ...
@@ -37,12 +40,26 @@ data("scExample", package = "singleCellTK")
37 40
 sce <- subsetSCECols(sce, colData = "type != 'EmptyDroplet'")
38 41
 rowData(sce)$Symbol <- rowData(sce)$feature_name
39 42
 rownames(sce) <- rowData(sce)$Symbol
40
-sce <- runNormalization(sce, normalizationMethod = "LogNormalize", useAssay = "counts", outAssayName = "logcounts")
41
-sce <- runDimReduce(inSCE = sce, method = "scaterPCA", useAssay = "logcounts", reducedDimName = "PCA")
42
-sce <- runDimReduce(inSCE = sce, method = "rTSNE", useReducedDim = "PCA", reducedDimName = "TSNE")
43
-sce <- runTSCAN (inSCE = sce, useReducedDim = "PCA", seed = NULL)
44
-sce <- runTSCANDEG(inSCE = sce, pathIndex = 4)
45
-plotTSCANPseudotimeGenes(inSCE = sce, pathIndex = 4, 
43
+sce <- runNormalization(sce, 
44
+                        normalizationMethod = "LogNormalize", 
45
+                        useAssay = "counts", 
46
+                        outAssayName = "logcounts")
47
+sce <- runDimReduce(inSCE = sce, 
48
+                    method = "scaterPCA", 
49
+                    useAssay = "logcounts", 
50
+                    reducedDimName = "PCA")
51
+sce <- runDimReduce(inSCE = sce, 
52
+                    method = "rTSNE", 
53
+                    useReducedDim = "PCA", 
54
+                    reducedDimName = "TSNE")
55
+sce <- runTSCAN (inSCE = sce, 
56
+                 useReducedDim = "PCA", 
57
+                 seed = NULL)
58
+terminalNodes <- listTSCANTerminalNodes(sce, analysisName = "Pseudotime")
59
+sce <- runTSCANDEG(inSCE = sce, 
60
+                   pathIndex = terminalNodes[1])
61
+plotTSCANPseudotimeGenes(inSCE = sce, 
62
+                         pathIndex = terminalNodes[1], 
46 63
                          direction = "increasing")
47 64
 }
48 65
 \author{
Browse code

documentation edits

nida pervaiz authored on 29/04/2022 10:21:57
Showing 1 changed files
... ...
@@ -15,16 +15,13 @@ plotTSCANPseudotimeGenes(
15 15
 \arguments{
16 16
 \item{inSCE}{Input \linkS4class{SingleCellExperiment} object.}
17 17
 
18
-
19 18
 \item{pathIndex}{Path index for which the pseudotime values should be used. PathIndex 
20 19
 corresponds to the terminal node of specific path from the root node to the terminal node.}
21 20
 
22
-\item{direction}{Which direction to use. Choices are increasing or 
23
-decreasing.}
21
+\item{direction}{Which direction to use. Choices are increasing or decreasing.}
24 22
 
25
-\item{n}{An integer. Only to plot this number of top genes that are 
26
-increasing/decreasing in expression with increasing pseudotime along the path
27
-in the MST. Default \code{10}.}
23
+\item{n}{An integer. Only to plot this number of top genes that are increasing/decreasing 
24
+in expression with increasing pseudotime along the path in the MST. Default 10}
28 25
 }
29 26
 \value{
30 27
 A plot with the top genes that increase/decrease in expression with 
Browse code

Merge branch 'devel' into devel

npervaiz authored on 29/04/2022 09:24:35 • GitHub committed on 29/04/2022 09:24:35
Showing 0 changed files
Browse code

documentation of TSCAN added

nida pervaiz authored on 29/04/2022 06:58:27
Showing 1 changed files
... ...
@@ -14,8 +14,8 @@ plotTSCANPseudotimeGenes(
14 14
 \arguments{
15 15
 \item{inSCE}{Input \linkS4class{SingleCellExperiment} object.}
16 16
 
17
-\item{pathIndex}{Path number for which the pseudotime values should be used. PathIndex 
18
-corresponds to one path from the root node to one of the terminal nodes.}
17
+\item{pathIndex}{Path index for which the pseudotime values should be used. PathIndex 
18
+corresponds to the terminal node of specific path from the root node to the terminal node.}
19 19
 
20 20
 \item{direction}{Which direction to use. Choices are increasing or decreasing.}
21 21
 
... ...
@@ -36,7 +36,7 @@ data("scExample", package = "singleCellTK")
36 36
 sce <- subsetSCECols(sce, colData = "type != 'EmptyDroplet'")
37 37
 rowData(sce)$Symbol <- rowData(sce)$feature_name
38 38
 rownames(sce) <- rowData(sce)$Symbol
39
-sce <- scaterlogNormCounts(sce, assayName = "logcounts")
39
+sce <- runNormalization(sce, normalizationMethod = "LogNormalize", useAssay = "counts", outAssayName = "logcounts")
40 40
 sce <- runDimReduce(inSCE = sce, method = "scaterPCA", useAssay = "logcounts", reducedDimName = "PCA")
41 41
 sce <- runDimReduce(inSCE = sce, method = "rTSNE", useReducedDim = "PCA", reducedDimName = "TSNE")
42 42
 sce <- runTSCAN (inSCE = sce, useReducedDim = "PCA", seed = NULL)
Browse code

Resolve R CMD check and BiocCheck issues

Yichen Wang authored on 22/04/2022 04:17:48
Showing 1 changed files
... ...
@@ -2,7 +2,8 @@
2 2
 % Please edit documentation in R/runTSCAN.R
3 3
 \name{plotTSCANPseudotimeGenes}
4 4
 \alias{plotTSCANPseudotimeGenes}
5
-\title{Run plotTSCANPseudotimeGenes function to plot genes with significant changes}
5
+\title{Run plotTSCANPseudotimeGenes function to plot genes with significant 
6
+changes}
6 7
 \usage{
7 8
 plotTSCANPseudotimeGenes(
8 9
   inSCE,
... ...
@@ -14,22 +15,25 @@ plotTSCANPseudotimeGenes(
14 15
 \arguments{
15 16
 \item{inSCE}{Input \linkS4class{SingleCellExperiment} object.}
16 17
 
17
-\item{pathIndex}{Path number for which the pseudotime values should be used. PathIndex 
18
-corresponds to one path from the root node to one of the terminal nodes.}
18
+\item{pathIndex}{Path number for which the pseudotime values should be used. 
19
+PathIndex corresponds to one path from the root node to one of the terminal 
20
+nodes.}
19 21
 
20
-\item{direction}{Which direction to use. Choices are increasing or decreasing.}
22
+\item{direction}{Which direction to use. Choices are increasing or 
23
+decreasing.}
21 24
 
22
-\item{n}{An integer. Only to plot this number of top genes that are increasing/decreasing 
23
-in expression with increasing pseudotime along the path in the MST. Default 10}
25
+\item{n}{An integer. Only to plot this number of top genes that are 
26
+increasing/decreasing in expression with increasing pseudotime along the path
27
+in the MST. Default \code{10}.}
24 28
 }
25 29
 \value{
26
-A plot with the top genes that increase/decrease in expression with increasing 
27
-pseudotime along the path in the MST
30
+A plot with the top genes that increase/decrease in expression with 
31
+increasing pseudotime along the path in the MST
28 32
 }
29 33
 \description{
30
-A wrapper function which visualizes outputs from the runTSCANDEG function. 
31
-Plots the genes that increase or decrease in expression with increasing pseudotime along 
32
-the path in the MST
34
+A wrapper function which visualizes outputs from the 
35
+\code{\link{runTSCANDEG}} function. Plots the genes that increase or decrease
36
+in expression with increasing pseudotime along the path in the MST.
33 37
 }
34 38
 \examples{
35 39
 data("scExample", package = "singleCellTK")
... ...
@@ -37,11 +41,14 @@ sce <- subsetSCECols(sce, colData = "type != 'EmptyDroplet'")
37 41
 rowData(sce)$Symbol <- rowData(sce)$feature_name
38 42
 rownames(sce) <- rowData(sce)$Symbol
39 43
 sce <- scaterlogNormCounts(sce, assayName = "logcounts")
40
-sce <- runDimReduce(inSCE = sce, method = "scaterPCA", useAssay = "logcounts", reducedDimName = "PCA")
41
-sce <- runDimReduce(inSCE = sce, method = "rTSNE", useReducedDim = "PCA", reducedDimName = "TSNE")
44
+sce <- runDimReduce(inSCE = sce, method = "scaterPCA", 
45
+                    useAssay = "logcounts", reducedDimName = "PCA")
46
+sce <- runDimReduce(inSCE = sce, method = "rTSNE", useReducedDim = "PCA", 
47
+                    reducedDimName = "TSNE")
42 48
 sce <- runTSCAN (inSCE = sce, useReducedDim = "PCA", seed = NULL)
43 49
 sce <- runTSCANDEG(inSCE = sce, pathIndex = 4)
44
-plotTSCANPseudotimeGenes(inSCE = sce, pathIndex = 4, direction = "increasing")
50
+plotTSCANPseudotimeGenes(inSCE = sce, pathIndex = 4, 
51
+                         direction = "increasing")
45 52
 }
46 53
 \author{
47 54
 Nida Pervaiz
Browse code

bugs fixed within TSCAN UI

nida pervaiz authored on 20/04/2022 06:47:51
Showing 1 changed files
... ...
@@ -34,14 +34,14 @@ the path in the MST
34 34
 \examples{
35 35
 data("scExample", package = "singleCellTK")
36 36
 sce <- subsetSCECols(sce, colData = "type != 'EmptyDroplet'")
37
-rowData(sce)$Symbol = rowData(sce)$feature_name
38
-rownames(sce) = rowData(sce)$Symbol
39
-sce <- singleCellTK::scaterlogNormCounts(sce, assayName = "logcounts")
40
-sce <- scater::runPCA(sce)
41
-sce <- scater::runTSNE (sce, dimred = "PCA")
37
+rowData(sce)$Symbol <- rowData(sce)$feature_name
38
+rownames(sce) <- rowData(sce)$Symbol
39
+sce <- scaterlogNormCounts(sce, assayName = "logcounts")
40
+sce <- runDimReduce(inSCE = sce, method = "scaterPCA", useAssay = "logcounts", reducedDimName = "PCA")
41
+sce <- runDimReduce(inSCE = sce, method = "rTSNE", useReducedDim = "PCA", reducedDimName = "TSNE")
42 42
 sce <- runTSCAN (inSCE = sce, useReducedDim = "PCA", seed = NULL)
43
-sce <- runTSCANDEG(inSCE = sce, pathIndex = 3)
44
-plotTSCANPseudotimeGenes(inSCE = sce, pathIndex = 3, direction = "increasing")
43
+sce <- runTSCANDEG(inSCE = sce, pathIndex = 4)
44
+plotTSCANPseudotimeGenes(inSCE = sce, pathIndex = 4, direction = "increasing")
45 45
 }
46 46
 \author{
47 47
 Nida Pervaiz
Browse code

bug fixes done

nida pervaiz authored on 12/04/2022 18:12:44
Showing 1 changed files
... ...
@@ -14,17 +14,22 @@ plotTSCANPseudotimeGenes(
14 14
 \arguments{
15 15
 \item{inSCE}{Input \linkS4class{SingleCellExperiment} object.}
16 16
 
17
-\item{pathIndex}{Path number for which the pseudotime values should be used. PathIndex corresponds to one path from the root node to one of the terminal nodes.}
17
+\item{pathIndex}{Path number for which the pseudotime values should be used. PathIndex 
18
+corresponds to one path from the root node to one of the terminal nodes.}
18 19
 
19 20
 \item{direction}{Which direction to use. Choices are increasing or decreasing.}
20 21
 
21
-\item{n}{An integer. Only to plot this number of top genes that are increasing/decreasing in expression with increasing pseudotime along the path in the MST. Default 10}
22
+\item{n}{An integer. Only to plot this number of top genes that are increasing/decreasing 
23
+in expression with increasing pseudotime along the path in the MST. Default 10}
22 24
 }
23 25
 \value{
24
-A plot with the top genes that increase/decrease in expression with increasing pseudotime along the path in the MST
26
+A plot with the top genes that increase/decrease in expression with increasing 
27
+pseudotime along the path in the MST
25 28
 }
26 29
 \description{
27
-A wrapper function which visualizes outputs from the runTSCANDEG function. Plots the genes that increase or decrease in expression with increasing pseudotime along the path in the MST
30
+A wrapper function which visualizes outputs from the runTSCANDEG function. 
31
+Plots the genes that increase or decrease in expression with increasing pseudotime along 
32
+the path in the MST
28 33
 }
29 34
 \examples{
30 35
 data("scExample", package = "singleCellTK")
Browse code

changes in examples

nida pervaiz authored on 12/04/2022 16:50:31
Showing 1 changed files
... ...
@@ -27,9 +27,16 @@ A plot with the top genes that increase/decrease in expression with increasing p
27 27
 A wrapper function which visualizes outputs from the runTSCANDEG function. Plots the genes that increase or decrease in expression with increasing pseudotime along the path in the MST
28 28
 }
29 29
 \examples{
30
+data("scExample", package = "singleCellTK")
31
+sce <- subsetSCECols(sce, colData = "type != 'EmptyDroplet'")
32
+rowData(sce)$Symbol = rowData(sce)$feature_name
33
+rownames(sce) = rowData(sce)$Symbol
34
+sce <- singleCellTK::scaterlogNormCounts(sce, assayName = "logcounts")
35
+sce <- scater::runPCA(sce)
36
+sce <- scater::runTSNE (sce, dimred = "PCA")
30 37
 sce <- runTSCAN (inSCE = sce, useReducedDim = "PCA", seed = NULL)
31
-sce <- runTSCANDEG(inSCE = sce, pathIndex = 6, discardCluster = 8)
32
-plotTSCANPseudotimeGenes(inSCE = sce, pathIndex = 6, direction = "increasing")
38
+sce <- runTSCANDEG(inSCE = sce, pathIndex = 3)
39
+plotTSCANPseudotimeGenes(inSCE = sce, pathIndex = 3, direction = "increasing")
33 40
 }
34 41
 \author{
35 42
 Nida Pervaiz
Browse code

errors fixed

nida pervaiz authored on 12/04/2022 14:02:58
Showing 1 changed files
... ...
@@ -27,9 +27,9 @@ A plot with the top genes that increase/decrease in expression with increasing p
27 27
 A wrapper function which visualizes outputs from the runTSCANDEG function. Plots the genes that increase or decrease in expression with increasing pseudotime along the path in the MST
28 28
 }
29 29
 \examples{
30
-sce <- runTSCAN (sce, useReducedDim = "PCA", seed = NULL)
31
-sce <- runTSCANDEG(sce, pathIndex = 6, discardCluster = 8)
32
-plotTSCANPseudotimeGenes(sce, pathIndex = 6, direction = "increasing")
30
+sce <- runTSCAN (inSCE = sce, useReducedDim = "PCA", seed = NULL)
31
+sce <- runTSCANDEG(inSCE = sce, pathIndex = 6, discardCluster = 8)
32
+plotTSCANPseudotimeGenes(inSCE = sce, pathIndex = 6, direction = "increasing")
33 33
 }
34 34
 \author{
35 35
 Nida Pervaiz
Browse code

changes in TSCAN

nida pervaiz authored on 12/04/2022 12:23:39
Showing 1 changed files
... ...
@@ -27,9 +27,9 @@ A plot with the top genes that increase/decrease in expression with increasing p
27 27
 A wrapper function which visualizes outputs from the runTSCANDEG function. Plots the genes that increase or decrease in expression with increasing pseudotime along the path in the MST
28 28
 }
29 29
 \examples{
30
-sce <- runTSCAN (sce, reducedDimName = "PCA", seed = NULL)
31
-sce <- runTSCANDEG(sce, pathIndex = 1, discardCluster = 8)
32
-plotTSCANPseudotimeGenes(sce, pathIndex = 1, direction = "increasing")
30
+sce <- runTSCAN (sce, useReducedDim = "PCA", seed = NULL)
31
+sce <- runTSCANDEG(sce, pathIndex = 6, discardCluster = 8)
32
+plotTSCANPseudotimeGenes(sce, pathIndex = 6, direction = "increasing")
33 33
 }
34 34
 \author{
35 35
 Nida Pervaiz
Browse code

Trajectory analysis added

nida pervaiz authored on 06/04/2022 16:05:50
Showing 1 changed files
1 1
new file mode 100644
... ...
@@ -0,0 +1,36 @@
1
+% Generated by roxygen2: do not edit by hand
2
+% Please edit documentation in R/runTSCAN.R
3
+\name{plotTSCANPseudotimeGenes}
4
+\alias{plotTSCANPseudotimeGenes}
5
+\title{Run plotTSCANPseudotimeGenes function to plot genes with significant changes}
6
+\usage{
7
+plotTSCANPseudotimeGenes(
8
+  inSCE,
9
+  pathIndex,
10
+  direction = c("increasing", "decreasing"),
11
+  n = 10
12
+)
13
+}
14
+\arguments{
15
+\item{inSCE}{Input \linkS4class{SingleCellExperiment} object.}
16
+
17
+\item{pathIndex}{Path number for which the pseudotime values should be used. PathIndex corresponds to one path from the root node to one of the terminal nodes.}
18
+
19
+\item{direction}{Which direction to use. Choices are increasing or decreasing.}
20
+
21
+\item{n}{An integer. Only to plot this number of top genes that are increasing/decreasing in expression with increasing pseudotime along the path in the MST. Default 10}
22
+}
23
+\value{
24
+A plot with the top genes that increase/decrease in expression with increasing pseudotime along the path in the MST
25
+}
26
+\description{
27
+A wrapper function which visualizes outputs from the runTSCANDEG function. Plots the genes that increase or decrease in expression with increasing pseudotime along the path in the MST
28
+}
29
+\examples{
30
+sce <- runTSCAN (sce, reducedDimName = "PCA", seed = NULL)
31
+sce <- runTSCANDEG(sce, pathIndex = 1, discardCluster = 8)
32
+plotTSCANPseudotimeGenes(sce, pathIndex = 1, direction = "increasing")
33
+}
34
+\author{
35
+Nida Pervaiz
36
+}