library(sciplot) # load sciplot package ## hfcsDA<-read.csv("FSCV_summary_data.csv") hfcsDA$ID<-as.factor(hfcsDA$id) # change id into a factor hfcsDA$Hz<-as.factor(hfcsDA$Hz) # change Hz into a factor ## plot peak DA release ########################################### ## line plot par(mfrow=c(1,1),mar=c(5.8,8,1.4,2),lwd=4,mgp=c(4.5, 1, 0),las=1) lineplot.CI(Hz,ED,group=group,data=hfcsDA,cex.lab="2.5",cex.axis="2.25", ylim=c(0,150),xlab = "Frequency (Hz)", ylab = "Evoked [DA] release (nM)", type = "b",col = c("blue","red"), pch = c(16,16),lty=c(1,1),cex=4,lwd=8,err.lty = 1,bty='l', font.axis = 2,font.lab=2,legend=T)