richarddmorey

data(attitude)
o = lm(rating ~ complaints, data = attitude)
par(mfrow=c(1,2), bg="white")
plot(rating ~ complaints, data = attitude)
abline(o, col = "red")
plot(o, ask = FALSE)

richarddmorey

x = rnorm(1000)
y = rnorm(1000)
plot(x,y,pch=19,col=rgb(0,0,1,.2))

richarddmorey

---
runtime: shiny
---

```{r echo=FALSE, message=FALSE, warning=FALSE}
require(shinyjs)
useShinyjs(rmd=TRUE)

observe({
  updateTextInput(session, 'ID', value = "Updated.")

richarddmorey

## Richard D. Morey
## 3 Nov 2017

# Redirect plots to temp directory
# We just want the variables
wd = getwd()
setwd(tempdir())
source('https://osf.io/x73zq/download?version=1')
setwd(wd)

richarddmorey

## Given a (scaled Cauchy) Bayes factor for the null against the 
## alternative (Rouder et al 2009), yields the t statistic
## that would yield it. The ... arguments are passed to 
## the ttest.tstat function.
bf.inv = Vectorize(function(b10, ...){
  fn = Vectorize(function(t,...){
    BayesFactor::ttest.tstat(t,...)[["bf"]]
  }, "t")
  t0 = optimize(function(t0, ...){

richarddmorey

find.ncp = Vectorize(function(a=.05, b=.2, df=1){
  cr = qchisq(1-a,df)
  q = optimize(function(q){
    (pchisq(cr,df,ncp=q/(1 - q)) - b)^2
  },interval = c(0,1))$minimum
  q / (1 - q)
},"b")

find.eb = Vectorize(function(a = 0.05, b = .2, df = 1){

richarddmorey

## Explore the sampling distribution of eta^2, omega^2-hat, and epsilon^2
## Richard D. Morey
## Sept 30, 2017

## Settings

## Number of participants in group
N = 10
## Number of groups
J = 3

richarddmorey

richarddmorey

---
title: "Normal meta-analysis"
author: "Richard D. Morey"
date: "04/07/2017"
output: html_document
---


```{r}
# data here

richarddmorey

x = scan()
0.78
0.71
0.69
0.71
0.73
0.68
0.69
0.64
0.64