Last chance! 50% off unlimited learning
Sale ends in
data(grain.us)
cgrain = scale(grain.us, scale = FALSE) # Center the time series only
fit = vglm(cgrain ~ 1, rrar(Rank = c(4, 1)),
eps = 1e-3, step = 0.5, trace = TRUE, maxit = 40)
summary(fit)
Run the code above in your browser using DataLab