overview of experiment

To be filled in.

methods

To be filled in.

bias

I considered females to be exhibited bias, or otherwise not particpating in my experiment if they did one of two things: if they spent more than 75% of 10 minutes where male stimuli were present on one side of the tank, and if they spent more than 50% of the time in the neutral zone of the tank during that same 10 minute period.

reading in the data

Doing some basic checks for side bias, etc.

There are a couple different ways I could do this. For now, I’ll just use the .txt file saved by the tracking script.

I need to go back and do this more precisely using the actual tracking data.

rm(list = ls())
setwd("~/Documents/transitivities/")
source("~/Documents/graphics/plotting_functions.R")
## 
## Attaching package: 'dplyr'
## 
## The following objects are masked from 'package:stats':
## 
##     filter, lag
## 
## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union
files <- list.files(path="~/Desktop/transitivity_trials/",pattern = "[le].txt",full.names = T)

files[1:5]
## [1] "/Users/lukereding/Desktop/transitivity_trials//Adeline_Intermediate_vs_Large.avi_remade.txt"
## [2] "/Users/lukereding/Desktop/transitivity_trials//Adeline_Intermediate_vs_Large.txt"           
## [3] "/Users/lukereding/Desktop/transitivity_trials//Adeline_large_vs_small.avi_remade.txt"       
## [4] "/Users/lukereding/Desktop/transitivity_trials//Adeline_large_vs_small.txt"                  
## [5] "/Users/lukereding/Desktop/transitivity_trials//Adeline_Small_vs_Intermediate.avi_remade.txt"
# initialize vectors
name = type_of_trial=background1_left=background1_neutral=background1_right=background2_left=background2_neutral=background2_right=males1_left=males1_neutral=males1_right=males2_left=males2_neutral=males2_right <- NULL

# get rid of files that have "remade" in them; this is temporary
for(i in 1:length(files)){
  if(length(grep("remade",files[i])) > 0){
    files <- files[-i]
  }
}

