Search Results

Idiomatic R I did some simple stuff first: pulled all your functions to the start of the script, reformatted some code/comments. … This side-effect is problematic in R, so I didn't include a pop_element function. To achieve pop_element you have to peek and then delete_min. …

You should always strive to find ways to define pure functions in your R code, and try to use vectorisation instead of for loops. …

The code was OK as far as base R goes. First thing I did was run styler and then lintr on your code; these two things help clean up the coding style in your scripts. … I also renamed your alpha / beta vectors - in R, these correspond to the shape and rate parameters that are passed into dgamma: # ---- experimental data num_observations <- 10 lambda <- .2 x <- rpois …

if (r < r_escape) { x1 <- x0^2 - y0^2 + c_x[nx] y1 <- 2 * x0 * y0 + c_y[ny] s <- ifelse(r > 0, s + min(r, 1 / r), s) x0 <- x1 y0 <- y1 r <- x0^2 + y0^2 … (df, x0^2 - y0^2 + cx) y1 <- with(df, 2 * x0 * y0 + cy) df %>% dplyr::mutate( s = ifelse(r > 0, s + pmin(r, 1 / r), s), x0 = x1, y0 = y1, r = x0^2 + y0^2 …

There's a couple of things you can do to improve the R-style and the efficiency of your code. … } x0 <- x1 y0 <- y1 if (r > 0) { s <- s + min(r, 1 / r) } # this is what we measure instead of escape count } Mxy[nx, ny] <- s } } Mmax <- …