data(indices)
t.fit <- fit.tmv(indices)
gauss.fit <- fit.gaussmv(indices)
## Compare the fitted Student-t and Gaussian density.
par(mfrow = c(1, 2))
## Once on the real scale...
plot(t.fit[1], type = "l")
lines(gauss.fit[1], col = "red")
## ...and once scaled to expectation = 0, variance = 1
plot(scale(t.fit)[1], type = "l")
lines(scale(gauss.fit)[1], col = "red")
Run the code above in your browser using DataLab