# NOT RUN {
testData = createData(sampleSize = 400, family = gaussian())
fittedModel <- lm(observedResponse ~ Environment1 , data = testData)
# response that was used to fit the model
getObservedResponse(fittedModel)
# predictions of the model for these points
getFitted(fittedModel)
# extract simulations from the model as matrix
getSimulations(fittedModel, nsim = 2)
# extract simulations from the model for refit (often requires different structure)
x = getSimulations(fittedModel, nsim = 2, type = "refit")
getRefit(fittedModel, x[[1]])
getRefit(fittedModel, getObservedResponse(fittedModel))
# }
Run the code above in your browser using DataLab