In this paper, ordered categorical variables are used to compare between linear and nonlinear Bayesian structural equation models, Gibbs Sampling method is applied for estimation and model comparison. Statistical inferences, which involve estimation of parameters and their standard deviations, and residuals analyses for testing the posited model, are discussed. The proposed procedure is illustrated by a simulation data obtained from R program. Data results are obtained from WinBUGS program