Skip to content

Instantly share code, notes, and snippets.

Avatar

Richard Morey richarddmorey

View GitHub Profile
@richarddmorey
richarddmorey / X2.R
Created Jul 29, 2020
Two proportions, p value
View 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),
@richarddmorey
richarddmorey / eqt_pvalue.R
Last active Jul 17, 2020
Compute p value for equivalence test
View 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, ...)
{
@richarddmorey
richarddmorey / repro_subways.R
Created Jun 21, 2020
Attempt to reproduce Maglio & Polman (2014) Figure 1
View 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(
@richarddmorey
richarddmorey / planning_precision.Rmd
Created Jun 12, 2020
Figures for "Power and precision" blog post
View 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
View 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"
@richarddmorey
richarddmorey / precision.R
Created Aug 21, 2019
CI "precision" reversal
View 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)
}
@richarddmorey
richarddmorey / bf_one_way
Created Jun 23, 2019
Bayes factor for one-way table
View bf_one_way
# Prior setup
prior_concentration = 10
# null
p0 = c(0.3, 0.3, 0.4)
# fake data, null
@richarddmorey
richarddmorey / eurovision2019.csv
Last active May 19, 2019
Eurovision 2016-2018 analysis
View 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
@richarddmorey
richarddmorey / demo.R
Last active Apr 3, 2019
Demonstrate reading bibtex for Rmarkdown from the web (e.g., from github)
View 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)
@richarddmorey
richarddmorey / install_flexjags.R
Last active Mar 19, 2019
Install flexjags in R for Essex Summer School 3F (Bayesian data analysis)
View install_flexjags.R
## First install the devtools package if necessary
devtools::install_github("richarddmorey/EssexSummerSchoolBayes", subdir = "shiny_demo/flexjags/flexjags")
flexjags::run_flexjags()
You can’t perform that action at this time.