import numpy as np
raw_data_points = np.array([ 53, 48, 40, 30, 30 ]) #spitballed based on the chart in Fig 19
sample_mean = np.mean( raw_data_points )
sample_standard_deviation = np.std( raw_data_points )
sample_std_err_of_mean = sample_standard_deviation / np.sqrt( len( raw_data_points ) )
sample_95pcnt_interval = 1.96 * sample_std_err_of_mean
print "Confidence interval: ", sample_mean, "+-", sample_95pcnt_interval
aW1wb3J0IG51bXB5IGFzIG5wCgpyYXdfZGF0YV9wb2ludHMgPSBucC5hcnJheShbIDUzLCA0OCwgNDAsIDMwLCAzMCBdKSAjc3BpdGJhbGxlZCBiYXNlZCBvbiB0aGUgY2hhcnQgaW4gRmlnIDE5CgpzYW1wbGVfbWVhbiA9IG5wLm1lYW4oIHJhd19kYXRhX3BvaW50cyApCnNhbXBsZV9zdGFuZGFyZF9kZXZpYXRpb24gPSBucC5zdGQoIHJhd19kYXRhX3BvaW50cyApCnNhbXBsZV9zdGRfZXJyX29mX21lYW4gPSBzYW1wbGVfc3RhbmRhcmRfZGV2aWF0aW9uIC8gbnAuc3FydCggbGVuKCByYXdfZGF0YV9wb2ludHMgKSApCnNhbXBsZV85NXBjbnRfaW50ZXJ2YWwgPSAxLjk2ICogc2FtcGxlX3N0ZF9lcnJfb2ZfbWVhbgoKcHJpbnQgIkNvbmZpZGVuY2UgaW50ZXJ2YWw6ICIsIHNhbXBsZV9tZWFuLCAiKy0iLCBzYW1wbGVfOTVwY250X2ludGVydmFs