-ace <- function(x, phy, type = "continuous", method = "ML", CI = TRUE,
- model = if (type == "continuous") "BM" else "ER",
- scaled = TRUE, kappa = 1, corStruct = NULL, ip = 0.1)
+ace <-
+ function(x, phy, type = "continuous",
+ method = if (type == "continuous") "REML" else "ML",
+ CI = TRUE, model = if (type == "continuous") "BM" else "ER",
+ scaled = TRUE, kappa = 1, corStruct = NULL, ip = 0.1,
+ use.expm = FALSE)