Load required packages:
library(circular)
##
## Attaching package: 'circular'
##
## The following objects are masked from 'package:stats':
##
## sd, var
library(ggplot2)
Load custom R scripts:
path_lib <- "C:/Users/Christie/Documents/Research/Projects/AttentionVSTM/Library"
lib <- list.files(path_lib)
for (i in 1:length(lib)) {
source(file.path(path_lib, lib[i]))
}
Create data frame of all subject data for Experiment 1:
path_data <- "C:/Users/Christie/Documents/Research/Projects/AttentionVSTM/Experiments/Experiment 2/data"
d1 <- read.data(path_data)
Remove poorly performing subjects, and practice trials:
data <- d1[which(d1$userid != 17 & d1$userid != 22 & d1$practice == "FALSE"),
]
Define angular difference as circular data:
angdiff <- circular(data$angdiff, type = "angle", units = "degrees")
Generate polar co-ordinate plot:
p <- ggplot(data, aes(x = angdiff)) + geom_histogram(aes(y = ..density..), breaks = seq(-90,
90, by = 180/21), colour = "black", fill = "gray")
p <- p + coord_polar(start = pi)
p <- p + labs(title = "Experiment 1: Distribution of offset on all tirals",
x = "Offset (degrees)", y = "Density") + theme_bw()
p