1531-NFL Quarterback Passer Rating.sagews

NFL Quarterback Passer Rating

Formula from http://www.nfl.com/help/quarterbackratingformula

Idea from http://demonstrations.wolfram.com/NFLQuarterbackPasserRating/

%cython
cdef int passpercent, comppass, passyrds, tdpass, intercept

cdef float bound(float x):
    return max(0, min(2.375, x))

cpdef float qb(float passatt, float comppass, float passyrds, float tdpass, float intercept):
    """
    quarterback rating formula
    """
    ret = bound((comppass/passatt * 100 - 30) * 0.05)
    ret += bound((passyrds/passatt - 3) * 0.25)
    ret += bound((tdpass/passatt) * 100 * 0.2)
    ret += bound((2.375 - (intercept/passatt) * 100 * 0.25))
    return ret * 100 / 6
    
cpdef float qb_avg(float passatt, float avgpassatt, float comppass, float avgpassyards, float tdpass, float intercept):
    """
    quarterback rating formula, average value arguments
    """
    ret = bound((avgpassatt - 30) * 0.05)
    ret += bound((avgpassyards - 3) * 0.25)
    ret += bound((tdpass/passatt) * 100 * 0.2)
    ret += bound((2.375 - (intercept/passatt) * 100 * 0.25))
    return ret * 100 / 6

qb(passatt=44, comppass=33, passyrds=11, tdpass=2, intercept=2)

73.295455932617188

var('y p')
p1 = plot3d(lambda y,p: qb_avg(44,p,5, y, 2, 2),(y, 1,50), (p, 1,100))
p2 = point3d((11, 50, qb(44, 33, 11, 2, 2)), size=30, rgbcolor='orange')
show(p1+p2)

@interact
def calc_quarterback_rating(
    passatt = (1,100,1),\
    comppass = (0,20,1),\
    passyrds = (1,2000,1),\
    tdpass = (0,10,1),\
    intercept = (0,20,1)):
        
    passydadj = n(max(20, min(100, passyrds/passatt))) 
    passpercent = 100*comppass/passatt
    avgpassyards = passyrds/passatt
    qb_rating = qb(passatt, comppass, passyrds, tdpass, intercept)
    
    print 'Pass Percentage:    %3.2f%%' % (passpercent)
    print 'Avg. Yards/Pass:    %3.2f' %(avgpassyards)
    print 'Quarterback-Rating: %3.2f' % (qb_rating)
    
    p1 = plot3d(lambda p,y: qb_avg(passatt, p, comppass, y, tdpass, intercept), (p, 1, 100),(y, 1,passydadj))
    p2 = point3d((passpercent, avgpassyards, qb_rating), size=30, rgbcolor='orange')
    show(p1+p2)

passatt  comppass  passyrds  tdpass  intercept  [removed] [removed] [removed] [removed]