{ "cells": [ { "cell_type": "code", "execution_count": 1, "id": "4a89aebf", "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import pandas as pd\n", "import matplotlib.pyplot as plt\n", "import json\n", "import os" ] }, { "cell_type": "code", "execution_count": 2, "id": "ce41f624", "metadata": {}, "outputs": [], "source": [ "samples = []\n", "values = []\n", "\n", "files = os.listdir('out')\n", "files.sort()\n", "\n", "for filename in files:\n", " if filename.endswith(\".json\"):\n", " with open('out/'+filename) as output:\n", " output_json = json.load(output)\n", " samples += output_json['x']\n", " values += output_json['y']" ] }, { "cell_type": "code", "execution_count": 3, "id": "1e2edd6f", "metadata": {}, "outputs": [], "source": [ "df = pd.DataFrame(data={'samples': samples, 'avg': 0, 'std': 0, 'values': values})" ] }, { "cell_type": "code", "execution_count": 4, "id": "67cd6668", "metadata": {}, "outputs": [], "source": [ "df['avg'] = df.apply(lambda row: np.array(row['values']).sum() / len(row['values']), axis=1)\n", "df['std'] = df.apply(lambda row: np.sqrt(np.array(row['values']).var()), axis=1)" ] }, { "cell_type": "code", "execution_count": 5, "id": "3080eb16", "metadata": { "scrolled": false }, "outputs": [ { "data": { "text/html": [ "
| \n", " | samples | \n", "avg | \n", "std | \n", "values | \n", "
|---|---|---|---|---|
| 0 | \n", "1 | \n", "1.599000 | \n", "2.799410 | \n", "[0.0, 0.0, 0.0, 6.500000000000005, 0.0, 6.5000... | \n", "
| 1 | \n", "2 | \n", "1.631500 | \n", "2.070829 | \n", "[0.0, 0.0, 3.2500000000000027, 6.5000000000000... | \n", "
| 2 | \n", "5 | \n", "1.664000 | \n", "1.272142 | \n", "[2.6000000000000023, 3.900000000000003, 0.0, 2... | \n", "
| 3 | \n", "10 | \n", "1.654900 | \n", "0.908579 | \n", "[0.6500000000000006, 3.2500000000000027, 1.950... | \n", "
| 4 | \n", "20 | \n", "1.133600 | \n", "0.673453 | \n", "[0.0, 1.3000000000000012, 0.6500000000000006, ... | \n", "
| 5 | \n", "50 | \n", "1.220960 | \n", "0.431168 | \n", "[1.1700000000000008, 1.1700000000000008, 1.040... | \n", "
| 6 | \n", "100 | \n", "1.106820 | \n", "0.310923 | \n", "[1.1700000000000008, 1.235000000000001, 1.4950... | \n", "
| 7 | \n", "200 | \n", "1.097460 | \n", "0.212839 | \n", "[0.9100000000000008, 1.4625000000000012, 0.975... | \n", "
| 8 | \n", "500 | \n", "1.102712 | \n", "0.136280 | \n", "[0.9880000000000008, 1.0790000000000008, 1.001... | \n", "
| 9 | \n", "1000 | \n", "1.106060 | \n", "0.100732 | \n", "[1.0075000000000007, 1.209000000000001, 0.9555... | \n", "
| 10 | \n", "2000 | \n", "1.112072 | \n", "0.070502 | \n", "[1.277250000000001, 1.098500000000001, 1.08875... | \n", "
| 11 | \n", "5000 | \n", "1.108662 | \n", "0.043625 | \n", "[1.090700000000001, 1.1089000000000009, 1.1232... | \n", "
| 12 | \n", "10000 | \n", "1.105114 | \n", "0.032244 | \n", "[1.1180000000000008, 1.0796500000000009, 1.122... | \n", "
| 13 | \n", "20000 | \n", "1.108188 | \n", "0.022374 | \n", "[1.0741250000000009, 1.123525000000001, 1.1371... | \n", "
| 14 | \n", "50000 | \n", "1.106677 | \n", "0.014118 | \n", "[1.1052600000000008, 1.1278800000000009, 1.115... | \n", "
| 15 | \n", "100000 | \n", "1.107418 | \n", "0.010126 | \n", "[1.1035700000000008, 1.0917400000000008, 1.111... | \n", "
| 16 | \n", "200000 | \n", "1.106957 | \n", "0.007035 | \n", "[1.108770000000001, 1.1150425000000008, 1.1036... | \n", "
| 17 | \n", "500000 | \n", "1.106960 | \n", "0.004477 | \n", "[1.102660000000001, 1.099852000000001, 1.11181... | \n", "
| 18 | \n", "1000000 | \n", "1.107122 | \n", "0.003119 | \n", "[1.1097190000000008, 1.106612000000001, 1.1135... | \n", "
| 19 | \n", "2000000 | \n", "1.107152 | \n", "0.002248 | \n", "[1.1059262500000009, 1.1094655000000009, 1.107... | \n", "
| 20 | \n", "5000000 | \n", "1.107034 | \n", "0.001477 | \n", "[1.105048100000001, 1.106229800000001, 1.11045... | \n", "
| 21 | \n", "10000000 | \n", "1.107038 | \n", "0.000980 | \n", "[1.107393300000001, 1.107971150000001, 1.10760... | \n", "
| 22 | \n", "20000000 | \n", "1.107093 | \n", "0.000682 | \n", "[1.107080650000001, 1.1058697000000008, 1.1065... | \n", "
| 23 | \n", "50000000 | \n", "1.107103 | \n", "0.000438 | \n", "[1.106985230000001, 1.106864330000001, 1.10673... | \n", "
| 24 | \n", "100000000 | \n", "1.107097 | \n", "0.000300 | \n", "[1.107135120000001, 1.106441765000001, 1.10693... | \n", "
| 25 | \n", "200000000 | \n", "1.107100 | \n", "0.000229 | \n", "[1.107198267500001, 1.1066414125000008, 1.1074... | \n", "
| 26 | \n", "500000000 | \n", "1.107083 | \n", "0.000136 | \n", "[1.107150577000001, 1.1068575310000008, 1.1071... | \n", "
| 27 | \n", "1000000000 | \n", "1.107091 | \n", "0.000095 | \n", "[1.107090458500001, 1.107013505000001, 1.10705... | \n", "
| 28 | \n", "2000000000 | \n", "1.107100 | \n", "0.000067 | \n", "[1.106987911250001, 1.1070893340000008, 1.1071... | \n", "
| 29 | \n", "5000000000 | \n", "1.107091 | \n", "0.000044 | \n", "[1.107101422700001, 1.107028651300001, 1.10703... | \n", "
| 30 | \n", "10000000000 | \n", "1.107071 | \n", "0.000040 | \n", "[1.1070528943500009, 1.107091172850001, 1.1070... | \n", "
| 31 | \n", "20000000000 | \n", "1.107083 | \n", "0.000028 | \n", "[1.107116142600001, 1.107052413675001, 1.10708... | \n", "
| 32 | \n", "50000000000 | \n", "1.107088 | \n", "0.000010 | \n", "[1.107083811470001, 1.107081873690001, 1.10711... | \n", "
| 33 | \n", "100000000000 | \n", "1.107081 | \n", "0.000000 | \n", "[1.107080884000001] | \n", "
| 34 | \n", "200000000000 | \n", "1.107081 | \n", "0.000000 | \n", "[1.107080584805001] | \n", "
| 35 | \n", "500000000000 | \n", "1.107090 | \n", "0.000000 | \n", "[1.107089721517001] | \n", "