Bayes Factor functions/packages

Question

Sim 2024-7-31

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2141801-bayes-factor-functions-packages

编辑： Walter Roberson 2024-8-2

I am using the following two-sample tests for non-normal distributions:

and I would like to calculate the Bayes Factor as well.

I found the bayesFactor Version 1.0.0 (253 KB) by Bart Krekelberg. However, to the best of my understanfding, that package has a limited number of implemented tests:

One sample t-test (bf.ttest)
Two sample t-test (bf.ttest2)
N-Way Anova with fixed and random effects, including continuous co-variates (bf.anova)
Regression (bf.regression)
Pearson Correlation (bf.corr)
Binomial Test (bf.binom)
Experimental Design & Power Analysis (bf.designAnalysis)

and I am not sure if they can be used as additional analysis to the 4 initially listed ones that I am employing.

Does anyone know if there are other Matlab functions/packages to calculate the Bayes factor, in relation to the two-sample tests I am currently using (i.e. the chi2gof, the kstest2, the ranksum, and the kruskalwallis)?

7 个评论
显示 5更早的评论隐藏 5更早的评论

Umar 2024-8-1

编辑：Walter Roberson 2024-8-2

在 MATLAB Online 中打开

@Sim,

After doing some thorough research online, I was able to dig out some research materials and gathered all the information to share with you now. So, to address your query, you can implement Bayesian models using probabilistic programming libraries like Stan, or PyMC3, and then interface with these libraries from Matlab using appropriate tools like the MATLAB Stan interface. For more information regarding Stan,PyMC3 please refer to

https://mc-stan.org/docs/stan-users-guide/index.html

https://www.pymc.io/projects/docs/en/stable/learn.html

Now, addressing your query about, Does anyone know if there are other Matlab functions/packages to calculate the Bayes factor, in relation to the two-sample tests I am currently using (i.e. the chi2gof, the kstest2, the ranksum, and the kruskalwallis)?

So, in order to calculate the Bayes Factor using the MATLAB Stan interface in relation to two-sample tests, ensure that you have the necessary data for your two-sample tests. This data should be appropriately formatted and ready for analysis. Then, select an appropriate Bayesian model that represents the hypotheses you want to compare. In the context of two-sample tests, you can define prior distributions for the parameters of interest and specify the likelihood function based on the test statistics used (e.g., chi-square, Kolmogorov-Smirnov, Wilcoxon rank-sum, or Kruskal-Wallis). Afterwards, you have to write the Bayesian model in the Stan modeling language by referring to their documentation. Here is a simple example of a Bayesian model for a two-sample test using Stan:

% Example Bayesian model in Stan for two-sample test
model_code = '
data {
 int<lower=0> N1; // Sample size for group 1
 int<lower=0> N2; // Sample size for group 2
 real y1[N1]; // Data for group 1
 real y2[N2]; // Data for group 2
}
parameters {
 real mu1; // Mean for group 1
 real mu2; // Mean for group 2
 real<lower=0> sigma1; // Standard deviation for group 1
 real<lower=0> sigma2; // Standard deviation for group 2
}
model {
 y1 ~ normal(mu1, sigma1); // Likelihood for group 1
 y2 ~ normal(mu2, sigma2); // Likelihood for group 2
 // Prior distributions can be specified here
}
';

Compile the Stan model using the stan function in MATLAB and make sure to provide the necessary data inputs and specify the number of iterations for sampling. Once the Bayesian model has been sampled using your defined methods, you can calculate the Bayes Factor by comparing the marginal likelihoods of the competing hypotheses. Please bear in mind that the Bayes Factor will quantify the relative strength of evidence for one hypothesis over another based on the observed data.

If you still have any further questions, please let me know.

Sim 2024-8-2

thanks a lot to everyone @Umar @Walter Roberson @Jeff Miller for your comments!

Many thanks @Umar for your research and suggestion on the way to proceed... it looks a bit complicated at a first sight, I will try it.. :-)

Umar 2024-8-2

Hi @Sim,

No problem, glad to help and good luck 👍 with your project.

请先登录，再进行评论。

请先登录，再回答此问题。

Bayes Factor functions/packages

7 个评论
显示 5更早的评论隐藏 5更早的评论

回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

Bayes Factor functions/packages

7 个评论 显示 5更早的评论隐藏 5更早的评论

回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

7 个评论
显示 5更早的评论隐藏 5更早的评论