Christopher Peters statwonk

Pinned

1. A visualization of the beta distribu...

   1

   library(tidyverse)

   2

   plot.new()

   3

   par(bg = "black")

   4

   plot(density(rbeta(1e6, 50, 50)), col = "pink3", lwd = 1, ylim = c(0, 8),

   5

        xlab = "", ylab = "", main = "Beta distribution", axes = FALSE,

2. A comparison of std.errors of the tr...

   1

   library(tidyverse)

   2

   library(lme4)

   3

   library(broom)

   4

   5

   N <- 1e4

3. Time to 10% ending / converting / fa...

   1

   library(tidyverse)

   2

   library(survival)

   3

   4

   rweibull_cens <- function(n, shape, scale) {

   5

     # http://statwonk.com/weibull.html

4. experimentr

   A place for experiment simulations

   R 1

5. openNHTSA

   An R wrapper for the U.S. Department of Transportation National Highway Traffic Safety Administration's API.

   R 21 3

6. Risk adds up. This code piece answer...

   1

   library(tidyverse)

   2

   expand.grid(

   3

     risk = seq(0.001, 0.02, 0.001),

   4

     units_of_exposure = seq_len(24*7)

   5

   ) %>% as_tibble() %>%