Modeling the role of oxygen vacancy on ferroelectric properties in thin films