for(i in 1:length(files)){
  
  x <- read.table(files[i]) %>% .[,1] %>% as.vector %>% factor
  pathname = strsplit(files[i],"//")[[1]][2]
  y <- strsplit(pathname, "_")[[1]] %>% tolower
  type_of_triall <- paste(y[2:4],collapse="_")
  type_of_triall <- strsplit(type_of_triall,"[.]")[[1]][1]
  n <- y[1]
  name <- append(name,n)
  type_of_trial <- append(type_of_trial,type_of_triall)
  
  print(files[i])
  
  # figure out which frame define the different parts of the trial:
  frame_end_background_1 <- floor(length(x)*0.2381)
  frame_start_males_1 <- frame_end_background_1 + 1
  frame_end_males1 <- floor(length(x)*0.4762)
  frame_starts_males2 <- floor(length(x)*0.5238)
  frame_end_males2 <- floor(length(x)*0.7619)
  frame_start_backgroud2 <- frame_end_males2 + 1
  
  # subset the dataframe into parts
  background1 <- x[1:frame_end_background_1]
  males1 <- x[frame_start_males_1:frame_end_males1]
  males2 <- x[frame_starts_males2:frame_end_males2]
  background2 <- x[frame_start_backgroud2:length(x)]
  
  background1_left <- append(background1_left, summary(factor(background1))["left"][[1]])
  background1_neutral <- append(background1_neutral, summary(factor(background1))["neutral"][[1]])
  background1_right <- append(background1_right, summary(factor(background1))["right"][[1]])
  
  background2_left <- append(background2_left, summary(factor(background2))["left"][[1]])
  background2_neutral <- append(background2_neutral, summary(factor(background2))["neutral"][[1]])
  background2_right <- append(background2_right, summary(factor(background2))["right"][[1]])
  
  males1_left <- append(males1_left, summary(factor(males1))["left"][[1]])
  males1_neutral <- append(males1_neutral, summary(factor(males1))["neutral"][[1]])
  males1_right <- append(males1_right, summary(factor(males1))["right"][[1]])
  
  males2_left <- append(males2_left, summary(factor(males2))["left"][[1]])
  males2_neutral <- append(males2_neutral, summary(factor(males2))["neutral"][[1]])
  males2_right <- append(males2_right, summary(factor(males2))["right"][[1]])
  
}
## [1] "/Users/lukereding/Desktop/transitivity_trials//Adeline_Intermediate_vs_Large.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Adeline_large_vs_small.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Adeline_Small_vs_Intermediate.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Amie_intermediate_vs_large.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Amie_large_vs_small.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Amie_small_vs_intermediate.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Aniya_intermediate_vs_large.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Aniya_large_vs_small.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Aniya_small_vs_intermediate.avi_remade.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Aniya_small_vs_intermediate.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Arianna_Intermediate_vs_Large.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Arianna_large_vs_small.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Arianna_small_vs_intermediate.avi_remade.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Arianna_small_vs_intermediate.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Autumn_large_vs_small.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Beulah_small_vs_intermediate.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Briella_small_vs_intermediate.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Cherie_intermediate_vs_large.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Cherie_Large_vs_Small.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Cherie_Small_vs_Intermediate.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Deloris_intermediate_vs_large.avi_remade.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Deloris_intermediate_vs_large.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Deloris_large_vs_small.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Deloris_small_vs_intermediate.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Elise_intermediate_vs_large.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Eliza_intermediate_vs_large.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Eliza_Large_vs_Small.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Eliza_Small_vs_Intermediate.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Erika_Intermediate_vs_Large.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Erika_large_vs_small.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Erika_small_vs_intermediate.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Jade_intermediate_vs_large.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Kim_intermediate_vs_large.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Kim_Large_vs_Small.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Kim_small_vs_intermediate.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Lottie_Intermediate_vs_Large.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Lottie_Large_vs_Small.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Lottie_small_vs_intermediate.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Madeline_intermediate_vs_large.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Madeline_Large_vs_Small.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Madeline_small_vs_intermediate.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Madyson_Intermediate_vs_Large.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Madyson_large_vs_small.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Madyson_Small_vs_Intermediate.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//May_intermediate_vs_large.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//May_large_vs_small.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//May_small_vs_intermediate.avi_remade.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//May_small_vs_intermediate.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Nettie_intermediate_vs_large.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Nettie_Large_vs_Small.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Nettie_small_vs_intermediate.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Sierra_intermediate_vs_large.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Sierra_Large_vs_Small.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Sierra_Small_vs_Intermediate.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Sophia_intermediate_vs_large.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Tamara_intermediate_vs_large.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Tamara_small_Vs_intermediate.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Tamra_Large_vs_Small.avi_remade.txt"
## [1] "/Users/lukereding/Desktop/transitivity_trials//Tamra_Large_vs_Small.txt"
# put everything together into a data frame
results <- data.frame(files,name, type_of_trial,background1_left,background1_neutral,background1_right,background2_left,background2_neutral,background2_right,males1_left,males1_neutral,males1_right,males2_left,males2_neutral,males2_right)

# replace NAs with zeros
results[is.na(results)] <- 0

results$trial_id <- with(results, paste(name,type_of_trial, sep="_"))

Let’s check to see if there’s an overall side bias:

prop_left_background1 <- with(results, background1_left / (background1_left + background1_right))
# if 0 / some number, return 0:
prop_left_background1[is.na(prop_left_background1)] <- 0

prop_left_background2 <- with(results, background2_left / (background2_left + background2_right))
prop_left_background2[is.na(prop_left_background2)] <- 0

# plot it
beeStripBox(list(prop_left_background1,prop_left_background2),lab = c("first background period", "second background period"),ylab= "proportion of time on left side of tank",main="no overall bias towards one side of the tank when viewing background",cex.main=0.9)
## [1] "data are a list"
abline(h=0.5)

Looks good.

What about during the male stimulus trials?

prop_left_males1 <- with(results, males1_left / (males1_left + males1_right))
prop_left_males1[is.na(prop_left_males1)] <- 0

prop_left_males2 <- with(results, males2_left / (males2_left + males2_right))
prop_left_males2[is.na(prop_left_males2)] <- 0

# plot
beeStripBox(list(prop_left_males1,prop_left_males2),lab = c("first male period", "second male period"),ylab= "proportion of time on left side of tank",main="no overall bias towards one side of the tank when viewing males",cex.main=0.9)
## [1] "data are a list"
abline(h=0.5)

