Am 27.03.2010 10:45, schrieb Philip Wong: > >> mah=mahalanobis(data,apply(data,2,mean),var(data)) >> > >> mah=mahalanobis(data,apply(data,2,mean),var(no_fisher)) >> > qqplot(qchisq(ppoints(data),ncol(data),data) >
As Joshua already pointed out: you are trying mathematical functions on names. Your Data containes e.g. Village names. You can exclude those with subset, generating a new data.frame. If you want averages by categories e.g. by a village have a look at the doby package, it has a function summaryby that is perfect for generating aggregated data with multiple functions. hth Stefan ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.