richarddmorey

library(tictoc)
df = data.frame(y = rnorm(10000000), x = rbinom(10000000,1,.5))

tic()
z = df[df$x == 1,]$y
toc()

rm(z)

tic()

richarddmorey

N = c(16, 16)
y = c(3, 9)
actual_result = prop.test(y, N, correct = FALSE)


compute_p_value = function(p0){
  ## Probability of all outcomes,
  ## assuming independence
  pr_X2s = 
    outer( dbinom(0:N[1], N[1], p0), 

richarddmorey

data(attitude)

library(cocor)
## EQtest at alpha = .05 needs a 90% CI
cocor(~ rating + complaints | rating + learning, data = attitude, alternative = "less", 
      null.value = .2, test = "zou2007", conf.level = .9)

zou_ci_func = function(conf.level, ...)
{

richarddmorey

## Get file off OSF
tf = tempfile(fileext = ".xlsx")
download.file(url = "https://osf.io/846cb/download", destfile = tf)

## Create summaries
library(tidyverse)

readxl::read_xlsx(tf) %>%
  group_by(STN_NUMBER, DIRECTION) %>%
  summarise(

richarddmorey

---
title: "Joint and conditional, binomial/beta example"
author: "Richard D. Morey"
date: "4 June 2016"
output: html_document
---

```{r,echo=FALSE,warning=FALSE}
library(rgl)
library(knitr)

richarddmorey

---
title             : "The title"
shorttitle        : "Title"

author: 
  - name          : "First Author"
    affiliation   : "1"
    corresponding : yes    # Define only one corresponding author
    address       : "Postal address"
    email         : "my@email.com"

richarddmorey

desired_precision = .6

N1 = 5
N2 = 6

# Function to compute CDFs of "precisions" (CI widths)
SS_func = function(W, n, alpha = .05){
  n * (n - 1) * W^2 / (4 * qt(1-alpha/2, n - 1)^2)
}

richarddmorey

# Prior setup
prior_concentration = 10
# null 
p0 = c(0.3, 0.3, 0.4)


# fake data, null

richarddmorey

richarddmorey

# This file is to demonstrate the output. 
# It is not needed.

testRmd = "https://gist.githubusercontent.com/richarddmorey/c4dd7d1fba2aa055a6f71286baedac25/raw/c256dd9f1ca0b22dfbce12443e1cb98ed5bb9d86/demo.Rmd"

tf = tempfile(fileext = ".Rmd")
tf_out = tempfile(fileext = ".html")
cat(RCurl::getURL(testRmd), file = tf)