In [4]:
import tables
with tables.open_file('../Data/SingleSourceGaussian/_PDUCB_kappa=70 gamma=-1e-5_exp_.0.h5') as data:
positions = data.root.positions.read()[0]
values = data.root.plume_measurements.read().T
In [5]:
print(positions.shape)
print(values.shape)
(1000, 3)
(1000, 1)
In [9]:
import plume.prediction as pred
In [28]:
reload(pred)
m = pred.Test(positions, values)
m.optimize()
[ 1. 15.]
(array([[-4032.96883807]]), array([ 478.64804255, -652.43576013]))
[ 0.99513519 15.99998817]
(array([[-4356.1536558]]), array([ 473.56976192, -630.46718979]))
[ 0.97567596 19.99994083]
(array([[-5512.81664291]]), array([ 433.73094331, -512.52807766]))
[ 0.1 38.06807734]
(array([[ 122592.7516661]]), array([-773570.97129183, 94343.70346406]))
[ 0.68640011 25.96867333]
(array([[-6844.7749728]]), array([ 422.57713208, -257.98629098]))
[ 0.1 32.01913417]
(array([[ 1974.97591319]]), array([-67084.36535201, 9873.93450678]))
[ 0.488315 28.01251028]
(array([[-7022.00338742]]), array([ 121.18049883, 125.08913229]))
[ 0.1 27.05938265]
(array([[-7079.78976924]]), array([-1936.62421231, 510.95954106]))
[ 0.33860266 27.64503807]
(array([[-7065.91070291]]), array([ 125.92511209, 148.74316045]))
[ 0.19394854 27.28998139]
(array([[-7106.43494645]]), array([-145.81948372, 257.23999976]))
[ 0.15489207 27.19411645]
(array([[-7109.09258225]]), array([-473.41024697, 326.52955136]))
[ 0.1 0.1]
(array([[-145.09617853]]), array([ 4191.18413781, -17.99949944]))
[ 0.14774728 23.66752822]
(array([[-6982.85676459]]), array([ 2053.49898417, -295.59954723]))
[ 0.15200632 25.76973922]
(array([[-7186.06531494]]), array([ 1094.4926673 , -64.86783161]))
[ 0.1 25.79554801]
(array([[-7241.04842648]]), array([ 914.80181891, 50.84626507]))
[ 0.1 25.74877581]
(array([[-7242.11111707]]), array([ 985.82558218, 39.22522748]))
[ 0.1 25.59090264]
(array([[-7243.75782669]]), array([ 1211.7751299 , 2.19519748]))
[ 0.1 25.58154368]
(array([[-7243.76517459]]), array([ 1.22452876e+03, 1.02287283e-01]))
[ 0.1 25.58108627]
(array([[-7243.76528137]]), array([ 1.22515031e+03, 2.78384210e-04]))
[ 0.1 25.58108503]
(array([[-7243.76528427]]), array([ 1.22515201e+03, 5.72500912e-06]))
[ 0.1 25.58108503]
In [30]:
for i in range(5, 150):
gp = pred.OnlineGP(pred.RBFKernel(200, 1.0), 1e-10)
gp.fit(positions[:i], values[:i])
gp.calc_neg_log_likelihood()
(array([[ 666.92450715]]), array([ 2.48753546, 7.94972867]))
(array([[ 659.05250673]]), array([ 2.94148687, 7.93197352]))
(array([[ 651.38937928]]), array([ 3.37036181, 7.91512255]))
(array([[ 643.32456722]]), array([ 3.48845966, 7.90996925]))
(array([[ 634.23683109]]), array([ 3.67415638, 7.90165163]))
(array([[ 625.01691153]]), array([ 3.78469509, 7.89626907]))
(array([[ 616.28646562]]), array([ 4.10536558, 7.87794627]))
(array([[ 606.83190524]]), array([ 4.27070726, 7.86895525]))
(array([[ 597.72512094]]), array([ 4.44427912, 7.86086966]))
(array([[ 588.76235737]]), array([ 4.53285413, 7.85586771]))
(array([[ 578.99990728]]), array([ 4.55868017, 7.85428077]))
(array([[ 569.5996939]]), array([ 4.61594834, 7.85117142]))
(array([[ 560.27767587]]), array([ 4.78495773, 7.84129916]))
(array([[ 550.4133697]]), array([ 4.82713417, 7.83850416]))
(array([[ 540.72332877]]), array([ 4.89693311, 7.83386429]))
(array([[ 530.9398232]]), array([ 5.04056355, 7.82664215]))
(array([[ 520.75413333]]), array([ 5.08128747, 7.82412106]))
(array([[ 511.75483855]]), array([ 5.14050223, 7.82132126]))
(array([[ 502.25335828]]), array([ 5.17652585, 7.81955516]))
(array([[ 492.59935793]]), array([ 5.36377025, 7.80819625]))
(array([[ 483.07806801]]), array([ 5.55688524, 7.79666769]))
(array([[ 473.2797989]]), array([ 5.65781892, 7.79083901]))
(array([[ 463.75008479]]), array([ 5.73299055, 7.78655785]))
(array([[ 453.89148345]]), array([ 5.78972062, 7.78293183]))
(array([[ 444.35538095]]), array([ 5.86458787, 7.77849972]))
(array([[ 434.56304965]]), array([ 5.91815948, 7.77443604]))
(array([[ 424.82156688]]), array([ 5.98593806, 7.76997721]))
(array([[ 415.0441472]]), array([ 6.05104991, 7.76518086]))
(array([[ 405.25850099]]), array([ 6.12786796, 7.75989745]))
(array([[ 395.47026118]]), array([ 6.21236257, 7.75405017]))
(array([[ 385.63605027]]), array([ 6.2861066 , 7.74885399]))
(array([[ 375.81358306]]), array([ 6.34908624, 7.74451606]))
(array([[ 365.93738554]]), array([ 6.41681992, 7.73947785]))
(array([[ 356.07517752]]), array([ 6.46488853, 7.73598365]))
(array([[ 345.91533166]]), array([ 6.50633578, 7.73280564]))
(array([[ 337.63751429]]), array([ 6.75654974, 7.71693515]))
(array([[ 329.10072617]]), array([ 6.96480374, 7.70639118]))
(array([[ 319.42373764]]), array([ 7.13440725, 7.69444878]))
(array([[ 310.05225607]]), array([ 7.32891571, 7.68240667]))
(array([[ 301.09711258]]), array([ 7.45730193, 7.67404167]))
(array([[ 291.28392583]]), array([ 7.50681935, 7.67090522]))
(array([[ 281.66332293]]), array([ 7.59273355, 7.66618255]))
(array([[ 272.53633632]]), array([ 7.71455156, 7.65692776]))
(array([[ 262.55295894]]), array([ 7.75245989, 7.65390367]))
(array([[ 252.68630671]]), array([ 7.80495028, 7.64969015]))
(array([[ 242.63294252]]), array([ 7.85947953, 7.64494254]))
(array([[ 232.97493606]]), array([ 7.9991447 , 7.63430865]))
(array([[ 223.06390541]]), array([ 8.13204536, 7.62360791]))
(array([[ 213.35377841]]), array([ 8.20329697, 7.61820163]))
(array([[ 203.42371481]]), array([ 8.27161684, 7.61256348]))
(array([[ 193.57845611]]), array([ 8.32446033, 7.60879174]))
(array([[ 183.55466831]]), array([ 8.37082973, 7.60485314]))
(array([[ 174.95921919]]), array([ 8.55596094, 7.59123662]))
(array([[ 166.93788983]]), array([ 8.91235508, 7.57477632]))
(array([[ 158.79157892]]), array([ 8.96777694, 7.57181997]))
(array([[ 149.59584174]]), array([ 9.04008135, 7.56763402]))
(array([[ 141.66186564]]), array([ 9.23524726, 7.55737963]))
(array([[ 132.38700092]]), array([ 9.47042699, 7.54283851]))
(array([[ 123.37601619]]), array([ 9.58641478, 7.53568048]))
(array([[ 113.95105959]]), array([ 9.68669788, 7.52954149]))
(array([[ 104.94941505]]), array([ 9.92735758, 7.51462028]))
(array([[ 95.41147825]]), array([ 10.00809535, 7.50910318]))
(array([[ 86.18233559]]), array([ 10.11391462, 7.50237804]))
(array([[ 76.52762158]]), array([ 10.15244465, 7.49981917]))
(array([[ 67.33235649]]), array([ 10.28781492, 7.49064229]))
(array([[ 57.50212403]]), array([ 10.42332069, 7.48044637]))
(array([[ 48.17773418]]), array([ 10.59061895, 7.46873045]))
(array([[ 38.42155475]]), array([ 10.65567416, 7.46435312]))
(array([[ 29.35859308]]), array([ 10.76135121, 7.45728814]))
(array([[ 20.01769569]]), array([ 10.9578228 , 7.44505251]))
(array([[ 10.57554175]]), array([ 11.07056766, 7.43801355]))
(array([[ 0.75261369]]), array([ 11.15749668, 7.43169608]))
(array([[-9.09160254]]), array([ 11.21756991, 7.42765372]))
(array([[-18.95401102]]), array([ 11.29089161, 7.4226755 ]))
(array([[-29.28201216]]), array([ 11.30578072, 7.42163708]))
(array([[-39.30662423]]), array([ 11.47512712, 7.40879657]))
(array([[-48.14744033]]), array([ 11.53934956, 7.40424035]))
(array([[-57.12288463]]), array([ 11.47974521, 7.4074263 ]))
(array([[-66.3607443]]), array([ 11.53220303, 7.40352362]))
(array([[-75.86473823]]), array([ 11.56819964, 7.40081723]))
(array([[-81.83679535]]), array([ 9.32049824, 7.54773468]))
(array([[-90.09264074]]), array([ 8.44687851, 7.60416208]))
(array([[-99.57827991]]), array([ 8.15776066, 7.62475112]))
(array([[-108.40023049]]), array([ 8.28709296, 7.61220032]))
(array([[-111.88215476]]), array([ 5.68876225, 7.77786993]))
(array([[-118.58721231]]), array([ 2.89919807, 7.95713847]))
(array([[-117.66508819]]), array([ 3.75579752, 7.87566675]))
(array([[-87.29471285]]), array([-16.20086598, 9.20293223]))
(array([[-83.78706061]]), array([-16.83224553, 9.26577831]))
(array([[-29.44997701]]), array([-36.19991339, 10.78012035]))
(array([[ 26.65221721]]), array([-44.59756798, 11.58573625]))
(array([[ 128.95259464]]), array([-98.08782071, 15.82537157]))
(array([[ 187.12148264]]), array([-124.56442308, 17.70856952]))
(array([[ 180.39864553]]), array([-120.32461274, 17.44298158]))
(array([[ 174.32204878]]), array([-124.41638017, 17.71184052]))
(array([[ 202.03476686]]), array([-147.88582578, 19.90689706]))
(array([[ 272.62073807]]), array([-186.93913884, 23.69702497]))
(array([[ 269.06335917]]), array([-191.44697753, 23.9012247 ]))
(array([[ 267.93883214]]), array([-211.05987484, 25.36397139]))
(array([[ 266.89111012]]), array([-204.42374176, 24.94585737]))
(array([[ 268.52760497]]), array([-223.97700662, 26.65905616]))
(array([[ 689.45347504]]), array([-289.96191955, 31.91948322]))
(array([[ 1404.37355629]]), array([-565.83333325, 53.669094 ]))
(array([[ 2323.33819235]]), array([-1036.14702034, 90.24054517]))
(array([[ 3307.51782033]]), array([-1889.7921164 , 157.62525954]))
(array([[ 3544.40731346]]), array([-2246.96882331, 187.00040087]))
(array([[ 3567.09771799]]), array([-2409.61451364, 199.59206047]))
(array([[ 3563.97367422]]), array([-2435.87282327, 201.54905244]))
(array([[ 3556.4332833]]), array([-2431.05439234, 201.23049821]))
(array([[ 3680.86154916]]), array([-2423.0943175 , 200.92586218]))
(array([[ 5687.8685826]]), array([-2189.67171478, 185.77958727]))
(array([[ 5826.67972432]]), array([-2065.12057334, 176.62879923]))
(array([[ 6614.64148167]]), array([-2069.29447556, 180.46221815]))
(array([[ 8034.68240485]]), array([-1882.5055244 , 164.19649087]))
(array([[ 10932.56146844]]), array([-3127.85507894, 261.86439891]))
(array([[ 11732.17233198]]), array([-4485.26750946, 367.12740118]))
(array([[ 11722.96033523]]), array([-4497.60752806, 368.06372715]))
(array([[ 11938.02542868]]), array([-4207.04344177, 347.94612662]))
(array([[ 13715.85745274]]), array([-4287.89152527, 361.52470494]))
(array([[ 13970.58952156]]), array([-4336.62154102, 362.8339204 ]))
(array([[ 15642.79881082]]), array([-4801.7411108 , 397.61365589]))
(array([[ 70508.63697127]]), array([-15640.45103455, 1374.73602715]))
(array([[ 147617.84505464]]), array([-82625.66300011, 6834.40121736]))
(array([[ 277044.73205632]]), array([-161980.28527832, 13841.37705427]))
(array([[ 299211.75589685]]), array([-193849.2277832 , 16342.47552109]))
(array([[ 394755.32945883]]), array([-243031.55744934, 20978.80837762]))
(array([[ 396305.54343813]]), array([-240797.37425995, 20860.27446252]))
(array([[ 433856.02060177]]), array([-262808.64431763, 23273.16009551]))
(array([[ 562910.16570153]]), array([-319371.68752289, 29525.15380931]))
(array([[ 604634.54454479]]), array([-341072.18295288, 32128.41426427]))
(array([[ 604785.65939884]]), array([-341825.81078911, 32217.32430971]))
(array([[ 655065.0655221]]), array([-377017.47119141, 35939.40427583]))
(array([[ 914320.1105354]]), array([-557692.69532013, 54199.39376681]))
(array([[ 921513.7134526]]), array([-544037.72674561, 52929.79643311]))
(array([[ 1086919.88239315]]), array([-568759.62106323, 55989.31088007]))
(array([[ 1090422.6594606]]), array([-570671.43305206, 56167.79889825]))
(array([[ 1295718.84297097]]), array([-633724.75268555, 62974.06821187]))
(array([[ 1317719.95262056]]), array([-641603.78643799, 63805.36361906]))
(array([[ 1317783.97634149]]), array([-642955.07075119, 63932.67639801]))
(array([[ 1317978.2123303]]), array([-643553.79808807, 63988.57819696]))
(array([[ 1690970.92295515]]), array([-699036.55419922, 70288.96598196]))
(array([[ 1708665.27685125]]), array([-684414.11611938, 68912.0253752 ]))
(array([[ 1716069.92554294]]), array([-699000.14810944, 70339.97539761]))
(array([[ 1729124.00550577]]), array([-727825.08906555, 73114.86306094]))
(array([[ 1928589.58966133]]), array([-763291.14865112, 77120.5305357 ]))
In [ ]:
Content source: jgosmann/plume
Similar notebooks: