A Monte Carlo method is used to calculate the profiles and the polarization of the Raman scattered O VI lines (λλ6827, 7088) in symbiotic stars, which are believed to be a binary system of a cool giant and a hot star with an emission nebula around it. A point-like isotropic UV radiation source is assumed and a simple spherical wind model is adopted for the kinematics of the scattering material from the cool giant.We first investigate the case where the incident line photons are described by a Gaussian profile having a width of 10 4 K. We subsequently investigate the effects of the extended ionized region and non-spherical wind models including a disk-type wind and a bipolar wind. The cases where the emission source is described by non-Gaussian profiles are briefly studied.Finally as an additional component for the kinematics of symbiotic stars the orbital motion of the hot component around the cool giant is included and the effect on the spectropolarimetry is investigated. In this case the polarization direction changes around the red part of the Raman-scattered emission lines, when the observer's line of sight is perpendicular to the orbital plane, and no such effect is seen when the line of sight lies in the orbital plane. Furthermore, complex peak structures are seen in the degree of polarization and the polarized flux, which have often been observed in several symbiotic systems including RR Tel.Brief observational consequences and predictions are discussed in relation to the present and future spectropolarimetry for symbiotic stars. It is concluded that spectropolarimetry may provide a powerful diagnostic of the physical conditions of symbiotic stars.