Fix all indentation errors

.lintr

@@ -6,6 +6,7 @@ linters: linters_with_defaults(
     object_length_linter = NULL,
     indentation_linter(
       indent = 2L,
+      hanging_indent_style = "never",
       assignment_as_infix = FALSE
     )
   )

R/Log_Likelihood.R

@@ -250,15 +250,12 @@ create.cov.upper.flex <- function(P, marg.var, marg.range, marg.smooth,
       j <- combs[iter, 2]
 
       smoothness.mat[i, j] <- (marg.smooth[i] + marg.smooth[j]) / 2
-      range.mat[i, j] <- 1 / sqrt(((1 / marg.range[i])^2 +
-        (1 / marg.range[j])^2) / 2)
+      range.mat[i, j] <- 1 / sqrt(((1 / marg.range[i])^2 + (1 / marg.range[j])^2) / 2)
 
       s1 <- sqrt(marg.var[i] * marg.var[j])
-      s2 <- ((1 / marg.range[i])^marg.smooth[i] *
-        (1 / marg.range[j])^marg.smooth[j]) /
+      s2 <- ((1 / marg.range[i])^marg.smooth[i] * (1 / marg.range[j])^marg.smooth[j]) /
         ((1 / range.mat[i, j])^(2 * smoothness.mat[i, j]))
-      s3 <- gamma(smoothness.mat[i, j]) / (sqrt(gamma(marg.smooth[i])) *
-        sqrt(gamma(marg.smooth[j])))
+      s3 <- gamma(smoothness.mat[i, j]) / (sqrt(gamma(marg.smooth[i])) * sqrt(gamma(marg.smooth[j])))
       s4 <- R.corr[iter]
       sig2.mat[i, j] <- s1 * s2 * s3 * s4
     }

R/PrestoGP_CreateU_Multivariate.R

