rasmusab

library(ggplot2)
library(plyr)

n_games <- readRDS("n_games.Rdata")
n_positions <- readRDS("n_positions.Rdata")
fullmoves_white <- readRDS("fullmoves_white.Rdata")

fits <- readRDS("chess_fits.Rdata")
fit <- fits$fit
fit <- fit 

rasmusab

# From https://aschinchon.wordpress.com/2015/06/23/the-2d-harmonograph-in-shiny/
# All code by Antonio S. Chinchón (@aschinchon), just rearranged so that you can directly
# run it by copy and pasting it into an R console.

library(shiny)
CreateDS = function () {
  f=jitter(sample(c(2,3),4, replace = TRUE))
  d=runif(4,0,1e-02)
  p=runif(4,0,pi)
  xt = function(t) exp(-d[1]*t)*sin(t*f[1]+p[1])+exp(-d[2]*t)*sin(t*f[2]+p[2])

rasmusab

---
title: "Easy Bayesian Bootstrap in R"
author: "Rasmus"
date: "04/06/2015"
output:
  html_document: 
    keep_md: yes
    self_contained: no
---

rasmusab

0   CLS
10  REM SETTING UP THE VARIABLES AND DATA
20 
30  REM NO OF PUPS PER DEN FROM A SAMPLE OF 16 WOLF DENS
40  DATA 5, 8, 7, 5, 3, 4, 3, 9, 5, 8, 5, 6, 5, 6, 4, 7 
50 
60  NDATA = 16
70  NSAMPLES = 1000
75  NBURNIN = 250
80  NPARAMS = 2

rasmusab

# A Bayesian model that calculates a probability that a couple is fertile 
# and pregnant. Please use this for fun only, not for any serious purpose 
# like *actually* trying to figure out whether you are pregnant.

# Enter your own period onsets here:
period_onset <- as.Date(c("2014-07-02", "2014-08-02", "2014-08-29", "2014-09-25", 
                          "2014-10-24", "2014-11-20", "2014-12-22", "2015-01-19"))
# If you have no dates you can just set days_between_periods to c() instead like:
# days_between_periods <- c()
days_between_periods <- as.numeric(diff(period_onset))

rasmusab

---
title: "A Bayesian Model to Calculate Whether My Wife is Pregnant or Not"
author: "Rasmus Bååth"
output: html_document
---

```{r echo = FALSE}
library(knitr)
options(digits=2)
```

rasmusab

### An example of building a TensorFlow model from R using rPython ###
# For this script you need to 
# 1. Have python 2.7 installed.
# 2. Install the rPython package in R.
# 3. Install Google's TensorFlow library as per these instructions:
#     http://www.tensorflow.org/get_started/os_setup.md#binary_installation

### Here is how to setup and run a trivial TensorFlow model ###

# Load TensorFlow (I couldn't get this to work without setting sys.argv... )

rasmusab

# Prior to the tutorial make sure that the script below runs without error on your R installation.
# What you need is a working installation of Stan: http://mc-stan.org/ .
# For installation instructions, see here: 
# https://github.com/stan-dev/rstan/wiki/RStan-Getting-Started

# After installation you should be able to run this script which should output
# some summary statistics and some pretty plots, :)

# Generating some fake data
set.seed(123)

rasmusab

# Prior to the tutorial make sure that the script below runs without error on your python installation.
# What you need is a working installation of Stan: http://mc-stan.org/ .
# For installation instructions, see here: 
# http://mc-stan.org/interfaces/pystan.html

# After installation you should be able to run this script which should output
# some summary statistics and some pretty plots, :)

# Fitting a simple binomial model using Stan

rasmusab

# fuzzycolor() takes a vector of color names and matches it against the
# xkcdcolors list of color names using edit distance. fuzzycolor() always
# returns a vector of hex color strings, perhaps the ones you wanted...
fuzzycolor <- function(color_names) {
  library(xkcdcolors)
  names_distance <- adist(color_names, xcolors(), ignore.case = TRUE, partial = TRUE)
  xkcd_colors <- xcolors()[ apply(names_distance, 1, which.min) ]
  hex_colors <- name2color(xkcd_colors)
  names(hex_colors) <- xkcd_colors
  hex_colors