daijiajie1@1: makeDriftGraph <- function(data, singermode) {
daijiajie1@1: 	
daijiajie1@1: 	subdata <- data[data$singermode == singermode, ]
daijiajie1@1: 	
daijiajie1@1: 	meanpitch <- median(subdata$pitchmedian-subdata$semitone)
daijiajie1@1: 	print(meanpitch)
daijiajie1@1: 	
daijiajie1@1: 	plot(I(subdata$pitchmedian - subdata$semitone) ~ I(subdata$run * 25 + subdata$noteid), 
daijiajie1@1: 	     pch = "+",
daijiajie1@1:          ylim = meanpitch + c(-1.5,1.5),
daijiajie1@1:          xlab = "note number",
daijiajie1@1:          ylab = "pitch after score subtraction")
daijiajie1@1: 	
daijiajie1@1: 	fit <- lm(I(subdata$pitchmedian - subdata$semitone) ~ I(subdata$run * 25 + subdata$noteid))
daijiajie1@1: 	abline(fit)
daijiajie1@1: }