mygls function

Our generalized least squares ls function

Internal function to adapt generalized least squares (gls) model with more details in output.

mygls(
  model,
  data = sys.frame(sys.parent()),
  correlation = NULL,
  weights = NULL,
  subset,
  method = c("REML", "ML"),
  na.action = na.fail,
  control = list(),
  verbose = FALSE
)

Arguments

• model: a formula.
• data: a data.frame containing variables called in model, correlation, weights, subset.
• correlation: a corCompSymm object. Default is NULL.
• weights: a varIdent object. Default is NULL.
• subset: an optional expression indicating which subset of the rows of data should be used in the fit. By default, all observations are included.
• method: a character string to choose the maximization method. Default is "REML".
• na.action: a function that indicates what should happen when the data contain NAs. Default is na.fail.
• control: a list of control values. Default is an empty list.
• verbose: an optional logical value. If TRUE information on the evolution of the iterative algorithm is printed. Default is FALSE.

Returns

a gls object