June 21, 2008

EURO-DNA-CALC 1.1 released

UPDATE (20 Feb 2011): EURO-DNA-CALC has been superseded by the Dodecad Ancestry Project.

The current version is 1.1.1 Download it if you obtained previous version (1.1) which had a bug for deCODEme data.

A new version of the EURO-DNA-CALC has been released. The previous version (1.0) was described here.

Version 1.1 also uses deCODEme/23andMe data for an individual but outputs a maximum likelihood estimate of his admixture proportions from the three groups NW European, SE European, and Ashkenazi Jewish.

More comments on this dna-forums thread (registration required).

Below are the admixture estimates for Greg and Lily Mendel of 23andMe. An examination of the confidence intervals output by the program suggests that non-NW European admixture in Greg is a strong possibility, while Lily could conceivably have no non-NW European admixture. On the other hand, Lily's 10% Ashkenazi and only 1% SE European admixture imbalance is interesting. As with all genetic tools, the results should be interpreted in the light of other information about a person's ethnic origin and genealogy.

37 comments:

  1. Thanks for the new version. I get 9 NW, 91 SW and 0 AS. Being from Italy (tuscany), that seems on the mark. Which algorithm are you using now?

    cacio

    ReplyDelete
  2. Just an exhaustive search over all possible admixture estimates (x,y,z) where x,y,z>=0 and x+y+z=100.

    If NW,SE,AJ are the frequencies of an allele in the three populations, then an (x,y,z) admixed individual will possess the allele with probability (x*NW+y*SE+z*AJ)/100.

    ReplyDelete
  3. Dinekes, I've added your tool to the genomics secondary market database. What is your full name so that I may credit you with authorship, or if you would prefer an alias, under what name may I list you?

    -Andrew Yates
    ThinkGene.com

    ReplyDelete
  4. Does this apply to 23andme v2, or does it only work with v1?

    ReplyDelete
  5. Does this apply to 23andme v2, or does it only work with v1?

    It should work fine with v2.

    ReplyDelete
  6. My results:

    > EuroDNACalc("23andme")
    [1] "NORTHWEST EURO: Maximum Likelihood Estimate=100% Interval=[86, 100]"
    [1] "SOUTHEAST EURO: Maximum Likelihood Estimate=0% Interval=[0, 14]"
    [1] "ASHKENAZI JEWISH: Maximum Likelihood Estimate=0% Interval=[0, 9]"
    >

    ReplyDelete
  7. This Ponto's EuroDNACalc.

    SE Euro 87% [53,100]

    NW Euro 8% [0,32]

    Ashken 5% [0,28]

    I have no know Jews, Ashkenazim or otherwise in my family going back as far as I have checked,which is 1500s. Speaking of Jews, in my country of origin, Sephardic Jews were present not Ashkenazim.

    Who are the SE Euros and the NW Euros? I am not of Greek descent.

    ReplyDelete
  8. NWE and AJ confidence intervals interesect 0% so I wouldn't give them any particular importance.

    ReplyDelete
  9. I ran my dna file from 23andme thru your calc program.
    I got:
    NW Euro 65%
    SE Euro 15%
    Ashkenazi 20%
    My father is from Germany and and far as I know he is not Jewish.
    My mothers mother is Romanian (both parents from Romania)
    and my mothers father ancestors have been here for hundreds of years with a little native american and african in the mix (1/2 % each).
    My question is, is it accurate that I am really 20% Jewish?
    Thanks!

    ReplyDelete
  10. You should look at the confidence intervals and not just the central estimates. Both the minor components have probably intervals that intersect zero.

    ReplyDelete
  11. I didn't notice the intervals only the pie chart. I re ran it and here are the results:
    "NORTHWEST EURO: Maximum Likelihood Estimate=65% Interval=[38, 89]"
    [1] "SOUTHEAST EURO: Maximum Likelihood Estimate=15% Interval=[0, 52]"
    [1] "ASHKENAZI JEWISH: Maximum Likelihood Estimate=20% Interval=[0, 42]"
    As I said, my father is from Germany and my mother is half Romanian and the other half mostly European from the UK.
    Can you explain what the Intervals mean? Am I 20% Jewish?
    Thanks

    ReplyDelete
  12. It doesn't seem to work for me whenever. I have followed every instruction, up until the part where you open the 23andme txt file on the command prompt.

    I get this:
    "Error in file(file, "rt") : cannot open the connection
    In addition: Warning message:
    In file(file, "rt") :
    cannot open file '23andme.txt': No such file or directory"

    How do i fix this

    ReplyDelete
  13. Is there any way to reverse-engineer the Frankish-Byzantine-Levantine components from one's "Global Similarity" scores in 23andMe? I'm trying to figure out a multiple-regression set up with the use of many people's results.

    ReplyDelete
  14. I get
    Error: unexpected symbol in "23andme"

    Any instructions?

    Thanks,

    Steven

    ReplyDelete
  15. Genespool2009 :

    I'm also African American and no it does not work. It shows me to be 64% Southeast European, 0% Northwest European, and 36% Ash. Jew. LOL This does not at all coincide with 23&me, they say I'm at least 62% African, but likely more because they have serious issues overestimating "Asian" (really Native American) ancestry in New World Blacks.

    My Central European girlfriend is shown as having 21% Ash Jew, 40% NW Euro, 39% SE Euro. Based on where she is from and her families name and known time in the area I doubt she has any Jewish ancestry that is in the last 200-300 years anyway...the location of her country makes it typical that she would come out as an even split between NW and SE...so with that it appears accurate and matches up with 23&me fairly well as she is almost equally related to Southern and Northern Europeans.

    ReplyDelete
  16. On second thought, in regard to my results, it does pare up with 23&me if you make the SE Euro AFRICAN and the Ash Jew "Other" (equally Western Euro and Asian) lol

    ReplyDelete
  17. My results:
    > EuroDNACalc("23andme")
    [1] "NORTHWEST EURO: Maximum Likelihood Estimate=45% Interval=[17, 67]"
    [1] "SOUTHEAST EURO: Maximum Likelihood Estimate=7% Interval=[0, 45]"
    [1] "ASHKENAZI JEWISH: Maximum Likelihood Estimate=48% Interval=[23, 69]"

    On ne peut rien vous cacher.
    Merci et bonne poursuite de votre superbe travail.

    Jacques Beaugrand

    ReplyDelete
  18. I can report two results for the Euro DNA Calculator, both for two users on 23andme (images attached):

    Dan Draghici (myself)
    52% NWE, 48% SEE and 0% Ashkenazi Jewish

    Some surprise given that most of my known DNA comes from Romania, in spite of some German blood; however, my Y-DNA haplogroup is U106+ L48+ L47- L148-.

    David Blais
    84% NWE, 16% SEE and 0% Ashkenazi Jewish

    No surprise here, with all known ancestry being overwhelmingly French, with a little English, as well.

    ReplyDelete
  19. I have run two dna files (23andME) through the EuroDNACalc program. I have a new computer with Windows 7 and I can't seem to run a third DNA file on my new computer. Is there a new version for Windows 7 I should use
    ?
    Thank you.

    ReplyDelete
  20. ctuemmler, I've ran EURO-DNA-CALC on both Win7 and XP without a problem. The only thing I can think of is to perhaps install the 64bit version of R if you are running a 64bit operating system.

    ReplyDelete
  21. 23andme for my twin and me

    58% European
    34% African
    8% Asian

    EuroDNACalc for my twin and me

    [1] "NORTHWEST EURO: Maximum Likelihood Estimate=0% Interval=[0, 10]"
    [1] "SOUTHEAST EURO: Maximum Likelihood Estimate=56% Interval=[32, 80]"
    [1] "ASHKENAZI JEWISH: Maximum Likelihood Estimate=44% Interval=[20, 68]"

    23andme for our mother

    57% European
    34% African
    9% Asian

    EuroDNACalc for our mother

    [1] "NORTHWEST EURO: Maximum Likelihood Estimate=0% Interval=[0, 18]"
    [1] "SOUTHEAST EURO: Maximum Likelihood Estimate=84% Interval=[57, 100]"
    [1] "ASHKENAZI JEWISH: Maximum Likelihood Estimate=16% Interval=[0, 41]"

    Considering that we only about 57-58% European, do this mean that my sister are I are about 25% Ashkenazi Jewish (58 x .44), whereas our mother is about 9% Ashkenazi Jewish (57 x .16)?

    ReplyDelete
  22. Considering that we only about 57-58% European, do this mean that my sister are I are about 25% Ashkenazi Jewish (58 x .44), whereas our mother is about 9% Ashkenazi Jewish (57 x .16)?

    This test is meant for Europeans only; non-Europeans will still get 3 numbers that add up to 100%, but they are meaningless.

    ReplyDelete
  23. I tried to use the calculator on my 23andme results and it is giving me an error. Any thoughts as to the cause?
    > EuroDNACalc("23andme")
    Error in if (allele1 == "-") { : missing value where TRUE/FALSE needed

    ReplyDelete
  24. hi,

    my wife is of iraqi Jewish descent on both sides.

    her EuroDNACalc results were as follows: "NORTHWEST EURO: Maximum Likelihood Estimate=0% Interval=[0, 22]"
    [1] "SOUTHEAST EURO: Maximum Likelihood Estimate=42% Interval=[13, 64]"
    [1] "ASHKENAZI JEWISH: Maximum Likelihood Estimate=58% Interval=[36, 80]"

    Does this mean that she is likely at least 36% Ashkenazi, or does it mean that of her European DNA (excluding her Iraqi/Near Eastern DNA), she is at least 36% Ashkenazi?

    Put another way (1) does EURODNACALC only count your European DNA strands; and (2) is Near Eastern/Middle Eastern DNA considered "European" for these purposes?

    ReplyDelete
  25. I'm Ashkenazi. I just ran the test and it came out

    Ashkenazi 92%
    SE Euro 8%

    with intervals

    > EuroDNACalc("23andme")
    [1] "NORTHWEST EURO: Maximum Likelihood Estimate=0% Interval=[0, 8]"
    [1] "SOUTHEAST EURO: Maximum Likelihood Estimate=8% Interval=[0, 30]"
    [1] "ASHKENAZI JEWISH: Maximum Likelihood Estimate=92% Interval=[70, 100]"

    I'd like to submit to DODECAD, but I understand that its closed.

    Thanks, Dienekes

    ReplyDelete
  26. Is the calculator working with 23andme V3 chip results?

    I have both V2 and V3 data and I am getting different results.
    V2: 83% NWE
    0% SEE
    27% ASK
    V3: 85% NWE
    65 SEE
    9% ASK

    Differences are not huge but still I think that they are significant

    ReplyDelete
  27. this euro dna is shit ,sorry for the term but how can my father be 54 % SE euro and me only 5 % ? we are both v3 data

    ReplyDelete
  28. I must be really cumputer stupid, but I cant seem to get this to work for me, I removed all of the comment lines from my 23and me file, but cant get past that. I go to change directory and cant find my data file. I changed the data file to a .R file and still no luck, any help would be appreciated

    ReplyDelete
  29. In change dir you are not going to find your data file, you will go to the directory where you saved the '23andme.txt' file.

    ReplyDelete
  30. "Error in row.names(B) : object 'X' not found"

    What's gone wrong there? Any help is appreciated!

    ReplyDelete
  31. Ok I figured it out - all I needed was to change the 23andme data file to 23andme.txt instead of the default name.

    Thanks for the application Dienekes

    ReplyDelete
  32. My known ancestry on my paternal side is Croat and Italian. Mother was adopted but my guess is NW Europe. My results came back as:

    [1] "NORTHWEST EURO: Maximum Likelihood Estimate=72% Interval=[52, 88]"
    [1] "SOUTHEAST EURO: Maximum Likelihood Estimate=0% Interval=[0, 23]"
    [1] "ASHKENAZI JEWISH: Maximum Likelihood Estimate=28% Interval=[11, 45]"

    Does this mean I have Jewish ancestry or do I take the results with a grain of salt? Thanks in advance.

    ReplyDelete
  33. >> Does this mean I have Jewish ancestry

    It does not.

    ReplyDelete
  34. Using R 2.13.1 on Linux (Fedora) I had to first load the RData file.


    load("/path/to/EuroDNACalc.RData")
    EuroDNACalc("23andme")

    ReplyDelete
  35. [1] "NORTHWEST EURO: Maximum Likelihood Estimate=15% Interval=[0, 42]"
    [1] "SOUTHEAST EURO: Maximum Likelihood Estimate=73% Interval=[38, 100]"
    [1] "ASHKENAZI JEWISH: Maximum Likelihood Estimate=12% Interval=[0, 34]"

    Sounds about right! Though, neither I nor my family knew we had Ashkenazi Jewish in our heritage before 23andMe. Now my dad's mom is less of a mystery :)

    ReplyDelete
  36. Is it still possible to download this?

    I've tried several times but it says forbidden.

    ReplyDelete

Stay on topic. Be polite. Use facts and arguments. Be Brief. Do not post back to back comments in the same thread, unless you absolutely have to. Don't quote excessively. Google before you ask.