alberthrocks

Gene Probability Python Code, "direct" conversion

Mar 21st, 2011
122
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
Python 0.81 KB | None | 0 0
  1. from random import *
  2. from sys import stdout as o
  3. genotype = [0,0]
  4. homoDom,heteroDom,homoRec,dom,rec=0,0,0,0,0
  5. for x in range(0,100):
  6.     genotype[0]=int(random()*2+1)
  7.     genotype[1]=int(random()*2+1)
  8.     if genotype[0] == 1:
  9.         o.write("R")
  10.         if genotype[1] == 1:
  11.             o.write("R\t|\thomozygous dominant\t|\tnormal skin\n")
  12.             homoDom+=1
  13.         else:
  14.             o.write("r\t|\theterozygous dominant\t|\tnormal skin\n")
  15.             heteroDom+=1
  16.     else:
  17.         o.write("r")
  18.         if genotype[1] == 1:
  19.             o.write("R\t|\theterozygous dominant\t|\tnormal skin\n")
  20.             heteroDom+=1
  21.         else:
  22.             o.write("r\t|\thomozygous recessive\t|\talbino\n")
  23.             homoRec+=1
  24. print "RR: "+str(homoDom)+"/100 Rr or rR: "+str(heteroDom)+"/100 rr: "+str(homoRec)+"/100."
  25. rec=homoRec
  26. dom=homoDom+heteroDom
  27. print str(dom)+"/100 were dominant. "+str(rec)+"/100 were recessive."
Advertisement
Add Comment
Please, Sign In to add comment