###     from package abd  and 
###     Broadbent, N.J., L.R. Squire, and R.E. Clark. 2004. Spatial
###     memory, recognition memory, and the hippocampus. _Proceedings of
###     the National Academy of Sciences (USA)_ 101: 14515-14520.

rm(list=ls())


D = read.csv("HippocampusLesions.csv",head=TRUE)
lesion = D[,1]
memory = D[,2]

o = order(lesion)
memory = memory[o]
lesion = lesion[o]


pdf("Hippocampus.pdf")

plot(lesion,memory,pch=21,cex=2,bg="pink",xlab="",ylab="",cex.axis=1.4)
out = lm(memory ~ lesion)
abline(out,lwd=4,col="purple")
title(ylab="Memory", line=2.4, cex.lab=2)
title(xlab="Lesions",line=3, cex.lab=2)

dev.off()