dmi3kno

library(magrittr)

logit <- function(u) log(u)-log(1-u)

raise_to_pow <- function(fun, pname_pow=".pow"){
  eval(parse(text=gsub(pattern=".pow", replacement=pname_pow, 
                  "function(u, .pow=1, ...) fun(u,...)^(.pow)")))
}

powlogit <- logit %>% raise_to_pow(pname_pow="e")

dmi3kno

library(rstan)
stanmodelcode <- "
data {
  int<lower=0> N;
  int<lower=0> y[N];
}
parameters {
  real<lower=0.00001> Theta;
}
model {

dmi3kno

library(minisvg)
lightpurple <- "#c0c0f2"
darkpurple <- "#9596ce"
x <- 5
s <- x*10
a <- x*sqrt(2)
chkrs_pat <- stag$pattern(width=2*s, height=2*s, id = 'chkrs_pattern', 
                    patternUnits = 'userSpaceOnUse',
                    stag$rect(x=0, y=0, width =s, height = s, fill=lightpurple),
                    stag$rect(x=s, y=0, width =s, height = s, fill=darkpurple),

dmi3kno

library(tidybayes)
library(ggplot2)
library(hrbrthemes)
m <- 100
s <- 10

dist_df <- tibble::tibble(y=seq.int(0.01,0.99, by=0.05), 
			  arg1_mean=qlnorm(p=y,meanlog=log(m), sdlog=log(s)),
			  arg2_sd=abs(0.5-y+0.1)*400+100)

dmi3kno

is.unset<-function(x)
  x==as.list(alist(a=))[[1]]

is.nameempty <- function(x){
  nx <- names(x)
  if(is.null(nx))
    return(rep(TRUE, length(x)))
  is.na(nx) | nx==""
}

dmi3kno

---
title: 'Tables with references'
author: dmi3kno
date: '2018-01-25'
abstract: |
  My goal was to prepare a minimum reproducible example.
output:
  pdf_document:
    citation_package: natbib
    keep_md: yes

dmi3kno

library(magick)
library(hocr)
library(bunny)
library(bbx)
library(dplyr)

img <- image_read("https://pbs.twimg.com/media/EXKzsiAWoAYH9Pv?format=png")

img_bw <- img %>% 
  image_threshold("white", threshold = "5%") %>% 

dmi3kno

# tips and ideas from
# https://stackoverflow.com/questions/5546346/how-to-place-and-center-text-in-an-svg-rectangle

library(minisvg)

doc <- svg_doc(width = 600, height = 600)

#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
# Centered text in circle
#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

dmi3kno

library(magick)
library(bunny)

# function for calculating mode
mode_n <- function(x, n=1){ 
  res <- with(rle(sort(x)), values[order(lengths, decreasing = TRUE)])
  res[n]
}

img <- image_read("EWOwbehXkAgxeHO.jpeg") %>% 

dmi3kno

library(tidyverse)
library(brms)
library(tidybayes)

mpg_model <- brm(
  bf(hwy ~ 0 + drv), 
  data = mpg,
  prior = c(
    set_prior("normal(20, 3)", class = "b", lb = 0, ub = 40),
    set_prior("student_t(3, 0, 10)", class = "sigma")