Richard Morey richarddmorey

## X2.R
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),

## eqt_pvalue.R

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, ...)
{

## repro_subways.R
## 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(

## planning_precision.Rmd
---
title: "Power and precision"
author: "Richard D. Morey"
date: "11/06/2020"
output:
  html_document:
    dev: png
    self_contained: no
editor_options:
  chunk_output_type: console

## papaja_example.Rmd
---
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"

## precision.R
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)
}

## bf_one_way


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


# fake data, null

## eurovision2019.csv

          
            To.country
            Jury
            Televoting

            
              Netherlands
              231
              261

            
              Italy
              212
              253

            
              Russia
              125
              244

            
              Switzerland
              148
              212

            
              Norway
              47
              291

            
              Sweden
              239
              93

            
              Azerbaijan
              197
              100

            
              North Macedonia
              237
              58

            
              Australia
              154
              131

## demo.R
# 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)


## install_flexjags.R

## First install the devtools package if necessary
devtools::install_github("richarddmorey/EssexSummerSchoolBayes", subdir = "shiny_demo/flexjags/flexjags")

flexjags::run_flexjags()
	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),

	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, ...)
	{
	## 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(
	---
	title: "Power and precision"
	author: "Richard D. Morey"
	date: "11/06/2020"
	output:
	html_document:
	dev: png
	self_contained: no
	editor_options:
	chunk_output_type: console
	---
	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"
	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)
	}



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


	# fake data, null
To.country	Jury	Televoting
Netherlands	231	261
Italy	212	253
Russia	125	244
Switzerland	148	212
Norway	47	291
Sweden	239	93
Azerbaijan	197	100
North Macedonia	237	58
Australia	154	131
	# 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)

	## First install the devtools package if necessary
	devtools::install_github("richarddmorey/EssexSummerSchoolBayes", subdir = "shiny_demo/flexjags/flexjags")

	flexjags::run_flexjags()