Tensorflow 黑客的贝叶斯方法第1章使用Edward
第一章示例:从文本消息数据推断行为 我们想在一系列每日短信统计中找到一个切换点。 假设序列可以分成两段,具有两个泊松分布 使用pymc3:(基于theano) 我使用edward:(基于tensorflow) 抛出此异常:Tensorflow 黑客的贝叶斯方法第1章使用Edward,tensorflow,pymc3,Tensorflow,Pymc3,第一章示例:从文本消息数据推断行为 我们想在一系列每日短信统计中找到一个切换点。 假设序列可以分成两段,具有两个泊松分布 使用pymc3:(基于theano) 我使用edward:(基于tensorflow) 抛出此异常: InvalidArgumentError回溯(最近一次调用上次) c:\python35\lib\site packages\tensorflow\python\client\session.py in\u do\u call(self,fn,*args) 1138尝试:
InvalidArgumentError回溯(最近一次调用上次)
c:\python35\lib\site packages\tensorflow\python\client\session.py in\u do\u call(self,fn,*args)
1138尝试:
->1139返回fn(*args)
1140错误除外。操作错误为e:
c:\python35\lib\site packages\tensorflow\python\client\session.py in\u run\fn(会话、提要、获取列表、目标列表、选项、运行元数据)
1120进货目录、进货目录、目标目录、,
->1121状态,运行(元数据)
1122
c:\python35\lib\contextlib.py in\uuuuu exit\uuuuuu(self、type、value、traceback)
65尝试:
--->66下一个(self.gen)
67除停止迭代外:
c:\python35\lib\site packages\tensorflow\python\framework\errors\u impl.py in raise\u exception\u on\u not\u ok\u status()
465兼容as_文本(pywrap_tensorflow.TF_消息(状态)),
-->466 pywrap_tensorflow.TF_GetCode(状态))
467最后:
InvalidArgumentError:索引=5000不在[0,5000]中
[[Node:ScatterUpdate\u 26=ScatterUpdate[T=DT\u FLOAT,Tindices=DT\u INT32,\u class=[“lock:@Variable\u 77”],use\u locking=true,\u device=“/job:localhost/replica:0/task:0/cpu:0”](变量\u 77,迭代\u 17/read,cond\u 17/Merge\u 3)]]
在处理上述异常期间,发生了另一个异常:
InvalidArgumentError回溯(最后一次最近调用)
在()
30
31适用于范围内的t(推断n_iter):
--->32 info_dict=推断。更新()
33打印(信息记录)
34
更新中的c:\python35\lib\site packages\edward\expressions\monte\u carlo.py(self,feed\u dict)
135
136 sess=获取会话()
-->137,接受率=系统运行([self.train,self.n\u accept\u over\t],feed\u dict)
138 t=自运行(自增量)
139
c:\python35\lib\site packages\tensorflow\python\client\session.py in run(self、fetches、feed\u dict、options、run\u元数据)
787尝试:
788结果=self.\u运行(无、取数、输入、选项、,
-->789运行(元数据)
790如果运行\u元数据:
791 proto_data=tf_session.tf_GetBuffer(run_metadata_ptr)
c:\python35\lib\site packages\tensorflow\python\client\session.py in\u run(self、handle、fetches、feed\u dict、options、run\u元数据)
995如果最终_获取或最终_目标:
996 results=self.\u do\u run(句柄、最终目标、最终获取、,
-->997提要内容(字符串、选项、运行元数据)
998其他:
999结果=[]
c:\python35\lib\site packages\tensorflow\python\client\session.py in\u do\u run(self、handle、target\u list、fetch\u list、feed\u dict、options、run\u元数据)
1130如果句柄为“无”:
1131返回self.\u do\u call(\u run\u fn,self.\u session,feed\u dict,fetch\u list,
->1132目标\u列表、选项、运行\u元数据)
1133其他:
1134返回self.\u do.\u调用(\u prun\u fn,self.\u会话,句柄,提要\u dict,
c:\python35\lib\site packages\tensorflow\python\client\session.py in\u do\u call(self,fn,*args)
1150键错误除外:
1151通行证
->1152提升类型(e)(节点定义、操作、消息)
1153
1154定义扩展图(自):
InvalidArgumentError:索引=5000不在[0,5000]中
[[Node:ScatterUpdate\u 26=ScatterUpdate[T=DT\u FLOAT,Tindices=DT\u INT32,\u class=[“lock:@Variable\u 77”],use\u locking=true,\u device=“/job:localhost/replica:0/task:0/cpu:0”](变量\u 77,迭代\u 17/read,cond\u 17/Merge\u 3)]]
由op'ScatterUpdate_26'引起,定义为:
文件“c:\python35\lib\runpy.py”,第193行,在作为主模块的运行模块中
“\uuuuuuuuuuuuuuuuuuuuuuuuuuuuuuuuuuuuuuuuuuuuuuuuuuuuuuuuuuuuuuuuuu
文件“c:\python35\lib\runpy.py”,第85行,在运行代码中
exec(代码、运行\全局)
文件“c:\python35\lib\site packages\ipykernel\u launcher.py”,第16行,在
app.launch_new_instance()
文件“c:\python35\lib\site packages\traitlets\config\application.py”,第658行,在launch\u实例中
app.start()
文件“c:\python35\lib\site packages\ipykernel\kernelapp.py”,第477行,在开始处
ioloop.ioloop.instance().start()
文件“c:\python35\lib\site packages\zmq\eventloop\ioloop.py”,第177行,在开始处
super(ZMQIOLoop,self).start()
文件“c:\python35\lib\site packages\tornado\ioloop.py”,第888行,在开始处
handler_func(fd_obj,事件)
文件“c:\python35\lib\site packages\tornado\stack\u context.py”,第277行,空包装
返回fn(*args,**kwargs)
文件“c:\python35\lib\site packages\zmq\eventloop\zmqstream.py”,第440行,在事件句柄中
self.\u handle\u recv()
文件“c:\python35\lib\site packages\zmq\eventloop\zmqstream.py”,第472行,在\u handle\u recv中
self.\u运行\u回调(回调,消息)
文件“c:\python35\lib\site packages\zmq\eventloop\zmqstream.py”,第414行,在运行回调中
回调(*args,**kwargs)
文件“c:\python35\lib\site packages\tornado\stack\u context.py”,第277行,空包装
返回fn(*args,**kwargs)
dispatcher中第283行的文件“c:\python35\lib\site packages\ipykernel\kernelbase.py”
返回self.dispatch\u shell(流,消息)
文件“c:\python35\lib\site packages\ipykernel\kernelbase.py”,第235行,在dispatch_shell中
处理程序(流、标识、消息)
文件“c:\python35\lib\site packages\ipykernel\kernelbase.py”,第399行,在执行请求中
用户\u表达式,允许\u stdin)
文件“c:\python35\lib\site packages\ipykernel\ipkernel.py”,第196行,在do_execute中
res=shell.run\u单元格(代码,store\u history=store\u history,silent=silent)
文件“c:\python35\lib\site packages\ipykernel\zmqshell.py”,第533行,在运行单元中
返回超级(ZMQInteractiveSh)
import pymc3 as pm
import theano.tensor as tt
with pm.Model() as model:
alpha = 1.0/count_data.mean() # Recall count_data is the
# variable that holds our txt counts
lambda_1 = pm.Exponential("lambda_1", alpha)
lambda_2 = pm.Exponential("lambda_2", alpha)
tau = pm.DiscreteUniform("tau", lower=0, upper=n_count_data - 1)
with model:
idx = np.arange(n_count_data) # Index
lambda_ = pm.math.switch(tau >= idx, lambda_1, lambda_2)
observation = pm.Poisson("obs", lambda_, observed=count_data)
step = pm.Metropolis()
trace = pm.sample(10000, tune=5000,step=step)
import edward as ed
import tensorflow as tf
from edward.models import Exponential,Uniform,Poisson,Empirical
from edward.models import PointMass
alpha_f = 1.0/count_data.mean()
alpha = tf.Variable(alpha_f, name="alpha", dtype=tf.float32)
# init
lambda_1 = Exponential(alpha)
lambda_2 = Exponential(alpha)
tau = Uniform(low=0.0,high=float(n_count_data - 1))
idx = np.arange(n_count_data)
lambda_ = tf.where(tau>=idx,tf.ones(shape=[n_count_data],dtype=tf.float32)*lambda_1,tf.ones(shape=[n_count_data],dtype=tf.float32)*lambda_2)
z = Poisson(lambda_,value=count_data)
T = 5000
qlambda_1 = Empirical(params=tf.Variable(tf.zeros([T])))
qlambda_2 = Empirical(params=tf.Variable(tf.zeros([T])))
qtau = Empirical(params=tf.Variable(tf.zeros([T])))
qz = Empirical(params=tf.Variable(tf.zeros([T,n_count_data])))
inference = ed.HMC({z:qz,lambda_1:qlambda_1,lambda_2:qlambda_2,tau:qtau})
inference.run()
for t in range(inference.n_iter):
info_dict = inference.update()
print(info_dict)
InvalidArgumentError Traceback (most recent call last)
c:\python35\lib\site-packages\tensorflow\python\client\session.py in _do_call(self, fn, *args)
1138 try:
-> 1139 return fn(*args)
1140 except errors.OpError as e:
c:\python35\lib\site-packages\tensorflow\python\client\session.py in _run_fn(session, feed_dict, fetch_list, target_list, options, run_metadata)
1120 feed_dict, fetch_list, target_list,
-> 1121 status, run_metadata)
1122
c:\python35\lib\contextlib.py in __exit__(self, type, value, traceback)
65 try:
---> 66 next(self.gen)
67 except StopIteration:
c:\python35\lib\site-packages\tensorflow\python\framework\errors_impl.py in raise_exception_on_not_ok_status()
465 compat.as_text(pywrap_tensorflow.TF_Message(status)),
--> 466 pywrap_tensorflow.TF_GetCode(status))
467 finally:
InvalidArgumentError: indices = 5000 is not in [0, 5000)
[[Node: ScatterUpdate_26 = ScatterUpdate[T=DT_FLOAT, Tindices=DT_INT32, _class=["loc:@Variable_77"], use_locking=true, _device="/job:localhost/replica:0/task:0/cpu:0"](Variable_77, iteration_17/read, cond_17/Merge_3)]]
During handling of the above exception, another exception occurred:
InvalidArgumentError Traceback (most recent call last)
<ipython-input-107-2c42d8ad11c6> in <module>()
30
31 for t in range(inference.n_iter):
---> 32 info_dict = inference.update()
33 print(info_dict)
34
c:\python35\lib\site-packages\edward\inferences\monte_carlo.py in update(self, feed_dict)
135
136 sess = get_session()
--> 137 _, accept_rate = sess.run([self.train, self.n_accept_over_t], feed_dict)
138 t = sess.run(self.increment_t)
139
c:\python35\lib\site-packages\tensorflow\python\client\session.py in run(self, fetches, feed_dict, options, run_metadata)
787 try:
788 result = self._run(None, fetches, feed_dict, options_ptr,
--> 789 run_metadata_ptr)
790 if run_metadata:
791 proto_data = tf_session.TF_GetBuffer(run_metadata_ptr)
c:\python35\lib\site-packages\tensorflow\python\client\session.py in _run(self, handle, fetches, feed_dict, options, run_metadata)
995 if final_fetches or final_targets:
996 results = self._do_run(handle, final_targets, final_fetches,
--> 997 feed_dict_string, options, run_metadata)
998 else:
999 results = []
c:\python35\lib\site-packages\tensorflow\python\client\session.py in _do_run(self, handle, target_list, fetch_list, feed_dict, options, run_metadata)
1130 if handle is None:
1131 return self._do_call(_run_fn, self._session, feed_dict, fetch_list,
-> 1132 target_list, options, run_metadata)
1133 else:
1134 return self._do_call(_prun_fn, self._session, handle, feed_dict,
c:\python35\lib\site-packages\tensorflow\python\client\session.py in _do_call(self, fn, *args)
1150 except KeyError:
1151 pass
-> 1152 raise type(e)(node_def, op, message)
1153
1154 def _extend_graph(self):
InvalidArgumentError: indices = 5000 is not in [0, 5000)
[[Node: ScatterUpdate_26 = ScatterUpdate[T=DT_FLOAT, Tindices=DT_INT32, _class=["loc:@Variable_77"], use_locking=true, _device="/job:localhost/replica:0/task:0/cpu:0"](Variable_77, iteration_17/read, cond_17/Merge_3)]]
Caused by op 'ScatterUpdate_26', defined at:
File "c:\python35\lib\runpy.py", line 193, in _run_module_as_main
"__main__", mod_spec)
File "c:\python35\lib\runpy.py", line 85, in _run_code
exec(code, run_globals)
File "c:\python35\lib\site-packages\ipykernel_launcher.py", line 16, in <module>
app.launch_new_instance()
File "c:\python35\lib\site-packages\traitlets\config\application.py", line 658, in launch_instance
app.start()
File "c:\python35\lib\site-packages\ipykernel\kernelapp.py", line 477, in start
ioloop.IOLoop.instance().start()
File "c:\python35\lib\site-packages\zmq\eventloop\ioloop.py", line 177, in start
super(ZMQIOLoop, self).start()
File "c:\python35\lib\site-packages\tornado\ioloop.py", line 888, in start
handler_func(fd_obj, events)
File "c:\python35\lib\site-packages\tornado\stack_context.py", line 277, in null_wrapper
return fn(*args, **kwargs)
File "c:\python35\lib\site-packages\zmq\eventloop\zmqstream.py", line 440, in _handle_events
self._handle_recv()
File "c:\python35\lib\site-packages\zmq\eventloop\zmqstream.py", line 472, in _handle_recv
self._run_callback(callback, msg)
File "c:\python35\lib\site-packages\zmq\eventloop\zmqstream.py", line 414, in _run_callback
callback(*args, **kwargs)
File "c:\python35\lib\site-packages\tornado\stack_context.py", line 277, in null_wrapper
return fn(*args, **kwargs)
File "c:\python35\lib\site-packages\ipykernel\kernelbase.py", line 283, in dispatcher
return self.dispatch_shell(stream, msg)
File "c:\python35\lib\site-packages\ipykernel\kernelbase.py", line 235, in dispatch_shell
handler(stream, idents, msg)
File "c:\python35\lib\site-packages\ipykernel\kernelbase.py", line 399, in execute_request
user_expressions, allow_stdin)
File "c:\python35\lib\site-packages\ipykernel\ipkernel.py", line 196, in do_execute
res = shell.run_cell(code, store_history=store_history, silent=silent)
File "c:\python35\lib\site-packages\ipykernel\zmqshell.py", line 533, in run_cell
return super(ZMQInteractiveShell, self).run_cell(*args, **kwargs)
File "c:\python35\lib\site-packages\IPython\core\interactiveshell.py", line 2718, in run_cell
interactivity=interactivity, compiler=compiler, result=result)
File "c:\python35\lib\site-packages\IPython\core\interactiveshell.py", line 2822, in run_ast_nodes
if self.run_code(code, result):
File "c:\python35\lib\site-packages\IPython\core\interactiveshell.py", line 2882, in run_code
exec(code_obj, self.user_global_ns, self.user_ns)
File "<ipython-input-107-2c42d8ad11c6>", line 27, in <module>
inference.run()
File "c:\python35\lib\site-packages\edward\inferences\inference.py", line 123, in run
self.initialize(*args, **kwargs)
File "c:\python35\lib\site-packages\edward\inferences\hmc.py", line 64, in initialize
return super(HMC, self).initialize(*args, **kwargs)
File "c:\python35\lib\site-packages\edward\inferences\monte_carlo.py", line 98, in initialize
self.train = self.build_update()
File "c:\python35\lib\site-packages\edward\inferences\hmc.py", line 116, in build_update
assign_ops.append(tf.scatter_update(variable, self.t, sample[z]))
File "c:\python35\lib\site-packages\tensorflow\python\ops\gen_state_ops.py", line 599, in scatter_update
name=name)
File "c:\python35\lib\site-packages\tensorflow\python\framework\op_def_library.py", line 767, in apply_op
op_def=op_def)
File "c:\python35\lib\site-packages\tensorflow\python\framework\ops.py", line 2506, in create_op
original_op=self._default_original_op, op_def=op_def)
File "c:\python35\lib\site-packages\tensorflow\python\framework\ops.py", line 1269, in __init__
self._traceback = _extract_stack()
InvalidArgumentError (see above for traceback): indices = 5000 is not in [0, 5000)
[[Node: ScatterUpdate_26 = ScatterUpdate[T=DT_FLOAT, Tindices=DT_INT32, _class=["loc:@Variable_77"], use_locking=true, _device="/job:localhost/replica:0/task:0/cpu:0"](Variable_77, iteration_17/read, cond_17/Merge_3)]]
--------------------------------------------------------------------------------