# RT test script vitality.fit(RT[,1],RT[,2],rc.data=T,se=75,gfit=75,pplot=1.3)