library(ggplot2)
+library(data.table)
gs <- read.delim("student_information.tsv")
d <- gs[,c(2,5)]
colnames(d) <- c("student.num", "discord.name")
-call.list <- do.call("rbind", lapply(list.files(".", pattern="^call_list-.*tsv$"), function (x) {read.delim(x)}))
+call.list <- do.call("rbind", lapply(list.files(".", pattern="^call_list-.*tsv$"), function (x) {read.delim(x)[,1:4]}))
colnames(call.list) <- gsub("_", ".", colnames(call.list))
+call.list$day <- as.Date(call.list$timestamp)
+
+## drop calls where the person wasn't present
+call.list.full <- call.list
+call.list[!call.list$answered,]
+call.list <- call.list[call.list$answered,]
+
call.counts <- data.frame(table(call.list$discord.name))
colnames(call.counts) <- c("discord.name", "num.calls")
attendance <- unlist(lapply(list.files(".", pattern="^attendance-.*tsv$"), function (x) {d <- read.delim(x); strsplit(d[[2]], ",")}))
-attendance.counts <- data.frame(table(attendance))
+file.to.attendance.list <- function (x) {
+ tmp <- read.delim(x)
+ d.out <- data.frame(discord.name=unlist(strsplit(tmp[[2]], ",")))
+ d.out$day <- rep(as.Date(tmp[[1]][1]), nrow(d.out))
+ return(d.out)
+}
+
+attendance <- do.call("rbind",
+ lapply(list.files(".", pattern="^attendance-.*tsv$"),
+ file.to.attendance.list))
+
+## create list of folks who are missing in class
+missing.in.class <- call.list.full[is.na(call.list.full$answered) |
+ (!is.na(call.list.full$answered) & !call.list.full$answered),
+ c("discord.name", "day")]
+
+missing.in.class <- unique(missing.in.class)
+
+setDT(attendance)
+setkey(attendance, discord.name, day)
+setDT(missing.in.class)
+setkey(missing.in.class, discord.name, day)
+
+## drop presence for people on missing days
+attendance[missing.in.class,]
+attendance <- as.data.frame(attendance[!missing.in.class,])
+
+attendance.counts <- data.frame(table(attendance$discord.name))
colnames(attendance.counts) <- c("discord.name", "num.present")
-d <- merge(d, attendance.counts, all.x=TRUE, all.y=TRUE, by="discord.name"); d
+d <- merge(d, attendance.counts,
+ all.x=TRUE, all.y=TRUE,
+ by="discord.name")
+
+days.list <- lapply(unique(attendance$day), function (day) {
+ day.total <- table(call.list.full$day == day)[["TRUE"]]
+ lapply(d$discord.name, function (discord.name) {
+ num.present <- nrow(attendance[attendance$day == day & attendance$discord.name == discord.name,])
+ if (num.present/day.total > 1) {print(day)}
+ data.frame(discord.name=discord.name,
+ days.present=(num.present/day.total))
+ })
+})
+
+days.tmp <- do.call("rbind", lapply(days.list, function (x) do.call("rbind", x)))
+
+days.tbl <- tapply(days.tmp$days.present, days.tmp$discord.name, sum)
+
+attendance.days <- data.frame(discord.name=names(days.tbl),
+ days.present=days.tbl,
+ days.absent=length(list.files(".", pattern="^attendance-.*tsv$"))-days.tbl)
+
+d <- merge(d, attendance.days,
+ all.x=TRUE, all.y=TRUE, by="discord.name")
+
+d[sort.list(d$days.absent), c("discord.name", "num.calls", "days.absent")]
+
+## make some visualizations of whose here/not here
+#######################################################
+
+png("questions_absence_histogram_combined.png", units="px", width=800, height=600)
+
+ggplot(d) +
+ aes(x=as.factor(num.calls), fill=days.absent, group=days.absent) +
+ geom_bar(color="black") +
+ scale_x_discrete("Number of questions asked") +
+ scale_y_continuous("Number of students") +
+ scale_fill_continuous("Days absent", low="red", high="blue")+
+ theme_bw()
+
+dev.off()
+
+png("questions_absenses_boxplots.png", units="px", width=800, height=600)
+
+ggplot(data=d) +
+ aes(x=as.factor(num.calls), y=days.absent) +
+ geom_boxplot() +
+ scale_x_discrete("Number of questions asked") +
+ scale_y_continuous("Days absent")
+
+dev.off()