... | ... |
@@ -16,7 +16,7 @@ output: github_document |
16 | 16 |
| _BioC_ branch | Status | Version | Dependencies | Rank | |
17 | 17 |
|- |- |- |- |- | |
18 | 18 |
| [Release](http://bioconductor.org/packages/release/bioc/html/POMA.html) | [](https://bioconductor.org/checkResults/release/bioc-LATEST/POMA/) | [](https://www.bioconductor.org/packages/POMA) | [](http://bioconductor.org/packages/release/bioc/html/POMA.html#since) | [](https://bioconductor.org/packages/stats/bioc/POMA) | |
19 |
-| [Devel](http://bioconductor.org/packages/devel/bioc/html/POMA.html) | [](https://bioconductor.org/checkResults/devel/bioc-LATEST/POMA/) | [](https://bioconductor.org/packages/devel/bioc/html/POMA.html) | [](http://bioconductor.org/packages/devel/bioc/html/POMA.html#since) | [](https://bioconductor.org/packages/stats/bioc/POMA) | |
|
19 |
+| [Devel](http://bioconductor.org/packages/devel/bioc/html/POMA.html) | [](https://bioconductor.org/checkResults/devel/bioc-LATEST/POMA/) | [](https://bioconductor.org/packages/devel/bioc/html/POMA.html) | [](http://bioconductor.org/packages/devel/bioc/html/POMA.html#since) | [](https://bioconductor.org/packages/stats/bioc/POMA) | |
|
20 | 20 |
|
21 | 21 |
<!-- badges: end --> |
22 | 22 |
|
... | ... |
@@ -17,7 +17,7 @@ v3](https://img.shields.io/badge/License-GPLv3-blue.svg)](https://www.gnu.org/li |
17 | 17 |
| *BioC* branch | Status | Version | Dependencies | Rank | |
18 | 18 |
|-------------------------------------------------------------------------|---------------------------------------------------------------------------------------------------------------------------------------------------------|----------------------------------------------------------------------------------------------------------------------------------------------------|------------------------------------------------------------------------------------------------------------------------------------------------------|------------------------------------------------------------------------------------------------------------------------------| |
19 | 19 |
| [Release](http://bioconductor.org/packages/release/bioc/html/POMA.html) | [](https://bioconductor.org/checkResults/release/bioc-LATEST/POMA/) | [](https://www.bioconductor.org/packages/POMA) | [](http://bioconductor.org/packages/release/bioc/html/POMA.html#since) | [](https://bioconductor.org/packages/stats/bioc/POMA) | |
20 |
-| [Devel](http://bioconductor.org/packages/devel/bioc/html/POMA.html) | [](https://bioconductor.org/checkResults/devel/bioc-LATEST/POMA/) | [](https://bioconductor.org/packages/devel/bioc/html/POMA.html) | [](http://bioconductor.org/packages/devel/bioc/html/POMA.html#since) | [](https://bioconductor.org/packages/stats/bioc/POMA) | |
|
20 |
+| [Devel](http://bioconductor.org/packages/devel/bioc/html/POMA.html) | [](https://bioconductor.org/checkResults/devel/bioc-LATEST/POMA/) | [](https://bioconductor.org/packages/devel/bioc/html/POMA.html) | [](http://bioconductor.org/packages/devel/bioc/html/POMA.html#since) | [](https://bioconductor.org/packages/stats/bioc/POMA) | |
|
21 | 21 |
|
22 | 22 |
<!-- badges: end --> |
23 | 23 |
|
... | ... |
@@ -19,7 +19,7 @@ link-citations: true |
19 | 19 |
|
20 | 20 |
**Compiled date**: `r Sys.Date()` |
21 | 21 |
|
22 |
-**Last edited**: 2023-12-06 |
|
22 |
+**Last edited**: 2023-12-07 |
|
23 | 23 |
|
24 | 24 |
**License**: `r packageDescription("POMA")[["License"]]` |
25 | 25 |
|
... | ... |
@@ -97,37 +97,37 @@ a <- PomaBoxplots(none, |
97 | 97 |
|
98 | 98 |
b <- PomaBoxplots(auto_scaling, |
99 | 99 |
x = "samples", |
100 |
- legend_position = "none") + |
|
100 |
+ theme_params = list(legend_title = FALSE, legend_position = "none")) + |
|
101 | 101 |
ggplot2::ggtitle("Auto Scaling") + |
102 | 102 |
ggplot2::theme(axis.text.x = ggplot2::element_blank()) |
103 | 103 |
|
104 | 104 |
c <- PomaBoxplots(level_scaling, |
105 | 105 |
x = "samples", |
106 |
- legend_position = "none") + |
|
106 |
+ theme_params = list(legend_title = FALSE, legend_position = "none")) + |
|
107 | 107 |
ggplot2::ggtitle("Level Scaling") + |
108 | 108 |
ggplot2::theme(axis.text.x = ggplot2::element_blank()) |
109 | 109 |
|
110 | 110 |
d <- PomaBoxplots(log_scaling, |
111 | 111 |
x = "samples", |
112 |
- legend_position = "none") + |
|
112 |
+ theme_params = list(legend_title = FALSE, legend_position = "none")) + |
|
113 | 113 |
ggplot2::ggtitle("Log Scaling") + |
114 | 114 |
ggplot2::theme(axis.text.x = ggplot2::element_blank()) |
115 | 115 |
|
116 | 116 |
e <- PomaBoxplots(log_transformation, |
117 | 117 |
x = "samples", |
118 |
- legend_position = "none") + |
|
118 |
+ theme_params = list(legend_title = FALSE, legend_position = "none")) + |
|
119 | 119 |
ggplot2::ggtitle("Log Transformation") + |
120 | 120 |
ggplot2::theme(axis.text.x = ggplot2::element_blank()) |
121 | 121 |
|
122 | 122 |
f <- PomaBoxplots(vast_scaling, |
123 | 123 |
x = "samples", |
124 |
- legend_position = "none") + |
|
124 |
+ theme_params = list(legend_title = FALSE, legend_position = "none")) + |
|
125 | 125 |
ggplot2::ggtitle("Vast Scaling") + |
126 | 126 |
ggplot2::theme(axis.text.x = ggplot2::element_blank()) |
127 | 127 |
|
128 | 128 |
g <- PomaBoxplots(log_pareto, |
129 | 129 |
x = "samples", |
130 |
- legend_position = "none") + |
|
130 |
+ theme_params = list(legend_title = FALSE, legend_position = "none")) + |
|
131 | 131 |
ggplot2::ggtitle("Log Pareto") + |
132 | 132 |
ggplot2::theme(axis.text.x = ggplot2::element_blank()) |
133 | 133 |
|
... | ... |
@@ -142,47 +142,47 @@ Here we can evaluate the normalization effects on features. |
142 | 142 |
```{r, message = FALSE, warning = FALSE} |
143 | 143 |
h <- PomaDensity(none, |
144 | 144 |
x = "features", |
145 |
- legend_position = "none") + |
|
145 |
+ theme_params = list(legend_title = FALSE, legend_position = "none")) + |
|
146 | 146 |
ggplot2::ggtitle("Not Normalized") |
147 | 147 |
|
148 | 148 |
i <- PomaDensity(auto_scaling, |
149 | 149 |
x = "features", |
150 |
- legend_position = "none") + |
|
150 |
+ theme_params = list(legend_title = FALSE, legend_position = "none")) + |
|
151 | 151 |
ggplot2::ggtitle("Auto Scaling") + |
152 | 152 |
ggplot2::theme(axis.title.x = ggplot2::element_blank(), |
153 | 153 |
axis.title.y = ggplot2::element_blank()) |
154 | 154 |
|
155 | 155 |
j <- PomaDensity(level_scaling, |
156 | 156 |
x = "features", |
157 |
- legend_position = "none") + |
|
157 |
+ theme_params = list(legend_title = FALSE, legend_position = "none")) + |
|
158 | 158 |
ggplot2::ggtitle("Level Scaling") + |
159 | 159 |
ggplot2::theme(axis.title.x = ggplot2::element_blank(), |
160 | 160 |
axis.title.y = ggplot2::element_blank()) |
161 | 161 |
|
162 | 162 |
k <- PomaDensity(log_scaling, |
163 | 163 |
x = "features", |
164 |
- legend_position = "none") + |
|
164 |
+ theme_params = list(legend_title = FALSE, legend_position = "none")) + |
|
165 | 165 |
ggplot2::ggtitle("Log Scaling") + |
166 | 166 |
ggplot2::theme(axis.title.x = ggplot2::element_blank(), |
167 | 167 |
axis.title.y = ggplot2::element_blank()) |
168 | 168 |
|
169 | 169 |
l <- PomaDensity(log_transformation, |
170 | 170 |
x = "features", |
171 |
- legend_position = "none") + |
|
171 |
+ theme_params = list(legend_title = FALSE, legend_position = "none")) + |
|
172 | 172 |
ggplot2::ggtitle("Log Transformation") + |
173 | 173 |
ggplot2::theme(axis.title.x = ggplot2::element_blank(), |
174 | 174 |
axis.title.y = ggplot2::element_blank()) |
175 | 175 |
|
176 | 176 |
m <- PomaDensity(vast_scaling, |
177 | 177 |
x = "features", |
178 |
- legend_position = "none") + |
|
178 |
+ theme_params = list(legend_title = FALSE, legend_position = "none")) + |
|
179 | 179 |
ggplot2::ggtitle("Vast Scaling") + |
180 | 180 |
ggplot2::theme(axis.title.x = ggplot2::element_blank(), |
181 | 181 |
axis.title.y = ggplot2::element_blank()) |
182 | 182 |
|
183 | 183 |
n <- PomaDensity(log_pareto, |
184 | 184 |
x = "features", |
185 |
- legend_position = "none") + |
|
185 |
+ theme_params = list(legend_title = FALSE, legend_position = "none")) + |
|
186 | 186 |
ggplot2::ggtitle("Log Pareto") + |
187 | 187 |
ggplot2::theme(axis.title.x = ggplot2::element_blank(), |
188 | 188 |
axis.title.y = ggplot2::element_blank()) |
... | ... |
@@ -19,7 +19,7 @@ link-citations: true |
19 | 19 |
|
20 | 20 |
**Compiled date**: `r Sys.Date()` |
21 | 21 |
|
22 |
-**Last edited**: 2023-12-06 |
|
22 |
+**Last edited**: 2023-12-07 |
|
23 | 23 |
|
24 | 24 |
**License**: `r packageDescription("POMA")[["License"]]` |
25 | 25 |
|
... | ... |
@@ -46,13 +46,6 @@ BiocManager::install("POMA") |
46 | 46 |
library(POMA) |
47 | 47 |
``` |
48 | 48 |
|
49 |
-Additionally, you can load other useful packages for this vignette: |
|
50 |
- |
|
51 |
-```{r, warning = FALSE, message = FALSE} |
|
52 |
-library(ggplot2) |
|
53 |
-library(plotly) |
|
54 |
-``` |
|
55 |
- |
|
56 | 49 |
# The POMA Workflow |
57 | 50 |
|
58 | 51 |
The `POMA` package functions are organized into three sequential, distinct blocks: Data Preparation, Pre-processing, and Statistical Analysis. |
... | ... |
@@ -122,30 +115,24 @@ normalized |
122 | 115 |
|
123 | 116 |
```{r, message = FALSE} |
124 | 117 |
PomaBoxplots(imputed, |
125 |
- x = "samples", |
|
126 |
- legend_position = "none") + |
|
127 |
- ggplot2::ggtitle("Not Normalized") # data before normalization |
|
118 |
+ x = "samples") # data before normalization |
|
128 | 119 |
``` |
129 | 120 |
|
130 | 121 |
```{r, message = FALSE} |
131 | 122 |
PomaBoxplots(normalized, |
132 |
- x = "samples", |
|
133 |
- legend_position = "none") + |
|
134 |
- ggplot2::ggtitle("Normalized") # data after normalization |
|
123 |
+ x = "samples") # data after normalization |
|
135 | 124 |
``` |
136 | 125 |
|
137 | 126 |
<!-- On the other hand, `PomaDensity` shows the distribution of all features before and after the normalization process. --> |
138 | 127 |
|
139 | 128 |
```{r, message = FALSE} |
140 | 129 |
PomaDensity(imputed, |
141 |
- x = "features") + |
|
142 |
- ggplot2::ggtitle("Not Normalized") # data before normalization |
|
130 |
+ x = "features") # data before normalization |
|
143 | 131 |
``` |
144 | 132 |
|
145 | 133 |
```{r, message = FALSE} |
146 | 134 |
PomaDensity(normalized, |
147 |
- x = "features") + |
|
148 |
- ggplot2::ggtitle("Normalized") # data after normalization |
|
135 |
+ x = "features") # data after normalization |
|
149 | 136 |
``` |
150 | 137 |
|
151 | 138 |
### Outlier Detection |
... | ... |
@@ -171,19 +158,22 @@ pre_processed |
171 | 158 |
#### T-test |
172 | 159 |
|
173 | 160 |
```{r} |
174 |
-PomaUnivariate(pre_processed, method = "ttest") |
|
161 |
+pre_processed %>% |
|
162 |
+ PomaUnivariate(method = "ttest") |
|
175 | 163 |
``` |
176 | 164 |
|
177 | 165 |
<!-- You can also compute a volcano plot using the T-test results. --> |
178 | 166 |
|
179 | 167 |
```{r} |
180 |
-PomaVolcano(imputed, pval = "adjusted", labels = TRUE) |
|
168 |
+imputed %>% |
|
169 |
+ PomaVolcano(pval = "adjusted", labels = TRUE) |
|
181 | 170 |
``` |
182 | 171 |
|
183 | 172 |
#### Wilcoxon Test |
184 | 173 |
|
185 | 174 |
```{r, warning = FALSE} |
186 |
-PomaUnivariate(pre_processed, method = "mann") |
|
175 |
+pre_processed %>% |
|
176 |
+ PomaUnivariate(method = "mann") |
|
187 | 177 |
``` |
188 | 178 |
|
189 | 179 |
<!-- ### Limma --> |