Pokročilé ladění

Paralelní programování v R

Nabeel Imam

Data Scientist

Výchozí chování při chybě

Frustrovaný programátor hledí na obrazovku a shodil všechny papíry ze stolu.

Ukázkový příklad

var_list <- list(1:10,
                 1:10,
                 c("a", "b", "c"))


lapply(var_list, sqrt)

Error in FUN(X[[i]], ...) : non-numeric
argument to mathematical function

Numerický
Numerický
Znakový

Nejednoznačnost

Antropomorfní počítač je zmaten, protože na jeho obrazovce se zobrazuje chyba.

¹ Designed by Freepik

Zachycení chyby

sqrt_custom <- function(var) {

  tryCatch(
    # Expression to evaluate
    sqrt(var),

    # What to do with an error
    error = function (e) return(e)
    )

}

Pracovník zachycuje a opravuje chybu pomocí klíče.

¹ Designed by Freepik

Zachycení chyby

lapply(var_list, sqrt_custom)

[[1]]
 [1] 1.000000 1.414214 1.732051 2.000000 2.236068 2.449490 2.645751 2.828427 3.000000
[10] 3.162278

[[2]]
 [1] 1.000000 1.414214 1.732051 2.000000 2.236068 2.449490 2.645751 2.828427 3.000000
[10] 3.162278

[[3]]
<simpleError in sqrt(var): non-numeric argument to mathematical function>

Zachycení chyb při paralelním zpracování

cl <- makeCluster(3)
parLapply(cl, var_list, sqrt_custom)
stopCluster(cl)

[[1]]
 [1] 1.000000 1.414214 1.732051 2.000000 2.236068 2.449490 2.645751 2.828427 3.000000
[10] 3.162278

[[2]]
 [1] 1.000000 1.414214 1.732051 2.000000 2.236068 2.449490 2.645751 2.828427 3.000000
[10] 3.162278

[[3]]
<simpleError in sqrt(var): non-numeric argument to mathematical function>

Data o porodech

print(ls_births)

$AK
 month  plurality
     1          1
     2          1
 ...
$AL
   month  plurality
      10          1
       9          1
 ...
...

Sumarizační funkce

summarise_births <- function(df) {


  tryCatch({

    df %>% 
      group_by(month) %>% 
      summarise(total = sum(plurality))

  },

    error = function (e) "Error! Check data")


}

Paralelní apply

cl <- makeCluster(4)
clusterEvalQ(cl, library(dplyr))

parLapply(cl, ls_births, summarise_births)

stopCluster(cl)

$AK
    month total
  1     1     7
  2     2    13
  ...

$AL
[1] "Error! Check data"

$AR
    month total
  1     1    28
  2     2    28
  ...
...

Prozkoumání zdroje

head(ls_births[["AL"]], n = 10)

     month plurality
263     10         1
335      9         1
473     12         1
474      6         1
475      9         1
839      9         1
1291    11     Twins
1369     4         1
1609     5         1
1610     5  Triplets

Future map

plan(multisession, workers = 4)
config <- furrr_options(packages = "dplyr")
future_map(ls_births, summarise_births,
           .options = config)
plan(sequential)

$AK
    month total
  1     1     7
  2     2    13
  ...

$AL
[1] "Error! Check data"

$AR
    month total
  1     1    28
  2     2    28
  ...
...

Případ foreach

cl <- makeCluster(4)
registerDoParallel(cl)

foreach(df = ls_births,
        .packages = "dplyr"
       ) %dopar% {
  summarise_births(df)
}

stopCluster(cl)

$AK
    month total
  1     1     7
  2     2    13
  ...

$AL
[1] "Error! Check data"

$AR
    month total
  1     1    28
  2     2    28
  ...
...

Lass uns üben!

Paralelní programování v R