Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Concise model/info / model.result for relational models #931

Open
Jammy2211 opened this issue Feb 8, 2024 · 0 comments
Open

Concise model/info / model.result for relational models #931

Jammy2211 opened this issue Feb 8, 2024 · 0 comments
Assignees

Comments

@Jammy2211
Copy link
Collaborator

You get the idea from these issues:

#925
#926
#930

Heres a script which runs a relational model:

https://github.com/Jammy2211/autolens_workspace/blob/release/scripts/multi/modeling/features/wavelength_dependence.py

Here is what the model.info looks like:

Total Free Parameters = 23

model                                                                           Collection (N=23)
    0                                                                           Collection (N=23)
        galaxies                                                                Collection (N=23)
            lens                                                                Galaxy (N=15)
                bulge                                                           Sersic (N=8)
                    intensity                                                   SumPrior (N=2)
                        self                                                    MultiplePrior (N=1)
                mass                                                            Isothermal (N=5)
                shear                                                           ExternalShear (N=2)
            source                                                              Galaxy (N=8)
                bulge                                                           Sersic (N=8)
                    intensity                                                   SumPrior (N=2)
                        self                                                    MultiplePrior (N=1)
    1                                                                           Collection (N=23)
        galaxies                                                                Collection (N=23)
            lens                                                                Galaxy (N=15)
                bulge                                                           Sersic (N=8)
                    intensity                                                   SumPrior (N=2)
                        self                                                    MultiplePrior (N=1)
                mass                                                            Isothermal (N=5)
                shear                                                           ExternalShear (N=2)
            source                                                              Galaxy (N=8)
                bulge                                                           Sersic (N=8)
                    intensity                                                   SumPrior (N=2)
                        self                                                    MultiplePrior (N=1)

0
    galaxies
        lens
            redshift                                                            0.5
            bulge
                centre
                    centre_0                                                    GaussianPrior [0], mean = 0.0, sigma = 0.3
                    centre_1                                                    GaussianPrior [1], mean = 0.0, sigma = 0.3
                ell_comps
                    ell_comps_0                                                 GaussianPrior [2], mean = 0.0, sigma = 0.3
                    ell_comps_1                                                 GaussianPrior [3], mean = 0.0, sigma = 0.3
                intensity
                    lens_c                                                      UniformPrior [22], lower_limit = -10.0, upper_limit = 10.0
                    self
                        wavelength                                              464
                        lens_m                                                  UniformPrior [21], lower_limit = -0.1, upper_limit = 0.1
                effective_radius                                                UniformPrior [5], lower_limit = 0.0, upper_limit = 30.0
                sersic_index                                                    UniformPrior [6], lower_limit = 0.8, upper_limit = 5.0
            mass
                centre
                    centre_0                                                    GaussianPrior [7], mean = 0.0, sigma = 0.1
                    centre_1                                                    GaussianPrior [8], mean = 0.0, sigma = 0.1
                ell_comps
                    ell_comps_0                                                 GaussianPrior [9], mean = 0.0, sigma = 0.3
                    ell_comps_1                                                 GaussianPrior [10], mean = 0.0, sigma = 0.3
                einstein_radius                                                 UniformPrior [11], lower_limit = 0.0, upper_limit = 8.0
            shear
                gamma_1                                                         UniformPrior [12], lower_limit = -0.2, upper_limit = 0.2
                gamma_2                                                         UniformPrior [13], lower_limit = -0.2, upper_limit = 0.2
        source
            redshift                                                            1.0
            bulge
                centre
                    centre_0                                                    GaussianPrior [14], mean = 0.0, sigma = 0.3
                    centre_1                                                    GaussianPrior [15], mean = 0.0, sigma = 0.3
                ell_comps
                    ell_comps_0                                                 GaussianPrior [16], mean = 0.0, sigma = 0.3
                    ell_comps_1                                                 GaussianPrior [17], mean = 0.0, sigma = 0.3
                intensity
                    source_c                                                    UniformPrior [24], lower_limit = -10.0, upper_limit = 10.0
                    self
                        wavelength                                              464
                        source_m                                                UniformPrior [23], lower_limit = -0.1, upper_limit = 0.1
                effective_radius                                                UniformPrior [19], lower_limit = 0.0, upper_limit = 30.0
                sersic_index                                                    UniformPrior [20], lower_limit = 0.8, upper_limit = 5.0
