{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "ANOVA Pizza Case:\n", "1 factor 3 level\n", "In a good pizza restaurant, the three pizzaiolos argue about who bakes the best pizza. They bake 5 pizzas each, and invite 15 random gourmets to sample the pizza. The grade the taste numerically, such that higher numbers mean better taste.\n", "\n", "First, lets input the data. We will use a dataframe." ] }, { "cell_type": "code", "execution_count": 51, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "'data.frame':\t15 obs. of 2 variables:\n", " $ Taste: num 51 45 33 45 67 23 43 23 43 45 ...\n", " $ Names: Factor w/ 3 levels \"Alberto\",\"Francesco\",..: 1 1 1 1 1 2 2 2 2 2 ...\n" ] } ], "source": [ "mydataframe <- data.frame(\"Taste\" = c(51,45,33,45,67,23,43,23,43,45,56,76,74,87,56), \"Names\" = c(rep(\"Alberto\",5),rep(\"Francesco\",5),rep(\"Paolo\",5)))\n", "str(mydataframe)" ] }, { "cell_type": "code", "execution_count": 52, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA0gAAANICAMAAADKOT/pAAADAFBMVEUAAAABAQECAgIDAwME\nBAQFBQUGBgYHBwcICAgJCQkKCgoLCwsMDAwNDQ0ODg4PDw8QEBARERESEhITExMUFBQVFRUW\nFhYXFxcYGBgZGRkaGhobGxscHBwdHR0eHh4fHx8gICAhISEiIiIjIyMkJCQlJSUmJiYnJyco\nKCgpKSkqKiorKyssLCwtLS0uLi4vLy8wMDAxMTEyMjIzMzM0NDQ1NTU2NjY3Nzc4ODg5OTk6\nOjo7Ozs8PDw9PT0+Pj4/Pz9AQEBBQUFCQkJDQ0NERERFRUVGRkZHR0dISEhJSUlKSkpLS0tM\nTExNTU1OTk5PT09QUFBRUVFSUlJTU1NUVFRVVVVWVlZXV1dYWFhZWVlaWlpbW1tcXFxdXV1e\nXl5fX19gYGBhYWFiYmJjY2NkZGRlZWVmZmZnZ2doaGhpaWlqampra2tsbGxtbW1ubm5vb29w\ncHBxcXFycnJzc3N0dHR1dXV2dnZ3d3d4eHh5eXl6enp7e3t8fHx9fX1+fn5/f3+AgICBgYGC\ngoKDg4OEhISFhYWGhoaHh4eIiIiJiYmKioqLi4uMjIyNjY2Ojo6Pj4+QkJCRkZGSkpKTk5OU\nlJSVlZWWlpaXl5eYmJiZmZmampqbm5ucnJydnZ2enp6fn5+goKChoaGioqKjo6OkpKSlpaWm\npqanp6eoqKipqamqqqqrq6usrKytra2urq6vr6+wsLCxsbGysrKzs7O0tLS1tbW2tra3t7e4\nuLi5ubm6urq7u7u8vLy9vb2+vr6/v7/AwMDBwcHCwsLDw8PExMTFxcXGxsbHx8fIyMjJycnK\nysrLy8vMzMzNzc3Ozs7Pz8/Q0NDR0dHS0tLT09PU1NTV1dXW1tbX19fY2NjZ2dna2trb29vc\n3Nzd3d3e3t7f39/g4ODh4eHi4uLj4+Pk5OTl5eXm5ubn5+fo6Ojp6enq6urr6+vs7Ozt7e3u\n7u7v7+/w8PDx8fHy8vLz8/P09PT19fX29vb39/f4+Pj5+fn6+vr7+/v8/Pz9/f3+/v7////i\nsF19AAAACXBIWXMAABJ0AAASdAHeZh94AAAgAElEQVR4nO3dC1xUZd7A8WeGm4oGpoiGBmFZ\nbzdB3pDCLM2yTNM1wdRatvVaumplGWprqWWrvdvuttVWm2K2ZVeK6m27qLWV7mbUFtZqXvJN\nTMtME0VTmfPOjUtJB/D8z5l5ht/382nOQR6eeT4MvxhmzpxRBgDLVKgXAEQCQgIEEBIggJAA\nAYQECCAkQAAhAQIICRBASIAAQgIEEBIggJAAAYQECCAkQAAhAQIICRBASIAAQgIEEBIggJAA\nAYQECCAkQAAhAQIICRBASIAAQgIEEBIggJAAAYQECCAkQAAhAQIICRBASIAAQgIEEBIggJAA\nAYQECCAkQAAhAQIICRBASIAAQgIEEBIggJAAAYQECCAkQAAhAQIICRBASIAAQgIEEBIggJAA\nAYQECCAkQAAhAQIICRBASIAAQgIEEBIggJAAAYQECCAkQAAhAQIICRBASIAAQgIEEBIggJAA\nAYQECCAkQAAhAQIICRBASIAAQgIEEBIggJAAAYQECCAkQAAhAQIICRBASIAAQgIEEBIggJAA\nAYQECCAkQAAhAQIICRBASIAAQgIEEBIggJAAAYQECCAkQAAhAQIICRDgQEj//gDQyr+b/lNu\nf0hrFKCZNU3+Mbc/pPfUD7ZfByDoB/Vek7+GkICfICRAACEBAggJEEBIgABCAgQQEiCAkAAB\nhAQIICRAACEBAggJEEBIgABCAgQQEiCAkAABhITm7Jvf5fe97kWP9YkICc3YynanTPjtlS0G\nVFqeiZDQfH113JQj3s2Gk0ZbnoqQ0HzNOLvKv13u/srqVISE5uv83wa2nrbPWJ2KkNB8df9D\ncCf9UatTERKar4GTAtv9LV63OhUhofl65Piv/dv/aXfA6lSEhObrUPZZawzjwP/ELLY8FSGh\nGds1zHX8f8Ucv8j6TISEZm3zsw+u2C8wDyEBAggJEEBIgABCAgQQEiCAkAABhAQIICRAACEB\nAggJEEBIgABCAgQQEiCAkAABhAQIICRAACEBAggJEEBIgABCAgQQEiCAkAABhAQIICRAACEB\nAggJEEBIgABCAgQQEiCAkAABhAQIICRAACEBAggJEEBIgABCAgQQEiCAkAABhAQIICRAACEB\nAggJEEBIgABCAgQQEiDA6ZA860qWPlayzmM+ipCgGWdDqpybovw6z600G0dI0IyjIe3rqdyZ\n+WPH5We4Vc5+k4GEBM04GtIMNWpbYK98hJppMpCQoBlHQ0rPqqrerepxsslAQoJmHA0pdmrt\n/pQ4k4GEBM04GlLS4Nr9QckmAwkJmnE0pBHuJdW7i10jTQYSEjTjaEgbE1RmYVFxcVFhhkrc\naDKQkKAZZ59HKstWQdllZuMICZpx+siG0oVj8vLGLCw1H0VI0Ez4HGu3Z/K4GgMICXoJn5B2\njsqrkaUqbLkOwCYhCWnLC6/uMR3wF0KCXpwNqahLyyE7jVuilGr1oNk4QoJmHA1plUtFq8sf\nV13yervUCpOBhATNOBrSsKiSqpeiu11aaRjF6gqTgYQEzTgaUtrl3ovL1Ue+/YvMDhEiJGjG\n0ZDifAetTlH+1/T9JtpkICFBM46GdMI13our1Rbf/vDjTAYSEjTjaEh9Ercb2xOPK/Tubo3v\naTKQkKAZR0N6WiUNSFJPuEY+PK+j+qPJQEKCZhwNyTNBqej5xm2+o1b7HTIZSEjQjMNHNmx6\ns9x7+ffJE544YjaMkKCZ8DnWri5CgmYICRBASIAAQgIEEBIggJAAAYQECCAkQAAhAQIICRBA\nSIAAQgIEEBIggJAAAYQECCAkQAAhAQIICRBASAh/Bz7aWBXqNTSAkBDuNg2MUqrNtMpQr8MU\nISHMrW/fb8XurX9L7R3Wbz5HSAhzl/T3n3GqPOkPoV6JGUJCeNvu+ldgZ/Z/h3Yh5ggJ4e0d\nV/BUoi8mhHYh5ggJ4e1fal9g5+n2oV2IOUJCeKto8WJgZ/wloV2IOUJCmBvf7WvfZkXsi6Fe\niRlCQpjbm5M8u3jpuJibQ70QU4SEcPfDPbmJnQe8HOplmCMkQAAhAQIICRBASIAAQgIEEBIg\ngJAAAYSECHNk/SufOP/SJUJCZFmaolqptgudfmk6ISGi3Bc7b5vx7V+Om+rw9RISIsmOVn/1\nb1e4S529YkJCJHkkxRPYOf9WZ6+YkBBJZlS/aOm6fGevmJAQSebkBnd+WeDsFRMSIslrcd/4\ntz90vt/ZKyYkRJIjZw084N1UXd9hj7NXTEiIKOtP7DqzaE5m23ccvl5CQmTZPeeitPNvLnf6\nagkJEEBIgABCAgQQEiCAkAABhAQIICRAACEBAggJEEBIgABCAgQQEiCAkAABhAQIICRAACEB\nAggJEEBIgABCAgQ4HZJnXcnSx0rWecxHERI042xIlXNTlF/nuZVm4wgJmnE0pH09lTszf+y4\n/Ay3ytlvMpCQoBlHQ5qhRm0L7JWPUDNNBhISNONoSOlZNW//VNXjZJOBhATNOBpSbJ13f5oS\nZzKQkKAZR0NKGly7PyjZZCAhQTOOhjTCvaR6d7FrpMlAQoJmHA1pY4LKLCwqLi4qzFCJG00G\nEhI04+zzSGXZKii7zGwcIUEzTh/ZULpwTF7emIUNvFMuIUEz4XOs3aZYVcdeW64DsEn4hOR5\n540aU/iNBL04GtKRxg7krh0042hIqs+yHxo1kJCgGWdDUirplg2NGEhI0IyzIWVf5laufs8c\namggIUEzzoZUYGyZ0VGp5MLN5gMJCZpxOiTDOPzsxS7l6v/8YZOBhATNOB+S16bpHZTqZDKQ\nkKCZkIRkGIee6usyGUhI0EyIQvJabzKQkKCZ0IVkhpCgGUdDWrqqkQMJCZoJn2Pt6iIkaIaQ\nAAGEBAggJEAAIQECCAkQQEiAAEICBBASIICQHFD+8A0znz4Q6lXAToRkvz/GpQ7ul9il6d9o\n6IOQbLc0pshjGPvHJnwR6pXAPoRkN0+XeYFt7rgQrwQ2IiS7faq2BnYeSgvtQmAnQrLbW67g\n2xS+HB/ahcBOhGS3T1V5YIffSJGMkOzmOZG/kZoBQrIdj9o1B4RkP55HagYIyQHbOLIh4hES\nIICQAAGEBAgQCGnLqj1Ci6lBSNCM5ZBWn63UG4bx5Blvia2JkKAdqyF9Ft96sC+kiviJcosi\nJOjGakgjYz/Z6QvJGNhdbE2EBO1YDSl5uBEIaVo7sTURErRjNaTo6cGQpseKrYmQoB2rISX9\nOhjSJalSSzIICdqxGtKQ5IP+kJa7CsTWREjQjtWQ3nVf9rYqef/GmJiP5RZFSNCN5eeRHoxW\nPjFLxJZkEBK0Y/3IhrWTstK6j18rtSA/QoJmONYOEGA1pKWbgztlS0XWE0BI0IzVkFR1P3Ml\nf08REjQjFtLtLpH1BBASNCMWUj6HCKEZsxTS8OHDVc5wn2HZ6grBVRESNGMpJFVHzibBVRES\nNGMppA0bNqh7Nvhs/l50VYQEzVj9G2m+7DOxQYQEzcg8Iftl0dOVEqupRkjQjNWQftftO8P4\nR2ulzpK8c0dI0IzVkHpe6LuILRyj5outiZCgHcsv7JtkGF+piYbRh3M2oBmzGlLMbw3jGfWa\nYdzSVm5RhATdWD75yXjDmOj2/n10Uyu5RRESdGM1pIs6bvu6fa53Z2i62JoICdqxGlKJiopV\nTxqG54RfyC2KkKAby88jLT7vvD97N2+1e1BqSQYhQTu8QhYQQEiAAOsh7XrxgXv9pJZkEBK0\nYzmku1pUv5BCbE2EBO1YDelJdc6d6qZ5fVUeJz9BM2Y1pF7JldvVq4bxeNRyuUUREnRjNaQ2\no40d6hXvzsC+YmsiJGjHakhxhcZ36gnvzqwEsTURErRjNaTUsYandaF352pCQjNmKaRXPzGu\nyPXeqzv+zYrnYs8XXBUhQTPWziJUYDzk2mqs8T0CHrVScFWEBM1YDcnvg1G517wvtSIfQoJm\nREISR0jQjNMhedaVLH2sZJ3HfBQhQTPWQsrfXkcjvrJybkrgcKLOc01P30VI0IzYKYsbc6zd\nvp7KnZk/dlx+hlvl7DcZSEjQjLWQ4lPraPgLZ6hR2wJ75SPUTJOBhATNOPo3UnpWVfVuVY+T\nTQYSEjTjaEixU2v3p8SZDCQkaMbRkJIG1+4PSjYZSEjQjKMhjXAvqd5d7BppMpCQoBlHQ9qY\noDILi4qLiwozVOJGk4GEBM1YCmm32UPY9SnLrn6sPLvMbBwhQTMyZxHyNPrdkUoXjsnLG7Ow\n1HwUIUEz1kL6quyQ9/LgdfEq7TGrKzn46EM1RhES9GItpEtTfU8MTVIJPVuplRZXsjU7q8aJ\nhAS9WAupve/whF2xKTuMjYkXN/rr1zx47wv7TEdw1w6asRTSEZfvbA1Pq7u8lzec0PAXrpi5\nyzB2XOB7sKH9S2YDCQmasRJSbo46Mzc3N02d7r3s6vJeNPCFlydVGZ4clfKrKX1VrNnjDYQE\nzVgJafW7rgWrV6/OiH7XezmrxbvvrWrgCzv1M4w31KW+B81fcA01GUhI0Iy1v5FSrjWMbS3O\n8+3OOaXhL4zJN4x5aq1/f0B7k4GEBM1YC+kmNWx2N/Wobzd38M+Or5F0gWHMDDYyMdZkICFB\nM9ZC2n2eUmqc73XjHwZyMndF3DbjCfW2fz8nzWQgIUEzFo9s8PzrhfX+nQ+f3NPwF76pzt1R\nefJ/rTOMQ7PUZJOBhATNOPtGY9NV/KjfREWf1au9SttpMo6QoBmBkLasasQvo6BFHQPHrLqG\nbjMbRkjQjOWQVp+t1BuG8eQZbzXu+l6dM/H6GUvKzUcREjRjNaTP4lsP9oVUET9RblGEBN1Y\nDWlk7Cc7fSEZA7uLrYmQoB2rISUPNwIhTWsntiZCgnashhQ9PRjSdLMnWJuKkKAZqyEl/ToY\n0iWpUksyCAnasRrSkOSD/pCWuwrE1kRI0I7VkN51X/a2Knn/xpiYj+UWRUjQjeXnkR6M9j/D\nGrPk50YfC0KCZqwf2bB2UlZa9/FrpRbkR0jQjLPH2jUWIUEzhAQIsBzS7ruv6tXTT2xNhATt\nWA1pTWIT3rGv0QgJmrEaUq6atfHAYT+5RRESdGM1pLiBcmupRUjQjOVDhK6XW0stQoJmrIZ0\ndaZHbjE1CAmasRpSecqNjX5Ll8YjJGjG8sPfnye27XW5n9iaCKkZqDoY6hUck0OHfuYTVkNa\n34mHv9FUnkXntIw6ZfreUK+jiQ4tODMm5swF9bZkNaRL1fj3vtjqd6zLqwchRTTPr+JnvPbO\nfV1PNzslW/g50LfD3StWzO/Q90A9n7QaUut+x7wuE4QU0f7Wyv9WJN93vzrUK2mS20/40rf5\nstMd9XzSakhtJU8eVIOQItoFUwPb12J2h3YhTeJJuT+w8+eUeh6pthrSL3KObVXmCCmiHf9c\nYFup/hnahTTJt+qTwM7H6tujP2s1pA3t7zhybOsyQ0gRLfGFwPagq+k/fKHzTfANiYy16puj\nP2s1pII+KnVIgd8xra5+hBTRzi0MbP8RVc//2sNWVYfgG6482sGGu3ZK8fA3mujhhM99m4O5\nQ0K9kia5Od3f/bfpN9fzSashfVTrGJdXH0KKaEcGtf/jh58vy+oi+ZSJ/fb2SF+0du2j6Vn1\nPf/FK2ThvMN3pymVeO3XoV5HE+2blqxU8rR99X2OkBASe74K9QqOyc6few7Zeki7XnzgXr9j\nWNbPISRoxnJId7XgwQbAakhPqnPuVDfN66vylsotipCgG6sh9Uqu3K5eNYzHo5bLLYqQoBur\nIbUZbexQr3h3BvYVWxMhQTuWT35SaHynnvDuzEoQWxMhQTtWQ0oda3ha+474uJqQ0IxZDemK\nXO+9uuPfrHgu9ny5RRESdGM1pIdcW401vkfAo1aKrYmQoB2RIxs+GJV7zfsiywkiJGjGakir\nJY9VrUFI0IzVkFxXyq2lFiFBM1ZDan+N3FpqERI0YzWkYafa8EpzQoJurIb0ebtJ++VWU42Q\noBmBcza07/dLztmAZo5zNgACOGcDIMBKSPeVyq6lFiFBM1ZCUvMNI+Ep2fUEEBI0YyWkuDu8\n+5IvjK1BSNCMlZC65ewkJMDHSkh3KRWl3FHVBFdFSNCMlZCO3Nf/LNX5jGqCqyIkaMby80jc\ntQOsh3STLQ+BExI0wymLAQGcshgQwCmLAQGcshgQwCmLAQGcshgQwCmLAQGanbJ43/vv1/vG\ng0BoaXXK4q+GuV0u9zA93zQREU2nUxbvSMtZXlGxPCdNtzfxReRz+pTFnnUlSx8rWecxH1V/\nSGMzKn2byoxxjbsuwDHOHiJUOTcl8ORt57mVZuPqDelIwrLAzrIEO86lB1hgNaSlm4M7ZY14\nQnZfT+XOzB87Lj/DrXLMTodXb0g71H8CO5+pHY1cKeAQsZdRzG3E76kZatS2wF75CDXTZGC9\nIe1WHwZ2StWeRqwScJBYSLe7Gv7C9Kyq6t2qHiebDKz/b6ST7w5s55t9KRAKYiHlt2v4C2On\n1u5PiTMZWH9If0r0nzvvo4Q/NXxVgKMshTR8+HCVM9xnWLa6ouEvTBpcuz8o2WRg/SFVXdNy\nwqJFE1peU1XPJ4FQshRSnfMVq5xNDX/hCPeS6t3FrpEmA3/uyIanB5500sCnG79SwCGWQtqw\nYYO6Z4PP5u8b84UbE1RmYVFxcVFhhkrcaDKQY+2gGat/I81f25SvLMuu/gWWXWY2jpCgGafP\n2VC6cExe3piFDZwyhZCgmfA5Z8PW7KwaJxIS9BI+52w4+OhDNUYREvTCORsAASE4Z8M/LmvX\nuvs9h82GEBI04+g5G5Iney+eiPLfERxs9koKQoJmHD1ngyowjG9bu2dt/u75Tupxk4GEBM04\nes4GX0h/Uf4D7v6pLjYZSEjQjKPnbPCFdL361L+fkWQykJCgGUfP2eAL6VoVeG3s0BiTgYQE\nzTh6zgZfSHcGX97ap73JQEKCZhw9REi54+Ji1Jv+/bRMk4GEBM04GtKpfnf5dkvVBJOBhATN\nhOqNxt5f+KHJZwkJmrESUuqPCa6KkKAZKyG180lUSsV7/0tsxDkbGo2QoBmrd+0qevV4pcKo\neCWzl+SPPiFBM1ZDmpoeONHj/vSp9Q8+JoQEzVgNKWVacGdaZ5H1BBASNGM1pNibgjs3mZ2n\nrqkICZqxGlK3tMAbf+1LPU1oRT6EBM1YDen3qnvxLmNXcXdl9ZwNdRESNGM1pKqxSqlo73/j\nJE9/SkjQjPUjG1YUnJ16dsFKofUEEBI0E6pDhMwREjRjNaTv5JZSByFBM1ZDalmwSm4xNQgJ\nmrEa0ilKnX1/o86g3xSEBM1YDcmzPD9WxY9u5HuaNxYhQTMCDzZ8fXdXpXo8xEGraMZEHrXz\nvH5ljGozwfSdWpqEkKAZmYe/t8zqoJRyXSX1buOEBM0IhHTkxQFu1WVO+f9eqEYJrYqQoBnL\nIX05O0W5Li4+4t31DGortCpCgmashjQwSh1/44bgB3dKHfNASNCM1ZDUOYsP1HxQep/AinwI\nCZqxGtIHckupg5CgGQ5aBQRYDmn33Vf16ukntiZCgnashrQmUQm9GXNdhATNWA0pV83aeOCw\nn9yiCAm6sfzWlwPl1lJLJqQdT93x4D8F5qlH1Yo/3FksftC71tb+dfbSzaFeROhYDSnperm1\n1BIJ6Xct2vc+Lap3ufWZjlJ2Rmz3cxMS/2bD1Jr6Ps+VfkHnqN9I3i/RitWQrs40e3fyYyUR\n0r0tH68yjM25px9oeGwTfdXhym+837qF0S+KT60pzyWn+t5d5PUOtvx/VQdWQypPubFSbjXV\nBEKqaPOIf7un458sL+enJvcI/I93+iniU2vq1bhN/u1K9/oQryRULD/8/Xli216X+4mtSSSk\nV1odDOxMvdTqVEfpen9gu1k11x+bn5pU/afyqX8M6TpCx2pI6zuF6cPff+0a3PlDd6tTHSX+\n5cC2yvW2+Nx6GjYpuNP/1pCuI3SshnSpGv/eF1v95BYlEdLzCcEzVhb2sTrVUbosCmy3q0/E\n59bT2LzgTuaCkK4jdKyG1Lqf3FpqCYS0M+Yl//ZQtzmWl/NTv7w4sP198hHxufX0ZMK3/u2n\n7jUhXkmoWA2p7US5tdSSeNTuho6+A2r3j+qwy/p6fuKzloW+gl5p9YD41Jo6dFafnd7N5tMH\nhXoloWI1pF/kyK2llkRIhwrc50+4Mim11PpyjvK/bVOvGpvlvs2GqTX1f93bDLz+4th+Uicb\n0I7VkDa0v8OGuzcyRzasmn3VpMX7BSY62q77x4+661NbptbU4Wem5Rf+3Y4nFfVgNaSCPip1\nSIGf2Jo41g7asfwKWRWmD38DTrIa0ke15BZFSNANr5AFBBASIICQAAGEBAggJL0deWnOpPtt\nOgZ93f2T5rzEQVCNQ0haW39mq95DT3NPk3xL+aCqae7ThvZudSavFGkUQtLZ3tQBvoNFX0u4\nXX7u2xNf815+OyB1r/zcEYiQdLYgNfDy5Cdb7JaeeneLJ/3bytTm+sKIpiEknV10S2B7qHWJ\n9NQvtjkU2LnZlhfKRBxC0ln36hd2d31Uemo7X2EciQhJZ/3s/I0UPLMWv5EahZB0xt9IYYOQ\ndMajdmGDkLTG80jhgpD0ZueRDes5sqHxCAkQQEiAAEICBBASIICQAAGEBAggJEAAIQECCAkQ\n4HRInnUlSx8rWdfAOaIJCZpxNqTKuSmB0xt3nmv6zrOEBM04GtK+nsqdmT92XH6GW+WYvU0E\nIUEzjoY0Q43aFtgrH6FmmgwkJGjG0ZDSs2oO96/qcbLJQEKCZhwNKXZq7f6UOJOBhATNOBpS\n0uDa/UHJJgMJCZpxNKQR7iXVu4tdI00GEhI042hIGxNUZmFRcXFRYYZK3GgykJCgGWefRyrL\nrn6fzOwys3GEBM04fWRD6cIxeXljFpaajyIkaCZ8jrXzvPNGjSmEBL2ET0ibYuu8Q7riZGrQ\nSmhCGl1k/nnu2kEzoQlJjTb/PCFBM46GNLOayvBemAwkJGjG0ZDUj5gMJCRoxtmQWt92r5/K\n8V6YDCQkaMbRkEo6dHo5MAN/IyGyOPtgwzdD1LXfG4SEiOP0o3aL2nR5nZAQcRx/+PuL3mpC\nBSEhwjj/PFLVgrh0QkKECcUTsp+cTUiIMCE5ssFzuIG3aiQkaCZ8Dlqti5CgGUICBBASIICQ\nAAGEBAggJEAAIQECCAkQQEiAAEICBBASIICQAAGEBAggJEAAIQECCAkQQEiAAEICBBASIICQ\nAAGEBAggJEAAIQECCAkQQEiAAEICBBASIICQAAGEBAggJEAAIQECCAkQQEiAAEICBBASIICQ\nAAGEBAggJEAAIQECCAkQQEiAAEICBBASIICQAAGEVKsqBNeJCEFIQR/lneju+uvNTl8tIgQh\nBTwbO/ixtx7pddxqh68XEYKQ/La3vsu38YxOO+DsFSNCEJLfglMCfyDtjS929ooRIQjJb9To\n4M4Ftzt7xYgQhOR31YTgzkWznL1iRAhC8pudFdgear/U2StGhCAkv/Uxy/zbecfvcfaKESEI\nKWBBTOEHO94bE/2Mw9eLCEFIQc+coVRUz7ecvlpECEKqsffTSuevFBGCkAABhAQIICRAACEB\nAggJEEBIgABCAgQQEiCAkAABTofkWVey9LGSdR7zUYQEzTgbUuXcFOXXea7p4TiEBM04GtK+\nnsqdmT92XH6GW+XsNxlISNCMoyHNUKO2BfbKR6iZJgMJCZpxNKT0rJpzMFb1ONlkICFBM46G\nFDu1dn9KnMlAQoJmHA0paXDt/qBkk4GEBM04GtII95Lq3cWukSYDCQmacTSkjQkqs7CouLio\nMEMlbjQZSEjQjLPPI5Vlq6DsMrNxhATNOH1kQ+nCMXl5YxaWmo8iJGgmfI612zkqr0YWIUEv\n4RPSnsnjagxQP9hyHYBNQhPSd3vNP/8eIUEvzob0xdgLp+401pylXL3XmY0jJGjG0ZB2dlRK\nZe5IVp2i1Am7TQYSEjTjaEi3qavfmqQuTVtr7B2i5pgMJCRoxtGQunc4bHjS1VPe3W9b9jQZ\nSEjQjKMhJQ7wXuSrHb79XokmAwkJmnE0pBZ53ovxgS8fFm0ykJCgGUdDOrGv92JKvH+/r9nR\n34QEzTgaUv/Otfvp/I2ECOJoSLerL6t3P1K3mAwkJGjG0ZCOHKg5Ddfq+Z+ZDCQkaCZ8jrWr\ni5CgGUICBBASIICQAAGEBAggJEAAIQECCAkQQEiAAEICBBASIICQAAGEBAggJEAAIQECCAkQ\nQEiAAEICBBASIICQAAGEBAggJEAAIQECCAkQQEiAAEICBBASIICQAAGEBAggJEAAIQECCAkQ\nQEiAAEICBBASIICQAAGEBAggJEAAIQECCAkQQEiAAEICBBASIICQAAGEBAggpLCz67X7X/rK\nnqmrPigq+qDKnrmbOUIKM565rVqc3iZ64gEb5i49U6WlqTNLbZi62SOkMDPjuMePGMbfU4bL\nT70uYeR2w9g+MmG9/NzNHiGFly3RJf7tx9Fvic899BKPb1N1yZXiU4OQwsufuwZ3LrlBeupD\nLV4O7LzU4pD03CCk8HJr/+DOxGHSU3+lgnfp1imbHstozggpvNyVHdwZMVp66gq1OrCzyrVP\nem4QUnh5J/oL/3Zv+8Xic/e4ObCd1kN8ahBSePH07vmNd7N/SFf5x7+fiX3et3k+9lnxqUFI\nYWZHj8Srbx/d6aTPbJh7flTvm2/uHTXfhqmbPUIKN4eWjLngmvsrbJn74+kDBkz/2JapmztC\nAgQQEiCAkAABhAQIICRAACEBAggJEEBIgABCAgQQEiCAkAABhAQIICRAACEBAggJEEBIgABC\nAgSEZ0hrFKCZNU3+Mbc/JOPfH4hInLjULnndbJv6YTXXtrl797Zt6rnqYdvm7pZn29QTE2V+\n1v7d9J9yB0ISkrzMtqnn5do29ffKvtPc/+pXtk1dqr63be7cebZNvSzZtqkbQkgGIR2NkJqK\nkAxCOhohNRUhGYR0NEJqKoGfDhwAAAbVSURBVEIyCOlohNRUhGQQ0tEIqakIySCkoxFSUxGS\nQUhHI6SmIiSDkI5GSE1FSAYhHY2QmkqfkLo8b9vUC/rYNvV+d5ltc48bZ9vUZe79ts3dZ4Ft\nUz/fxbapG6JPSFsO2zb1/u22TW1ssm/q776zb24bl73dvkYPb7Ft6oboExIQxggJEEBIgABC\nAgQQEiCAkAABhAQIICRAACEBAggJEEBIgABCAgQQEiCAkAABhAQIICRAgBYhVSy76rSWx+U+\nUmXXFZQoNdOWid8c3CG28xUrbZjZ83zflBYnDVslOedzk86LV8OrP9o4Mjnu5JlSr8KrO7fw\n7fnjZRs23pwmtAjpXhWbk9c7Wl1hU0nfJLe25zt/q4q7IL9POzvmvl4lXD3lMrerSHDOLHVc\nt5qfyLJE16ApPVROpfzcwrfnj5Zt2HhzmtEipGcf2OO9/LSDesKe+Yd0us2W7/widW65d1P1\nrfzUm1T7bd7NC0ryLAUrN3heqvmJzFaLvUsfoebKzy18e/5o2YZ9N6cpLUIKmq/G2zLvIvXy\nvXZ853/oGL9DftaAN9UA36YquqXsvDU/kaUqw7cpd3f2iM8dJHh71p3arpvTnE4hPaAm2zHt\nF22uNWz5zr+qRh1YNuvON8V+Eusoj0rynbHlJTVEdt6an8iFqtC/zVDrxOcOErw960xt281p\nTqOQPDnqDRumrerdZY893/k5avIpvvdRPNeO30vzVOI1Uy+Pvnyn7LQ1P5FjVOCvr3xVIj53\ngOTtWTu1fTenOY1Cmq2G2jHtAvW6Yc93fpKKOnVlxScXqwvl5zaMJ47zNnqq9F+NNT+RearY\nvx2nHhOfO0Dy9qyd2r6b05w+Id2nethx9s9P4iYYNn3nr1PR//Fu9p1wDO/t26DbXbd8sb/0\nkuAdMDFHhTRWLRWf20/09qyZ2sab05w2Id2jsuw4H6Kn+0kVhk3f+RnqTP+2QD0oPvdraoRv\nU9klSvaciI7dtZO9PauntvPmNKdLSLPVuXvsmPdw7VvCj5aee4nq5d9OUfdKT21MVg/7t3nq\nBdF5j3qwIdOeBxuEb8/qqe28Oc1pEtIN6sIKWyauGu2XozJGSz616Vfuan/It+0r/MPuM0EF\nTkXfW70qOm+dh78zfZtt7hQ7Hv6Wvj2rp7bz5jSnRUhVY1V/oSfYf4Y99wWGqtmG70Zuv098\n6r+pjlu9mxJXK9lf1HWfkF3i/c6PknpCtu7c8rfnTx7H4K5d/RYo94gCn3vsugZ7vvPb0tS5\nEwe6Y+R/IRlH+qj44ZMvVqJ/fj1XUHCRSisouMn3QVmCe/DULNVT6Ce+7tzCt+ePlu1HSPWb\nXn2/t79d12DTd37nb1Jj2v3ChsfsDOOH32e3jkoatFxyzpnBb3Oq/6ONI5Ji02dI/TKtO7fw\n7fnjZfsQEqAnQgIEEBIggJAAAYQECCAkQAAhAQIICRBASIAAQgIEEBIggJAAAYQECCAkQAAh\nAQIICRBASIAAQgIEEBIggJAAAYQECCAkQAAhAQIICRBASIAAQgIEEBIggJAAAYQECCAkQAAh\nAQIICRBASIAAQtJYu9RQrwDVCCncHVAJP/cpQgofhBTuCEkLhBTuCEkLhBTuAiF9pAq+HNGu\nxX+/4v+3qntPi+s8tSIY0qqhyTGdRv3HuzdY/cn3D7PU6BAttvkipHBXHVLf5B7XXRnl/ofv\n38ap1JumpfdKTPV98LA76drp+bHx/zSMXSfGfWgYb7pP3x/KFTdLhBTuqkNSszyGsVQN8n6w\nUnXfZxj7M1Wq94PPYvpXejcftz7be/le9CkVX3dsuTaUC26eCCncVYd04mHvxpOQ7L0sUMW+\nz7ziD2mSenunz2C1xfvRfDXyYvVI6FbbbBFSuKsOabD/ozNivRdnq12+/Qp/SFmq2mrvR57+\nSo0I1VKbM0IKdzUPNvg/6h7lvUiNDnwqPtV7kaZK3gjY4/u3vyj1r1Ass7kjpHBXT0g/+o3U\nXb1fZ/Tnrdu6zzrg8BJBSOGvnpB+9DfSeHVj7eCDma7XZqrxTq8RhBT26glpReBRux7+kMqi\nY5b7PlOxzPA98jDdOJKrngrRWpsxQgp39YRkjFVptc8jPRrt6n/rzYPizzCMYtXzsGF8efxx\nm0K12maLkMJdfSFV/b5bbErNkQ0fXdMltu0ZE1Ya/9c2YbPvH15Q5/wQkrU2Y4QECCAkQAAh\nAQIICRBASIAAQgIEEBIggJAAAYQECCAkQAAhAQIICRBASIAAQgIEEBIggJAAAYQECCAkQAAh\nAQIICRBASIAAQgIEEBIggJAAAYQECCAkQAAhAQIICRBASIAAQgIEEBIg4P8BkS4QHQeRyMkA\nAAAASUVORK5CYII=", "text/plain": [ "plot without title" ] }, "metadata": { "image/png": { "height": 420, "width": 420 }, "text/plain": { "height": 420, "width": 420 } }, "output_type": "display_data" } ], "source": [ "plot(mydataframe$Taste)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We have 3 groups with n_i=5 points each, so the degrees of freedom for SST:3-1=2 and for SSE:15-3=12. Indeed, we find 12+2=14=n-1. \n", "Now we need to compute the statistical parameters of the groups." ] }, { "cell_type": "code", "execution_count": 57, "metadata": {}, "outputs": [], "source": [ "# Compute all variances = standard deviation ^2\n", "s1sq=sd(mydataframe$Taste[1:5])^2\n", "s2sq=sd(mydataframe$Taste[6:10])^2\n", "s3sq=sd(mydataframe$Taste[11:15])^2\n", "\n", "# Compute all group averages\n", "a1=mean(mydataframe$Taste[1:5])\n", "a2=mean(mydataframe$Taste[6:10])\n", "a3=mean(mydataframe$Taste[11:15])\n", "\n", "# Compute global average\n", "a=mean(mydataframe$Taste)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Next we can compute SS, SST, SSE:" ] }, { "cell_type": "code", "execution_count": 60, "metadata": {}, "outputs": [ { "data": { "text/html": [ "4883.73333333333" ], "text/latex": [ "4883.73333333333" ], "text/markdown": [ "4883.73333333333" ], "text/plain": [ "[1] 4883.733" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/html": [ "3022.93333333333" ], "text/latex": [ "3022.93333333333" ], "text/markdown": [ "3022.93333333333" ], "text/plain": [ "[1] 3022.933" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/html": [ "1860.8" ], "text/latex": [ "1860.8" ], "text/markdown": [ "1860.8" ], "text/plain": [ "[1] 1860.8" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "SS=sum((data-a)^2)\n", "SST= 5*(a1-a)^2 + 5*(a2-a)^2 + 5*(a3-a)^2\n", "SSE= (5-1)*s1sq + (5-1)*s2sq + (5-1)*s3sq\n", "SS\n", "SST\n", "SSE" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Always test that SST+SSE=SS!!" ] }, { "cell_type": "code", "execution_count": 61, "metadata": {}, "outputs": [ { "data": { "text/html": [ "4883.73333333333" ], "text/latex": [ "4883.73333333333" ], "text/markdown": [ "4883.73333333333" ], "text/plain": [ "[1] 4883.733" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "SST+SSE" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Now we compute the mean squares by normalizing the sum of squares by their respective degrees of freedom" ] }, { "cell_type": "code", "execution_count": 62, "metadata": {}, "outputs": [ { "data": { "text/html": [ "1511.46666666667" ], "text/latex": [ "1511.46666666667" ], "text/markdown": [ "1511.46666666667" ], "text/plain": [ "[1] 1511.467" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/html": [ "155.066666666667" ], "text/latex": [ "155.066666666667" ], "text/markdown": [ "155.066666666667" ], "text/plain": [ "[1] 155.0667" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "MST=SST/2\n", "MSE=SSE/12\n", "MST\n", "MSE" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Now we compute the F-statistic:" ] }, { "cell_type": "code", "execution_count": 63, "metadata": {}, "outputs": [ { "data": { "text/html": [ "9.74720550300946" ], "text/latex": [ "9.74720550300946" ], "text/markdown": [ "9.74720550300946" ], "text/plain": [ "[1] 9.747206" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "F=MST/MSE\n", "F" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "This is a large F, but to make a statistical statement we need to compare it to the quantile of the F-distribution. Large values of F disprove H0, so we find the maximal F at which we still accept H0 at a level of significance 0.05 as:" ] }, { "cell_type": "code", "execution_count": 64, "metadata": {}, "outputs": [ { "data": { "text/html": [ "3.88529383465239" ], "text/latex": [ "3.88529383465239" ], "text/markdown": [ "3.88529383465239" ], "text/plain": [ "[1] 3.885294" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "F_crit=qf(0.95,2,12) #(q)uantile of the (f)-distribution at alpha, with df1 and df2 degrees of freedom.\n", "F_crit" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Indeed, our F is significantly above the critical value, so we can reject the null hypothesis that all pizza are equivalent. Clearly there are better and worse bakers in the restaurant." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Let us redo all with internal R commands:" ] }, { "cell_type": "code", "execution_count": 65, "metadata": {}, "outputs": [], "source": [ "aov.out=aov(Taste~Names, data=mydataframe)" ] }, { "cell_type": "code", "execution_count": 66, "metadata": {}, "outputs": [ { "data": { "text/plain": [ " Df Sum Sq Mean Sq F value Pr(>F) \n", "Names 2 3023 1511.5 9.747 0.00306 **\n", "Residuals 12 1861 155.1 \n", "---\n", "Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "summary(aov.out)" ] } ], "metadata": { "kernelspec": { "display_name": "R", "language": "R", "name": "ir" }, "language_info": { "codemirror_mode": "r", "file_extension": ".r", "mimetype": "text/x-r-source", "name": "R", "pygments_lexer": "r", "version": "3.6.3" } }, "nbformat": 4, "nbformat_minor": 4 }