{"id":9209,"date":"2022-12-22T18:48:08","date_gmt":"2022-12-22T18:48:08","guid":{"rendered":"https:\/\/www.blopig.com\/blog\/?p=9209"},"modified":"2022-12-22T18:48:10","modified_gmt":"2022-12-22T18:48:10","slug":"quality-stats","status":"publish","type":"post","link":"https:\/\/www.blopig.com\/blog\/2022\/12\/quality-stats\/","title":{"rendered":"Quality Stats"},"content":{"rendered":"\n

Disclaimer – the title is a Quality Street pun only and bears no relation to the quality of the data or analysis presented below. This whole blog post is basically to discredit the personal chocolate preferences of a group member who shall remain nameless. Safe to say though, they Vostly overestimated people’s love for the Toffee Finger. Long live the Orange Creme.<\/em><\/p>\n\n\n\n

\"\"<\/a><\/figure>\n\n\n\n\n\n\n\n

In the run-up to Christmas, I was arguing with another member of the group about which are the best and worst Quality Street chocolates. This is clearly an important topic, with YouGov<\/a> previously dedicating vast resources (I assume) attempting to answer this very question.<\/p>\n\n\n\n

However, as the YouGov poll did not perfectly align with my personal and very accurate preferences, I decided to run another, better experiment. For this experiment, I bought a tub of Quality Street, counted all the chocolates, and then left the tub out in the common area for hungry opiglets to consume. I then recounted the chocolates at various points over the next two days to find out which flavours disappeared first, and perhaps more importantly, which sad chocolates were taken only after all other options were exhausted.<\/p>\n\n\n\n

As expected, crowd favourites The Purple One and The Green Triangle were quick to go, along with the Fudge and Milk Choc Block. The cremes, controversially my personal favourites, sadly performed only averagely. However, to my great delight, the Toffee Finger comprehensively beat all other competition (including the Coconut Eclair!) to take the wooden spoon and provide me with a moderate degree of smugness in the end.<\/p>\n\n\n\n

In an attempt to make this blog post somewhat useful, I’ve included the code I used to make the results plot below. This code should allow you to sort a DataFrame using a custom list, pivot the data when you’re an idiot and type it up the wrong way round, and make a DIY colour palette for your plots.<\/p>\n\n\n\n

And remember, even though you may be appalled by others’ chocolate preferences, this actually makes them the perfect person to sit down and share a box with during these festive times.<\/p>\n\n\n\n

Happy holidays!<\/p>\n\n\n\n

import os\nimport pandas as pd\nimport numpy as np\nimport matplotlib.pyplot as plt<\/pre>\n\n\n\n
Read in data\u00b6<\/a><\/h3>\n\n\n\n
# data I manually recorded in a csv\nquality_df = pd.read_csv(\"Quality_data.csv\", names=[\"flavour\", \"t1\", \"t2\", \"t3\", \"t4\", \"t5\", \"t6\", \"t7\", \"t8\", \"t9\"])\n\n# YouGov ranking - https:\/\/twitter.com\/yougov\/status\/940868550700527616?lang=en-GB\nyougov_ranking = [\"The Purple One\", \"The Green Triangle\", \"Caramel Swirl\", \"Strawberry Delight\", \"Orange Creme\", \"Milk Choc Block\", \"Fudge\", \"Toffee Finger\", \"Orange Chocolate Crunch\", \"Toffee Penny\", \"Coconut Eclair\"]\n\n# sort data (reverse so that plot is in sensible order later)\nquality_df[\"flavour\"] = quality_df[\"flavour\"].astype(\"category\")\nquality_df[\"flavour\"] = quality_df[\"flavour\"].cat.set_categories(yougov_ranking)\nquality_df = quality_df.sort_values([\"flavour\"], ascending=False).reset_index(drop=True)\n\nquality_df<\/pre>\n\n\n\n
<\/th>flavour<\/th>t1<\/th>t2<\/th>t3<\/th>t4<\/th>t5<\/th>t6<\/th>t7<\/th>t8<\/th>t9<\/th><\/tr><\/thead>
0<\/th>Coconut Eclair<\/td>5<\/td>5<\/td>5<\/td>2<\/td>2<\/td>2<\/td>2<\/td>2<\/td>0<\/td><\/tr>
1<\/th>Toffee Penny<\/td>6<\/td>6<\/td>6<\/td>3<\/td>3<\/td>1<\/td>1<\/td>1<\/td>0<\/td><\/tr>
2<\/th>Orange Chocolate Crunch<\/td>5<\/td>5<\/td>5<\/td>5<\/td>4<\/td>2<\/td>1<\/td>0<\/td>0<\/td><\/tr>
3<\/th>Toffee Finger<\/td>7<\/td>6<\/td>6<\/td>5<\/td>5<\/td>5<\/td>5<\/td>4<\/td>2<\/td><\/tr>
4<\/th>Fudge<\/td>8<\/td>6<\/td>4<\/td>0<\/td>0<\/td>0<\/td>0<\/td>0<\/td>0<\/td><\/tr>
5<\/th>Milk Choc Block<\/td>4<\/td>3<\/td>1<\/td>0<\/td>0<\/td>0<\/td>0<\/td>0<\/td>0<\/td><\/tr>
6<\/th>Orange Creme<\/td>6<\/td>6<\/td>6<\/td>3<\/td>3<\/td>2<\/td>2<\/td>0<\/td>0<\/td><\/tr>
7<\/th>Strawberry Delight<\/td>7<\/td>5<\/td>5<\/td>3<\/td>2<\/td>0<\/td>0<\/td>0<\/td>0<\/td><\/tr>
8<\/th>Caramel Swirl<\/td>7<\/td>7<\/td>5<\/td>2<\/td>1<\/td>1<\/td>1<\/td>1<\/td>1<\/td><\/tr>
9<\/th>The Green Triangle<\/td>4<\/td>4<\/td>3<\/td>1<\/td>0<\/td>0<\/td>0<\/td>0<\/td>0<\/td><\/tr>
10<\/th>The Purple One<\/td>5<\/td>4<\/td>4<\/td>0<\/td>0<\/td>0<\/td>0<\/td>0<\/td>0<\/td><\/tr><\/tbody><\/table><\/figure>\n\n\n\n
Pivot data for plotting\u00b6<\/a><\/h3>\n\n\n\n
quality_df = quality_df.T\n\nnew_header = quality_df.iloc[0]   # grab the first row for the header\nquality_df = quality_df[1:]       # take the data less the header row\nquality_df.columns = new_header   # set the header row as the df header\n\nquality_df = quality_df.reset_index(drop=True)\nquality_df.index.names = [\"time\"]\nquality_df<\/pre>\n\n\n\n
flavour<\/th>Coconut Eclair<\/th>Toffee Penny<\/th>Orange Chocolate Crunch<\/th>Toffee Finger<\/th>Fudge<\/th>Milk Choc Block<\/th>Orange Creme<\/th>Strawberry Delight<\/th>Caramel Swirl<\/th>The Green Triangle<\/th>The Purple One<\/th><\/tr>
time<\/th><\/th><\/th><\/th><\/th><\/th><\/th><\/th><\/th><\/th><\/th><\/th><\/tr><\/thead>
0<\/th>5<\/td>6<\/td>5<\/td>7<\/td>8<\/td>4<\/td>6<\/td>7<\/td>7<\/td>4<\/td>5<\/td><\/tr>
1<\/th>5<\/td>6<\/td>5<\/td>6<\/td>6<\/td>3<\/td>6<\/td>5<\/td>7<\/td>4<\/td>4<\/td><\/tr>
2<\/th>5<\/td>6<\/td>5<\/td>6<\/td>4<\/td>1<\/td>6<\/td>5<\/td>5<\/td>3<\/td>4<\/td><\/tr>
3<\/th>2<\/td>3<\/td>5<\/td>5<\/td>0<\/td>0<\/td>3<\/td>3<\/td>2<\/td>1<\/td>0<\/td><\/tr>
4<\/th>2<\/td>3<\/td>4<\/td>5<\/td>0<\/td>0<\/td>3<\/td>2<\/td>1<\/td>0<\/td>0<\/td><\/tr>
5<\/th>2<\/td>1<\/td>2<\/td>5<\/td>0<\/td>0<\/td>2<\/td>0<\/td>1<\/td>0<\/td>0<\/td><\/tr>
6<\/th>2<\/td>1<\/td>1<\/td>5<\/td>0<\/td>0<\/td>2<\/td>0<\/td>1<\/td>0<\/td>0<\/td><\/tr>
7<\/th>2<\/td>1<\/td>0<\/td>4<\/td>0<\/td>0<\/td>0<\/td>0<\/td>1<\/td>0<\/td>0<\/td><\/tr>
8<\/th>0<\/td>0<\/td>0<\/td>2<\/td>0<\/td>0<\/td>0<\/td>0<\/td>1<\/td>0<\/td>0<\/td><\/tr><\/tbody><\/table><\/figure>\n\n\n\n
Plot data\u00b6<\/a><\/h3>\n\n\n\n
flavours = quality_df.columns.tolist()\ntime = quality_df.index.tolist()\n\ndata = [quality_df[flavour].tolist() for flavour in flavours]\nnormalised_data = np.zeros_like(data).astype(float)\nfor i in range(len(time)):\n    normalised_data[:,i] = (np.array(data)[:,i]) \/ np.array(data)[:,i].sum()<\/pre>\n\n\n\n
flavour_to_colour = {\"The Purple One\":          \"purple\",\n                     \"The Green Triangle\":      \"limegreen\",\n                     \"Caramel Swirl\":           \"gold\",\n                     \"Strawberry Delight\":      \"red\",\n                     \"Orange Creme\":            \"darkorange\",\n                     \"Milk Choc Block\":         \"darkgreen\",\n                     \"Fudge\":                   \"fuchsia\",\n                     \"Toffee Finger\":           \"chocolate\",\n                     \"Orange Chocolate Crunch\": \"orangered\",\n                     \"Toffee Penny\":            \"goldenrod\",\n                     \"Coconut Eclair\":          \"mediumblue\"}\n\npalette = [colour for colour in flavour_to_colour.values()]\npalette.reverse()<\/pre>\n\n\n\n
# stacked area plot\nplt.stackplot(time, normalised_data, labels=flavours, colors=palette)\n\nplt.legend(reversed(plt.legend().legendHandles), reversed(flavours), bbox_to_anchor=(1.04, 1), loc=\"upper left\")\n\nplt.xlim([time[0], time[-1]])\nplt.ylim([0,1])\nplt.xlabel(\"Random times I checked the tub\")\nplt.title(\"Quality Street tub composition over time\\n(chocolates ordered according to YouGov ranking)\")\nplt.show()<\/pre>\n\n\n\n
\"\"<\/a><\/figure>\n\n\n\n
<\/p>\n\n\n\n
<\/p>\n\n\n\n
<\/p>\n","protected":false},"excerpt":{"rendered":"
Disclaimer – the title is a Quality Street pun only and bears no relation to the quality of the data or analysis presented below. This whole blog post is basically to discredit the personal chocolate preferences of a group member who shall remain nameless. Safe to say though, they Vostly overestimated people’s love for the […]<\/p>\n","protected":false},"author":91,"featured_media":0,"comment_status":"closed","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"nf_dc_page":"","wikipediapreview_detectlinks":true,"_monsterinsights_skip_tracking":false,"_monsterinsights_sitenote_active":false,"_monsterinsights_sitenote_note":"","_monsterinsights_sitenote_category":0,"ngg_post_thumbnail":0,"_jetpack_memberships_contains_paid_content":false,"footnotes":""},"categories":[621,46],"tags":[],"ppma_author":[558],"class_list":["post-9209","post","type-post","status-publish","format-standard","hentry","category-data-visualization","category-humour"],"jetpack_featured_media_url":"","jetpack_sharing_enabled":true,"authors":[{"term_id":558,"user_id":91,"is_guest":0,"slug":"lewis","display_name":"Lewis Chinery","avatar_url":"https:\/\/secure.gravatar.com\/avatar\/29bddf38b6dd9db3c7161683ddd1a5fc6bad04b6f40e83334e5d43972380ed54?s=96&d=mm&r=g","0":null,"1":"","2":"","3":"","4":"","5":"","6":"","7":"","8":""}],"_links":{"self":[{"href":"https:\/\/www.blopig.com\/blog\/wp-json\/wp\/v2\/posts\/9209","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/www.blopig.com\/blog\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/www.blopig.com\/blog\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/www.blopig.com\/blog\/wp-json\/wp\/v2\/users\/91"}],"replies":[{"embeddable":true,"href":"https:\/\/www.blopig.com\/blog\/wp-json\/wp\/v2\/comments?post=9209"}],"version-history":[{"count":5,"href":"https:\/\/www.blopig.com\/blog\/wp-json\/wp\/v2\/posts\/9209\/revisions"}],"predecessor-version":[{"id":9237,"href":"https:\/\/www.blopig.com\/blog\/wp-json\/wp\/v2\/posts\/9209\/revisions\/9237"}],"wp:attachment":[{"href":"https:\/\/www.blopig.com\/blog\/wp-json\/wp\/v2\/media?parent=9209"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.blopig.com\/blog\/wp-json\/wp\/v2\/categories?post=9209"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.blopig.com\/blog\/wp-json\/wp\/v2\/tags?post=9209"},{"taxonomy":"author","embeddable":true,"href":"https:\/\/www.blopig.com\/blog\/wp-json\/wp\/v2\/ppma_author?post=9209"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}