[R] Best way to test for numeric digits?

Leonard Mada |eo@m@d@ @end|ng |rom @yon|c@eu
Wed Oct 18 20:35:38 CEST 2023


Dear Rui,

On 10/18/2023 8:45 PM, Rui Barradas wrote:
> split_chem_elements <- function(x, rm.digits = TRUE) {
>   regex <- "(?<=[A-Z])(?![a-z]|$)|(?<=.)(?=[A-Z])|(?<=[a-z])(?=[^a-z])"
>   if(rm.digits) {
>     stringr::str_replace_all(mol, regex, "#") |>
>       strsplit("#|[[:digit:]]") |>
>       lapply(\(x) x[nchar(x) > 0L])
>   } else {
>     strsplit(x, regex, perl = TRUE)
>   }
> }
>
> split.symbol.character = function(x, rm.digits = TRUE) {
>   # Perl is partly broken in R 4.3, but this works:
>   regex <- "(?<=[A-Z])(?![a-z]|$)|(?<=.)(?=[A-Z])|(?<=[a-z])(?=[^a-z])"
>   s <- strsplit(x, regex, perl = TRUE)
>   if(rm.digits) {
>     s <- lapply(s, \(x) x[grep("[[:digit:]]+", x, invert = TRUE)])
>   }
>   s
> }

You have a glitch (mol is hardcoded) in the code of the first function. 
The times are similar, after correcting for that glitch.

Note:
- grep("[[:digit:]]", ...) behaves almost twice as slow as grep("[0-9]", 
...)!
- corrected results below;

Sincerely,

Leonard
#######

split_chem_elements <- function(x, rm.digits = TRUE) {
   regex <- "(?<=[A-Z])(?![a-z]|$)|(?<=.)(?=[A-Z])|(?<=[a-z])(?=[^a-z])"
   if(rm.digits) {
     stringr::str_replace_all(x, regex, "#") |>
       strsplit("#|[[:digit:]]") |>
       lapply(\(x) x[nchar(x) > 0L])
   } else {
     strsplit(x, regex, perl = TRUE)
   }
}

split.symbol.character = function(x, rm.digits = TRUE) {
   # Perl is partly broken in R 4.3, but this works:
   regex <- "(?<=[A-Z])(?![a-z]|$)|(?<=.)(?=[A-Z])|(?<=[a-z])(?=[^a-z])"
   s <- strsplit(x, regex, perl = TRUE)
   if(rm.digits) {
     s <- lapply(s, \(x) x[grep("[0-9]", x, invert = TRUE)])
   }
   s
}

mol <- c("CCl3F", "Li4Al4H16", "CCl2CO2AlPO4SiO4Cl")
mol10000 <- rep(mol, 10000)

system.time(
   split_chem_elements(mol10000)
)
#   user  system elapsed
#   0.58    0.00    0.58

system.time(
   split.symbol.character(mol10000)
)
#   user  system elapsed
#   0.67    0.00    0.67



More information about the R-help mailing list