1
    galaxies
        lens
            redshift                                                            0.5
            bulge
                centre
                    centre_0                                                    GaussianPrior [0], mean = 0.0, sigma = 0.3
                    centre_1                                                    GaussianPrior [1], mean = 0.0, sigma = 0.3
                ell_comps
                    ell_comps_0                                                 GaussianPrior [2], mean = 0.0, sigma = 0.3
                    ell_comps_1                                                 GaussianPrior [3], mean = 0.0, sigma = 0.3
                intensity
                    lens_c                                                      UniformPrior [22], lower_limit = -10.0, upper_limit = 10.0
                    self
                        wavelength                                              658
                        lens_m                                                  UniformPrior [21], lower_limit = -0.1, upper_limit = 0.1
                effective_radius                                                UniformPrior [5], lower_limit = 0.0, upper_limit = 30.0
                sersic_index                                                    UniformPrior [6], lower_limit = 0.8, upper_limit = 5.0
            mass
                centre
                    centre_0                                                    GaussianPrior [7], mean = 0.0, sigma = 0.1
                    centre_1                                                    GaussianPrior [8], mean = 0.0, sigma = 0.1
                ell_comps
                    ell_comps_0                                                 GaussianPrior [9], mean = 0.0, sigma = 0.3
                    ell_comps_1                                                 GaussianPrior [10], mean = 0.0, sigma = 0.3
                einstein_radius                                                 UniformPrior [11], lower_limit = 0.0, upper_limit = 8.0
            shear
                gamma_1                                                         UniformPrior [12], lower_limit = -0.2, upper_limit = 0.2
                gamma_2                                                         UniformPrior [13], lower_limit = -0.2, upper_limit = 0.2
        source
            redshift                                                            1.0
            bulge
                centre
                    centre_0                                                    GaussianPrior [14], mean = 0.0, sigma = 0.3
                    centre_1                                                    GaussianPrior [15], mean = 0.0, sigma = 0.3
                ell_comps
                    ell_comps_0                                                 GaussianPrior [16], mean = 0.0, sigma = 0.3
                    ell_comps_1                                                 GaussianPrior [17], mean = 0.0, sigma = 0.3
                intensity
                    source_c                                                    UniformPrior [24], lower_limit = -10.0, upper_limit = 10.0
                    self
                        wavelength                                              658
                        source_m                                                UniformPrior [23], lower_limit = -0.1, upper_limit = 0.1
                effective_radius                                                UniformPrior [19], lower_limit = 0.0, upper_limit = 30.0
                sersic_index                                                    UniformPrior [20], lower_limit = 0.8, upper_limit = 5.0

Here is what the model.results looks like:

Bayesian Evidence                                                               -1624746.95432675
Maximum Log Likelihood                                                          -1624734.15071202
Maximum Log Posterior                                                           -1624731.23908322

model                                                                           Collection (N=23)
    0                                                                           Collection (N=23)
        galaxies                                                                Collection (N=23)
            lens                                                                Galaxy (N=15)
                bulge                                                           Sersic (N=8)
                    intensity                                                   SumPrior (N=2)
                        self                                                    MultiplePrior (N=1)
                mass                                                            Isothermal (N=5)
                shear                                                           ExternalShear (N=2)
            source                                                              Galaxy (N=8)
                bulge                                                           Sersic (N=8)
                    intensity                                                   SumPrior (N=2)
                        self                                                    MultiplePrior (N=1)
    1                                                                           Collection (N=23)
        galaxies                                                                Collection (N=23)
            lens                                                                Galaxy (N=15)
                bulge                                                           Sersic (N=8)
                    intensity                                                   SumPrior (N=2)
                        self                                                    MultiplePrior (N=1)
                mass                                                            Isothermal (N=5)
                shear                                                           ExternalShear (N=2)
            source                                                              Galaxy (N=8)
                bulge                                                           Sersic (N=8)
                    intensity                                                   SumPrior (N=2)
                        self                                                    MultiplePrior (N=1)

Maximum Log Likelihood Model:

1
    galaxies
        lens
            bulge
                centre
                    centre_0                                                    0.005
                    centre_1                                                    -0.106
                ell_comps
                    ell_comps_0                                                 0.041
                    ell_comps_1                                                 0.410
                effective_radius                                                0.926
                sersic_index                                                    1.327
                intensity
                    self
                        lens_m                                                  0.005
                    lens_c                                                      1.804
            mass
                centre
                    centre_0                                                    -0.040
                    centre_1                                                    0.033
                ell_comps
                    ell_comps_0                                                 -0.009
                    ell_comps_1                                                 0.131
                einstein_radius                                                 5.786
            shear
                gamma_1                                                         0.101
                gamma_2                                                         -0.069
        source
            bulge
                centre
                    centre_0                                                    -0.259
                    centre_1                                                    -0.368
                ell_comps
                    ell_comps_0                                                 0.248
                    ell_comps_1                                                 -0.194
                effective_radius                                                0.697
                sersic_index                                                    2.363
                intensity
                    self
                        source_m                                                -0.030
                    source_c                                                    2.754

 WARNING: The samples have not converged enough to compute a PDF and model errors. 
 The model below over estimates errors. 



Summary (1.0 sigma limits):

1
    galaxies
        lens
            bulge
                centre
                    centre_0                                                    0.0054 (-0.6725, 0.5170)
                    centre_1                                                    -0.1059 (-0.6539, 0.8056)
                ell_comps
                    ell_comps_0                                                 0.0413 (-0.8034, 0.7618)
                    ell_comps_1                                                 0.4105 (-0.5712, 0.9085)
                effective_radius                                                0.9263 (0.0215, 21.6129)
                sersic_index                                                    1.3269 (0.8286, 4.5588)
                intensity
                    self
                        lens_m                                                  0.0 (-0.1, 0.1)
                    lens_c                                                      1.8042 (-9.8177, 9.7379)
            mass
                centre
                    centre_0                                                    -0.0396 (-0.1951, 0.1702)
                    centre_1                                                    0.0331 (-0.1350, 0.2387)
                ell_comps
                    ell_comps_0                                                 -0.0093 (-0.6904, 0.8435)
                    ell_comps_1                                                 0.1311 (-0.5822, 0.6198)
                einstein_radius                                                 5.7859 (0.1773, 7.9742)
            shear
                gamma_1                                                         0.1015 (-0.1945, 0.1886)
                gamma_2                                                         -0.0689 (-0.1960, 0.1896)
        source
            bulge
                centre
                    centre_0                                                    -0.2588 (-0.7563, 0.9348)
                    centre_1                                                    -0.3681 (-0.7200, 0.9345)
                ell_comps
                    ell_comps_0                                                 0.2477 (-0.5623, 0.5847)
                    ell_comps_1                                                 -0.1940 (-0.6462, 0.6469)
                effective_radius                                                0.6968 (0.0055, 18.2304)
                sersic_index                                                    2.3634 (0.8211, 4.9060)
                intensity
                    self
                        source_m                                                -0.0 (-0.1, 0.1)
                    source_c                                                    2.7541 (-9.1619, 9.3511)

instances

0
    galaxies
        lens
            redshift                                                            0.5
        source
            redshift                                                            1.0
1
    galaxies
        lens
            redshift                                                            0.5
        source
            redshift                                                            1.0
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants