First of all, the function obtains a vector with the SNPs sorted using the data frame with the genomic
positions (see 'dataSNPs.pos' argument).
Then, the columns which indicate where the information about the
genotypes is in our data frame, are sorted using this vector.
This information is useful when WGassociation
function is called since it
allow the user to summaryze the results with the SNPs sorted by genomic position