A method to predict the creep behavior of fiber-reinforced ceramic composites at high temperatures was suggested based on finite element modeling using constituent creep equations of fiber and matrix and showed good agreement with the experimental results. The effects of matrix creep behavior, fiber volume fraction, and residual stresses on the composite creep behavior were also investigated. The results showed that the primary behavior of composites was greatly affected by that of matrix but post-primary behavior was governed by fiber creep characteristics. The increase of fiber volume fraction from 15 vol% to 30 vol% caused the 50% and 40% decrease of steady-state creep rates and total creep strains at $1200^{circ}C$, 180MPa, respectively. Feasible compressive residual stresses in the matrix caused by different thermal expansion coefficients between the fiber and the matrix could significantly reduce total creep strains of the composite. The creep deformation mechanism in the fiber-reinforced ceramic composites could be explained by the stress transfer and redistribution in the fiber and matrix due to different creep characteristics of its constituents.