powered by
spmHrf(RT, p = c(6, 16, 1, 1, 6, 0, 32))
RT<-1/16 can.hrf<-spmHrf(RT)[[1]] x<-seq(0,32,RT) plot(x,can.hrf,lwd=7,col="red",type="l", main="Canonical HRF",xlab="Time(seconds)", ylab="Intensity(A.U.)")
Run the code above in your browser using DataLab