Looks good.

Looking at the aggregate data

Now let’s read in the trial data. This .csv file contains the location of each male during the trial.

** Doublecheck to make sure all the data were entered correctly here. **

trial_data<-read.csv("trial_data.csv",header=T)
trial_data$female %<>% tolower
trial_data$trial_id <- with(trial_data, paste(female,test, sep="_"))
head(trial_data)
##     female                  test    date time.of.trial weight   sl
## 1 madeline intermediate_vs_large 9/28/15         11:00   0.62 25.6
## 2 madeline small_vs_intermediate 9/30/15         10:00   0.81 26.4
## 3 madeline        large_vs_small 10/2/15         10:30   0.59 26.1
## 4   lottie small_vs_intermediate 9/28/15         11:35   0.98 28.3
## 5   lottie intermediate_vs_large 9/30/15         10:40   0.67 29.5
## 6   lottie        large_vs_small 10/2/15         11:10   0.68 27.5
##   water_temp male_left_side_part_4 male_right_side_part_4 name bias
## 1       80.8                 large           intermediate luke   no
## 2       81.2          intermediate                  small  joe   no
## 3       80.1                 small                  large  joe   no
## 4       80.3                 small           intermediate luke   no
## 5       80.6          intermediate                  large  joe   no
## 6       80.0                 large                  small  joe   no
##                         trial_id
## 1 madeline_intermediate_vs_large
## 2 madeline_small_vs_intermediate
## 3        madeline_large_vs_small
## 4   lottie_small_vs_intermediate
## 5   lottie_intermediate_vs_large
## 6          lottie_large_vs_small
# merge dataframes together

trials <- merge(trial_data,results,by="trial_id")

For my purposes right now, I’m just going to look at any trial where there was not evidence of bias.

no_bias <- subset(trials,trials$bias=="no")
nrow(no_bias) # number of trials we're dealing with
## [1] 38
# make sure there are no typos
summary(no_bias$male_left_side_part_4) %>% length == 3
## [1] TRUE
summary(no_bias$male_right_side_part_4) %>% length == 3
## [1] TRUE

Now we need to associate a side of a tank with a given type of male.

To maintain my sanity, for now, I’m going to break the trial up into two parts: parts where the female sees males, and parts where she sees the background.

no_bias %<>% mutate(frames_male_left_part4 = males2_left + males1_right)
no_bias %<>% mutate(frames_male_right_part4 = males2_right + males1_left)

no_bias %<>% mutate(background_left = background1_left + background2_left)
no_bias %<>% mutate(background_right = background1_right + background2_right)


no_bias %<>% mutate(frames_large = ifelse(male_left_side_part_4 == "large", frames_male_right_part4, ifelse(male_right_side_part_4 == "large",  frames_male_left_part4, NA)))


no_bias %<>% mutate(frames_small = ifelse(male_left_side_part_4 == "small", frames_male_right_part4, ifelse(male_right_side_part_4 == "small",  frames_male_left_part4, NA)))

no_bias %<>% mutate(frames_intermediate = ifelse(male_left_side_part_4 == "intermediate", frames_male_right_part4, ifelse(male_right_side_part_4 == "intermediate",  frames_male_left_part4, NA)))

Let’s now only look at data from the large vs. small male experiments:

large_vs_small <- subset(no_bias, no_bias$type_of_trial == "large_vs_small")
nrow(large_vs_small)
## [1] 13
# plot
with(large_vs_small, beeStripBox(list(frames_large/10, frames_small/10),lab=c("large","small"),ylab="seconds",main= "large vs small male trials"))
## [1] "data are a list"

# plot differences
with(large_vs_small, beeStripBox(list(frames_large/10 - frames_small/10),lab=c("large - small"),ylab="difference in seconds", main= "large vs small male trials"))
## [1] "data are a list"
abline(h=0,col="#00000050",lwd=1.2)
with(large_vs_small, wilcox.test(frames_large/10 - frames_small/10))
## 
##  Wilcoxon signed rank test
## 
## data:  frames_large/10 - frames_small/10
## V = 71, p-value = 0.08032
## alternative hypothesis: true location is not equal to 0
legend("bottomleft", paste("p = ",with(large_vs_small, wilcox.test(frames_large/10 - frames_small/10))$p.value %>% round(3)," by Wilcoxon",sep=""), bty="n") 

