diff --git a/NEWS.md b/NEWS.md
index 6e8c9248f..8e7ababd0 100644
--- a/NEWS.md
+++ b/NEWS.md
@@ -72,6 +72,8 @@
    - Note: when interpolation is linear/midpoint for
      factors/characters it changes to locf with a warning (#759)
 
+   - Also note, that the default keep interpolation is `na`
+
 - Now you can specify the interpolation method per covariate in the model:
 
   - `linear(var1, var2)` says both `var1` and `var2` would use linear
diff --git a/R/rxsolve.R b/R/rxsolve.R
index 00db8f555..92fede70e 100644
--- a/R/rxsolve.R
+++ b/R/rxsolve.R
@@ -705,7 +705,7 @@ rxSolve <- function(object, params = NULL, events = NULL, inits = NULL,
                     cores,
                     covsInterpolation = c("locf", "linear", "nocb", "midpoint"),
                     naInterpolation = c("locf", "nocb"),
-                    keepInterpolation=c("locf", "nocb", "na"),
+                    keepInterpolation=c("na", "locf", "nocb"),
                     addCov = TRUE, sigma = NULL, sigmaDf = NULL,
                     sigmaLower = -Inf, sigmaUpper = Inf,
                     nCoresRV = 1L, sigmaIsChol = FALSE,