ax.text(5.5e-8, 1.2, '1st') ax.text(3.4e-8, 1.2, '15th') ax.text(1.1e-8, 1.2, '11th') # Mean over 2000-2008 #ax = plot_me(good_XXth_medians_smooth, years=(2000,2008)) #ax.text(2004, 6.2e-8, '1st') #ax.text(2005, 1.3e-8, '11th') #ax.text(2001, 3.3e-8, '15th') # Plot raw data # This is a total mess #plot_me(good_XXth_data) # Plot good medians ax = plot_me(good_XXth_medians_smooth) ax.plot(good_XXth_medians_smooth['year'], good_XXth_medians_smooth['11th'], 'b-', linewidth=4) ax.text(1935, 4e-7, '1st') ax.text(1925, 4.5e-8, '11th') ax.text(1920, 2.58e-7, '31st') # Plot good and xxth together ax = plot_me(good_XXth_medians_smooth, sum_bad_XXth_and_11th_medians_smooth) ax.plot(good_XXth_medians_smooth['year'], good_XXth_medians_smooth['11th'], 'b-', linewidth=4) ax.plot(good_XXth_medians_smooth['year'], sum_bad_XXth_and_11th_medians_smooth['allth'], 'b-', linewidth=4) ax.plot(good_XXth_medians_smooth['year'], sum_XXth_minus_bad_medians_smooth['diffth'], 'b-', linewidth=4) ax.text(1925, 4.5e-8, '11th') ax.text(1927, 10e-8, 'sum') ax.text(1895, 2.5e-8, 'xxth') # Plot good and nth together
ax.text(5.5e-8, 1.2, '1st') ax.text(3.4e-8, 1.2, '15th') ax.text(1.1e-8, 1.2, '11th') # Mean over 2000-2008 #ax = plot_me(good_XXth_medians_smooth, years=(2000,2008)) #ax.text(2004, 6.2e-8, '1st') #ax.text(2005, 1.3e-8, '11th') #ax.text(2001, 3.3e-8, '15th') # Plot raw data # This is a total mess #plot_me(good_XXth_data) # Plot good medians ax = plot_me(good_XXth_medians_smooth) ax.plot(good_XXth_medians_smooth['year'], good_XXth_medians_smooth['11th'], 'b-', linewidth=4) ax.text(1935, 4e-7, '1st') ax.text(1925, 4.5e-8, '11th') ax.text(1920, 2.58e-7, '31st') # Plot good and xxth together ax = plot_me(good_XXth_medians_smooth, sum_bad_XXth_and_11th_medians_smooth) ax.plot(good_XXth_medians_smooth['year'], good_XXth_medians_smooth['11th'], 'b-', linewidth=4) ax.plot(good_XXth_medians_smooth['year'],
only_22d = filter_for_specific_day(old_style_data, '22d') only_23d = filter_for_specific_day(old_style_data, '23d') # Sum old and new style sum_2nd_2d = sum_over_patterns(only_2nd.merge(only_2d, on='year'), new_name='2xx') sum_3rd_3d = sum_over_patterns(only_3rd.merge(only_3d, on='year'), new_name='3xx') sum_22nd_22d = sum_over_patterns(only_22nd.merge(only_22d, on='year'), new_name='22xx') sum_23rd_23d = sum_over_patterns(only_23rd.merge(only_23d, on='year'), new_name='23xx') # Merge data sum_new_old = merge_data(sum_2nd_2d, sum_3rd_3d, sum_22nd_22d, sum_23rd_23d) sum_new_old_medians = median_over_months(sum_new_old) sum_new_old_medians_smooth = smooth_data(sum_new_old_medians) # Plot new-style ax = plot_me(good_XXth_medians_smooth) ax.plot(good_XXth_medians_smooth['year'], good_XXth_medians_smooth['2nd'], 'g-', linewidth=4) ax.plot(good_XXth_medians_smooth['year'], good_XXth_medians_smooth['3rd'], 'r-', linewidth=4) ax.plot(good_XXth_medians_smooth['year'], good_XXth_medians_smooth['22nd'], 'c-', linewidth=4) ax.plot(good_XXth_medians_smooth['year'], good_XXth_medians_smooth['23rd'], 'm-', linewidth=4) ax.text(1878, 13.2e-8, '2nd', color='g') ax.text(1889, 11.6e-8, '3rd', color='r') ax.text(1902, 10.7e-8, '22nd', color='c') ax.text(1916, 9.4e-8, '23rd', color='m') ax.text(1935, 4e-7, '1st') ax.text(1890, 2.58e-7, '4th') ax.text(1933, 3.3e-8, '11th') ax.text(1920, 2.58e-7, '31st') # Plot old-style ax = plot_me(good_XXth_medians_smooth)
sum_2nd_2d = sum_over_patterns(only_2nd.merge(only_2d, on='year'), new_name='2xx') sum_3rd_3d = sum_over_patterns(only_3rd.merge(only_3d, on='year'), new_name='3xx') sum_22nd_22d = sum_over_patterns(only_22nd.merge(only_22d, on='year'), new_name='22xx') sum_23rd_23d = sum_over_patterns(only_23rd.merge(only_23d, on='year'), new_name='23xx') # Merge data sum_new_old = merge_data(sum_2nd_2d, sum_3rd_3d, sum_22nd_22d, sum_23rd_23d) sum_new_old_medians = median_over_months(sum_new_old) sum_new_old_medians_smooth = smooth_data(sum_new_old_medians) # Plot new-style ax = plot_me(good_XXth_medians_smooth) ax.plot(good_XXth_medians_smooth['year'], good_XXth_medians_smooth['2nd'], 'g-', linewidth=4) ax.plot(good_XXth_medians_smooth['year'], good_XXth_medians_smooth['3rd'], 'r-', linewidth=4) ax.plot(good_XXth_medians_smooth['year'], good_XXth_medians_smooth['22nd'], 'c-', linewidth=4) ax.plot(good_XXth_medians_smooth['year'], good_XXth_medians_smooth['23rd'], 'm-',