OK. Maybe 2 errors. I already found one: On 05/27/04 16:43, Jonathan Baron wrote: >Sumdiffs <- rowMeans(Mothers[,-(1:4)]-Daughters[,-(1:4)]) should be Sumdiffs <- rowMeans(Mothers[,-(1:4)])-rowMeans(Daughters[,-(1:4)])