R/DAMsleep3.R

In ueno-tr/flysleep: flysleep: Sleep analysis using DAM data.

# Analyze fly sleep based on DAM data
# 3日分のデータから1日ずつの睡眠のグラフをプロット

DAMsleep3 <- function(hrs,n) { #1-3日のうち、除外する日をnとして入れる
  library(readxl)
  summary <- read_excel("summary.xls")

  out <- meanSEM(hrs, summary)
  out2 <- transform(out, hr=c(1:24), day=c(rep("day1",24),rep("day2",24),rep("day3",24)))
  out3 <- subset(out2, day!=paste("day",n, sep="")) #n日目を除外する

  ##グラフに出力
  library(ggplot2)
  p <- ggplot(out3, aes_string(x = "hr", y = "mean", "group=interaction(names(out3)[2], day)", colour=names(out3)[2], linetype="day") ) + geom_line() + ylab("sleep (min/hr)")

  errors <- aes(ymax = mean + SEM, ymin = mean - SEM)
  p <- p + geom_errorbar(errors, width = 0.2) + geom_point(size = 2)
  p

  ##pngファイルに出力
  ggsave(file = "sleep3graph.png", dpi = 200, width =15 , height =10 )
}