This study shows how a Gibbs sampling approach can be used for Bayesian inference of inbreeding depression. The method presented is mainly concerned with organisms that can be both selfed and outcrossed. Tests performed on simulated data with unequal variances and missing observations show that the method works well. Real data from the plant Scabiosa canescens is also analyzed.