Let’s now only look at data from the large vs. intermediate male experiments:

large_vs_intermediate <- subset(no_bias, no_bias$type_of_trial == "intermediate_vs_large")
nrow(large_vs_intermediate)
## [1] 14
# plot
with(large_vs_intermediate, beeStripBox(list(frames_large/10, frames_intermediate/10),lab=c("large","intermediate"),ylab="seconds",main= "large vs intermediate male trials"))
## [1] "data are a list"

# plot differences
with(large_vs_intermediate, beeStripBox(list(frames_large/10 - frames_intermediate/10),lab=c("large - intermediate"),ylab="difference in seconds", main= "large vs intermediate male trials"))
## [1] "data are a list"
abline(h=0,col="#00000050",lwd=1.2)
with(large_vs_intermediate, wilcox.test(frames_large/10 - frames_intermediate/10))
## 
##  Wilcoxon signed rank test
## 
## data:  frames_large/10 - frames_intermediate/10
## V = 73, p-value = 0.2166
## alternative hypothesis: true location is not equal to 0
legend("bottomleft", paste("p = ",with(large_vs_intermediate, wilcox.test(frames_large/10 - frames_intermediate/10))$p.value %>% round(3)," by Wilcoxon",sep=""), bty="n")

Let’s now only look at data from the intermediate vs. small male experiments:

intermediate_vs_small <- subset(no_bias, no_bias$type_of_trial == "small_vs_intermediate")
nrow(intermediate_vs_small)
## [1] 11
# plot
with(intermediate_vs_small, beeStripBox(list(frames_intermediate/10, frames_small/10),lab=c("intermediate","small"),ylab="seconds",main= "intermediate vs small male trials"))
## [1] "data are a list"

# plot differences
with(intermediate_vs_small, beeStripBox(list(frames_intermediate/10 - frames_small/10),lab=c("intermediate - small"),ylab="difference in seconds", main= "intermediate vs small male trials"))
## [1] "data are a list"
abline(h=0,col="#00000050",lwd=1.2)
with(intermediate_vs_small, wilcox.test(frames_intermediate/10 - frames_small/10))
## 
##  Wilcoxon signed rank test
## 
## data:  frames_intermediate/10 - frames_small/10
## V = 44, p-value = 0.3652
## alternative hypothesis: true location is not equal to 0
legend("bottomleft", paste("p = ",with(intermediate_vs_small, wilcox.test(frames_intermediate/10 - frames_small/10))$p.value %>% round(3)," by Wilcoxon",sep=""), bty="n")

analyzing patterns of preference in each trial

The idea: If females are really expressing a preference for one of the fish, we should see a bias towards one side of the tank when males are initially shown, and a bias towards the opposite side of the tank in the next part of the trial when the males are switched sides.

source("./stats/generate_null_distributions.R")
files <- list.files(path="~/Desktop/transitivity_trials/",pattern = "[ltg][el].txt",full.names = T)
pass<-c()

# set up a dataframe
result <- data.frame(rep(NA,length(files)),rep(NA,length(files)),rep(NA,length(files)),rep(NA,length(files)))
names(result)<-c("observed","cutoff","verdict","name")

