Created
February 8, 2016 21:50
-
-
Save CJ-Davies/e69fa23d8d9d20ba6471 to your computer and use it in GitHub Desktop.
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| # usage | |
| # $ R CMD BATCH script.r | |
| library(ggplot2) | |
| library(reshape2) | |
| library(gridExtra) | |
| ################################################################################ | |
| # INPUT | |
| ################################################################################ | |
| # traditional every 0.5 seconds | |
| resStat <- read.csv("stationary_elapsed.log",header=T,sep="\t") | |
| # 1-3 every 0.5 seconds | |
| res13 <- read.csv("1-3_log_elapsed_0.5.log",header=T,sep="\t") | |
| # 1-4 every 0.1 seconds | |
| res14 <- read.csv("1-4_log_elapsed_0.1.log",header=T,sep="\t") | |
| resStat_combined_data <- data.frame(time=resStat$timestamp, delta_x=resStat$delta_x, delta_z=resStat$delta_z, yaw=resStat$left_y, pitch=resStat$left_x, frames=resStat$framerate) | |
| #resStat_combined_data <- subset(resStat_combined_data_,time>0&&time<=175) | |
| res13_combined_data <- data.frame(time=res13$timestamp, delta_x=res13$delta_x, delta_z=res13$delta_z, yaw=res13$left_y, pitch=res13$left_x, environment=res13$environment, opacity=res13$left_opacity, frames=res13$framerate) | |
| res14_combined_data <- data.frame(time=res14$timestamp, delta_x=res14$delta_x, delta_z=res14$delta_z, yaw=res14$left_y, pitch=res14$left_x, environment=res14$environment, opacity=res14$left_opacity, frames=res14$framerate) | |
| merged_combined_data <- merge(res13_combined_data,res14_combined_data,by="time",all=T) | |
| ################################################################################ | |
| # PLOTS | |
| ################################################################################ | |
| # black & red for lines | |
| linePal <- c("#000000", "#D55E00") | |
| # blue, pink, green, yellow & dark blue for background | |
| rectPal <- c("#56B4E9","#CC79A7","#009E73","#F0E442", "#0072B2") | |
| ################################################################################ | |
| # 1-3 2up | |
| p131 <- ggplot(res13_combined_data,aes(time))+geom_rect(data=res13_combined_data,aes(xmin=(round(time*2))/2,xmax=(((round(time*2))/2)+0.5),fill=environment,ymin=-180,ymax=Inf),alpha=0.4)+geom_line(aes(y=yaw,col="yaw"),alpha=1,na.rm=T)+geom_line(aes(y=pitch,col="pitch"),alpha=1)+xlab("Time (seconds)")+coord_cartesian(ylim = c(-200, 200))+ylab("Rotation (degrees)")+ggtitle("Pitch & yaw against time for PR scenario\nwith transitions 1-3")+theme(legend.position="top",legend.title=element_blank(), axis.text=element_text(size=16),axis.title=element_text(size=18),legend.text=element_text(size=16),plot.title=element_text(size=20))+scale_fill_manual(values=rectPal,breaks=c("real","transition 1","transition 2","transition 3"))+scale_colour_manual(values=linePal)+scale_y_continuous(breaks = c(-180,-135,-90,-45,0,45,90,135,180))+scale_x_continuous(breaks=c(0,50,100,150,200,250),limits=c(0,285)) | |
| p132 <- ggplot(res13_combined_data,aes(time))+geom_rect(data=res13_combined_data,aes(xmin=(round(time*2))/2,xmax=(((round(time*2))/2)+0.5),fill=environment,ymin=0,ymax=Inf),alpha=0.4)+geom_line(aes(y=delta_x,col="East/West"),alpha=1)+geom_line(aes(y=delta_z,col="North/South"),alpha=1)+xlab("Time (seconds)")+ylab("Distance (meters)")+ggtitle("Distance moved against time for PR scenario\nwith transitions 1-3")+theme(legend.position="top",legend.title=element_blank(), axis.text=element_text(size=16),axis.title=element_text(size=18),legend.text=element_text(size=16),plot.title=element_text(size=20))+scale_fill_manual(values=rectPal)+scale_colour_manual(values=linePal)+scale_x_continuous(breaks=c(0,50,100,150,200,250),limits=c(0,285)) | |
| gA <- ggplotGrob(p131) | |
| gB <- ggplotGrob(p132) | |
| maxWidth = grid::unit.pmax(gA$widths[2:5], gB$widths[2:5]) | |
| gA$widths[2:5] <- as.list(maxWidth) | |
| gB$widths[2:5] <- as.list(maxWidth) | |
| png("1-3_2up.png",width=1920,height=1080) | |
| grid.arrange(gA, gB, ncol=1) | |
| dev.off() | |
| ################################################################################ | |
| # 1-3 opacity against time | |
| p131 <- ggplot(res13_combined_data,aes(time))+geom_rect(data=res13_combined_data,aes(xmin=time,xmax=time+0.5,fill=environment,ymin=-180,ymax=Inf),alpha=0.4)+geom_line(aes(y=opacity,col="opacity"),alpha=1,na.rm=T)+xlab("Time (seconds)")+coord_cartesian(ylim = c(0, 1))+ylab("Opacity")+ggtitle("Opacity against time for PR scenario\nwith transitions 1-3")+theme(legend.position="top",legend.title=element_blank(), axis.text=element_text(size=16),axis.title=element_text(size=18),legend.text=element_text(size=16),plot.title=element_text(size=20))+scale_fill_manual(values=rectPal,breaks=c("real","transition 1","transition 2","transition 3","transition 4"))+scale_colour_manual(values=linePal)+scale_x_continuous(breaks=c(0,50,100,150,200,250),limits=c(0,285)) | |
| gA <- ggplotGrob(p131) | |
| png("1-3_opacity.png",width=1920,height=540) | |
| plot(gA) | |
| dev.off() | |
| ################################################################################ | |
| # 1-3 2up opacity | |
| p131 <- ggplot(res13_combined_data,aes(time))+geom_rect(data=res13_combined_data,aes(xmin=(round(time*2))/2,xmax=(((round(time*2))/2)+0.5),fill=environment,ymin=-180,ymax=Inf),alpha=0.4)+geom_line(aes(y=yaw,col="yaw"),alpha=1,na.rm=T)+geom_line(aes(y=pitch,col="pitch"),alpha=1)+xlab("Time (seconds)")+coord_cartesian(ylim = c(-200, 200))+ylab("Rotation (degrees)")+ggtitle("Pitch & yaw against time for PR scenario\nwith transitions 1-3")+theme(legend.position="top",legend.title=element_blank(), axis.text=element_text(size=16),axis.title=element_text(size=18),legend.text=element_text(size=16),plot.title=element_text(size=20))+scale_fill_manual(values=rectPal,breaks=c("real","transition 1","transition 2","transition 3"))+scale_colour_manual(values=linePal)+scale_y_continuous(breaks = c(-180,-135,-90,-45,0,45,90,135,180))+scale_x_continuous(breaks=c(0,50,100,150,200,250),limits=c(0,285)) | |
| p132 <- ggplot(res13_combined_data,aes(time))+geom_rect(data=res13_combined_data,aes(xmin=time,xmax=time+0.5,fill=environment,ymin=-180,ymax=Inf),alpha=0.4)+geom_line(aes(y=opacity,col="opacity"),alpha=1,na.rm=T)+xlab("Time (seconds)")+coord_cartesian(ylim = c(0, 1))+ylab("Opacity")+ggtitle("Opacity against time for PR scenario\nwith transitions 1-3")+theme(legend.position="top",legend.title=element_blank(), axis.text=element_text(size=16),axis.title=element_text(size=18),legend.text=element_text(size=16),plot.title=element_text(size=20))+scale_fill_manual(values=rectPal,breaks=c("real","transition 1","transition 2","transition 3","transition 4"))+scale_colour_manual(values=linePal)+scale_x_continuous(breaks=c(0,50,100,150,200,250),limits=c(0,285)) | |
| gA <- ggplotGrob(p131) | |
| gB <- ggplotGrob(p132) | |
| maxWidth = grid::unit.pmax(gA$widths[2:5], gB$widths[2:5]) | |
| gA$widths[2:5] <- as.list(maxWidth) | |
| gB$widths[2:5] <- as.list(maxWidth) | |
| png("1-3_2up_opacity.png",width=1920,height=1080) | |
| grid.arrange(gA, gB, ncol=1) | |
| dev.off() | |
| ################################################################################ | |
| # 1-3 3up | |
| p131 <- ggplot(res13_combined_data,aes(time))+geom_rect(data=res13_combined_data,aes(xmin=(round(time*2))/2,xmax=(((round(time*2))/2)+0.5),fill=environment,ymin=-180,ymax=Inf),alpha=0.4)+geom_line(aes(y=yaw,col="yaw"),alpha=1,na.rm=T)+geom_line(aes(y=pitch,col="pitch"),alpha=1)+xlab("Time (seconds)")+coord_cartesian(ylim = c(-200, 200))+ylab("Rotation (degrees)")+ggtitle("Pitch & yaw against time for PR scenario\nwith transitions 1-3")+theme(legend.position="top",legend.title=element_blank(), axis.text=element_text(size=16),axis.title=element_text(size=18),legend.text=element_text(size=16),plot.title=element_text(size=20))+scale_fill_manual(values=rectPal,breaks=c("real","transition 1","transition 2","transition 3"))+scale_colour_manual(values=linePal)+scale_y_continuous(breaks = c(-180,-135,-90,-45,0,45,90,135,180))+scale_x_continuous(breaks=c(0,50,100,150,200,250),limits=c(0,285)) | |
| p132 <- ggplot(res13_combined_data,aes(time))+geom_rect(data=res13_combined_data,aes(xmin=(round(time*2))/2,xmax=(((round(time*2))/2)+0.5),fill=environment,ymin=0,ymax=Inf),alpha=0.4)+geom_line(aes(y=delta_x,col="East/West"),alpha=1)+geom_line(aes(y=delta_z,col="North/South"),alpha=1)+xlab("Time (seconds)")+ylab("Distance (meters)")+ggtitle("Distance moved against time for PR scenario\nwith transitions 1-3")+theme(legend.position="top",legend.title=element_blank(), axis.text=element_text(size=16),axis.title=element_text(size=18),legend.text=element_text(size=16),plot.title=element_text(size=20))+scale_fill_manual(values=rectPal)+scale_colour_manual(values=linePal)+scale_x_continuous(breaks=c(0,50,100,150,200,250),limits=c(0,285)) | |
| p133 <- ggplot(res13_combined_data,aes(time))+geom_rect(data=res13_combined_data,aes(xmin=time,xmax=time+0.5,fill=environment,ymin=-180,ymax=Inf),alpha=0.4)+geom_line(aes(y=opacity,col="opacity"),alpha=1,na.rm=T)+xlab("Time (seconds)")+coord_cartesian(ylim = c(0, 1))+ylab("Opacity")+ggtitle("Opacity against time for PR scenario\nwith transitions 1-3")+theme(legend.position="top",legend.title=element_blank(), axis.text=element_text(size=16),axis.title=element_text(size=18),legend.text=element_text(size=16),plot.title=element_text(size=20))+scale_fill_manual(values=rectPal,breaks=c("real","transition 1","transition 2","transition 3","transition 4"))+scale_colour_manual(values=linePal)+scale_x_continuous(breaks=c(0,50,100,150,200,250),limits=c(0,285)) | |
| gA <- ggplotGrob(p131) | |
| gB <- ggplotGrob(p132) | |
| gC <- ggplotGrob(p133) | |
| maxWidth = grid::unit.pmax(gA$widths[2:5], gB$widths[2:5], gC$widths[2:5]) | |
| gA$widths[2:5] <- as.list(maxWidth) | |
| gB$widths[2:5] <- as.list(maxWidth) | |
| gC$widths[2:5] <- as.list(maxWidth) | |
| png("1-3_3up.png",width=1920,height=1357) | |
| grid.arrange(gA, gB, gC, ncol=1) | |
| dev.off() | |
| ################################################################################ | |
| # 1-4 2up | |
| p141 <- ggplot(res14_combined_data,aes(time))+geom_rect(data=res14_combined_data,aes(xmin=time,xmax=time+0.1,fill=environment,ymin=-180,ymax=Inf),alpha=0.4)+geom_line(aes(y=yaw,col="yaw"),alpha=1,na.rm=T)+geom_line(aes(y=pitch,col="pitch"),alpha=1)+xlim(0,225)+xlab("Time (seconds)")+coord_cartesian(ylim = c(-200, 200))+ylab("Rotation (degrees)")+ggtitle("Pitch & yaw against time for PR scenario\nwith transitions 1-4")+theme(legend.position="top",legend.title=element_blank(), axis.text=element_text(size=16),axis.title=element_text(size=18),legend.text=element_text(size=16),plot.title=element_text(size=20))+scale_fill_manual(values=rectPal,breaks=c("real","transition 1","transition 2","transition 3","transition 4"))+scale_colour_manual(values=linePal)+scale_y_continuous(breaks = c(-180,-135,-90,-45,0,45,90,135,180)) | |
| p142 <- ggplot(res14_combined_data,aes(time))+geom_rect(data=res14_combined_data,aes(xmin=time,xmax=time+0.1,fill=environment,ymin=0,ymax=Inf),alpha=0.4)+geom_line(aes(y=delta_x,col="East/West"),alpha=1)+geom_line(aes(y=delta_z,col="North/South"),alpha=1)+xlim(0,225)+xlab("Time (seconds)")+ylab("Distance (meters)")+ggtitle("Distance moved against time for PR scenario\nwith transitions 1-4")+theme(legend.position="top",legend.title=element_blank(), axis.text=element_text(size=16),axis.title=element_text(size=18),legend.text=element_text(size=16),plot.title=element_text(size=20))+scale_fill_manual(values=rectPal,breaks=c("real","transition 1","transition 2","transition 3","transition 4"))+scale_colour_manual(values=linePal)+ylim(0,30) | |
| gA <- ggplotGrob(p141) | |
| gB <- ggplotGrob(p142) | |
| maxWidth = grid::unit.pmax(gA$widths[2:5], gB$widths[2:5]) | |
| gA$widths[2:5] <- as.list(maxWidth) | |
| gB$widths[2:5] <- as.list(maxWidth) | |
| png("1-4_2up.png",width=1920,height=1080) | |
| grid.arrange(gA, gB, ncol=1) | |
| dev.off() | |
| ################################################################################ | |
| # 1-4 opacity against time | |
| p141 <- ggplot(res14_combined_data,aes(time))+geom_rect(data=res14_combined_data,aes(xmin=time,xmax=time+0.1,fill=environment,ymin=-180,ymax=Inf),alpha=0.4)+geom_line(aes(y=opacity,col="opacity"),alpha=1,na.rm=T)+xlab("Time (seconds)")+coord_cartesian(ylim = c(0, 1))+ylab("Opacity")+ggtitle("Opacity against time for PR scenario\nwith transitions 1-4")+theme(legend.position="top",legend.title=element_blank(), axis.text=element_text(size=16),axis.title=element_text(size=18),legend.text=element_text(size=16),plot.title=element_text(size=20))+scale_fill_manual(values=rectPal,breaks=c("real","transition 1","transition 2","transition 3","transition 4"))+scale_colour_manual(values=linePal)+scale_x_continuous(breaks=c(0,50,100,150,200,250),limits=c(0,225)) | |
| gA <- ggplotGrob(p141) | |
| png("1-4_opacity.png",width=1920,height=540) | |
| plot(gA) | |
| dev.off() | |
| ################################################################################ | |
| # 1-4 2up opacity | |
| p141 <- ggplot(res14_combined_data,aes(time))+geom_rect(data=res14_combined_data,aes(xmin=time,xmax=time+0.1,fill=environment,ymin=-180,ymax=Inf),alpha=0.4)+geom_line(aes(y=yaw,col="yaw"),alpha=1,na.rm=T)+geom_line(aes(y=pitch,col="pitch"),alpha=1)+xlim(0,225)+xlab("Time (seconds)")+coord_cartesian(ylim = c(-200, 200))+ylab("Rotation (degrees)")+ggtitle("Pitch & yaw against time for PR scenario\nwith transitions 1-4")+theme(legend.position="top",legend.title=element_blank(), axis.text=element_text(size=16),axis.title=element_text(size=18),legend.text=element_text(size=16),plot.title=element_text(size=20))+scale_fill_manual(values=rectPal,breaks=c("real","transition 1","transition 2","transition 3","transition 4"))+scale_colour_manual(values=linePal)+scale_y_continuous(breaks = c(-180,-135,-90,-45,0,45,90,135,180)) | |
| p142 <- ggplot(res14_combined_data,aes(time))+geom_rect(data=res14_combined_data,aes(xmin=time,xmax=time+0.1,fill=environment,ymin=-180,ymax=Inf),alpha=0.4)+geom_line(aes(y=opacity,col="opacity"),alpha=1,na.rm=T)+xlab("Time (seconds)")+coord_cartesian(ylim = c(0, 1))+ylab("Opacity")+ggtitle("Opacity against time for PR scenario\nwith transitions 1-4")+theme(legend.position="top",legend.title=element_blank(), axis.text=element_text(size=16),axis.title=element_text(size=18),legend.text=element_text(size=16),plot.title=element_text(size=20))+scale_fill_manual(values=rectPal,breaks=c("real","transition 1","transition 2","transition 3","transition 4"))+scale_colour_manual(values=linePal)+scale_x_continuous(breaks=c(0,50,100,150,200,250),limits=c(0,225)) | |
| gA <- ggplotGrob(p141) | |
| gB <- ggplotGrob(p142) | |
| maxWidth = grid::unit.pmax(gA$widths[2:5], gB$widths[2:5]) | |
| gA$widths[2:5] <- as.list(maxWidth) | |
| gB$widths[2:5] <- as.list(maxWidth) | |
| png("1-4_2up_opacity.png",width=1920,height=1080) | |
| grid.arrange(gA, gB, ncol=1) | |
| dev.off() | |
| ################################################################################ | |
| # 1-4 3up | |
| p141 <- ggplot(res14_combined_data,aes(time))+geom_rect(data=res14_combined_data,aes(xmin=time,xmax=time+0.1,fill=environment,ymin=-180,ymax=Inf),alpha=0.4)+geom_line(aes(y=yaw,col="yaw"),alpha=1,na.rm=T)+geom_line(aes(y=pitch,col="pitch"),alpha=1)+xlim(0,225)+xlab("Time (seconds)")+coord_cartesian(ylim = c(-200, 200))+ylab("Rotation (degrees)")+ggtitle("Pitch & yaw against time for PR scenario\nwith transitions 1-4")+theme(legend.position="top",legend.title=element_blank(), axis.text=element_text(size=16),axis.title=element_text(size=18),legend.text=element_text(size=16),plot.title=element_text(size=20))+scale_fill_manual(values=rectPal,breaks=c("real","transition 1","transition 2","transition 3","transition 4"))+scale_colour_manual(values=linePal)+scale_y_continuous(breaks = c(-180,-135,-90,-45,0,45,90,135,180)) | |
| p142 <- ggplot(res14_combined_data,aes(time))+geom_rect(data=res14_combined_data,aes(xmin=time,xmax=time+0.1,fill=environment,ymin=0,ymax=Inf),alpha=0.4)+geom_line(aes(y=delta_x,col="East/West"),alpha=1)+geom_line(aes(y=delta_z,col="North/South"),alpha=1)+xlim(0,225)+xlab("Time (seconds)")+ylab("Distance (meters)")+ggtitle("Distance moved against time for PR scenario\nwith transitions 1-4")+theme(legend.position="top",legend.title=element_blank(), axis.text=element_text(size=16),axis.title=element_text(size=18),legend.text=element_text(size=16),plot.title=element_text(size=20))+scale_fill_manual(values=rectPal,breaks=c("real","transition 1","transition 2","transition 3","transition 4"))+scale_colour_manual(values=linePal)+ylim(0,30) | |
| p143 <- ggplot(res14_combined_data,aes(time))+geom_rect(data=res14_combined_data,aes(xmin=time,xmax=time+0.1,fill=environment,ymin=-180,ymax=Inf),alpha=0.4)+geom_line(aes(y=opacity,col="opacity"),alpha=1,na.rm=T)+xlab("Time (seconds)")+coord_cartesian(ylim = c(0, 1))+ylab("Opacity")+ggtitle("Opacity against time for PR scenario\nwith transitions 1-4")+theme(legend.position="top",legend.title=element_blank(), axis.text=element_text(size=16),axis.title=element_text(size=18),legend.text=element_text(size=16),plot.title=element_text(size=20))+scale_fill_manual(values=rectPal,breaks=c("real","transition 1","transition 2","transition 3","transition 4"))+scale_colour_manual(values=linePal)+scale_x_continuous(breaks=c(0,50,100,150,200,250),limits=c(0,225)) | |
| gA <- ggplotGrob(p141) | |
| gB <- ggplotGrob(p142) | |
| gC <- ggplotGrob(p143) | |
| maxWidth = grid::unit.pmax(gA$widths[2:5], gB$widths[2:5], gC$widths[2:5]) | |
| gA$widths[2:5] <- as.list(maxWidth) | |
| gB$widths[2:5] <- as.list(maxWidth) | |
| gC$widths[2:5] <- as.list(maxWidth) | |
| png("1-4_3up.png",width=1920,height=1357) | |
| grid.arrange(gA, gB, gC, ncol=1) | |
| dev.off() | |
| ################################################################################ | |
| # traditional 2-up | |
| pStat1 <- ggplot(resStat_combined_data,aes(time))+geom_line(aes(y=yaw,col="yaw"),alpha=1,na.rm=T)+geom_line(aes(y=pitch,col="pitch"),alpha=1)+xlab("Time (seconds)")+coord_cartesian(ylim = c(-200, 200))+ylab("Rotation (degrees)")+ggtitle("Pitch & yaw against time for traditional scenario")+theme(legend.position="top",legend.title=element_blank(), axis.text=element_text(size=16),axis.title=element_text(size=18),legend.text=element_text(size=16),plot.title=element_text(size=20))+scale_colour_manual(values=linePal)+scale_x_continuous(breaks=c(0,25,50,75,100,125,150,175),limits=c(0,175))+scale_y_continuous(breaks = c(-180,-135,-90,-45,0,45,90,135,180)) | |
| pStat2 <- ggplot(resStat_combined_data,aes(time))+geom_line(aes(y=delta_x,col="East/West"),alpha=1,na.rm=T)+geom_line(aes(y=delta_z,col="North/South"),alpha=1)+xlab("Time (seconds)")+ylab("Distance (meters)")+ggtitle("Distance moved against time for traditional scenario")+theme(legend.position="top",legend.title=element_blank(), axis.text=element_text(size=16),axis.title=element_text(size=18),legend.text=element_text(size=16),plot.title=element_text(size=20))+scale_colour_manual(values=linePal)+scale_x_continuous(breaks=c(0,25,50,75,100,125,150,175),limits=c(0,175)) | |
| gA <- ggplotGrob(pStat1) | |
| gB <- ggplotGrob(pStat2) | |
| maxWidth = grid::unit.pmax(gA$widths[2:5], gB$widths[2:5]) | |
| gA$widths[2:5] <- as.list(maxWidth) | |
| gB$widths[2:5] <- as.list(maxWidth) | |
| png("stationary_2up.png",width=1920,height=1080) | |
| grid.arrange(gA, gB, ncol=1) | |
| dev.off() | |
| ################################################################################ | |
| # pitch & yaw against time for traditional + 1-3 + 1-4 | |
| pStat1 <- ggplot(resStat_combined_data,aes(time))+geom_line(aes(y=yaw,col="yaw"),alpha=1,na.rm=T)+geom_line(aes(y=pitch,col="pitch"),alpha=1)+xlab("Time (seconds)")+coord_cartesian(ylim = c(-200, 200))+ylab("Rotation (degrees)")+ggtitle("Pitch and yaw against time for\nseated VR scenario\nfor participant 8")+theme(legend.position="top",legend.title=element_blank(), axis.text=element_text(size=24),axis.title=element_text(size=26),legend.text=element_text(size=24),plot.title=element_text(size=28))+scale_x_continuous(breaks=c(0,25,50,75,100,125,150,175),limits=c(0,175))+scale_y_continuous(breaks = c(-180,-135,-90,-45,0,45,90,135,180)) | |
| p131 <- ggplot(res13_combined_data,aes(time))+geom_line(aes(y=yaw,col="yaw"),alpha=1,na.rm=T)+geom_line(aes(y=pitch,col="pitch"),alpha=1)+xlab("Time (seconds)")+coord_cartesian(ylim = c(-200, 200))+ylab("Rotation (degrees)")+ggtitle("Pitch and yaw against time for\nparallel reality scenario with\ntransitions 1-3 for participant 8")+theme(legend.position="top",legend.title=element_blank(), axis.text=element_text(size=24),axis.title=element_text(size=26),legend.text=element_text(size=24),plot.title=element_text(size=28))+scale_fill_manual(values=rectPal,breaks=c("real","transition 1","transition 2","transition 3"))+scale_y_continuous(breaks = c(-180,-135,-90,-45,0,45,90,135,180))+scale_x_continuous(breaks=c(0,50,100,150,200,250),limits=c(0,285)) | |
| p141 <- ggplot(res14_combined_data,aes(time))+geom_line(aes(y=yaw,col="yaw"),alpha=1,na.rm=T)+geom_line(aes(y=pitch,col="pitch"),alpha=1)+xlim(0,225)+xlab("Time (seconds)")+coord_cartesian(ylim = c(-200, 200))+ylab("Rotation (degrees)")+ggtitle("Pitch and yaw against time for\nparallel reality scenario with\ntransitions 1-4 for participant 8")+theme(legend.position="top",legend.title=element_blank(), axis.text=element_text(size=24),axis.title=element_text(size=26),legend.text=element_text(size=24),plot.title=element_text(size=28))+scale_fill_manual(values=rectPal,breaks=c("real","transition 1","transition 2","transition 3","transition 4"))+scale_y_continuous(breaks = c(-180,-135,-90,-45,0,45,90,135,180)) | |
| gA <- ggplotGrob(pStat1) | |
| gB <- ggplotGrob(p131) | |
| gC <- ggplotGrob(p141) | |
| maxWidth = grid::unit.pmax(gA$widths[2:5], gB$widths[2:5], gC$widths[2:5]) | |
| gA$widths[2:5] <- as.list(maxWidth) | |
| gB$widths[2:5] <- as.list(maxWidth) | |
| gC$widths[2:5] <- as.list(maxWidth) | |
| png("8-pitch-yaw-trad-1-3-1-4.png",width=1920,height=678) | |
| grid.arrange(gA, gB, gC, ncol=3) | |
| dev.off() | |
| ################################################################################ | |
| p131 <- ggplot(res13_combined_data,aes(time))+geom_rect(data=res13_combined_data,aes(xmin=(round(time*2))/2,xmax=(((round(time*2))/2)+0.5),fill=environment,ymin=-180,ymax=Inf),alpha=0.4)+geom_line(aes(y=yaw,col="yaw"),alpha=1,na.rm=T)+geom_line(aes(y=pitch,col="pitch"),alpha=1)+xlab("Time (seconds)")+coord_cartesian(ylim = c(-200, 200))+ylab("Rotation (degrees)")+ggtitle("Pitch and yaw against time for parallel reality\nscenario with transitions 1-3 for participant 8")+theme(legend.position="top",legend.title=element_blank(), axis.text=element_text(size=24),axis.title=element_text(size=26),legend.text=element_text(size=24),plot.title=element_text(size=28))+scale_fill_manual(values=rectPal,breaks=c("real","transition 1","transition 2","transition 3"))+scale_colour_manual(values=linePal)+scale_y_continuous(breaks = c(-180,-135,-90,-45,0,45,90,135,180))+scale_x_continuous(breaks=c(0,50,100,150,200,250),limits=c(0,285)) | |
| gA <- ggplotGrob(p131) | |
| maxWidth = grid::unit.pmax(gA$widths[2:5]) | |
| gA$widths[2:5] <- as.list(maxWidth) | |
| png("8-1-3-pitch-yaw.png",width=1920,height=678) | |
| grid.arrange(gA, ncol=1) | |
| dev.off() | |
| ################################################################################ | |
| p141 <- ggplot(res14_combined_data,aes(time))+geom_rect(data=res14_combined_data,aes(xmin=time,xmax=time+0.1,fill=environment,ymin=-180,ymax=Inf),alpha=0.4)+geom_line(aes(y=yaw,col="yaw"),alpha=1,na.rm=T)+geom_line(aes(y=pitch,col="pitch"),alpha=1)+xlim(0,225)+xlab("Time (seconds)")+coord_cartesian(ylim = c(-200, 200))+ylab("Rotation (degrees)")+ggtitle("Pitch & yaw against time for PR scenario\nwith transitions 1-4 for participant 8")+theme(legend.position="top",legend.title=element_blank(), axis.text=element_text(size=24),axis.title=element_text(size=26),legend.text=element_text(size=24),plot.title=element_text(size=28))+scale_fill_manual(values=rectPal,breaks=c("real","transition 1","transition 2","transition 3","transition 4"))+scale_colour_manual(values=linePal)+scale_y_continuous(breaks = c(-180,-135,-90,-45,0,45,90,135,180)) | |
| gA <- ggplotGrob(p141) | |
| maxWidth = grid::unit.pmax(gA$widths[2:5]) | |
| gA$widths[2:5] <- as.list(maxWidth) | |
| png("8-1-4-pitch-yaw.png",width=1920,height=678) | |
| grid.arrange(gA, ncol=1) | |
| dev.off() | |
| ################################################################################ | |
| p131 <- ggplot(res13_combined_data,aes(time))+geom_rect(data=res13_combined_data,aes(xmin=(round(time*2))/2,xmax=(((round(time*2))/2)+0.5),fill=environment,ymin=-180,ymax=Inf),alpha=0.4)+geom_line(aes(y=yaw,col="yaw"),alpha=1,na.rm=T)+geom_line(aes(y=pitch,col="pitch"),alpha=1)+xlab("Time (seconds)")+coord_cartesian(ylim = c(-200, 200))+ylab("Rotation (degrees)")+ggtitle("Pitch and yaw against time for parallel reality\nscenario with transitions 1-3 for participant 8")+theme(legend.position="top",legend.title=element_blank(), axis.text=element_text(size=24),axis.title=element_text(size=26),legend.text=element_text(size=24),plot.title=element_text(size=28))+scale_fill_manual(values=rectPal,breaks=c("real","transition 1","transition 2","transition 3"))+scale_colour_manual(values=linePal)+scale_y_continuous(breaks = c(-180,-135,-90,-45,0,45,90,135,180))+scale_x_continuous(breaks=c(0,50,100,150,200,250),limits=c(0,285)) | |
| p132 <- ggplot(res13_combined_data,aes(time))+geom_rect(data=res13_combined_data,aes(xmin=(round(time*2))/2,xmax=(((round(time*2))/2)+0.5),fill=environment,ymin=0,ymax=Inf),alpha=0.4)+geom_line(aes(y=delta_x,col="East/West"),alpha=1)+geom_line(aes(y=delta_z,col="North/South"),alpha=1)+xlab("Time (seconds)")+ylab("Distance (meters)")+ggtitle("Distance moved against time for parallel reality\nscenario with transitions 1-3 for participant 8")+theme(legend.position="top",legend.title=element_blank(), axis.text=element_text(size=24),axis.title=element_text(size=26),legend.text=element_text(size=24),plot.title=element_text(size=28))+scale_fill_manual(values=rectPal)+scale_colour_manual(values=linePal)+scale_x_continuous(breaks=c(0,50,100,150,200,250),limits=c(0,285)) | |
| gA <- ggplotGrob(p131) | |
| gB <- ggplotGrob(p132) | |
| maxWidth = grid::unit.pmax(gA$widths[2:5], gB$widths[2:5]) | |
| gA$widths[2:5] <- as.list(maxWidth) | |
| gB$widths[2:5] <- as.list(maxWidth) | |
| png("8-1-3_2up.png",width=1920,height=1356) | |
| grid.arrange(gA, gB, ncol=1) | |
| dev.off() | |
| ################################################################################ | |
| p1 <- ggplot(resStat_combined_data,aes(time))+geom_line(aes(y=frames,col="framerate"),alpha=1,na.rm=T)+xlab("Time (seconds)")+ylim(0,75)+ylab("Framerate (frames per second)")+ggtitle("Framerate against time for seated VR\nscenario for participant 8")+theme(legend.position="top",legend.title=element_blank(), axis.text=element_text(size=24),axis.title=element_text(size=26),legend.text=element_text(size=24),plot.title=element_text(size=28))+scale_x_continuous(breaks=c(0,25,50,75,100,125,150,175),limits=c(0,175))+scale_colour_manual(values=linePal) | |
| p2 <- ggplot(res13_combined_data,aes(time))+geom_rect(data=res13_combined_data,aes(xmin=(round(time*2))/2,xmax=(((round(time*2))/2)+0.5),fill=environment,ymin=0,ymax=Inf),alpha=0.4)+ylim(0,75)+geom_line(aes(y=frames,col="framerate"),alpha=1)+xlab("Time (seconds)")+ylab("Framerate (frames per second)")+ggtitle("Framerate against time for parallel reality\nscenario with transitions 1-3 for participant 8")+theme(legend.position="top",legend.title=element_blank(), axis.text=element_text(size=24),axis.title=element_text(size=26),legend.text=element_text(size=24),plot.title=element_text(size=28))+scale_fill_manual(values=rectPal)+scale_colour_manual(values=linePal)+scale_x_continuous(breaks=c(0,50,100,150,200,250),limits=c(0,285)) | |
| p3 <- ggplot(res14_combined_data,aes(time))+geom_rect(data=res14_combined_data,aes(xmin=time,xmax=time+0.1,fill=environment,ymin=0,ymax=Inf),alpha=0.4)+geom_line(aes(y=frames,col="framerate"),alpha=1,na.rm=T)+ylim(0,75)+xlim(0,225)+xlab("Time (seconds)")+ylab("Framerate (frames per second)")+ggtitle("Framerate against time for PR scenario\nwith transitions 1-4 for participant 8")+theme(legend.position="top",legend.title=element_blank(), axis.text=element_text(size=24),axis.title=element_text(size=26),legend.text=element_text(size=24),plot.title=element_text(size=28))+scale_fill_manual(values=rectPal,breaks=c("real","transition 1","transition 2","transition 3","transition 4"))+scale_colour_manual(values=linePal) | |
| gp1 <- ggplotGrob(p1) | |
| gp2 <- ggplotGrob(p2) | |
| gp3 <- ggplotGrob(p3) | |
| maxWidth = grid::unit.pmax(gp1$widths[2:5], gp2$widths[2:5], gp3$widths[2:5]) | |
| maxHeight = grid::unit.pmax(gp1$heights[2:5], gp2$heights[2:5], gp3$heights[2:5]) | |
| gp1$widths[2:5] <- as.list(maxWidth) | |
| gp2$widths[2:5] <- as.list(maxWidth) | |
| gp3$widths[2:5] <- as.list(maxWidth) | |
| gp1$heights[2:5] <- as.list(maxHeight) | |
| gp2$heights[2:5] <- as.list(maxHeight) | |
| gp3$heights[2:5] <- as.list(maxHeight) | |
| png("8-frames-3-up.png",width=1920,height=1920) | |
| grid.arrange(gp1, gp2, gp3, ncol=1) | |
| dev.off() | |
| ################################################################################ | |
| # STATS | |
| ################################################################################ | |
| # only want to calculate stats over the time frames that we've plotted | |
| res13_combined_data <- subset(res13_combined_data,time>0&time<=285) | |
| res14_combined_data <- subset(res14_combined_data,time>0&time<=225) | |
| resStat_combined_data <- subset(resStat_combined_data,time>0&time<=175) | |
| ################################################################################ | |
| res13_combined_data$chunks <- c(NA, res13_combined_data$environment[-1] != head(res13_combined_data$environment,-1)) | |
| res13_combined_data$chunks <- cumsum(c(1, res13_combined_data$chunks[-1])) | |
| res14_combined_data$chunks <- c(NA, res14_combined_data$environment[-1] != head(res14_combined_data$environment,-1)) | |
| res14_combined_data$chunks <- cumsum(c(1, res14_combined_data$chunks[-1])) | |
| ################################################################################ | |
| mean(resStat_combined_data$frames) | |
| sd(resStat_combined_data$frames) | |
| mean(res13_combined_data$frames) | |
| sd(res13_combined_data$frames) | |
| mean(res14_combined_data$frames) | |
| sd(res14_combined_data$frames) | |
| ################################################################################ | |
| # standard deviation in pitch & yaw overall for stationary | |
| sd(resStat_combined_data$pitch) | |
| sd(resStat_combined_data$yaw) | |
| # standard deviation in pitch & yaw overall for 1-3 | |
| sd(res13_combined_data$pitch) | |
| sd(res13_combined_data$yaw) | |
| # standard deviation in pitch & yaw overall for 1-4 | |
| sd(res14_combined_data$pitch) | |
| sd(res14_combined_data$yaw) | |
| ################################################################################ | |
| # comparison of deviation in pitch & yaw overall between traditional & 1-3 & 1-4 | |
| var.test(resStat_combined_data$pitch, res13_combined_data$pitch) | |
| var.test(resStat_combined_data$pitch, res14_combined_data$pitch) | |
| var.test(res13_combined_data$pitch, res14_combined_data$pitch) | |
| var.test(resStat_combined_data$yaw, res13_combined_data$yaw) | |
| var.test(resStat_combined_data$yaw, res14_combined_data$yaw) | |
| var.test(res13_combined_data$yaw, res14_combined_data$yaw) | |
| ################################################################################ | |
| # comparison of deviation in pitch & deviation in yaw within each scenario | |
| var.test(resStat_combined_data$pitch, resStat_combined_data$yaw) | |
| var.test(res13_combined_data$pitch, res13_combined_data$yaw) | |
| var.test(res14_combined_data$pitch, res14_combined_data$yaw) | |
| ################################################################################ | |
| # weighted mean of standard deviation in yaw for each environment (1-3) | |
| weighted.mean(by((subset(res13_combined_data,environment=="real"))$yaw,(subset(res13_combined_data,environment=="real"))$chunks,sd), (by((subset(res13_combined_data,environment=="real")), (subset(res13_combined_data,environment=="real"))$chunks, nrow)),na.rm=T) | |
| weighted.mean(by((subset(res13_combined_data,environment!="real"))$yaw,(subset(res13_combined_data,environment!="real"))$chunks,sd), (by((subset(res13_combined_data,environment!="real")), (subset(res13_combined_data,environment!="real"))$chunks, nrow)),na.rm=T) | |
| weighted.mean(by((subset(res13_combined_data,environment=="transition 1"))$yaw,(subset(res13_combined_data,environment=="transition 1"))$chunks,sd), (by((subset(res13_combined_data,environment=="transition 1")), (subset(res13_combined_data,environment=="transition 1"))$chunks, nrow)),na.rm=T) | |
| weighted.mean(by((subset(res13_combined_data,environment=="transition 2"))$yaw,(subset(res13_combined_data,environment=="transition 2"))$chunks,sd), (by((subset(res13_combined_data,environment=="transition 2")), (subset(res13_combined_data,environment=="transition 2"))$chunks, nrow)),na.rm=T) | |
| weighted.mean(by((subset(res13_combined_data,environment=="transition 3"))$yaw,(subset(res13_combined_data,environment=="transition 3"))$chunks,sd), (by((subset(res13_combined_data,environment=="transition 3")), (subset(res13_combined_data,environment=="transition 3"))$chunks, nrow)),na.rm=T) | |
| ################################################################################ | |
| # weighted mean of standard deviation in yaw for each environment (1-4) | |
| weighted.mean(by((subset(res14_combined_data,environment=="real"))$yaw,(subset(res14_combined_data,environment=="real"))$chunks,sd), (by((subset(res14_combined_data,environment=="real")), (subset(res14_combined_data,environment=="real"))$chunks, nrow)),na.rm=T) | |
| weighted.mean(by((subset(res14_combined_data,environment!="real"))$yaw,(subset(res14_combined_data,environment!="real"))$chunks,sd), (by((subset(res14_combined_data,environment!="real")), (subset(res14_combined_data,environment!="real"))$chunks, nrow)),na.rm=T) | |
| weighted.mean(by((subset(res14_combined_data,environment=="transition 1"))$yaw,(subset(res14_combined_data,environment=="transition 1"))$chunks,sd), (by((subset(res14_combined_data,environment=="transition 1")), (subset(res14_combined_data,environment=="transition 1"))$chunks, nrow)),na.rm=T) | |
| weighted.mean(by((subset(res14_combined_data,environment=="transition 2"))$yaw,(subset(res14_combined_data,environment=="transition 2"))$chunks,sd), (by((subset(res14_combined_data,environment=="transition 2")), (subset(res14_combined_data,environment=="transition 2"))$chunks, nrow)),na.rm=T) | |
| weighted.mean(by((subset(res14_combined_data,environment=="transition 3"))$yaw,(subset(res14_combined_data,environment=="transition 3"))$chunks,sd), (by((subset(res14_combined_data,environment=="transition 3")), (subset(res14_combined_data,environment=="transition 3"))$chunks, nrow)),na.rm=T) | |
| weighted.mean(by((subset(res14_combined_data,environment=="transition 4"))$yaw,(subset(res14_combined_data,environment=="transition 4"))$chunks,sd), (by((subset(res14_combined_data,environment=="transition 4")), (subset(res14_combined_data,environment=="transition 4"))$chunks, nrow)),na.rm=T) | |
| ################################################################################ | |
| # weighted mean deviation in pitch & yaw between real & virtual parts of 1-3 | |
| weighted.mean(by((subset(res13_combined_data,environment=="real"))$pitch,(subset(res13_combined_data,environment=="real"))$chunks,sd), (by((subset(res13_combined_data,environment=="real")), (subset(res13_combined_data,environment=="real"))$chunks, nrow)),na.rm=T) | |
| weighted.mean(by((subset(res13_combined_data,environment!="real"))$pitch,(subset(res13_combined_data,environment!="real"))$chunks,sd), (by((subset(res13_combined_data,environment!="real")), (subset(res13_combined_data,environment!="real"))$chunks, nrow)),na.rm=T) | |
| weighted.mean(by((subset(res13_combined_data,environment=="real"))$yaw,(subset(res13_combined_data,environment=="real"))$chunks,sd), (by((subset(res13_combined_data,environment=="real")), (subset(res13_combined_data,environment=="real"))$chunks, nrow)),na.rm=T) | |
| weighted.mean(by((subset(res13_combined_data,environment!="real"))$yaw,(subset(res13_combined_data,environment!="real"))$chunks,sd), (by((subset(res13_combined_data,environment!="real")), (subset(res13_combined_data,environment!="real"))$chunks, nrow)),na.rm=T) | |
| # weighted mean deviation in pitch & yaw between real & virtual parts of 1-4 | |
| weighted.mean(by((subset(res14_combined_data,environment=="real"))$pitch,(subset(res14_combined_data,environment=="real"))$chunks,sd), (by((subset(res14_combined_data,environment=="real")), (subset(res14_combined_data,environment=="real"))$chunks, nrow)),na.rm=T) | |
| weighted.mean(by((subset(res14_combined_data,environment!="real"))$pitch,(subset(res14_combined_data,environment!="real"))$chunks,sd), (by((subset(res14_combined_data,environment!="real")), (subset(res14_combined_data,environment!="real"))$chunks, nrow)),na.rm=T) | |
| weighted.mean(by((subset(res14_combined_data,environment=="real"))$yaw,(subset(res14_combined_data,environment=="real"))$chunks,sd), (by((subset(res14_combined_data,environment=="real")), (subset(res14_combined_data,environment=="real"))$chunks, nrow)),na.rm=T) | |
| weighted.mean(by((subset(res14_combined_data,environment!="real"))$yaw,(subset(res14_combined_data,environment!="real"))$chunks,sd), (by((subset(res14_combined_data,environment!="real")), (subset(res14_combined_data,environment!="real"))$chunks, nrow)),na.rm=T) | |
| ################################################################################ | |
| # weighted mean deviation in pitch & yaw between traditional & virtual parts of 1-3 | |
| sd(resStat_combined_data$pitch) | |
| weighted.mean(by((subset(res13_combined_data,environment!="real"))$pitch,(subset(res13_combined_data,environment!="real"))$chunks,sd), (by((subset(res13_combined_data,environment!="real")), (subset(res13_combined_data,environment!="real"))$chunks, nrow)),na.rm=T) | |
| weighted.mean(by((subset(res14_combined_data,environment!="real"))$pitch,(subset(res14_combined_data,environment!="real"))$chunks,sd), (by((subset(res14_combined_data,environment!="real")), (subset(res14_combined_data,environment!="real"))$chunks, nrow)),na.rm=T) | |
| sd(resStat_combined_data$yaw) | |
| weighted.mean(by((subset(res13_combined_data,environment!="real"))$yaw,(subset(res13_combined_data,environment!="real"))$chunks,sd), (by((subset(res13_combined_data,environment!="real")), (subset(res13_combined_data,environment!="real"))$chunks, nrow)),na.rm=T) | |
| weighted.mean(by((subset(res14_combined_data,environment!="real"))$yaw,(subset(res14_combined_data,environment!="real"))$chunks,sd), (by((subset(res14_combined_data,environment!="real")), (subset(res14_combined_data,environment!="real"))$chunks, nrow)),na.rm=T) | |
| ################################################################################ | |
| # mean, sd & total for length of time spent in each environment | |
| ################################################################################ | |
| # 1-3 | |
| # real | |
| mean((by((subset(res13_combined_data,environment=="real"))$time, (subset(res13_combined_data,environment=="real"))$chunks, tail, n=1)) - (by((subset(res13_combined_data,environment=="real"))$time, (subset(res13_combined_data,environment=="real"))$chunks, head, n=1))) | |
| sd((by((subset(res13_combined_data,environment=="real"))$time, (subset(res13_combined_data,environment=="real"))$chunks, tail, n=1)) - (by((subset(res13_combined_data,environment=="real"))$time, (subset(res13_combined_data,environment=="real"))$chunks, head, n=1))) | |
| sum((by((subset(res13_combined_data,environment=="real"))$time, (subset(res13_combined_data,environment=="real"))$chunks, tail, n=1)) - (by((subset(res13_combined_data,environment=="real"))$time, (subset(res13_combined_data,environment=="real"))$chunks, head, n=1))) | |
| # !real | |
| mean((by((subset(res13_combined_data,environment!="real"))$time, (subset(res13_combined_data,environment!="real"))$chunks, tail, n=1)) - (by((subset(res13_combined_data,environment!="real"))$time, (subset(res13_combined_data,environment!="real"))$chunks, head, n=1))) | |
| sd((by((subset(res13_combined_data,environment!="real"))$time, (subset(res13_combined_data,environment!="real"))$chunks, tail, n=1)) - (by((subset(res13_combined_data,environment!="real"))$time, (subset(res13_combined_data,environment!="real"))$chunks, head, n=1))) | |
| sum((by((subset(res13_combined_data,environment!="real"))$time, (subset(res13_combined_data,environment!="real"))$chunks, tail, n=1)) - (by((subset(res13_combined_data,environment!="real"))$time, (subset(res13_combined_data,environment!="real"))$chunks, head, n=1))) | |
| # transition 1 | |
| mean((by((subset(res13_combined_data,environment=="transition 1"))$time, (subset(res13_combined_data,environment=="transition 1"))$chunks, tail, n=1)) - (by((subset(res13_combined_data,environment=="transition 1"))$time, (subset(res13_combined_data,environment=="transition 1"))$chunks, head, n=1))) | |
| sd((by((subset(res13_combined_data,environment=="transition 1"))$time, (subset(res13_combined_data,environment=="transition 1"))$chunks, tail, n=1)) - (by((subset(res13_combined_data,environment=="transition 1"))$time, (subset(res13_combined_data,environment=="transition 1"))$chunks, head, n=1))) | |
| sum((by((subset(res13_combined_data,environment=="transition 1"))$time, (subset(res13_combined_data,environment=="transition 1"))$chunks, tail, n=1)) - (by((subset(res13_combined_data,environment=="transition 1"))$time, (subset(res13_combined_data,environment=="transition 1"))$chunks, head, n=1))) | |
| # transition 2 | |
| mean((by((subset(res13_combined_data,environment=="transition 2"))$time, (subset(res13_combined_data,environment=="transition 2"))$chunks, tail, n=1)) - (by((subset(res13_combined_data,environment=="transition 2"))$time, (subset(res13_combined_data,environment=="transition 2"))$chunks, head, n=1))) | |
| sd((by((subset(res13_combined_data,environment=="transition 2"))$time, (subset(res13_combined_data,environment=="transition 2"))$chunks, tail, n=1)) - (by((subset(res13_combined_data,environment=="transition 2"))$time, (subset(res13_combined_data,environment=="transition 2"))$chunks, head, n=1))) | |
| sum((by((subset(res13_combined_data,environment=="transition 2"))$time, (subset(res13_combined_data,environment=="transition 2"))$chunks, tail, n=1)) - (by((subset(res13_combined_data,environment=="transition 2"))$time, (subset(res13_combined_data,environment=="transition 2"))$chunks, head, n=1))) | |
| # transition 3 | |
| mean((by((subset(res13_combined_data,environment=="transition 3"))$time, (subset(res13_combined_data,environment=="transition 3"))$chunks, tail, n=1)) - (by((subset(res13_combined_data,environment=="transition 3"))$time, (subset(res13_combined_data,environment=="transition 3"))$chunks, head, n=1))) | |
| sd((by((subset(res13_combined_data,environment=="transition 3"))$time, (subset(res13_combined_data,environment=="transition 3"))$chunks, tail, n=1)) - (by((subset(res13_combined_data,environment=="transition 3"))$time, (subset(res13_combined_data,environment=="transition 3"))$chunks, head, n=1))) | |
| sum((by((subset(res13_combined_data,environment=="transition 3"))$time, (subset(res13_combined_data,environment=="transition 3"))$chunks, tail, n=1)) - (by((subset(res13_combined_data,environment=="transition 3"))$time, (subset(res13_combined_data,environment=="transition 3"))$chunks, head, n=1))) | |
| ################################################################################ | |
| # 1-4 | |
| # real | |
| mean((by((subset(res14_combined_data,environment=="real"))$time, (subset(res14_combined_data,environment=="real"))$chunks, tail, n=1)) - (by((subset(res14_combined_data,environment=="real"))$time, (subset(res14_combined_data,environment=="real"))$chunks, head, n=1))) | |
| sd((by((subset(res14_combined_data,environment=="real"))$time, (subset(res14_combined_data,environment=="real"))$chunks, tail, n=1)) - (by((subset(res14_combined_data,environment=="real"))$time, (subset(res14_combined_data,environment=="real"))$chunks, head, n=1))) | |
| sum((by((subset(res14_combined_data,environment=="real"))$time, (subset(res14_combined_data,environment=="real"))$chunks, tail, n=1)) - (by((subset(res14_combined_data,environment=="real"))$time, (subset(res14_combined_data,environment=="real"))$chunks, head, n=1))) | |
| # !real | |
| mean((by((subset(res14_combined_data,environment!="real"))$time, (subset(res14_combined_data,environment!="real"))$chunks, tail, n=1)) - (by((subset(res14_combined_data,environment!="real"))$time, (subset(res14_combined_data,environment!="real"))$chunks, head, n=1))) | |
| sd((by((subset(res14_combined_data,environment!="real"))$time, (subset(res14_combined_data,environment!="real"))$chunks, tail, n=1)) - (by((subset(res14_combined_data,environment!="real"))$time, (subset(res14_combined_data,environment!="real"))$chunks, head, n=1))) | |
| sum((by((subset(res14_combined_data,environment!="real"))$time, (subset(res14_combined_data,environment!="real"))$chunks, tail, n=1)) - (by((subset(res14_combined_data,environment!="real"))$time, (subset(res14_combined_data,environment!="real"))$chunks, head, n=1))) | |
| # !real & !transition 4 | |
| mean((by((subset(res14_combined_data,environment!="real"&environment!="transition 4"))$time, (subset(res14_combined_data,environment!="real"&environment!="transition 4"))$chunks, tail, n=1)) - (by((subset(res14_combined_data,environment!="real"&environment!="transition 4"))$time, (subset(res14_combined_data,environment!="real"&environment!="transition 4"))$chunks, head, n=1))) | |
| sd((by((subset(res14_combined_data,environment!="real"&environment!="transition 4"))$time, (subset(res14_combined_data,environment!="real"&environment!="transition 4"))$chunks, tail, n=1)) - (by((subset(res14_combined_data,environment!="real"&environment!="transition 4"))$time, (subset(res14_combined_data,environment!="real"&environment!="transition 4"))$chunks, head, n=1))) | |
| sum((by((subset(res14_combined_data,environment!="real"&environment!="transition 4"))$time, (subset(res14_combined_data,environment!="real"&environment!="transition 4"))$chunks, tail, n=1)) - (by((subset(res14_combined_data,environment!="real"&environment!="transition 4"))$time, (subset(res14_combined_data,environment!="real"&environment!="transition 4"))$chunks, head, n=1))) | |
| # transition 1 | |
| mean((by((subset(res14_combined_data,environment=="transition 1"))$time, (subset(res14_combined_data,environment=="transition 1"))$chunks, tail, n=1)) - (by((subset(res14_combined_data,environment=="transition 1"))$time, (subset(res14_combined_data,environment=="transition 1"))$chunks, head, n=1))) | |
| sd((by((subset(res14_combined_data,environment=="transition 1"))$time, (subset(res14_combined_data,environment=="transition 1"))$chunks, tail, n=1)) - (by((subset(res14_combined_data,environment=="transition 1"))$time, (subset(res14_combined_data,environment=="transition 1"))$chunks, head, n=1))) | |
| sum((by((subset(res14_combined_data,environment=="transition 1"))$time, (subset(res14_combined_data,environment=="transition 1"))$chunks, tail, n=1)) - (by((subset(res14_combined_data,environment=="transition 1"))$time, (subset(res14_combined_data,environment=="transition 1"))$chunks, head, n=1))) | |
| # transition 2 | |
| mean((by((subset(res14_combined_data,environment=="transition 2"))$time, (subset(res14_combined_data,environment=="transition 2"))$chunks, tail, n=1)) - (by((subset(res14_combined_data,environment=="transition 2"))$time, (subset(res14_combined_data,environment=="transition 2"))$chunks, head, n=1))) | |
| sd((by((subset(res14_combined_data,environment=="transition 2"))$time, (subset(res14_combined_data,environment=="transition 2"))$chunks, tail, n=1)) - (by((subset(res14_combined_data,environment=="transition 2"))$time, (subset(res14_combined_data,environment=="transition 2"))$chunks, head, n=1))) | |
| sum((by((subset(res14_combined_data,environment=="transition 2"))$time, (subset(res14_combined_data,environment=="transition 2"))$chunks, tail, n=1)) - (by((subset(res14_combined_data,environment=="transition 2"))$time, (subset(res14_combined_data,environment=="transition 2"))$chunks, head, n=1))) | |
| # transition 3 | |
| mean((by((subset(res14_combined_data,environment=="transition 3"))$time, (subset(res14_combined_data,environment=="transition 3"))$chunks, tail, n=1)) - (by((subset(res14_combined_data,environment=="transition 3"))$time, (subset(res14_combined_data,environment=="transition 3"))$chunks, head, n=1))) | |
| sd((by((subset(res14_combined_data,environment=="transition 3"))$time, (subset(res14_combined_data,environment=="transition 3"))$chunks, tail, n=1)) - (by((subset(res14_combined_data,environment=="transition 3"))$time, (subset(res14_combined_data,environment=="transition 3"))$chunks, head, n=1))) | |
| sum((by((subset(res14_combined_data,environment=="transition 3"))$time, (subset(res14_combined_data,environment=="transition 3"))$chunks, tail, n=1)) - (by((subset(res14_combined_data,environment=="transition 3"))$time, (subset(res14_combined_data,environment=="transition 3"))$chunks, head, n=1))) | |
| # transition 4 | |
| mean((by((subset(res14_combined_data,environment=="transition 4"))$time, (subset(res14_combined_data,environment=="transition 4"))$chunks, tail, n=1)) - (by((subset(res14_combined_data,environment=="transition 4"))$time, (subset(res14_combined_data,environment=="transition 4"))$chunks, head, n=1))) | |
| sd((by((subset(res14_combined_data,environment=="transition 4"))$time, (subset(res14_combined_data,environment=="transition 4"))$chunks, tail, n=1)) - (by((subset(res14_combined_data,environment=="transition 4"))$time, (subset(res14_combined_data,environment=="transition 4"))$chunks, head, n=1))) | |
| sum((by((subset(res14_combined_data,environment=="transition 4"))$time, (subset(res14_combined_data,environment=="transition 4"))$chunks, tail, n=1)) - (by((subset(res14_combined_data,environment=="transition 4"))$time, (subset(res14_combined_data,environment=="transition 4"))$chunks, head, n=1))) | |
| ################################################################################ | |
| # number of transitions into virtual | |
| # 1-3 | |
| length(unique(subset(res13_combined_data, res13_combined_data$environment=="real")$chunks)) | |
| length(unique(subset(res13_combined_data, res13_combined_data$environment!="real")$chunks)) | |
| length(unique(subset(res13_combined_data, res13_combined_data$environment=="transition 1")$chunks)) | |
| length(unique(subset(res13_combined_data, res13_combined_data$environment=="transition 2")$chunks)) | |
| length(unique(subset(res13_combined_data, res13_combined_data$environment=="transition 3")$chunks)) | |
| # 1-4 | |
| length(unique(subset(res14_combined_data, res14_combined_data$environment=="real")$chunks)) | |
| length(unique(subset(res14_combined_data, res14_combined_data$environment!="real")$chunks)) | |
| length(unique(subset(res14_combined_data, res14_combined_data$environment=="transition 1")$chunks)) | |
| length(unique(subset(res14_combined_data, res14_combined_data$environment=="transition 2")$chunks)) | |
| length(unique(subset(res14_combined_data, res14_combined_data$environment=="transition 3")$chunks)) | |
| length(unique(subset(res14_combined_data, res14_combined_data$environment=="transition 4")$chunks)) | |
| ################################################################################ |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment