#R code for generating plot of mds


#Reads data in from plink.mds and stores it in data

read.table("plink.mds",header=T)->data


#Plots the 1st two axis of variation
plot(data$C1,data$C2,col=data$SOL+3,xlab="MDS component 1", ylab="MDS component 2",main="First two MDS components")