#1. create data x = rnorm(1000, 0, 4) y = rnorm(1000, 0, 1) #2. calculate means and euclidean distance meanx = mean(x) meany = mean(y) eucliddist = sqrt((x-meanx)*(x-meanx) + (y-meany)*(y-meany)) #3.calculate mahalanobis distance data = cbind(x,y) datamean = cbind(meanx, meany) datacov = cov(data) mahaldist = mahalanobis(data, datamean, datacov) #4. normalize the data x_norm = (x - min(x))/(max(x)-min(x)) y_norm = (y - min(y))/(max(y)-min(y)) #5. calculate euclidean distance of normalized data meanx_norm = mean(x_norm) meany_norm = mean(y_norm) eucliddist_norm = sqrt((x_norm-meanx_norm)*(x_norm-meanx_norm) + (y_norm-meany_norm)*(y_norm-meany_norm)) #6. calculate mahalanobis distance of normalized data data_norm = cbind(x_norm,y_norm) datamean_norm = cbind(meanx_norm, meany_norm) datacov_norm = cov(data_norm) mahaldist_norm = mahalanobis(data_norm, datamean_norm, datacov_norm) #7. Exportieren der Daten im csv Format exportdata = cbind(x, y, eucliddist, mahaldist, x_norm, y_norm, eucliddist_norm, mahaldist_norm) write.table(exportdata, "mahalanobis1.csv", row.names=FALSE, col.names=c("x", "y", "eucliddist", "mahaldist", "x_norm", "y_norm", "eucliddist_norm", "mahaldist_norm"), sep=";")