OTOH looking at your model code again my change will allow prob>1 so rather than using exp(cvec) it appears you will need to transform prob with something like a logistic to make sense. -- David A. Fournier P.O. Box 2040, Sidney, B.C. V8l 3S3 Canada Phone/FAX 250-655-3364 http://otter-rsch.com