We present our spectroscopic data around Hβ of the symbiotic star V1016 Cyg obtained with the Bohyunsan Echelle Spectrograph, in order to secure the broad emission feature at around λ 4850, which is formed through Raman scattering of He II λ 972. The total cross-section around Lyγ is given approximately by σ (λ) 1.7 × 10 −28 [λ Lyγ /(λ − λ Lyγ )] 2 cm 2 , with λ Lyγ being the line centre wavelength of Lyγ . We find a centre shift redward by an amount λ = +0.64 Å in the Raman scattered He II λ 4850. This redward centre shift is exactly analogous to the effect for the Raman scattered He II λ 6545 blueward of Hα discussed in our previous study. We compute the branching ratios of Raman scattering into the level 2s and levels 3s and 3d, which are subsequently incorporated in our Monte Carlo code. Using this code, we present the centre shift of the 4850 feature as a function of the neutral hydrogen column density in the scattering region. Assuming that He II λ 972 emission is characterized by a single Gaussian profile, the redward peak shift observed in the Raman scattered He II λ 4850 feature corresponds to the neutral column density N H I = 1.2 × 10 21 cm −2 . Assuming that the covering factor ∼ 0.1 of the scattering region with respect to the He II emission region and adopting a simple spherical stellar wind model, we may place an upper boundṀ 3.6 × 10 −7 M yr −1 for the massloss rate of the giant component of V1016 Cyg. Our estimate can be severely affected by the kinematics of the scattering and He II emission regions and the exact atomic physics, about which brief discussions are presented.