Question for the mess slip distribution result in Scenario II

[vic-chou]

Hi, when I complete the running Scenario II, I plot the slip distribution (the following figure I attached), but the slip distribution is mess and different from the slip distribution in tutorial in the website. Can you give some suggestion to solve this problem?

[hvasbath]

What regularization and initialization did you use? How many chains in the SMC sampler did you use? Looks like it is not converged. For example in the file you sent me for debugging the plot there was no regularization used ...
If you did not use least-squares-solution as the starting points of the Markov Chains, but used random init you will have to use MANY Markov Chains with a large number of steps. E.g.: n_chains = 8000 and n_steps=500. This will produce huge sampling results-please keep that in mind of your disk-space!

[vic-chou]

I attached the configure file, in this runing, I set the regularization to laplacian based on the tutorial. The regularization and initialization is baed on the tutorial (I don't know if my understanding of the tutorial is right or not). The chains in the SMC sampler is set to 1000. How can I use least-squares-solution as the starting points of the Markov Chains? I don't see the procedure in the tutorial.
--- !beat.BEATconfig
name: Laquila
date: '2009-04-06'
event: !pf.Event
lat: 42.29
lon: 13.35
time: 2009-04-06 01:32:49.190000
name: 200904060132A
depth: 12000.0
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
project_dir: /data4/gyxu/beat/test_beat/Laquila
problem_config: !beat.ProblemConfig
mode: ffo
mode_config: !beat.FFOConfig
regularization: laplacian
npatches: 132
initialization: random
source_type: RectangularSource
stf_type: HalfSinusoid
decimation_factors:
geodetic: 4
n_sources: 1
datatypes: [geodetic]
hyperparameters:
h_SAR: !beat.heart.Parameter
name: h_SAR
form: Uniform
lower: [0.0]
upper: [6.0]
testvalue: [3.0]
h_laplacian: !beat.heart.Parameter
name: h_laplacian
form: Uniform
lower: [-3.0]
upper: [2.0]
testvalue: [-0.5]
priors:
uparr: !beat.heart.Parameter
name: uparr
form: Uniform
lower: [-0.1]
upper: [2.0]
testvalue: [1.15]
uperp: !beat.heart.Parameter
name: uperp
form: Uniform
lower: [-1.0]
upper: [1.0]
testvalue: [0.5]
geodetic_config: !beat.GeodeticConfig
datadir: /data4/gyxu/beat/test_beat/laquila_kite_container-master
names: [Laquila_dscxn, Laquila_ascxn]
types: [SAR]
calc_data_cov: true
interpolation: multilinear
fit_plane: false
dataset_specific_residual_noise_estimation: true
gf_config: !beat.GeodeticLinearGFConfig
store_superdir: /data4/gyxu/beat/test_beat/GF/Laquila
reference_model_idx: 0
n_variations: [0, 1]
earth_model_name: ak135-f-continental.m
nworkers: 1
reference_sources:
- !beat.sources.RectangularSource
lat: 42.29
lon: 13.35
north_shift: 5320.626993768318
east_shift: 10834.953621207327
elevation: 0.0
depth: 2783.512543671438
time: 2009-04-06 01:32:49.190000
stf: !pf.HalfSinusoidSTF
duration: 0.0
anchor: -1.0
stf_mode: post
strike: 145.5721440593029
dip: 54.984442677677144
rake: -114.77908006170307
length: 13046.578107041581
width: 9770.686350202324
velocity: 3500.0
slip: 0.5359876408166948
patch_widths: [2.0]
patch_lengths: [2.0]
extension_widths: [0.6]
extension_lengths: [0.4]
sample_rate: 1.1574074074074073e-05
sampler_config: !beat.SamplerConfig
name: SMC
progressbar: true
buffer_size: 5000
parameters: !beat.SMCConfig
tune_interval: 50
check_bnd: true
rm_flag: true
n_jobs: 5
n_steps: 100
n_chains: 1000
coef_variation: 1.0
stage: 0
proposal_dist: MultivariateNormal
update_covariances: false
hyper_sampler_config: !beat.SamplerConfig
name: Metropolis
progressbar: true
buffer_size: 5000
parameters: !beat.MetropolisConfig
tune_interval: 50
proposal_dist: Normal
check_bnd: true
rm_flag: false
n_jobs: 1
n_steps: 25000
n_chains: 20
thin: 2
burn: 0.5

[hvasbath]

It is written here under Markov Chain initialization:
https://hvasbath.github.io/beat/examples/FFO_static.html#sample-the-solution-space

But you are right, I did not write explicitely to set the initialization flag (under mode_config) from "random" to "lsq". Please do so! Will update the tutorial accordingly. Thanks!

[vic-chou]

I set the initialization to "lsq", but the slip distribution is still mess and the synthetic displacements is bad. Is there a problem with my understanding of the tutorial? Can you give me some suggestion? I attached the slip distribution and configure file.

the config_ffo.yaml file
--- !beat.BEATconfig
name: Laquila
date: '2009-04-06'
event: !pf.Event
lat: 42.29
lon: 13.35
time: 2009-04-06 01:32:49.190000
name: 200904060132A
depth: 12000.0
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
project_dir: /data4/gyxu/beat/test_beat/Laquila
problem_config: !beat.ProblemConfig
mode: ffo
mode_config: !beat.FFOConfig
regularization: laplacian
npatches: 132
initialization: lsq
source_type: RectangularSource
stf_type: HalfSinusoid
decimation_factors:
geodetic: 4
n_sources: 1
datatypes: [geodetic]
hyperparameters:
h_SAR: !beat.heart.Parameter
name: h_SAR
form: Uniform
lower: [0.0]
upper: [7.0]
testvalue: [3.5]
h_laplacian: !beat.heart.Parameter
name: h_laplacian
form: Uniform
lower: [-3.0]
upper: [2.0]
testvalue: [-0.5]
priors:
uparr: !beat.heart.Parameter
name: uparr
form: Uniform
lower: [-0.1]
upper: [2.0]
testvalue: [1.15]
uperp: !beat.heart.Parameter
name: uperp
form: Uniform
lower: [-1.0]
upper: [1.0]
testvalue: [0.5]
geodetic_config: !beat.GeodeticConfig
datadir: /data4/gyxu/beat/test_beat/laquila_kite_container-master
names: [Laquila_dscxn, Laquila_ascxn]
types: [SAR]
calc_data_cov: true
interpolation: multilinear
fit_plane: false
dataset_specific_residual_noise_estimation: true
gf_config: !beat.GeodeticLinearGFConfig
store_superdir: /data4/gyxu/beat/test_beat/GF/Laquila
reference_model_idx: 0
n_variations: [0, 1]
earth_model_name: ak135-f-continental.m
nworkers: 1
reference_sources:
- !beat.sources.RectangularSource
lat: 42.29
lon: 13.35
north_shift: 5320.626993768318
east_shift: 10834.953621207327
elevation: 0.0
depth: 2783.512543671438
time: 2009-04-06 01:32:49.190000
stf: !pf.HalfSinusoidSTF
duration: 0.0
anchor: -1.0
stf_mode: post
strike: 145.5721440593029
dip: 54.984442677677144
rake: -114.77908006170307
length: 13046.578107041581
width: 9770.686350202324
velocity: 3500.0
slip: 0.5359876408166948
patch_widths: [2.0]
patch_lengths: [2.0]
extension_widths: [0.6]
extension_lengths: [0.4]
sample_rate: 1.1574074074074073e-05
sampler_config: !beat.SamplerConfig
name: SMC
progressbar: true
buffer_size: 5000
parameters: !beat.SMCConfig
tune_interval: 50
check_bnd: true
rm_flag: true
n_jobs: 10
n_steps: 100
n_chains: 1000
coef_variation: 1.0
stage: 0
proposal_dist: MultivariateNormal
update_covariances: false
hyper_sampler_config: !beat.SamplerConfig
name: Metropolis
progressbar: true
buffer_size: 5000
parameters: !beat.MetropolisConfig
tune_interval: 50
proposal_dist: Normal
check_bnd: true
rm_flag: false
n_jobs: 1
n_steps: 25000
n_chains: 20
thin: 2
burn: 0.5

[hvasbath]

Of course after setting it to lsq you will have to rerun the sampling after deleting all your old sampling results. As mentioned earlier you should set the n_chains parameters and n_steps parameters to higher values. For the tutorial I used 3000 n_chains and 300 nsteps. Did you also do the hyper parameter estimation including the laplacian? Your bounds there are rather narrow. Try with -5, 5.

Please also do not use any decimation factor in the ffo= set it to 1 and rerun your linear GF calculations!

[vic-chou]

Yes, I run the "beat sample Laquila --hypers --mode=ffo" to do the hyper parameter estimation including the laplacian. The bounds -3~2 is updated after running the "beat sample Laquila --hypers --mode=ffo". Should I set the bounds to -5,5 without running the command "beat sample Laquila --hypers --mode=ffo"? Does the parameters "n_chains" and "nsteps" refer to theparameters in "sampler_config" section not in "hyper_sampler_config" section?

[hvasbath]

Yes set the parameters by hand without running the hypers-otherwise they will be overwritten again.
Of course this refers to the sampler_config.

[vic-chou]

OK, now I can obtain the same slip distribution in the tutorial under your suggestion. Thank you!

[hvasbath]

Great!