# R code for analyzing the properly designed "experiments" in Li (1998)

#### Mixed-effects logistic regression

library(lme4)

# Claim 01

claim01 = read.table("claim01_log.txt",T)
attach(claim01)
claim01.glmm = lmer(Judgment~You+(1|Speaker),family=binomial,data=claim01)
claim01_log = round(summary(claim01.glmm)@coefs,4)
detach(claim01)

# Claim 02

claim02 = read.table("claim02_log.txt",T)
attach(claim02)
claim02.glmm = lmer(Judgment~HowMany+(1|Speaker),family=binomial,data=claim02)
claim02_log = round(summary(claim02.glmm)@coefs,4)
detach(claim02)

# Claim 03

claim03 = read.table("claim03_log.txt",T)
attach(claim03)
claim03.glmm = lmer(Judgment~You+(1|Speaker),family=binomial,data=claim03)
claim03_log = round(summary(claim03.glmm)@coefs,4)
detach(claim03)

# Claim 04

claim04 = read.table("claim04_log.txt",T)
attach(claim04)
claim04.glmm = lmer(Judgment~You+(1|Speaker),family=binomial,data=claim04)
claim04_log = round(summary(claim04.glmm)@coefs,4)
detach(claim04)

# Claim 05

claim05 = read.table("claim05_log.txt",T)
attach(claim05)
claim05.glmm = lmer(Judgment~Num*Bind+(1|Speaker),family=binomial,data=claim05)
claim05_log = round(summary(claim05.glmm)@coefs,4)
detach(claim05)

# Claim 06

claim06 = read.table("claim06_log.txt",T)
attach(claim06)
claim06.glmm = lmer(Judgment~You*Scope+(1|Speaker),family=binomial,data=claim06)
claim06_log = round(summary(claim06.glmm)@coefs,4)
detach(claim06)

############## Exact within-group tests

# Claim 01

claim01 = read.table("claim01_wg.txt",T)
attach(claim01)
claim01_wg = round(small.exp(Judgment,You,group=Speaker),4)
detach(claim01)

# Claim 02

claim02 = read.table("claim02_wg.txt",T)
attach(claim02)
claim02_wg = round(small.exp(Judgment,HowMany,group=Speaker),4)
detach(claim02)

# Claim 03

claim03 = read.table("claim03_wg.txt",T)
attach(claim03)
claim03_wg = small.exp(Judgment,You,group=Speaker)
detach(claim03)

# Claim 04

claim04 = read.table("claim04_wg.txt",T)
attach(claim04)
claim04_wg = round(small.exp(Judgment,You,group=Speaker),4)
detach(claim04)

# Claim 05

claim05 = read.table("claim05_wg.txt",T)
attach(claim05)
claim05_wg = round(small.exp(Judgment,Num,Bind,group=Speaker),4)
detach(claim05)

# Claim 06

claim06 = read.table("claim06_wg.txt",T)
attach(claim06)
claim06_wg = round(small.exp(Judgment,You,Scope,group=Speaker),4)
detach(claim06)

############## Exact between-group tests

# Claim 01

claim01 = read.table("claim01_bg.txt",T)
attach(claim01)
claim01_bg = round(small.exp(Judgment,You),4)
detach(claim01)

# Claim 02

claim02 = read.table("claim02_bg.txt",T)
attach(claim02)
claim02_bg = round(small.exp(Judgment,HowMany),4)
detach(claim02)

# Claim 03

claim03 = read.table("claim03_bg.txt",T)
attach(claim03)
claim03_bg = round(small.exp(Judgment,You),4)
detach(claim03)

# Claim 04

claim04 = read.table("claim04_bg.txt",T)
attach(claim04)
claim04_bg = round(small.exp(Judgment,You),4)
detach(claim04)

# Claim 05

claim05 = read.table("claim05_bg.txt",T)
attach(claim05)
claim05_bg = round(small.exp(Judgment,Num,Bind),4)
detach(claim05)

# Claim 06

claim06 = read.table("claim06_bg.txt",T)
attach(claim06)
claim06_bg = round(small.exp(Judgment,You,Scope),4)
detach(claim06)

############# Results

claim01_log
claim01_wg
claim01_bg

claim02_log
claim02_wg
claim02_bg

claim03_log
claim03_wg
claim03_bg

claim04_log
claim04_wg
claim04_bg

claim05_log
claim05_wg
claim05_bg

claim06_log
claim06_wg
claim06_bg