Hi, I'm not sure how use curve(dexGAUS(… None of the following four works: rt<- rexGAUS(100, mu=300, nu=100, sigma=35) m1<-gamlss(rt~1, family=exGAUS) curve(dexGAUS(rt=x, mu=300 ,sigma=35,nu=100), 100, 600, main = "The ex- GAUS density mu=300 ,sigma=35,nu=100") curve(dexGAUS(x=rt, mu=300 ,sigma=35,nu=100), 100, 600, main = "The ex- GAUS density mu=300 ,sigma=35,nu=100") curve(dexGAUS(y=rt, mu=300 ,sigma=35,nu=100), 100, 600, main = "The ex- GAUS density mu=300 ,sigma=35,nu=100") curve(dexGAUS(rt=y, mu=300 ,sigma=35,nu=100), 100, 600, main = "The ex- GAUS density mu=300 ,sigma=35,nu=100") But this does: y<- rexGAUS(100, mu=300, nu=100, sigma=35) m1<-gamlss(y~1, family=exGAUS) curve(dexGAUS(y=x, mu=300 ,sigma=35,nu=100), 100, 600, main = "The ex- GAUS density mu=300 ,sigma=35,nu=100") _____________________________ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400 Charlottesville, VA 22904-4400 Parcels: Room 102 Gilmer Hall McCormick Road Charlottesville, VA 22903 Office: B011 +1-434-982-4729 Lab: B019 +1-434-982-4751 Fax: +1-434-982-4766 WWW: http://www.people.virginia.edu/~mk9y/ [[alternative HTML version deleted]]