# loop
library("plyr")
## -------------------------------------------------------------------------
## You have loaded plyr after dplyr - this is likely to cause problems.
## If you need functions from both plyr and dplyr, please load plyr first, then dplyr:
## library(plyr); library(dplyr)
## -------------------------------------------------------------------------
## 
## Attaching package: 'plyr'
## 
## The following objects are masked from 'package:dplyr':
## 
##     arrange, count, desc, failwith, id, mutate, rename, summarise,
##     summarize
for(i in 1:length(files)){
  
  print(files[i])
  
  x<- male_monte(files[i])
  
  result[i,] <- c(x[[1]],x[[2]],x[[3]],files[i])
  
}
## [1] "/Users/lukereding/Desktop/transitivity_trials//Adeline_Intermediate_vs_Large.txt"
## [1] 12549
## [1] neutral neutral neutral neutral neutral
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Adeline_large_vs_small.txt"
## [1] 12577
## [1] neutral neutral neutral neutral neutral
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Adeline_Small_vs_Intermediate.txt"
## [1] 12548
## [1] neutral neutral neutral neutral neutral
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Amie_intermediate_vs_large.txt"
## [1] 12533
## [1] neutral neutral neutral neutral neutral
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Amie_large_vs_small.txt"
## [1] 12575
## [1] neutral neutral neutral neutral neutral
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Amie_small_vs_intermediate.txt"
## [1] 12542
## [1] neutral neutral neutral neutral neutral
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Aniya_intermediate_vs_large.txt"
## [1] 12582
## [1] neutral neutral neutral neutral neutral
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Aniya_large_vs_small.txt"
## [1] 12537
## [1] neutral neutral neutral neutral neutral
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Aniya_small_vs_intermediate.txt"
## [1] 12520
## [1] neutral neutral neutral neutral neutral
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Arianna_Intermediate_vs_Large.txt"
## [1] 12535
## [1] neutral neutral neutral neutral neutral
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Arianna_large_vs_small.txt"
## [1] 12537
## [1] neutral neutral neutral neutral neutral
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Arianna_small_vs_intermediate.txt"
## [1] 12582
## [1] neutral neutral neutral neutral right  
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Autumn_large_vs_small.txt"
## [1] 12553
## [1] neutral neutral neutral neutral neutral
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Beulah_small_vs_intermediate.txt"
## [1] 12600
## [1] neutral neutral neutral neutral neutral
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Briella_small_vs_intermediate.txt"
## [1] 12572
## [1] neutral neutral neutral neutral neutral
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Cherie_intermediate_vs_large.txt"
## [1] 12567
## [1] neutral neutral neutral neutral neutral
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Cherie_Large_vs_Small.txt"
## [1] 12519
## [1] neutral neutral neutral neutral neutral
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Cherie_Small_vs_Intermediate.txt"
## [1] 12533
## [1] neutral neutral neutral neutral neutral
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Deloris_intermediate_vs_large.txt"
## [1] 12566
## [1] neutral neutral neutral neutral neutral
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Deloris_large_vs_small.txt"
## [1] 12531
## [1] neutral neutral neutral neutral neutral
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Deloris_small_vs_intermediate.txt"
## [1] 12583
## [1] neutral neutral neutral neutral neutral
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Elise_intermediate_vs_large.txt"
## [1] 12600
## [1] neutral neutral neutral left    left   
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Eliza_intermediate_vs_large.txt"
## [1] 12569
## [1] neutral neutral neutral neutral neutral
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Eliza_Large_vs_Small.txt"
## [1] 12529
## [1] neutral neutral neutral neutral neutral
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Eliza_Small_vs_Intermediate.txt"
## [1] 12546
## [1] neutral neutral neutral neutral neutral
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Erika_Intermediate_vs_Large.txt"
## [1] 12534
## [1] neutral neutral neutral neutral neutral
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Erika_large_vs_small.txt"
## [1] 12547
## [1] neutral neutral neutral right   right  
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Erika_small_vs_intermediate.txt"
## [1] 12535
## [1] neutral neutral neutral neutral neutral
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Jade_intermediate_vs_large.txt"
## [1] 12599
## [1] neutral neutral neutral neutral neutral
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Kim_intermediate_vs_large.txt"
## [1] 12567
## [1] neutral neutral neutral neutral neutral
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Kim_Large_vs_Small.txt"
## [1] 12558
## [1] neutral neutral neutral neutral neutral
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Kim_small_vs_intermediate.txt"
## [1] 12558
## [1] neutral neutral neutral neutral neutral
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Lottie_Intermediate_vs_Large.txt"
## [1] 12532
## [1] neutral neutral neutral neutral neutral
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Lottie_Large_vs_Small.txt"
## [1] 12595
## [1] neutral neutral neutral neutral neutral
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Lottie_small_vs_intermediate.txt"
## [1] 12438
## [1] neutral neutral neutral neutral neutral
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Madeline_intermediate_vs_large.txt"
## [1] 12564
## [1] neutral neutral neutral neutral neutral
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Madeline_Large_vs_Small.txt"
## [1] 12587
## [1] neutral neutral left    left    left   
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Madeline_small_vs_intermediate.txt"
## [1] 12552
## [1] neutral neutral neutral neutral right  
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Madyson_Intermediate_vs_Large.txt"
## [1] 12562
## [1] neutral neutral neutral neutral neutral
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Madyson_large_vs_small.txt"
## [1] 12537
## [1] neutral neutral neutral neutral neutral
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Madyson_Small_vs_Intermediate.txt"
## [1] 12525
## [1] neutral neutral neutral neutral neutral
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//May_intermediate_vs_large.txt"
## [1] 12557
## [1] neutral neutral neutral neutral neutral
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//May_large_vs_small.txt"
## [1] 12554
## [1] neutral neutral neutral neutral neutral
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//May_small_vs_intermediate.txt"
## [1] 12564
## [1] neutral neutral neutral neutral neutral
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Nettie_intermediate_vs_large.txt"
## [1] 12589
## [1] neutral neutral neutral neutral left   
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Nettie_Large_vs_Small.txt"
## [1] 12579
## [1] neutral neutral neutral neutral neutral
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Nettie_small_vs_intermediate.txt"
## [1] 12593
## [1] neutral neutral neutral neutral neutral
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Sierra_intermediate_vs_large.txt"
## [1] 12549
## [1] neutral neutral neutral neutral right  
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Sierra_Large_vs_Small.txt"
## [1] 12544
## [1] neutral neutral neutral neutral neutral
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Sierra_Small_vs_Intermediate.txt"
## [1] 12570
## [1] neutral neutral neutral neutral neutral
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Sophia_intermediate_vs_large.txt"
## [1] 12593
## [1] neutral neutral neutral neutral neutral
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Tamara_intermediate_vs_large.txt"
## [1] 12564
## [1] neutral neutral neutral neutral neutral
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Tamara_small_Vs_intermediate.txt"
## [1] 12573
## [1] neutral neutral neutral neutral neutral
## Levels: left neutral right
## [1] "on trial 1000"

