Nothing
model.fit <-
function (y_obs, y_exp, plotVal = FALSE)
{
n <- length(y_obs)
sum_y_obs <- sum(y_obs)
resids <- y_obs-y_exp
resids_sq <- resids^2
resids_abs <- abs(resids)
resids_abs_sum <- sum(resids_abs)
resids_sq_sum <- sum(resids_sq)
pseudorsq <- (var(y_obs)-var(resids))/var(y_obs)
mape <- mean(resids_abs/y_obs)
globerr <- resids_abs_sum/sum_y_obs
if (plotVal == TRUE)
{
plot(c(0, max(y_obs)), c(0, max(y_obs)), type="l", main="Model fit", xlab = "Observed", ylab = "Expected")
points (y_obs, y_exp, col = "blue", pch = 20, cex = 1.5)
}
return (list(resids_sq_sum=resids_sq_sum, pseudorsq=pseudorsq, globerr=globerr, mape=mape))
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.