@@ -112,8 +112,7 @@ sparseNN <- function(ordered_locs, n_neighbors,
   for (row in 1:n_neighbors) {
     # for the locations from 1 to n_neighbors, use the entire locs list to find the neighbors
     nn <- knn_indices(
-      ordered_locs[1:
-      (n_neighbors + 1), , drop = FALSE][-row, ,
+      ordered_locs[1:(n_neighbors + 1), , drop = FALSE][-row, ,
         drop = FALSE
       ],
       ordered_locs[row, , drop = FALSE], n_neighbors,
@@ -365,8 +364,7 @@ createUMultivariate <- function(vec.approx, params, cov_func = NULL) {
     ajj <- 1 / rangep[ondx[2]]
     aij <- sqrt((aii^2 + ajj^2) / 2)
     K1 <- rho.mat[ondx[1], ondx[2]] * sqrt(sig2[ondx[1]]) * sqrt(sig2[ondx[2]]) *
-      aii^vii * ajj^vjj * gamma(vij) / (aij^(2 * vij) * sqrt(gamma(vii) *
-        gamma(vjj))) *
+      aii^vii * ajj^vjj * gamma(vij) / (aij^(2 * vij) * sqrt(gamma(vii) * gamma(vjj))) *
       cov_func(dist_func(olocs[1, , drop = FALSE], olocs[2, , drop = FALSE], ),
         smoothness = vij, alpha = aij
       )
@@ -421,8 +419,7 @@ createUMultivariate <- function(vec.approx, params, cov_func = NULL) {
     ajj <- 1 / rangep[ondx[2]]
     aij <- sqrt((aii^2 + ajj^2) / 2)
     K1 <- rho.mat[ondx[1], ondx[2]] * sqrt(sig2[ondx[1]]) * sqrt(sig2[ondx[2]]) *
-      aii^vii * ajj^vjj * gamma(vij) / (aij^(2 * vij) * sqrt(gamma(vii) *
-        gamma(vjj))) *
+      aii^vii * ajj^vjj * gamma(vij) / (aij^(2 * vij) * sqrt(gamma(vii) * gamma(vjj))) *
       cov_func(dist_func(olocs[1, , drop = FALSE], olocs[2, , drop = FALSE], ),
         smoothness = vij, alpha = aij
       )
@@ -461,8 +458,7 @@ createUMultivariate <- function(vec.approx, params, cov_func = NULL) {
         # positive definite. See equation (9) in Apanasovich (2011).
         K1[j] <- rho.mat[ondx[i], ondx[cur.q[j]]] *
           sqrt(sig2[ondx[i]]) * sqrt(sig2[ondx[cur.q[j]]]) *
-          aii^vii * ajj^vjj * gamma(vij) / (aij^(2 * vij) *
-            sqrt(gamma(vii) * gamma(vjj))) *
+          aii^vii * ajj^vjj * gamma(vij) / (aij^(2 * vij) * sqrt(gamma(vii) * gamma(vjj))) *
           cov_func(
             dist_func(
               olocs[i, , drop = FALSE],

R/PrestoGP_Model.R

@@ -521,8 +521,7 @@ setMethod("scale_locs", "PrestoGPModel", function(model, locs) {
       for (j in 1:model@nscale) {
         locs.out[[i]][, model@scaling == j] <-
           locs[[i]][, model@scaling == j] /
-            model@covparams[model@param_sequence[2, 1] +
-              model@nscale * (i - 1) + j - 1]
+            model@covparams[model@param_sequence[2, 1] + model@nscale * (i - 1) + j - 1]
       }
     }
     return(locs.out)
@@ -535,25 +534,20 @@ setMethod("transform_covariance_parameters", "PrestoGPModel", function(model) {
     model@covparams <- c(
       exp(model@logparams[1:model@param_sequence[2, 2]]),
       gtools::inv.logit(
-        model@logparams[model@param_sequence[3, 1]:
-        model@param_sequence[3, 2]],
+        model@logparams[model@param_sequence[3, 1]:model@param_sequence[3, 2]],
         0, 2.5
       ),
-      exp(model@logparams[model@param_sequence[4, 1]:
-      model@param_sequence[4, 2]]),
-      tanh(model@logparams[model@param_sequence[5, 1]:
-      model@param_sequence[5, 2]])
+      exp(model@logparams[model@param_sequence[4, 1]:model@param_sequence[4, 2]]),
+      tanh(model@logparams[model@param_sequence[5, 1]:model@param_sequence[5, 2]])
     )
   } else {
     model@covparams <- c(
       exp(model@logparams[1:model@param_sequence[2, 2]]),
       gtools::inv.logit(
-        model@logparams[model@param_sequence[3, 1]:
-        model@param_sequence[3, 2]],
+        model@logparams[model@param_sequence[3, 1]:model@param_sequence[3, 2]],
         0, 2.5
       ),
-      exp(model@logparams[model@param_sequence[4, 1]:
-      model@param_sequence[4, 2]]), 1
+      exp(model@logparams[model@param_sequence[4, 1]:model@param_sequence[4, 2]]), 1
     )
   }
   invisible(model)

R/PrestoGP_Multivariate_Vecchia.R

@@ -45,78 +45,76 @@ setMethod("initialize", "MultivariateVecchiaModel", function(.Object, n_neighbor
 #' prediction <- prestogp_predict(model, X.test, locs.test)
 #' Vec.mean <- prediction[[1]]
 #' Vec.sds <- prediction[[2]]
-setMethod("prestogp_predict",
-  "MultivariateVecchiaModel",
+setMethod("prestogp_predict", "MultivariateVecchiaModel",
   function(model, X, locs, m = NULL, ordering.pred = c("obspred", "general"), pred.cond = c("independent", "general"), return.values = c("mean", "meanvar")) {
-  # validate parameters
-  ordering.pred <- match.arg(ordering.pred)
-  pred.cond <- match.arg(pred.cond)
-  return.values <- match.arg(return.values)
-  X <- check_input_pred(model, X, locs)
-  if (is.null(m)) { # m defaults to the value used for training
-    m <- model@n_neighbors
-  }
-  if (m < model@min_m) {
-    stop(paste("m must be at least ", model@min_m, sep = ""))
-  }
-  if (m >= nrow(model@X_train)) {
-    warning("Conditioning set size m chosen to be >=n. Changing to m=n-1")
-    m <- nrow(model@X_train) - 1
-  }
-
-  # Vecchia prediction at new locations
-  Vecchia.Pred <- predict(model@linear_model, newx = X, s = model@linear_model$lambda[model@lambda_1se_idx])
-  # Vecchia trend prediction at observed data
-  Vecchia.hat <- predict(model@linear_model, newx = model@X_train, s = model@linear_model$lambda[model@lambda_1se_idx])
+    # validate parameters
+    ordering.pred <- match.arg(ordering.pred)
+    pred.cond <- match.arg(pred.cond)
+    return.values <- match.arg(return.values)
+    X <- check_input_pred(model, X, locs)
+    if (is.null(m)) { # m defaults to the value used for training
+      m <- model@n_neighbors
+    }
+    if (m < model@min_m) {
+      stop(paste("m must be at least ", model@min_m, sep = ""))
+    }
+    if (m >= nrow(model@X_train)) {
+      warning("Conditioning set size m chosen to be >=n. Changing to m=n-1")
+      m <- nrow(model@X_train) - 1
+    }
 
-  # Test set prediction
-  res <- model@Y_train - Vecchia.hat
+    # Vecchia prediction at new locations
+    Vecchia.Pred <- predict(model@linear_model, newx = X, s = model@linear_model$lambda[model@lambda_1se_idx])
+    # Vecchia trend prediction at observed data
+    Vecchia.hat <- predict(model@linear_model, newx = model@X_train, s = model@linear_model$lambda[model@lambda_1se_idx])
 
-  locs.train.scaled <- scale_locs(model, model@locs_train)
-  locs.scaled <- scale_locs(model, locs)
-  vec.approx.test <- vecchia_Mspecify(locs.train.scaled, m,
-    locs.list.pred = locs.scaled,
-    ordering.pred = ordering.pred,
-    pred.cond = pred.cond
-  )
+    # Test set prediction
+    res <- model@Y_train - Vecchia.hat
 
-  ## carry out prediction
-  if (!model@apanasovich) {
-    params <- model@covparams
-    param.seq <- model@param_sequence
-    pred <- vecchia_Mprediction(res, vec.approx.test,
-      c(
-        params[1:param.seq[1, 2]],
-        rep(1, param.seq[2, 2] - param.seq[2, 1] + 1),
-        params[param.seq[3, 1]:
-        param.seq[5, 2]]
-      ),
-      return.values = return.values
-    )
-  } else {
-    pred <- vecchia_Mprediction(res, vec.approx.test, model@covparams,
-      return.values = return.values
+    locs.train.scaled <- scale_locs(model, model@locs_train)
+    locs.scaled <- scale_locs(model, locs)
+    vec.approx.test <- vecchia_Mspecify(locs.train.scaled, m,
+      locs.list.pred = locs.scaled,
+      ordering.pred = ordering.pred,
+      pred.cond = pred.cond
     )
-  }
 
-  # prediction function can return both mean and sds
-  # returns a list with elements mu.pred,mu.obs,var.pred,var.obs,V.ord
-  Vec.mean <- pred$mu.pred + Vecchia.Pred # residual + mean trend
-  if (return.values == "mean") {
-    return.list <- list(means = Vec.mean)
-  } else {
-    warning("Variance estimates do not include model fitting variance and are anticonservative. Use with caution.")
-    vec.var <- pred$var.pred
-    ndx.out <- NULL
-    for (i in seq_along(length(locs))) {
-      vec.sds[ndx.out == i] <- sqrt(vec.var[ndx.out == i] +
-        model@covparams[model@param_sequence[4, i]])
+    ## carry out prediction
+    if (!model@apanasovich) {
+      params <- model@covparams
+      param.seq <- model@param_sequence
+      pred <- vecchia_Mprediction(res, vec.approx.test,
+        c(
+          params[1:param.seq[1, 2]],
+          rep(1, param.seq[2, 2] - param.seq[2, 1] + 1),
+          params[param.seq[3, 1]:param.seq[5, 2]]
+        ),
+        return.values = return.values
+      )
+    } else {
+      pred <- vecchia_Mprediction(res, vec.approx.test, model@covparams,
+        return.values = return.values
+      )
     }
-    return.list <- list(means = Vec.mean, sds = vec.sds)
-  }
 
-  return(return.list)
-})
+    # prediction function can return both mean and sds
+    # returns a list with elements mu.pred,mu.obs,var.pred,var.obs,V.ord
+    Vec.mean <- pred$mu.pred + Vecchia.Pred # residual + mean trend
+    if (return.values == "mean") {
+      return.list <- list(means = Vec.mean)
+    } else {
+      warning("Variance estimates do not include model fitting variance and are anticonservative. Use with caution.")
+      vec.var <- pred$var.pred
+      ndx.out <- NULL
+      for (i in seq_along(length(locs))) {
+        vec.sds[ndx.out == i] <- sqrt(vec.var[ndx.out == i] + model@covparams[model@param_sequence[4, i]])
+      }
+      return.list <- list(means = Vec.mean, sds = vec.sds)
+    }
+
+    return(return.list)
+  }
+)
 
 setMethod("check_input", "MultivariateVecchiaModel", function(model, Y, X, locs) {
   if (!is.list(locs)) {
@@ -222,8 +220,7 @@ setMethod("specify", "MultivariateVecchiaModel", function(model) {
       for (j in 1:model@nscale) {
         olocs.scaled[model@vecchia_approx$ondx == i, model@scaling == j] <-
           olocs.scaled[model@vecchia_approx$ondx == i, model@scaling == j] *
-            model@covparams[model@param_sequence[2, 1] +
-              model@nscale * (i - 1) + j - 1]
+            model@covparams[model@param_sequence[2, 1] + model@nscale * (i - 1) + j - 1]
       }
     }
     model@vecchia_approx$locsord <- olocs.scaled
@@ -294,12 +291,8 @@ setMethod("transform_data", "MultivariateVecchiaModel", function(model, Y, X) {
       vecchia.approx = vecchia.approx,
       c(
         params[1:param.seq[1, 2]],
-        rep(
-          1,
-          param.seq[2, 2] - param.seq[2, 1] + 1
-        ),
-        params[param.seq[3, 1]:
-        param.seq[5, 2]]
+        rep(1, param.seq[2, 2] - param.seq[2, 1] + 1),
+        params[param.seq[3, 1]:param.seq[5, 2]]
       )
     )
   } else {

R/PrestoGP_Util_Functions.R

@@ -232,8 +232,7 @@ revMat <- function(mat) {
 }
 
 #' @export
-vecchia_Mprediction <- function(z, vecchia.approx, covparms, var.exact = NULL,
-                                return.values = "mean") {
+vecchia_Mprediction <- function(z, vecchia.approx, covparms, var.exact = NULL, return.values = "mean") {
   GPvecchia:::removeNAs()
   U.obj <- createUMultivariate(vecchia.approx, covparms)
   V.ord <- U2V(U.obj)