##run lordif first
age.dif <- lordif(Anxiety[paste("R",1:29,sep="")],Anxiety$age)
plotDIF(age.dif,labels=c("Younger","Older"))
age.dif.MC <- montecarlo(age.dif,alpha=.05,nr=500)
plotMC(age.dif.MC,mfrow=c(1,1))
Run the code above in your browser using DataLab