aammd

m <- 0.086

o <- (1-m)/m
curve((sqrt(o) - sqrt(o^-1)) / sqrt(x^2 - x + 1), xlim = c(0,50))
abline(h = 1/3)

aammd

---
title: "Exploring observations"
author: "Andrew, based on Tim's lecture!"
date: "18/05/2020"
output: html_document
---

```{r setup, include=FALSE}
knitr::opts_chunk$set(echo = TRUE)
```

aammd

library(tidyverse)
#> Warning: package 'tibble' was built under R version 4.0.2
#> Warning: package 'tidyr' was built under R version 4.0.2
#> Warning: package 'dplyr' was built under R version 4.0.2
tibble(x = 0:42,
       y1 = pbinom(x, 42, 0.3),
       y2 = 1-pbeta(0.3, x+1, 42-x)) %>% 
  ggplot(aes(x = x, y = y1)) + 
 geom_point(col = "green", size = 3) + 

aammd

library(shiny)
library(bslib)
library(thematic)
library(ggplot2)


thematic::thematic_shiny(font = "auto")


# take a vector, make a histogram