# data(fmdata1)
# s <- c(12.9, 1.5, 18.3, 5.2, 5.0, 2.3)
# w <- c(1.4, 0.3, 0.5, 1.6, 1.7, 1.7)
# par <- c(0,0, s, w)
# Estimated model by rav:
# fit1 <- rav(fmdata1, lev=c(3,3))
# Fitted model by original parameters:
# fit2 <- rav.indexes(param=par, lev=c(3,3), data=fmdata1)Run the code above in your browser using DataLab