Exemplo n.º 1
0
def getLL(ys, d, dof=4, fn="tdf/Venture/posteriorDict"):
  global pd
  global fac
  global maxPd
  if dof == 4:
    if pd == None:
      #print "loading Dict"
      with open(fn,'r') as f:
        fac, pd = cPickle.load(f)
      maxPd = max(pd.values())
    try:
      return pd[round(d*fac)]
    except:
      return float("-inf")
  elif dof == 21:
    if pd == None:
      #print "loading Dict"
      with open(fn + "21",'r') as f:
        fac, pd = cPickle.load(f)
      maxPd = max(pd.values())
    try:
      return pd[round(d*fac)]
    except:
      return float("-inf")
  else:
    return cp.logLiks(ys, d, base=2)
Exemplo n.º 2
0
def getLL(ys,d,dof=4):
  global pd
  global fac
  if dof == 4:
    if pd == None:
      print "loading Dict"
      with open("posteriorDict",'r') as f:
        fac, pd = cPickle.load(f)
    return pd[round(d*fac)]
  elif dof == 21:
    if pd == None:
      print "loading Dict"
      with open("posteriorDict21",'r') as f:
        fac, pd = cPickle.load(f)
    return pd[round(d*fac)]
  else:
    return cp.logLiks(ys, d, base=2)