1
0
mirror of https://github.com/msberends/AMR.git synced 2024-12-27 11:26:11 +01:00
AMR/R/p.symbol.R

61 lines
2.4 KiB
R
Raw Normal View History

2018-07-01 21:40:37 +02:00
# ==================================================================== #
# TITLE #
# Antimicrobial Resistance (AMR) Analysis #
# #
2019-01-02 23:24:07 +01:00
# SOURCE #
# https://gitlab.com/msberends/AMR #
2018-07-01 21:40:37 +02:00
# #
# LICENCE #
2019-01-02 23:24:07 +01:00
# (c) 2019 Berends MS (m.s.berends@umcg.nl), Luz CF (c.f.luz@umcg.nl) #
2018-07-01 21:40:37 +02:00
# #
2019-01-02 23:24:07 +01:00
# This R package is free software; you can freely use and distribute #
# it for both personal and commercial purposes under the terms of the #
# GNU General Public License version 2.0 (GNU GPL-2), as published by #
# the Free Software Foundation. #
# #
# This R package was created for academic research and was publicly #
# released in the hope that it will be useful, but it comes WITHOUT #
# ANY WARRANTY OR LIABILITY. #
2019-04-05 18:47:39 +02:00
# Visit our website for more info: https://msberends.gitlab.io/AMR. #
2018-07-01 21:40:37 +02:00
# ==================================================================== #
#' Symbol of a p value
#'
#' Return the symbol related to the p value: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1. Values above \code{p = 1} will return \code{NA}.
2018-07-01 21:40:37 +02:00
#' @param p p value
#' @param emptychar text to show when \code{p > 0.1}
#' @return Text
2019-01-02 23:24:07 +01:00
#' @inheritSection AMR Read more on our website!
2018-07-01 21:40:37 +02:00
#' @export
p.symbol <- function(p, emptychar = " ") {
setting.bak <- options()$scipen
2018-07-01 21:40:37 +02:00
options(scipen = 999)
s <- vector(mode = "character", length = length(p))
2018-07-01 21:40:37 +02:00
for (i in 1:length(p)) {
if (is.na(p[i])) {
s[i] <- NA_character_
2018-07-01 21:40:37 +02:00
next
}
if (p[i] > 1) {
s[i] <- NA_character_
2018-07-01 21:40:37 +02:00
next
} else {
p_test <- p[i]
}
if (p_test > 0.1) {
s[i] <- emptychar
} else if (p_test > 0.05) {
s[i] <- '.'
} else if (p_test > 0.01) {
s[i] <- '*'
} else if (p_test > 0.001) {
s[i] <- '**'
} else if (p_test >= 0) {
s[i] <- '***'
}
}
options(scipen = setting.bak)
2018-07-01 21:40:37 +02:00
s
}