## [1] "/Users/lukereding/Desktop/transitivity_trials//Tamra_Large_vs_Small.txt"
## [1] 12575
## [1] neutral neutral neutral neutral neutral
## Levels: left neutral right
## [1] "on trial 1000"

summary(result$verdict %>% factor)
## 0.0199720556253186 0.0213799503912178 0.0218380551751172 
##                  1                  1                  1 
##  0.022005253498899 0.0228604519714872 0.0230017087619016 
##                  1                  1                  1 
## 0.0231974216651467 0.0233670440143904 0.0238083541742731 
##                  1                  1                  1 
## 0.0239138675720244 0.0249349790234518 0.0251662605391063 
##                  1                  1                  1 
## 0.0251686326718017 0.0253981006744662 0.0258823520406083 
##                  1                  1                  1 
## 0.0260725085857668 0.0261785781694999 0.0262198118023189 
##                  1                  1                  1 
## 0.0263256950132197 0.0270221749764198 0.0274427034428921 
##                  1                  1                  1 
## 0.0276387544819715 0.0278341297360253 0.0280195281875897 
##                  1                  1                  1 
## 0.0284309983377286 0.0286796851704885 0.0287074501870041 
##                  1                  1                  1 
## 0.0287157064951473 0.0289203834409998 0.0293449348904927 
##                  1                  1                  1 
## 0.0295136327151927 0.0298110841862939 0.0299606722031315 
##                  1                  1                  1 
##  0.029972594495238 0.0300362895429693 0.0300641709569598 
##                  1                  1                  1 
## 0.0302821360985733 0.0302979210226296 0.0303322030362621 
##                  1                  1                  1 
## 0.0304776594591911 0.0304808831541206 0.0306383684858692 
##                  1                  1                  1 
## 0.0308427444258213 0.0315013984521141 0.0321489891134357 
##                  1                  1                  1 
## 0.0329920855070339 0.0334568721307393 0.0340276328776179 
##                  1                  1                  1 
## 0.0346548417278248   0.03564494615326 0.0416894429069035 
##                  1                  1                  1 
## 0.0431479870823143  0.043759856705085 0.0477812251664059 
##                  1                  1                  1