phenotypic_data_scaled <- scale(as.matrix(AllData[,c(4:52)])) # scale all phenotypic data around 0, between [-4,+4] phenotypic_data_scaled[1:10,1:10] #double check ...