Advertisement
Peradam

MDLP K17

Aug 29th, 2014
992
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 4.21 KB | None | 0 0
  1. Scroll to the bottom for instructions for other calculators.
  2.  
  3. I'll assume that you're using Windows, but refer to DIYDodecad's README.txt file if you're running Linux. I don't think there's Mac OS support for DIYDodecad, but you might have luck using Wine.
  4.  
  5. Admixture:
  6.  
  7. 1. Download and install R: http://www.r-project.org/
  8. 2. Download DIYDodecad: http://dodecad.blogspot.co.uk/2011/09/do-it-yourself-dodecad-v-21.html
  9. 3. Download the MDLP K17 package: https://drive.google.com/file/d/0B6n7iMc2P-yQaWhVWTZkczJtLUk/edit?usp=sharing
  10. 4. Download your 23andMe raw data file: https://www.23andme.com/you/download/
  11. 5. Extract and place the DIYDodecad, K17 & 23andMe files into the same folder.
  12. 6. Open R and change the directory (File > Change Dir) to the folder where you placed your DIY, K17 & raw data files.
  13. 7. Enter this command in R:
  14. source('standardize.r')
  15.  
  16. 8. Enter this command but replace X with the name of your raw data file:
  17. standardize('X.txt', company='23andMe')
  18.  
  19. 9. Enter this command:
  20. system('DIYDodecadWin k17.par')
  21.  
  22. That's it. Just be patient, the results might take some time to generate.
  23.  
  24.  
  25. Oracle:
  26.  
  27. Once you've obtained your Admixture results, you can then calculate your Oracle results using the TOraclerelease.RData file and TOracle() function.
  28.  
  29. 1. Open the 'TOraclerelease.RData' file in R (File > Load Workspace, or double-clicking the file should do it).
  30. 2. Enter this command but replace X with your Admixture results in sequence, separated by commas and excluding the percentage signs:
  31. TOracle(c(X), k=20)
  32.  
  33. For example, my Admixture results are:
  34.  
  35. 25.22% Caucasian-Basal
  36. 0.37% Melano-Austronesian
  37. 1.28% Uralic
  38. 0.41% Ancestral_East_Siberian
  39. 0.86% Ancestral_West_Siberian
  40. 21.64% Ancestral_North_Indian
  41. 0.73% Ancestral_South_Indian
  42. 4.75% West_European_HG
  43. 0.01% Ancestral_Sami-Finnic
  44. 0.78% South_East_Asian
  45. 0.57% Archaic_African
  46. 22.55% Ancestral_Mediterranean_EEF
  47. 1.06% Circumpolar
  48. 17.65% Near-East-Basal
  49. 0.01% Ancestral_East_European_ANE
  50. 1.55% African_Sub_Saharian
  51. 0.55% Amerindian
  52.  
  53. So I would type:
  54.  
  55. TOracle(c(25.22, 0.37, 1.28, 0.41, 0.86, 21.64, 0.73, 4.75, 0.01, 0.78, 0.57, 22.55, 1.06, 17.65, 0.01, 1.55, 0.55), k=20)
  56.  
  57. That gives you Single Population Sharing results.
  58.  
  59. 3. For Mixed Mode Population Sharing, add 'mixedmode=T' to the command:
  60.  
  61. TOracle(c(X), k=20, mixedmode=T)
  62.  
  63. 4. If you want greater or less than 20 results, change k=20 to k=X, where X is the number of results you want to generate.
  64. For example:
  65.  
  66. TOracle(c(X), k=250, mixedmode=T)
  67.  
  68.  
  69. Using other MDLP, Dodecad, Eurogenes or Harappa calculators:
  70.  
  71. The same Admixture steps apply if you'd prefer to use a different calculator.
  72. Follow the steps in the Admixture section above, but ignore step 3 and download your chosen calculator package instead.
  73. In steps 5 & 6, "K17" refers to your chosen calculator.
  74. In step 9, replace k17.par with X.par, where X is the name of the par file included in your chosen calculator package.
  75.  
  76. If an Oracle RData file is included with your chosen calculator package, the Oracle instructions above will work if used with the correct function.
  77. The correct function differs with each calculator: MDLP K17 uses TOracle(), Dodecad K12b uses DodecadOracle() and HarappaWorld uses HarappaOracle()
  78. If the function to use isn't included in a README file with the package or on the project's website, open the RData file in R as normal and type the following command to locate the function:
  79. ls()
  80.  
  81.  
  82. ***Note:
  83. In step 8 of the Admixture section above, the "standardize" function creates the file "genotype.txt". Keep this file and DIYDodecadWin.exe for future use.
  84. The genotype.txt file allows you to skip a couple of steps and can be used for all DIY calculators released by MDLP, Dodecad, Eurogenes and Harappa projects.
  85.  
  86. 1. Download your chosen calculator package and extract it.
  87. 2. Copy the genotype.txt file, DIYDodecadWin.exe file and the package contents of your chosen calculator into the same folder.
  88. 3. Open R and change the directory (File > Change Dir) to the folder where you placed your genotype.txt, DIYDodecadWin.exe and calculator files.
  89. 4. Enter this command but replace X with the name of the par file included with your chosen calculator:
  90. system('DIYDodecadWin X.par')
  91.  
  92. That's it.
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement