12345678910111213141516171819202122 |
- library(tidyverse)
- library(R.matlab)
- res = readMat('test_RS04.mat')
- dat = as.tibble(res$trials)
- names(dat) = c('Duration', 'session','trlNo','blkNo','dur1','pdur','production','vrep','Reproduction')
- mrep = dat %>% group_by(session, Duration) %>%
- summarise(mRep = mean(Reproduction), n = n(), se = sd(Reproduction)/sqrt(n-1))
- mrep$session = factor(mrep$session, labels = c("RWalk", "Random"))
- mrep%>%
- ggplot(aes(Duration, mRep, color = session, group = session)) +
- geom_point() + geom_line() +
- geom_errorbar(aes(ymin = mRep - se, ymax = mRep +se), width = 0.05) +
- theme_classic() +
- geom_abline(slope = 1)
- ggplot(dat, aes(trlNo, Duration)) + geom_line() + geom_point()
|