Skip to content

Instantly share code, notes, and snippets.

@aakhmetz
Created October 19, 2020 05:09
Show Gist options
  • Save aakhmetz/a0180cb3d4684ae9d1cdc84108de56ce to your computer and use it in GitHub Desktop.
Save aakhmetz/a0180cb3d4684ae9d1cdc84108de56ce to your computer and use it in GitHub Desktop.
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [],
"source": [
"libraries = c(\"dplyr\", \"magrittr\")\n",
"for(x in libraries) { library(x, character.only=TRUE, warn.conflicts=FALSE) }"
]
},
{
"cell_type": "code",
"execution_count": 11,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"<table>\n",
"<caption>A data.frame: 23 × 3</caption>\n",
"<thead>\n",
"\t<tr><th scope=col>ID</th><th scope=col>onset</th><th scope=col>death</th></tr>\n",
"\t<tr><th scope=col>&lt;int&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;dbl&gt;</th></tr>\n",
"</thead>\n",
"<tbody>\n",
"\t<tr><td> 2</td><td>17</td><td>33</td></tr>\n",
"\t<tr><td> 3</td><td>19</td><td>48</td></tr>\n",
"\t<tr><td> 4</td><td>21</td><td>31</td></tr>\n",
"\t<tr><td> 5</td><td>21</td><td>47</td></tr>\n",
"\t<tr><td> 6</td><td>21</td><td>30</td></tr>\n",
"\t<tr><td> 7</td><td>22</td><td>41</td></tr>\n",
"\t<tr><td> 8</td><td>22</td><td>37</td></tr>\n",
"\t<tr><td> 9</td><td>23</td><td>32</td></tr>\n",
"\t<tr><td>10</td><td>23</td><td>31</td></tr>\n",
"\t<tr><td>11</td><td>23</td><td>35</td></tr>\n",
"\t<tr><td>12</td><td>23</td><td>54</td></tr>\n",
"\t<tr><td>13</td><td>23</td><td>30</td></tr>\n",
"\t<tr><td>14</td><td>25</td><td>31</td></tr>\n",
"\t<tr><td>15</td><td>26</td><td>33</td></tr>\n",
"\t<tr><td>16</td><td>26</td><td>30</td></tr>\n",
"\t<tr><td>17</td><td>26</td><td>33</td></tr>\n",
"\t<tr><td>18</td><td>27</td><td>48</td></tr>\n",
"\t<tr><td>19</td><td>31</td><td>39</td></tr>\n",
"\t<tr><td>20</td><td>39</td><td>51</td></tr>\n",
"\t<tr><td>21</td><td>41</td><td>48</td></tr>\n",
"\t<tr><td>22</td><td>41</td><td>55</td></tr>\n",
"\t<tr><td>23</td><td>42</td><td>60</td></tr>\n",
"\t<tr><td>24</td><td>44</td><td>62</td></tr>\n",
"</tbody>\n",
"</table>\n"
],
"text/latex": [
"A data.frame: 23 × 3\n",
"\\begin{tabular}{lll}\n",
" ID & onset & death\\\\\n",
" <int> & <dbl> & <dbl>\\\\\n",
"\\hline\n",
"\t 2 & 17 & 33\\\\\n",
"\t 3 & 19 & 48\\\\\n",
"\t 4 & 21 & 31\\\\\n",
"\t 5 & 21 & 47\\\\\n",
"\t 6 & 21 & 30\\\\\n",
"\t 7 & 22 & 41\\\\\n",
"\t 8 & 22 & 37\\\\\n",
"\t 9 & 23 & 32\\\\\n",
"\t 10 & 23 & 31\\\\\n",
"\t 11 & 23 & 35\\\\\n",
"\t 12 & 23 & 54\\\\\n",
"\t 13 & 23 & 30\\\\\n",
"\t 14 & 25 & 31\\\\\n",
"\t 15 & 26 & 33\\\\\n",
"\t 16 & 26 & 30\\\\\n",
"\t 17 & 26 & 33\\\\\n",
"\t 18 & 27 & 48\\\\\n",
"\t 19 & 31 & 39\\\\\n",
"\t 20 & 39 & 51\\\\\n",
"\t 21 & 41 & 48\\\\\n",
"\t 22 & 41 & 55\\\\\n",
"\t 23 & 42 & 60\\\\\n",
"\t 24 & 44 & 62\\\\\n",
"\\end{tabular}\n"
],
"text/markdown": [
"\n",
"A data.frame: 23 × 3\n",
"\n",
"| ID &lt;int&gt; | onset &lt;dbl&gt; | death &lt;dbl&gt; |\n",
"|---|---|---|\n",
"| 2 | 17 | 33 |\n",
"| 3 | 19 | 48 |\n",
"| 4 | 21 | 31 |\n",
"| 5 | 21 | 47 |\n",
"| 6 | 21 | 30 |\n",
"| 7 | 22 | 41 |\n",
"| 8 | 22 | 37 |\n",
"| 9 | 23 | 32 |\n",
"| 10 | 23 | 31 |\n",
"| 11 | 23 | 35 |\n",
"| 12 | 23 | 54 |\n",
"| 13 | 23 | 30 |\n",
"| 14 | 25 | 31 |\n",
"| 15 | 26 | 33 |\n",
"| 16 | 26 | 30 |\n",
"| 17 | 26 | 33 |\n",
"| 18 | 27 | 48 |\n",
"| 19 | 31 | 39 |\n",
"| 20 | 39 | 51 |\n",
"| 21 | 41 | 48 |\n",
"| 22 | 41 | 55 |\n",
"| 23 | 42 | 60 |\n",
"| 24 | 44 | 62 |\n",
"\n"
],
"text/plain": [
" ID onset death\n",
"1 2 17 33 \n",
"2 3 19 48 \n",
"3 4 21 31 \n",
"4 5 21 47 \n",
"5 6 21 30 \n",
"6 7 22 41 \n",
"7 8 22 37 \n",
"8 9 23 32 \n",
"9 10 23 31 \n",
"10 11 23 35 \n",
"11 12 23 54 \n",
"12 13 23 30 \n",
"13 14 25 31 \n",
"14 15 26 33 \n",
"15 16 26 30 \n",
"16 17 26 33 \n",
"17 18 27 48 \n",
"18 19 31 39 \n",
"19 20 39 51 \n",
"20 21 41 48 \n",
"21 22 41 55 \n",
"22 23 42 60 \n",
"23 24 44 62 "
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"# data\n",
"df = data.frame(ID = 2:24, \n",
" onset = c(17, 19, 21, 21, 21, 22, 22, 23, 23, 23, 23, 23, 25, 26, 26, 26, 27, 31, 39, 41, 41, 42, 44),\n",
" death = c(33, 48, 31, 47, 30, 41, 37, 32, 31, 35, 54, 30, 31, 33, 30, 33, 48, 39, 51, 48, 55, 60, 62))\n",
"df"
]
},
{
"cell_type": "code",
"execution_count": 15,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"<table>\n",
"<caption>A data.frame: 23 × 4</caption>\n",
"<thead>\n",
"\t<tr><th scope=col>ID</th><th scope=col>onset</th><th scope=col>death</th><th scope=col>time</th></tr>\n",
"\t<tr><th scope=col>&lt;int&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;dbl&gt;</th></tr>\n",
"</thead>\n",
"<tbody>\n",
"\t<tr><td> 2</td><td>17</td><td>33</td><td>16</td></tr>\n",
"\t<tr><td> 3</td><td>19</td><td>48</td><td>29</td></tr>\n",
"\t<tr><td> 4</td><td>21</td><td>31</td><td>10</td></tr>\n",
"\t<tr><td> 5</td><td>21</td><td>47</td><td>26</td></tr>\n",
"\t<tr><td> 6</td><td>21</td><td>30</td><td> 9</td></tr>\n",
"\t<tr><td> 7</td><td>22</td><td>41</td><td>19</td></tr>\n",
"\t<tr><td> 8</td><td>22</td><td>37</td><td>15</td></tr>\n",
"\t<tr><td> 9</td><td>23</td><td>32</td><td> 9</td></tr>\n",
"\t<tr><td>10</td><td>23</td><td>31</td><td> 8</td></tr>\n",
"\t<tr><td>11</td><td>23</td><td>35</td><td>12</td></tr>\n",
"\t<tr><td>12</td><td>23</td><td>54</td><td>31</td></tr>\n",
"\t<tr><td>13</td><td>23</td><td>30</td><td> 7</td></tr>\n",
"\t<tr><td>14</td><td>25</td><td>31</td><td> 6</td></tr>\n",
"\t<tr><td>15</td><td>26</td><td>33</td><td> 7</td></tr>\n",
"\t<tr><td>16</td><td>26</td><td>30</td><td> 4</td></tr>\n",
"\t<tr><td>17</td><td>26</td><td>33</td><td> 7</td></tr>\n",
"\t<tr><td>18</td><td>27</td><td>48</td><td>21</td></tr>\n",
"\t<tr><td>19</td><td>31</td><td>39</td><td> 8</td></tr>\n",
"\t<tr><td>20</td><td>39</td><td>51</td><td>12</td></tr>\n",
"\t<tr><td>21</td><td>41</td><td>48</td><td> 7</td></tr>\n",
"\t<tr><td>22</td><td>41</td><td>55</td><td>14</td></tr>\n",
"\t<tr><td>23</td><td>42</td><td>60</td><td>18</td></tr>\n",
"\t<tr><td>24</td><td>44</td><td>62</td><td>18</td></tr>\n",
"</tbody>\n",
"</table>\n"
],
"text/latex": [
"A data.frame: 23 × 4\n",
"\\begin{tabular}{llll}\n",
" ID & onset & death & time\\\\\n",
" <int> & <dbl> & <dbl> & <dbl>\\\\\n",
"\\hline\n",
"\t 2 & 17 & 33 & 16\\\\\n",
"\t 3 & 19 & 48 & 29\\\\\n",
"\t 4 & 21 & 31 & 10\\\\\n",
"\t 5 & 21 & 47 & 26\\\\\n",
"\t 6 & 21 & 30 & 9\\\\\n",
"\t 7 & 22 & 41 & 19\\\\\n",
"\t 8 & 22 & 37 & 15\\\\\n",
"\t 9 & 23 & 32 & 9\\\\\n",
"\t 10 & 23 & 31 & 8\\\\\n",
"\t 11 & 23 & 35 & 12\\\\\n",
"\t 12 & 23 & 54 & 31\\\\\n",
"\t 13 & 23 & 30 & 7\\\\\n",
"\t 14 & 25 & 31 & 6\\\\\n",
"\t 15 & 26 & 33 & 7\\\\\n",
"\t 16 & 26 & 30 & 4\\\\\n",
"\t 17 & 26 & 33 & 7\\\\\n",
"\t 18 & 27 & 48 & 21\\\\\n",
"\t 19 & 31 & 39 & 8\\\\\n",
"\t 20 & 39 & 51 & 12\\\\\n",
"\t 21 & 41 & 48 & 7\\\\\n",
"\t 22 & 41 & 55 & 14\\\\\n",
"\t 23 & 42 & 60 & 18\\\\\n",
"\t 24 & 44 & 62 & 18\\\\\n",
"\\end{tabular}\n"
],
"text/markdown": [
"\n",
"A data.frame: 23 × 4\n",
"\n",
"| ID &lt;int&gt; | onset &lt;dbl&gt; | death &lt;dbl&gt; | time &lt;dbl&gt; |\n",
"|---|---|---|---|\n",
"| 2 | 17 | 33 | 16 |\n",
"| 3 | 19 | 48 | 29 |\n",
"| 4 | 21 | 31 | 10 |\n",
"| 5 | 21 | 47 | 26 |\n",
"| 6 | 21 | 30 | 9 |\n",
"| 7 | 22 | 41 | 19 |\n",
"| 8 | 22 | 37 | 15 |\n",
"| 9 | 23 | 32 | 9 |\n",
"| 10 | 23 | 31 | 8 |\n",
"| 11 | 23 | 35 | 12 |\n",
"| 12 | 23 | 54 | 31 |\n",
"| 13 | 23 | 30 | 7 |\n",
"| 14 | 25 | 31 | 6 |\n",
"| 15 | 26 | 33 | 7 |\n",
"| 16 | 26 | 30 | 4 |\n",
"| 17 | 26 | 33 | 7 |\n",
"| 18 | 27 | 48 | 21 |\n",
"| 19 | 31 | 39 | 8 |\n",
"| 20 | 39 | 51 | 12 |\n",
"| 21 | 41 | 48 | 7 |\n",
"| 22 | 41 | 55 | 14 |\n",
"| 23 | 42 | 60 | 18 |\n",
"| 24 | 44 | 62 | 18 |\n",
"\n"
],
"text/plain": [
" ID onset death time\n",
"1 2 17 33 16 \n",
"2 3 19 48 29 \n",
"3 4 21 31 10 \n",
"4 5 21 47 26 \n",
"5 6 21 30 9 \n",
"6 7 22 41 19 \n",
"7 8 22 37 15 \n",
"8 9 23 32 9 \n",
"9 10 23 31 8 \n",
"10 11 23 35 12 \n",
"11 12 23 54 31 \n",
"12 13 23 30 7 \n",
"13 14 25 31 6 \n",
"14 15 26 33 7 \n",
"15 16 26 30 4 \n",
"16 17 26 33 7 \n",
"17 18 27 48 21 \n",
"18 19 31 39 8 \n",
"19 20 39 51 12 \n",
"20 21 41 48 7 \n",
"21 22 41 55 14 \n",
"22 23 42 60 18 \n",
"23 24 44 62 18 "
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"df %<>% mutate(time = death - onset)\n",
"df "
]
},
{
"cell_type": "code",
"execution_count": 18,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"<table>\n",
"<caption>A grouped_df: 16 × 2</caption>\n",
"<thead>\n",
"\t<tr><th scope=col>time</th><th scope=col>n</th></tr>\n",
"\t<tr><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;int&gt;</th></tr>\n",
"</thead>\n",
"<tbody>\n",
"\t<tr><td> 4</td><td>1</td></tr>\n",
"\t<tr><td> 6</td><td>1</td></tr>\n",
"\t<tr><td> 7</td><td>4</td></tr>\n",
"\t<tr><td> 8</td><td>2</td></tr>\n",
"\t<tr><td> 9</td><td>2</td></tr>\n",
"\t<tr><td>10</td><td>1</td></tr>\n",
"\t<tr><td>12</td><td>2</td></tr>\n",
"\t<tr><td>14</td><td>1</td></tr>\n",
"\t<tr><td>15</td><td>1</td></tr>\n",
"\t<tr><td>16</td><td>1</td></tr>\n",
"\t<tr><td>18</td><td>2</td></tr>\n",
"\t<tr><td>19</td><td>1</td></tr>\n",
"\t<tr><td>21</td><td>1</td></tr>\n",
"\t<tr><td>26</td><td>1</td></tr>\n",
"\t<tr><td>29</td><td>1</td></tr>\n",
"\t<tr><td>31</td><td>1</td></tr>\n",
"</tbody>\n",
"</table>\n"
],
"text/latex": [
"A grouped\\_df: 16 × 2\n",
"\\begin{tabular}{ll}\n",
" time & n\\\\\n",
" <dbl> & <int>\\\\\n",
"\\hline\n",
"\t 4 & 1\\\\\n",
"\t 6 & 1\\\\\n",
"\t 7 & 4\\\\\n",
"\t 8 & 2\\\\\n",
"\t 9 & 2\\\\\n",
"\t 10 & 1\\\\\n",
"\t 12 & 2\\\\\n",
"\t 14 & 1\\\\\n",
"\t 15 & 1\\\\\n",
"\t 16 & 1\\\\\n",
"\t 18 & 2\\\\\n",
"\t 19 & 1\\\\\n",
"\t 21 & 1\\\\\n",
"\t 26 & 1\\\\\n",
"\t 29 & 1\\\\\n",
"\t 31 & 1\\\\\n",
"\\end{tabular}\n"
],
"text/markdown": [
"\n",
"A grouped_df: 16 × 2\n",
"\n",
"| time &lt;dbl&gt; | n &lt;int&gt; |\n",
"|---|---|\n",
"| 4 | 1 |\n",
"| 6 | 1 |\n",
"| 7 | 4 |\n",
"| 8 | 2 |\n",
"| 9 | 2 |\n",
"| 10 | 1 |\n",
"| 12 | 2 |\n",
"| 14 | 1 |\n",
"| 15 | 1 |\n",
"| 16 | 1 |\n",
"| 18 | 2 |\n",
"| 19 | 1 |\n",
"| 21 | 1 |\n",
"| 26 | 1 |\n",
"| 29 | 1 |\n",
"| 31 | 1 |\n",
"\n"
],
"text/plain": [
" time n\n",
"1 4 1\n",
"2 6 1\n",
"3 7 4\n",
"4 8 2\n",
"5 9 2\n",
"6 10 1\n",
"7 12 2\n",
"8 14 1\n",
"9 15 1\n",
"10 16 1\n",
"11 18 2\n",
"12 19 1\n",
"13 21 1\n",
"14 26 1\n",
"15 29 1\n",
"16 31 1"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"df %>% group_by(time) %>% count"
]
},
{
"cell_type": "code",
"execution_count": 25,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"<table>\n",
"<caption>A data.frame: 23 × 2</caption>\n",
"<thead>\n",
"\t<tr><th scope=col>ID</th><th scope=col>time</th></tr>\n",
"\t<tr><th scope=col>&lt;int&gt;</th><th scope=col>&lt;dbl&gt;</th></tr>\n",
"</thead>\n",
"<tbody>\n",
"\t<tr><td>16</td><td> 4</td></tr>\n",
"\t<tr><td>14</td><td> 6</td></tr>\n",
"\t<tr><td>13</td><td> 7</td></tr>\n",
"\t<tr><td>15</td><td> 7</td></tr>\n",
"\t<tr><td>17</td><td> 7</td></tr>\n",
"\t<tr><td>21</td><td> 7</td></tr>\n",
"\t<tr><td>10</td><td> 8</td></tr>\n",
"\t<tr><td>19</td><td> 8</td></tr>\n",
"\t<tr><td> 6</td><td> 9</td></tr>\n",
"\t<tr><td> 9</td><td> 9</td></tr>\n",
"\t<tr><td> 4</td><td>10</td></tr>\n",
"\t<tr><td>11</td><td>12</td></tr>\n",
"\t<tr><td>20</td><td>12</td></tr>\n",
"\t<tr><td>22</td><td>14</td></tr>\n",
"\t<tr><td> 8</td><td>15</td></tr>\n",
"\t<tr><td> 2</td><td>16</td></tr>\n",
"\t<tr><td>23</td><td>18</td></tr>\n",
"\t<tr><td>24</td><td>18</td></tr>\n",
"\t<tr><td> 7</td><td>19</td></tr>\n",
"\t<tr><td>18</td><td>21</td></tr>\n",
"\t<tr><td> 5</td><td>26</td></tr>\n",
"\t<tr><td> 3</td><td>29</td></tr>\n",
"\t<tr><td>12</td><td>31</td></tr>\n",
"</tbody>\n",
"</table>\n"
],
"text/latex": [
"A data.frame: 23 × 2\n",
"\\begin{tabular}{ll}\n",
" ID & time\\\\\n",
" <int> & <dbl>\\\\\n",
"\\hline\n",
"\t 16 & 4\\\\\n",
"\t 14 & 6\\\\\n",
"\t 13 & 7\\\\\n",
"\t 15 & 7\\\\\n",
"\t 17 & 7\\\\\n",
"\t 21 & 7\\\\\n",
"\t 10 & 8\\\\\n",
"\t 19 & 8\\\\\n",
"\t 6 & 9\\\\\n",
"\t 9 & 9\\\\\n",
"\t 4 & 10\\\\\n",
"\t 11 & 12\\\\\n",
"\t 20 & 12\\\\\n",
"\t 22 & 14\\\\\n",
"\t 8 & 15\\\\\n",
"\t 2 & 16\\\\\n",
"\t 23 & 18\\\\\n",
"\t 24 & 18\\\\\n",
"\t 7 & 19\\\\\n",
"\t 18 & 21\\\\\n",
"\t 5 & 26\\\\\n",
"\t 3 & 29\\\\\n",
"\t 12 & 31\\\\\n",
"\\end{tabular}\n"
],
"text/markdown": [
"\n",
"A data.frame: 23 × 2\n",
"\n",
"| ID &lt;int&gt; | time &lt;dbl&gt; |\n",
"|---|---|\n",
"| 16 | 4 |\n",
"| 14 | 6 |\n",
"| 13 | 7 |\n",
"| 15 | 7 |\n",
"| 17 | 7 |\n",
"| 21 | 7 |\n",
"| 10 | 8 |\n",
"| 19 | 8 |\n",
"| 6 | 9 |\n",
"| 9 | 9 |\n",
"| 4 | 10 |\n",
"| 11 | 12 |\n",
"| 20 | 12 |\n",
"| 22 | 14 |\n",
"| 8 | 15 |\n",
"| 2 | 16 |\n",
"| 23 | 18 |\n",
"| 24 | 18 |\n",
"| 7 | 19 |\n",
"| 18 | 21 |\n",
"| 5 | 26 |\n",
"| 3 | 29 |\n",
"| 12 | 31 |\n",
"\n"
],
"text/plain": [
" ID time\n",
"1 16 4 \n",
"2 14 6 \n",
"3 13 7 \n",
"4 15 7 \n",
"5 17 7 \n",
"6 21 7 \n",
"7 10 8 \n",
"8 19 8 \n",
"9 6 9 \n",
"10 9 9 \n",
"11 4 10 \n",
"12 11 12 \n",
"13 20 12 \n",
"14 22 14 \n",
"15 8 15 \n",
"16 2 16 \n",
"17 23 18 \n",
"18 24 18 \n",
"19 7 19 \n",
"20 18 21 \n",
"21 5 26 \n",
"22 3 29 \n",
"23 12 31 "
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"df %>% arrange(time) %>% select(-onset, -death) "
]
},
{
"cell_type": "code",
"execution_count": 40,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"<table>\n",
"<caption>A data.frame: 23 × 3</caption>\n",
"<thead>\n",
"\t<tr><th scope=col>rank</th><th scope=col>ID</th><th scope=col>time</th></tr>\n",
"\t<tr><th scope=col>&lt;int&gt;</th><th scope=col>&lt;int&gt;</th><th scope=col>&lt;dbl&gt;</th></tr>\n",
"</thead>\n",
"<tbody>\n",
"\t<tr><td> 1</td><td>16</td><td> 4</td></tr>\n",
"\t<tr><td> 2</td><td>14</td><td> 6</td></tr>\n",
"\t<tr><td> 3</td><td>13</td><td> 7</td></tr>\n",
"\t<tr><td> 4</td><td>15</td><td> 7</td></tr>\n",
"\t<tr><td> 5</td><td>17</td><td> 7</td></tr>\n",
"\t<tr><td> 6</td><td>21</td><td> 7</td></tr>\n",
"\t<tr><td> 7</td><td>10</td><td> 8</td></tr>\n",
"\t<tr><td> 8</td><td>19</td><td> 8</td></tr>\n",
"\t<tr><td> 9</td><td> 6</td><td> 9</td></tr>\n",
"\t<tr><td>10</td><td> 9</td><td> 9</td></tr>\n",
"\t<tr><td>11</td><td> 4</td><td>10</td></tr>\n",
"\t<tr><td>12</td><td>11</td><td>12</td></tr>\n",
"\t<tr><td>13</td><td>20</td><td>12</td></tr>\n",
"\t<tr><td>14</td><td>22</td><td>14</td></tr>\n",
"\t<tr><td>15</td><td> 8</td><td>15</td></tr>\n",
"\t<tr><td>16</td><td> 2</td><td>16</td></tr>\n",
"\t<tr><td>17</td><td>23</td><td>18</td></tr>\n",
"\t<tr><td>18</td><td>24</td><td>18</td></tr>\n",
"\t<tr><td>19</td><td> 7</td><td>19</td></tr>\n",
"\t<tr><td>20</td><td>18</td><td>21</td></tr>\n",
"\t<tr><td>21</td><td> 5</td><td>26</td></tr>\n",
"\t<tr><td>22</td><td> 3</td><td>29</td></tr>\n",
"\t<tr><td>23</td><td>12</td><td>31</td></tr>\n",
"</tbody>\n",
"</table>\n"
],
"text/latex": [
"A data.frame: 23 × 3\n",
"\\begin{tabular}{lll}\n",
" rank & ID & time\\\\\n",
" <int> & <int> & <dbl>\\\\\n",
"\\hline\n",
"\t 1 & 16 & 4\\\\\n",
"\t 2 & 14 & 6\\\\\n",
"\t 3 & 13 & 7\\\\\n",
"\t 4 & 15 & 7\\\\\n",
"\t 5 & 17 & 7\\\\\n",
"\t 6 & 21 & 7\\\\\n",
"\t 7 & 10 & 8\\\\\n",
"\t 8 & 19 & 8\\\\\n",
"\t 9 & 6 & 9\\\\\n",
"\t 10 & 9 & 9\\\\\n",
"\t 11 & 4 & 10\\\\\n",
"\t 12 & 11 & 12\\\\\n",
"\t 13 & 20 & 12\\\\\n",
"\t 14 & 22 & 14\\\\\n",
"\t 15 & 8 & 15\\\\\n",
"\t 16 & 2 & 16\\\\\n",
"\t 17 & 23 & 18\\\\\n",
"\t 18 & 24 & 18\\\\\n",
"\t 19 & 7 & 19\\\\\n",
"\t 20 & 18 & 21\\\\\n",
"\t 21 & 5 & 26\\\\\n",
"\t 22 & 3 & 29\\\\\n",
"\t 23 & 12 & 31\\\\\n",
"\\end{tabular}\n"
],
"text/markdown": [
"\n",
"A data.frame: 23 × 3\n",
"\n",
"| rank &lt;int&gt; | ID &lt;int&gt; | time &lt;dbl&gt; |\n",
"|---|---|---|\n",
"| 1 | 16 | 4 |\n",
"| 2 | 14 | 6 |\n",
"| 3 | 13 | 7 |\n",
"| 4 | 15 | 7 |\n",
"| 5 | 17 | 7 |\n",
"| 6 | 21 | 7 |\n",
"| 7 | 10 | 8 |\n",
"| 8 | 19 | 8 |\n",
"| 9 | 6 | 9 |\n",
"| 10 | 9 | 9 |\n",
"| 11 | 4 | 10 |\n",
"| 12 | 11 | 12 |\n",
"| 13 | 20 | 12 |\n",
"| 14 | 22 | 14 |\n",
"| 15 | 8 | 15 |\n",
"| 16 | 2 | 16 |\n",
"| 17 | 23 | 18 |\n",
"| 18 | 24 | 18 |\n",
"| 19 | 7 | 19 |\n",
"| 20 | 18 | 21 |\n",
"| 21 | 5 | 26 |\n",
"| 22 | 3 | 29 |\n",
"| 23 | 12 | 31 |\n",
"\n"
],
"text/plain": [
" rank ID time\n",
"1 1 16 4 \n",
"2 2 14 6 \n",
"3 3 13 7 \n",
"4 4 15 7 \n",
"5 5 17 7 \n",
"6 6 21 7 \n",
"7 7 10 8 \n",
"8 8 19 8 \n",
"9 9 6 9 \n",
"10 10 9 9 \n",
"11 11 4 10 \n",
"12 12 11 12 \n",
"13 13 20 12 \n",
"14 14 22 14 \n",
"15 15 8 15 \n",
"16 16 2 16 \n",
"17 17 23 18 \n",
"18 18 24 18 \n",
"19 19 7 19 \n",
"20 20 18 21 \n",
"21 21 5 26 \n",
"22 22 3 29 \n",
"23 23 12 31 "
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"df %>% arrange(time) %>% mutate(rank = 1:n()) %>% \n",
" select(-onset, -death) %>% select(rank, everything())"
]
},
{
"cell_type": "code",
"execution_count": 46,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"<table>\n",
"<caption>A data.frame: 23 × 4</caption>\n",
"<thead>\n",
"\t<tr><th scope=col>ID</th><th scope=col>time</th><th scope=col>weight</th><th scope=col>percentile_below</th></tr>\n",
"\t<tr><th scope=col>&lt;int&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;dbl&gt;</th></tr>\n",
"</thead>\n",
"<tbody>\n",
"\t<tr><td>16</td><td> 4</td><td>0.04347826</td><td>0.00000000</td></tr>\n",
"\t<tr><td>14</td><td> 6</td><td>0.04347826</td><td>0.04347826</td></tr>\n",
"\t<tr><td>13</td><td> 7</td><td>0.04347826</td><td>0.08695652</td></tr>\n",
"\t<tr><td>15</td><td> 7</td><td>0.04347826</td><td>0.13043478</td></tr>\n",
"\t<tr><td>17</td><td> 7</td><td>0.04347826</td><td>0.17391304</td></tr>\n",
"\t<tr><td>21</td><td> 7</td><td>0.04347826</td><td>0.21739130</td></tr>\n",
"\t<tr><td>10</td><td> 8</td><td>0.04347826</td><td>0.26086957</td></tr>\n",
"\t<tr><td>19</td><td> 8</td><td>0.04347826</td><td>0.30434783</td></tr>\n",
"\t<tr><td> 6</td><td> 9</td><td>0.04347826</td><td>0.34782609</td></tr>\n",
"\t<tr><td> 9</td><td> 9</td><td>0.04347826</td><td>0.39130435</td></tr>\n",
"\t<tr><td> 4</td><td>10</td><td>0.04347826</td><td>0.43478261</td></tr>\n",
"\t<tr><td>11</td><td>12</td><td>0.04347826</td><td>0.47826087</td></tr>\n",
"\t<tr><td>20</td><td>12</td><td>0.04347826</td><td>0.52173913</td></tr>\n",
"\t<tr><td>22</td><td>14</td><td>0.04347826</td><td>0.56521739</td></tr>\n",
"\t<tr><td> 8</td><td>15</td><td>0.04347826</td><td>0.60869565</td></tr>\n",
"\t<tr><td> 2</td><td>16</td><td>0.04347826</td><td>0.65217391</td></tr>\n",
"\t<tr><td>23</td><td>18</td><td>0.04347826</td><td>0.69565217</td></tr>\n",
"\t<tr><td>24</td><td>18</td><td>0.04347826</td><td>0.73913043</td></tr>\n",
"\t<tr><td> 7</td><td>19</td><td>0.04347826</td><td>0.78260870</td></tr>\n",
"\t<tr><td>18</td><td>21</td><td>0.04347826</td><td>0.82608696</td></tr>\n",
"\t<tr><td> 5</td><td>26</td><td>0.04347826</td><td>0.86956522</td></tr>\n",
"\t<tr><td> 3</td><td>29</td><td>0.04347826</td><td>0.91304348</td></tr>\n",
"\t<tr><td>12</td><td>31</td><td>0.04347826</td><td>0.95652174</td></tr>\n",
"</tbody>\n",
"</table>\n"
],
"text/latex": [
"A data.frame: 23 × 4\n",
"\\begin{tabular}{llll}\n",
" ID & time & weight & percentile\\_below\\\\\n",
" <int> & <dbl> & <dbl> & <dbl>\\\\\n",
"\\hline\n",
"\t 16 & 4 & 0.04347826 & 0.00000000\\\\\n",
"\t 14 & 6 & 0.04347826 & 0.04347826\\\\\n",
"\t 13 & 7 & 0.04347826 & 0.08695652\\\\\n",
"\t 15 & 7 & 0.04347826 & 0.13043478\\\\\n",
"\t 17 & 7 & 0.04347826 & 0.17391304\\\\\n",
"\t 21 & 7 & 0.04347826 & 0.21739130\\\\\n",
"\t 10 & 8 & 0.04347826 & 0.26086957\\\\\n",
"\t 19 & 8 & 0.04347826 & 0.30434783\\\\\n",
"\t 6 & 9 & 0.04347826 & 0.34782609\\\\\n",
"\t 9 & 9 & 0.04347826 & 0.39130435\\\\\n",
"\t 4 & 10 & 0.04347826 & 0.43478261\\\\\n",
"\t 11 & 12 & 0.04347826 & 0.47826087\\\\\n",
"\t 20 & 12 & 0.04347826 & 0.52173913\\\\\n",
"\t 22 & 14 & 0.04347826 & 0.56521739\\\\\n",
"\t 8 & 15 & 0.04347826 & 0.60869565\\\\\n",
"\t 2 & 16 & 0.04347826 & 0.65217391\\\\\n",
"\t 23 & 18 & 0.04347826 & 0.69565217\\\\\n",
"\t 24 & 18 & 0.04347826 & 0.73913043\\\\\n",
"\t 7 & 19 & 0.04347826 & 0.78260870\\\\\n",
"\t 18 & 21 & 0.04347826 & 0.82608696\\\\\n",
"\t 5 & 26 & 0.04347826 & 0.86956522\\\\\n",
"\t 3 & 29 & 0.04347826 & 0.91304348\\\\\n",
"\t 12 & 31 & 0.04347826 & 0.95652174\\\\\n",
"\\end{tabular}\n"
],
"text/markdown": [
"\n",
"A data.frame: 23 × 4\n",
"\n",
"| ID &lt;int&gt; | time &lt;dbl&gt; | weight &lt;dbl&gt; | percentile_below &lt;dbl&gt; |\n",
"|---|---|---|---|\n",
"| 16 | 4 | 0.04347826 | 0.00000000 |\n",
"| 14 | 6 | 0.04347826 | 0.04347826 |\n",
"| 13 | 7 | 0.04347826 | 0.08695652 |\n",
"| 15 | 7 | 0.04347826 | 0.13043478 |\n",
"| 17 | 7 | 0.04347826 | 0.17391304 |\n",
"| 21 | 7 | 0.04347826 | 0.21739130 |\n",
"| 10 | 8 | 0.04347826 | 0.26086957 |\n",
"| 19 | 8 | 0.04347826 | 0.30434783 |\n",
"| 6 | 9 | 0.04347826 | 0.34782609 |\n",
"| 9 | 9 | 0.04347826 | 0.39130435 |\n",
"| 4 | 10 | 0.04347826 | 0.43478261 |\n",
"| 11 | 12 | 0.04347826 | 0.47826087 |\n",
"| 20 | 12 | 0.04347826 | 0.52173913 |\n",
"| 22 | 14 | 0.04347826 | 0.56521739 |\n",
"| 8 | 15 | 0.04347826 | 0.60869565 |\n",
"| 2 | 16 | 0.04347826 | 0.65217391 |\n",
"| 23 | 18 | 0.04347826 | 0.69565217 |\n",
"| 24 | 18 | 0.04347826 | 0.73913043 |\n",
"| 7 | 19 | 0.04347826 | 0.78260870 |\n",
"| 18 | 21 | 0.04347826 | 0.82608696 |\n",
"| 5 | 26 | 0.04347826 | 0.86956522 |\n",
"| 3 | 29 | 0.04347826 | 0.91304348 |\n",
"| 12 | 31 | 0.04347826 | 0.95652174 |\n",
"\n"
],
"text/plain": [
" ID time weight percentile_below\n",
"1 16 4 0.04347826 0.00000000 \n",
"2 14 6 0.04347826 0.04347826 \n",
"3 13 7 0.04347826 0.08695652 \n",
"4 15 7 0.04347826 0.13043478 \n",
"5 17 7 0.04347826 0.17391304 \n",
"6 21 7 0.04347826 0.21739130 \n",
"7 10 8 0.04347826 0.26086957 \n",
"8 19 8 0.04347826 0.30434783 \n",
"9 6 9 0.04347826 0.34782609 \n",
"10 9 9 0.04347826 0.39130435 \n",
"11 4 10 0.04347826 0.43478261 \n",
"12 11 12 0.04347826 0.47826087 \n",
"13 20 12 0.04347826 0.52173913 \n",
"14 22 14 0.04347826 0.56521739 \n",
"15 8 15 0.04347826 0.60869565 \n",
"16 2 16 0.04347826 0.65217391 \n",
"17 23 18 0.04347826 0.69565217 \n",
"18 24 18 0.04347826 0.73913043 \n",
"19 7 19 0.04347826 0.78260870 \n",
"20 18 21 0.04347826 0.82608696 \n",
"21 5 26 0.04347826 0.86956522 \n",
"22 3 29 0.04347826 0.91304348 \n",
"23 12 31 0.04347826 0.95652174 "
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"df %>% arrange(time) %>% select(-onset, -death) %>% \n",
" mutate(weight = 1/nrow(df), percentile_below = cumsum(weight), \n",
" percentile_below = lag(percentile_below, default=0))"
]
},
{
"cell_type": "code",
"execution_count": 41,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"<table>\n",
"<caption>A data.frame: 23 × 2</caption>\n",
"<thead>\n",
"\t<tr><th scope=col>ID</th><th scope=col>time</th></tr>\n",
"\t<tr><th scope=col>&lt;int&gt;</th><th scope=col>&lt;dbl&gt;</th></tr>\n",
"</thead>\n",
"<tbody>\n",
"\t<tr><td>16</td><td> 4</td></tr>\n",
"\t<tr><td>14</td><td> 6</td></tr>\n",
"\t<tr><td>13</td><td> 7</td></tr>\n",
"\t<tr><td>15</td><td> 7</td></tr>\n",
"\t<tr><td>17</td><td> 7</td></tr>\n",
"\t<tr><td>21</td><td> 7</td></tr>\n",
"\t<tr><td>10</td><td> 8</td></tr>\n",
"\t<tr><td>19</td><td> 8</td></tr>\n",
"\t<tr><td> 6</td><td> 9</td></tr>\n",
"\t<tr><td> 9</td><td> 9</td></tr>\n",
"\t<tr><td> 4</td><td>10</td></tr>\n",
"\t<tr><td>11</td><td>12</td></tr>\n",
"\t<tr><td>20</td><td>12</td></tr>\n",
"\t<tr><td>22</td><td>14</td></tr>\n",
"\t<tr><td> 8</td><td>15</td></tr>\n",
"\t<tr><td> 2</td><td>16</td></tr>\n",
"\t<tr><td>23</td><td>18</td></tr>\n",
"\t<tr><td>24</td><td>18</td></tr>\n",
"\t<tr><td> 7</td><td>19</td></tr>\n",
"\t<tr><td>18</td><td>21</td></tr>\n",
"\t<tr><td> 5</td><td>26</td></tr>\n",
"\t<tr><td> 3</td><td>29</td></tr>\n",
"\t<tr><td>12</td><td>31</td></tr>\n",
"</tbody>\n",
"</table>\n"
],
"text/latex": [
"A data.frame: 23 × 2\n",
"\\begin{tabular}{ll}\n",
" ID & time\\\\\n",
" <int> & <dbl>\\\\\n",
"\\hline\n",
"\t 16 & 4\\\\\n",
"\t 14 & 6\\\\\n",
"\t 13 & 7\\\\\n",
"\t 15 & 7\\\\\n",
"\t 17 & 7\\\\\n",
"\t 21 & 7\\\\\n",
"\t 10 & 8\\\\\n",
"\t 19 & 8\\\\\n",
"\t 6 & 9\\\\\n",
"\t 9 & 9\\\\\n",
"\t 4 & 10\\\\\n",
"\t 11 & 12\\\\\n",
"\t 20 & 12\\\\\n",
"\t 22 & 14\\\\\n",
"\t 8 & 15\\\\\n",
"\t 2 & 16\\\\\n",
"\t 23 & 18\\\\\n",
"\t 24 & 18\\\\\n",
"\t 7 & 19\\\\\n",
"\t 18 & 21\\\\\n",
"\t 5 & 26\\\\\n",
"\t 3 & 29\\\\\n",
"\t 12 & 31\\\\\n",
"\\end{tabular}\n"
],
"text/markdown": [
"\n",
"A data.frame: 23 × 2\n",
"\n",
"| ID &lt;int&gt; | time &lt;dbl&gt; |\n",
"|---|---|\n",
"| 16 | 4 |\n",
"| 14 | 6 |\n",
"| 13 | 7 |\n",
"| 15 | 7 |\n",
"| 17 | 7 |\n",
"| 21 | 7 |\n",
"| 10 | 8 |\n",
"| 19 | 8 |\n",
"| 6 | 9 |\n",
"| 9 | 9 |\n",
"| 4 | 10 |\n",
"| 11 | 12 |\n",
"| 20 | 12 |\n",
"| 22 | 14 |\n",
"| 8 | 15 |\n",
"| 2 | 16 |\n",
"| 23 | 18 |\n",
"| 24 | 18 |\n",
"| 7 | 19 |\n",
"| 18 | 21 |\n",
"| 5 | 26 |\n",
"| 3 | 29 |\n",
"| 12 | 31 |\n",
"\n"
],
"text/plain": [
" ID time\n",
"1 16 4 \n",
"2 14 6 \n",
"3 13 7 \n",
"4 15 7 \n",
"5 17 7 \n",
"6 21 7 \n",
"7 10 8 \n",
"8 19 8 \n",
"9 6 9 \n",
"10 9 9 \n",
"11 4 10 \n",
"12 11 12 \n",
"13 20 12 \n",
"14 22 14 \n",
"15 8 15 \n",
"16 2 16 \n",
"17 23 18 \n",
"18 24 18 \n",
"19 7 19 \n",
"20 18 21 \n",
"21 5 26 \n",
"22 3 29 \n",
"23 12 31 "
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"df %>% arrange(time) %>% select(-onset, -death)"
]
},
{
"cell_type": "code",
"execution_count": 49,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"13.6086956521739"
],
"text/latex": [
"13.6086956521739"
],
"text/markdown": [
"13.6086956521739"
],
"text/plain": [
"[1] 13.6087"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"mean(df$time)"
]
},
{
"cell_type": "code",
"execution_count": 54,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"5.1"
],
"text/latex": [
"5.1"
],
"text/markdown": [
"5.1"
],
"text/plain": [
"[1] 5.1"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"c(3, 4, 5, 5, 5, 5, 5, 6, 6, 7) %>% mean"
]
},
{
"cell_type": "code",
"execution_count": 61,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"5.6"
],
"text/latex": [
"5.6"
],
"text/markdown": [
"5.6"
],
"text/plain": [
"[1] 5.6"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"c(1, 2, 3, 4, 4, 5, 6, 8, 10, 13) %>% mean"
]
},
{
"cell_type": "code",
"execution_count": 62,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"<style>\n",
".list-inline {list-style: none; margin:0; padding: 0}\n",
".list-inline>li {display: inline-block}\n",
".list-inline>li:not(:last-child)::after {content: \"\\00b7\"; padding: 0 .5ex}\n",
"</style>\n",
"<ol class=list-inline><li>4</li><li>31</li></ol>\n"
],
"text/latex": [
"\\begin{enumerate*}\n",
"\\item 4\n",
"\\item 31\n",
"\\end{enumerate*}\n"
],
"text/markdown": [
"1. 4\n",
"2. 31\n",
"\n",
"\n"
],
"text/plain": [
"[1] 4 31"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"c(min(df$time), max(df$time))"
]
},
{
"cell_type": "code",
"execution_count": 64,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"<style>\n",
".dl-inline {width: auto; margin:0; padding: 0}\n",
".dl-inline>dt, .dl-inline>dd {float: none; width: auto; display: inline-block}\n",
".dl-inline>dt::after {content: \":\\0020\"; padding-right: .5ex}\n",
".dl-inline>dt:not(:first-of-type) {padding-left: .5ex}\n",
"</style><dl class=dl-inline><dt>25%</dt><dd>7.5</dd><dt>75%</dt><dd>18</dd></dl>\n"
],
"text/latex": [
"\\begin{description*}\n",
"\\item[25\\textbackslash{}\\%] 7.5\n",
"\\item[75\\textbackslash{}\\%] 18\n",
"\\end{description*}\n"
],
"text/markdown": [
"25%\n",
": 7.575%\n",
": 18\n",
"\n"
],
"text/plain": [
" 25% 75% \n",
" 7.5 18.0 "
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"quantile(df$time, c(.25, .75))"
]
},
{
"cell_type": "code",
"execution_count": 65,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"562949953421312"
],
"text/latex": [
"562949953421312"
],
"text/markdown": [
"562949953421312"
],
"text/plain": [
"[1] 5.6295e+14"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"2^50/2"
]
},
{
"cell_type": "code",
"execution_count": 66,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"1225"
],
"text/latex": [
"1225"
],
"text/markdown": [
"1225"
],
"text/plain": [
"[1] 1225"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"25*49"
]
},
{
"cell_type": "code",
"execution_count": 67,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"9"
],
"text/latex": [
"9"
],
"text/markdown": [
"9"
],
"text/plain": [
"[1] 9"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"mean(c(1,6,10,12,16))"
]
},
{
"cell_type": "code",
"execution_count": 68,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"33"
],
"text/latex": [
"33"
],
"text/markdown": [
"33"
],
"text/plain": [
"[1] 33"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"var(c(1,6,10,12,16))"
]
},
{
"cell_type": "code",
"execution_count": 70,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"5.74456264653803"
],
"text/latex": [
"5.74456264653803"
],
"text/markdown": [
"5.74456264653803"
],
"text/plain": [
"[1] 5.744563"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"sd(c(1,6,10,12,16))"
]
},
{
"cell_type": "code",
"execution_count": 78,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"3.65"
],
"text/latex": [
"3.65"
],
"text/markdown": [
"3.65"
],
"text/plain": [
"[1] 3.65"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"sd(c(1,3,7,9)) %>% round(2)"
]
},
{
"cell_type": "code",
"execution_count": 79,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"4.32"
],
"text/latex": [
"4.32"
],
"text/markdown": [
"4.32"
],
"text/plain": [
"[1] 4.32"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"sd(c(1,3,5,11)) %>% round(2)"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
}
],
"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": "4.0.2"
}
},
"nbformat": 4,
"nbformat_minor": 4
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment