Problem: H0 marginalization in the supernovae_Union2.f90
Posted: October 19 2015
Dear Prof. Antony Lewis,
I have some questions about the H0 marginalization in the supernovae_Union2.f90.
In the old version of cosmomc, such as the Oct. 2012, the chisq = AT-BT**2/sn_sumninv. This form is what we usually used formula: chi^2 = A-B^2/C, where the distance moduli mu=5 log10(dL) is independent of the nuisance parameter H0. However, in the latest version, the chisq = dot_product(diffs,matmul(this%sn_ninv,diffs)). I don't konw whether the formula: diffs=5 log10(dL)+25 is independent of H0. I don't know how to deal with it.
(1 ) diffs = 5 log10(dL) %independent of H0
chisq = dot_product(diffs,matmul(this%sn_ninv,diffs))
(2) diffs = 5 log10(dL) + \mu0 (H0) %depend on H0
chisq = dot_product(diffs,matmul(this%sn_ninv,diffs))
(3) diffs = 5 log10(dL) %independent of H0
AT = dot_product(diffs,matmul(sn_ninv,diffs))
BT = SUM(matmul(sn_ninv,diffs))
chisq = AT-BT**2/sn_sumninv
I want to know which is right for the above three formula when we use the Union2.1 data.
Thanks very much for your help.
Best wishes,
mingjian
I have some questions about the H0 marginalization in the supernovae_Union2.f90.
In the old version of cosmomc, such as the Oct. 2012, the chisq = AT-BT**2/sn_sumninv. This form is what we usually used formula: chi^2 = A-B^2/C, where the distance moduli mu=5 log10(dL) is independent of the nuisance parameter H0. However, in the latest version, the chisq = dot_product(diffs,matmul(this%sn_ninv,diffs)). I don't konw whether the formula: diffs=5 log10(dL)+25 is independent of H0. I don't know how to deal with it.
(1 ) diffs = 5 log10(dL) %independent of H0
chisq = dot_product(diffs,matmul(this%sn_ninv,diffs))
(2) diffs = 5 log10(dL) + \mu0 (H0) %depend on H0
chisq = dot_product(diffs,matmul(this%sn_ninv,diffs))
(3) diffs = 5 log10(dL) %independent of H0
AT = dot_product(diffs,matmul(sn_ninv,diffs))
BT = SUM(matmul(sn_ninv,diffs))
chisq = AT-BT**2/sn_sumninv
I want to know which is right for the above three formula when we use the Union2.1 data.
Thanks very much for your help.
Best wishes,
mingjian