sample error with pymc version 5.17.0 #134
Comments
|
Hi, I am happy you like the software! Please update to the latest version 2.0.3 that error was resolved with that update. Thank you! |
|
Thanks for your reply, I updated beat to version 2.0.3 but it still reports the same error，I uploaded my data and yaml file and don't know where the error was |
|
Whats your pymc version? |
|
|
|
Please try installing an earlier version 5.16.2? There were major changes in pymc with 5.17.0 (just released this week) maybe that broke things? |
|
It's starting to work properly Thank you very much for your guidance |
|
Great thanks for reporting! Will need to pin the version of pymc then in the next release. |
I used beat to complete the data processing of column 3 and got the same result as the case. However, I used my own data and the following error occurred. I only had insar data, so I made a new yaml file，How should I deal with this error(/opt/beat) root@ubuntu:~/Desktop# beat sample mydz models - INFO Analysing problem ... models - INFO Using RectangularSource for 1 sources for event --- !pf.Event geodetic - INFO Number of geodetic datasets: 1 config - INFO not solving for opening_fraction, got fixed at 0.0 |
|
Please see the FAQ, yours is point 5: https://pyrocko.org/beat/docs/current/faq.html |
|
Thank you, Professor. I have completed List 3 and List 4 according to the instructions and the issue on github. beat is very attractive to me and I think it is a great software. I want to know more about it |
Dear developers, thank you very much for providing such excellent software. I found that I could not solve this error when I processed column 3, but I handled it normally on another computer without any error. My current beat version is 2.0.1, how should I deal with this error？
(/opt/beat) root@ubuntu:~/Desktop# beat sample Laquila
WARNING (pytensor.tensor.blas): Using NumPy C-API based implementation for BLAS functions.
To enable 'bem' mode packages 'pygmsh' and 'cutde' need to be installed.
config - INFO All hierarchicals ok!
config - INFO All hyperparameters ok!
config - INFO All priors ok!
models - INFO ... Initialising Geometry Optimizer ...
models - INFO Analysing problem ...
models - INFO ---------------------
models - INFO Using RectangularSource for 1 sources for event --- !pf.Event
lat: 42.29
lon: 13.35
time: '2009-04-06 01:32:49.190000057'
depth: 12000.0
name: '200904060132A'
magnitude: 6.343080192483292
region: 'CENTRAL ITALY'
catalog: 'gCMT'
moment_tensor: !pf.MomentTensor
mnn: 1.43e+18
mee: 1.87e+18
mdd: -3.3e+18
mne: 1.77e+18
mnd: -1.43e+18
med: 2.6900000000000003e+17
strike1: 120.23408298515041
dip1: 54.240869089580485
rake1: -112.81739742081386
strike2: 335.98575923255856
dip2: 41.58440373860804
rake2: -61.69749587601104
moment: 3.6696131948749036e+18
magnitude: 6.343080192483292
duration: 7.0
geodetic - INFO Number of geodetic datasets: 2
geodetic - INFO Initialising corrections ...
heart - INFO Setting up Ramps correction for Laquila_dscxn
heart - INFO Masking data for Ramps estimation!
heart - INFO Not correcting Laquila_dscxn for Euler Pole
heart - INFO Not correcting Laquila_dscxn for Strain Rate
heart - INFO Setting up Ramps correction for Laquila_ascxn
heart - INFO Masking data for Ramps estimation!
heart - INFO Not correcting Laquila_ascxn for Euler Pole
heart - INFO Not correcting Laquila_ascxn for Strain Rate
geodetic - INFO Number of geodetic data points: 419
geodetic - INFO Initialising geometry geodetic composite ...
models - INFO ... Building model ...
/opt/beat/lib/python3.11/site-packages/pymc/model/core.py:1453: UserWarning: To disable default transform, please use default_transform=None instead of transform=None. Setting transform to None will not have any effect in future.
warnings.warn(
/opt/beat/lib/python3.11/site-packages/pymc/model/core.py:1453: UserWarning: To disable default transform, please use default_transform=None instead of transform=None. Setting transform to None will not have any effect in future.
warnings.warn(
/opt/beat/lib/python3.11/site-packages/pymc/model/core.py:1453: UserWarning: To disable default transform, please use default_transform=None instead of transform=None. Setting transform to None will not have any effect in future.
warnings.warn(
/opt/beat/lib/python3.11/site-packages/pymc/model/core.py:1453: UserWarning: To disable default transform, please use default_transform=None instead of transform=None. Setting transform to None will not have any effect in future.
warnings.warn(
/opt/beat/lib/python3.11/site-packages/pymc/model/core.py:1453: UserWarning: To disable default transform, please use default_transform=None instead of transform=None. Setting transform to None will not have any effect in future.
warnings.warn(
/opt/beat/lib/python3.11/site-packages/pymc/model/core.py:1453: UserWarning: To disable default transform, please use default_transform=None instead of transform=None. Setting transform to None will not have any effect in future.
warnings.warn(
/opt/beat/lib/python3.11/site-packages/pymc/model/core.py:1453: UserWarning: To disable default transform, please use default_transform=None instead of transform=None. Setting transform to None will not have any effect in future.
warnings.warn(
/opt/beat/lib/python3.11/site-packages/pymc/model/core.py:1453: UserWarning: To disable default transform, please use default_transform=None instead of transform=None. Setting transform to None will not have any effect in future.
warnings.warn(
/opt/beat/lib/python3.11/site-packages/pymc/model/core.py:1453: UserWarning: To disable default transform, please use default_transform=None instead of transform=None. Setting transform to None will not have any effect in future.
warnings.warn(
/opt/beat/lib/python3.11/site-packages/pymc/model/core.py:1453: UserWarning: To disable default transform, please use default_transform=None instead of transform=None. Setting transform to None will not have any effect in future.
warnings.warn(
/opt/beat/lib/python3.11/site-packages/pymc/model/core.py:1453: UserWarning: To disable default transform, please use default_transform=None instead of transform=None. Setting transform to None will not have any effect in future.
warnings.warn(
/opt/beat/lib/python3.11/site-packages/pymc/model/core.py:1453: UserWarning: To disable default transform, please use default_transform=None instead of transform=None. Setting transform to None will not have any effect in future.
warnings.warn(
/opt/beat/lib/python3.11/site-packages/pymc/model/core.py:1453: UserWarning: To disable default transform, please use default_transform=None instead of transform=None. Setting transform to None will not have any effect in future.
warnings.warn(
models - INFO Initialised hyperparameter h_SAR with size 2
models - INFO Initialized 2 hyperparameters in total!
geodetic - INFO Evaluating config for Ramps corrections for datasets...
/opt/beat/lib/python3.11/site-packages/pymc/model/core.py:1453: UserWarning: To disable default transform, please use default_transform=None instead of transform=None. Setting transform to None will not have any effect in future.
warnings.warn(
/opt/beat/lib/python3.11/site-packages/pymc/model/core.py:1453: UserWarning: To disable default transform, please use default_transform=None instead of transform=None. Setting transform to None will not have any effect in future.
warnings.warn(
/opt/beat/lib/python3.11/site-packages/pymc/model/core.py:1453: UserWarning: To disable default transform, please use default_transform=None instead of transform=None. Setting transform to None will not have any effect in future.
warnings.warn(
/opt/beat/lib/python3.11/site-packages/pymc/model/core.py:1453: UserWarning: To disable default transform, please use default_transform=None instead of transform=None. Setting transform to None will not have any effect in future.
warnings.warn(
/opt/beat/lib/python3.11/site-packages/pymc/model/core.py:1453: UserWarning: To disable default transform, please use default_transform=None instead of transform=None. Setting transform to None will not have any effect in future.
warnings.warn(
/opt/beat/lib/python3.11/site-packages/pymc/model/core.py:1453: UserWarning: To disable default transform, please use default_transform=None instead of transform=None. Setting transform to None will not have any effect in future.
warnings.warn(
geodetic - INFO Evaluating config for Euler Pole corrections for datasets...
geodetic - INFO No Euler Pole correction!
geodetic - INFO Evaluating config for Strain Rate corrections for datasets...
geodetic - INFO No Strain Rate correction!
geodetic - INFO Initialized 6 hierarchical parameters.
geodetic - INFO Geodetic optimization on:
depth, dip, east_shift, length, north_shift, rake, slip, strike, width
geodetic - INFO Retrieving geodetic data-covariances with structure "import" for Laquila_dscxn ...
geodetic - INFO Retrieving geodetic data-covariances with structure "import" for Laquila_ascxn ...
geodetic - INFO Applying corrections! ...
models - INFO Model building was successful!
models - INFO Using "bin" backend to store samples!
models - INFO ... Initiate Sequential Monte Carlo ...
n_chains=1000, tune_interval=10, n_jobs=4, proposal_distribution: MultivariateCauchy,
metropolis - INFO Creating initial population for 1000 chains ...
Sampling: [Laquila_ascxn_azimuth_ramp, Laquila_ascxn_offset, Laquila_ascxn_range_ramp, Laquila_dscxn_azimuth_ramp, Laquila_dscxn_offset, Laquila_dscxn_range_ramp, depth, dip, east_shift, h_SAR, length, north_shift, nucleation_x, nucleation_y, rake, slip, strike, time, width]
pymc.sampling.forward - INFO Sampling: [Laquila_ascxn_azimuth_ramp, Laquila_ascxn_offset, Laquila_ascxn_range_ramp, Laquila_dscxn_azimuth_ramp, Laquila_dscxn_offset, Laquila_dscxn_range_ramp, depth, dip, east_shift, h_SAR, length, north_shift, nucleation_x, nucleation_y, rake, slip, strike, time, width]
metropolis - INFO Compiling model graph ...
Traceback (most recent call last):
File "/opt/beat/bin/beat", line 8, in
sys.exit(main())
^^^^^^
File "/opt/beat/lib/python3.11/site-packages/beat/apps/beat.py", line 2442, in main
globals()"command_" + command
File "/opt/beat/lib/python3.11/site-packages/beat/apps/beat.py", line 1064, in command_sample
step = problem.init_sampler(hypers=options.hypers)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/opt/beat/lib/python3.11/site-packages/beat/models/problems.py", line 175, in init_sampler
step = sampler.SMC(
^^^^^^^^^^^^
File "/opt/beat/lib/python3.11/site-packages/beat/sampler/smc.py", line 96, in init
super(SMC, self).init(
File "/opt/beat/lib/python3.11/site-packages/beat/sampler/metropolis.py", line 125, in init
self.compile_model_graph(model)
File "/opt/beat/lib/python3.11/site-packages/beat/sampler/metropolis.py", line 178, in compile_model_graph
out_vars=[model.varlogp],
^^^^^^^^^^^^^
File "/opt/beat/lib/python3.11/site-packages/pymc/model/core.py", line 862, in varlogp
return self.logp(vars=self.free_RVs)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/opt/beat/lib/python3.11/site-packages/pymc/model/core.py", line 746, in logp
rv_logps = transformed_conditional_logp(
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/opt/beat/lib/python3.11/site-packages/pymc/logprob/basic.py", line 598, in transformed_conditional_logp
temp_logp_terms = conditional_logp(
^^^^^^^^^^^^^^^^^
File "/opt/beat/lib/python3.11/site-packages/pymc/logprob/basic.py", line 481, in conditional_logp
fgraph = construct_ir_fgraph(rv_values, ir_rewriter=ir_rewriter)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/opt/beat/lib/python3.11/site-packages/pymc/logprob/rewriting.py", line 243, in construct_ir_fgraph
replacements = tuple((rv, valued_rv(rv, value)) for rv, value in ir_rv_values.items())
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/opt/beat/lib/python3.11/site-packages/pymc/logprob/rewriting.py", line 243, in
replacements = tuple((rv, valued_rv(rv, value)) for rv, value in ir_rv_values.items())
^^^^^^^^^^^^^^^^^^^^
File "/opt/beat/lib/python3.11/site-packages/pytensor/graph/op.py", line 304, in call
compute_test_value(node)
File "/opt/beat/lib/python3.11/site-packages/pytensor/graph/op.py", line 125, in compute_test_value
thunk()
File "/opt/beat/lib/python3.11/site-packages/pytensor/graph/op.py", line 524, in rval
r = p(n, [x[0] for x in i], o)
^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/opt/beat/lib/python3.11/site-packages/pymc/logprob/abstract.py", line 239, in perform
raise NotImplementedError("ValuedVar should not be present in the final graph!")
NotImplementedError: ValuedVar should not be present in the final graph!
(/opt/beat) root@ubuntu:~/Desktop#
BEAT_log.txt
The text was updated successfully, but these errors were encountered: