diff --git a/fsl/wrappers/__init__.py b/fsl/wrappers/__init__.py
index 44812bce93bdae7386270d9963054360ceb7a2f6..87f410cb9c5475ab4542e96c074d0c04327763e1 100644
--- a/fsl/wrappers/__init__.py
+++ b/fsl/wrappers/__init__.py
@@ -40,9 +40,12 @@ the :func:`.bet` function uses ``mask`` instead of ``m``).
 Two exceptions to the above are :class:`.fslmaths` and :class:`.fslstats`,
 which provide a more object-oriented interface::
 
-    from fsl.wrappers import fslmaths
+    from fsl.wrappers import fslmaths, fslstats
+
     fslmaths('image.nii').mas('mask.nii').bin().run('output.nii')
 
+    imean, imin, imax = fslstats('image.nii').k('mask.nii').m.R.run()
+
 
 Wrapper functions for commands which accept NIfTI image or numeric text files
 will for the most part accept either in-memory ``nibabel`` images/Numpy arrays