data("dataHetIV")
het <- hetErrorsIV(y~X1+X2+P|P|IIV(X1, X2),
data = dataHetIV)
# returns the fitted values
predict(het)
# using the data used for fitting also for predicting,
# correctly results in fitted values
all.equal(predict(het, dataHetIV), fitted(het)) # TRUE
Run the code above in your browser using DataLab