|
8 | 8 | "outputs": [], |
9 | 9 | "source": [ |
10 | 10 | "import pytest\n", |
11 | | - "import lqrt\n", |
12 | 11 | "import numpy as np\n", |
13 | | - "import scipy as sp\n", |
14 | | - "import pandas as pd" |
| 12 | + "from math import gamma" |
15 | 13 | ] |
16 | 14 | }, |
17 | 15 | { |
|
21 | 19 | "metadata": {}, |
22 | 20 | "outputs": [], |
23 | 21 | "source": [ |
24 | | - "from dabest._stats_tools import effsize\n", |
25 | | - "from dabest import Dabest, PermutationTest" |
| 22 | + "from dabest import Dabest, PermutationTest\n", |
| 23 | + "from data.mocked_data_test_06 import df_test, df_test_control, df_test_treatment1, dabest_default_kwargs" |
26 | 24 | ] |
27 | 25 | }, |
28 | 26 | { |
|
32 | 30 | "metadata": {}, |
33 | 31 | "outputs": [], |
34 | 32 | "source": [ |
35 | | - "# Data for tests.\n", |
36 | | - "# See: Asheber Abebe. Introduction to Design and Analysis of Experiments \n", |
37 | | - "# with the SAS, from Example: Two-way RM Design Pg 137.\n", |
38 | | - "hr = [72, 78, 71, 72, 66, 74, 62, 69, 69, 66, 84, 80, 72, 65, 75, 71, \n", |
39 | | - " 86, 83, 82, 83, 79, 83, 73, 75, 73, 62, 90, 81, 72, 62, 69, 70]\n", |
40 | | - "\n", |
41 | | - "# Add experiment column\n", |
42 | | - "e1 = np.repeat('Treatment1', 8).tolist()\n", |
43 | | - "e2 = np.repeat('Control', 8).tolist()\n", |
44 | | - "experiment = e1 + e2 + e1 + e2\n", |
45 | | - "\n", |
46 | | - "# Add a `Drug` column as the first variable\n", |
47 | | - "d1 = np.repeat('AX23', 8).tolist()\n", |
48 | | - "d2 = np.repeat('CONTROL', 8).tolist()\n", |
49 | | - "drug = d1 + d2 + d1 + d2\n", |
50 | | - "\n", |
51 | | - "# Add a `Time` column as the second variable\n", |
52 | | - "t1 = np.repeat('T1', 16).tolist()\n", |
53 | | - "t2 = np.repeat('T2', 16).tolist()\n", |
54 | | - "time = t1 + t2\n", |
55 | | - "\n", |
56 | | - "# Add an `id` column for paired data plotting.\n", |
57 | | - "id_col = []\n", |
58 | | - "for i in range(1, 9):\n", |
59 | | - " id_col.append(str(i)+\"a\")\n", |
60 | | - "for i in range(1, 9):\n", |
61 | | - " id_col.append(str(i)+\"c\")\n", |
62 | | - "id_col.extend(id_col)\n", |
63 | | - "\n", |
64 | | - "# Combine samples and gender into a DataFrame.\n", |
65 | | - "df_test = pd.DataFrame({'ID' : id_col,\n", |
66 | | - " 'Drug' : drug,\n", |
67 | | - " 'Time' : time, \n", |
68 | | - " 'Experiment': experiment,\n", |
69 | | - " 'Heart Rate': hr\n", |
70 | | - " })\n", |
71 | | - "\n", |
72 | | - "\n", |
73 | | - "df_test_control = df_test[df_test[\"Experiment\"]==\"Control\"]\n", |
74 | | - "df_test_control = df_test_control.pivot(index=\"ID\", columns=\"Time\", values=\"Heart Rate\")\n", |
75 | | - "\n", |
76 | | - "\n", |
77 | | - "df_test_treatment1 = df_test[df_test[\"Experiment\"]==\"Treatment1\"]\n", |
78 | | - "df_test_treatment1 = df_test_treatment1.pivot(index=\"ID\", columns=\"Time\", values=\"Heart Rate\")\n", |
79 | | - "\n", |
80 | | - "\n", |
81 | | - "# kwargs for Dabest class init.\n", |
82 | | - "dabest_default_kwargs = dict(ci=95, \n", |
83 | | - " resamples=5000, random_seed=12345,\n", |
84 | | - " idx=None, proportional=False, mini_meta=False\n", |
85 | | - " )\n", |
86 | | - "\n", |
87 | 33 | "# example of unpaired delta-delta calculation\n", |
88 | 34 | "unpaired = Dabest(data = df_test, x = [\"Time\", \"Drug\"], y = \"Heart Rate\", \n", |
89 | 35 | " delta2 = True, experiment = \"Experiment\",\n", |
|
333 | 279 | "metadata": {}, |
334 | 280 | "outputs": [], |
335 | 281 | "source": [ |
336 | | - "from math import gamma\n", |
337 | 282 | "hedges_g = unpaired.hedges_g.results['difference'].to_list()\n", |
338 | 283 | "a = 8*2-2\n", |
339 | 284 | "fac = gamma(a/2)/(np.sqrt(a/2)*gamma((a-1)/2))\n", |
|
360 | 305 | "metadata": {}, |
361 | 306 | "outputs": [], |
362 | 307 | "source": [ |
363 | | - "from math import gamma\n", |
364 | 308 | "hedges_g = paired.hedges_g.results['difference'].to_list()\n", |
365 | 309 | "a = 8*2-2\n", |
366 | 310 | "fac = gamma(a/2)/(np.sqrt(a/2)*gamma((a-1)/2))\n", |
|
387 | 331 | "metadata": {}, |
388 | 332 | "outputs": [], |
389 | 333 | "source": [ |
390 | | - "from math import gamma\n", |
391 | 334 | "hedges_g = paired_specified_level.hedges_g.results['difference'].to_list()\n", |
392 | 335 | "a = 8*2-2\n", |
393 | 336 | "fac = gamma(a/2)/(np.sqrt(a/2)*gamma((a-1)/2))\n", |
|
0 commit comments