############## Esta función permite dibujar
############## el intervalo de confianza simultáneo
############## de la recta. Se llama confianza

############## Ejemplo: confianza(recta)

confianza<-
function(object, alpha = 0.05, plot.it = F, ...)
{
	f <- predict(object, se.fit = T)
	p <- length(coef(object))
	fit <- f$fit
	adjust <- (p * qf(1 - alpha, p, length(fit) - p))^0.5 * f$se.fit
	lower <- fit - adjust
	upper <- fit + adjust
	if(plot.it) {
		y <- fit + resid(object)
		plot(fit, y, type = "n")
		abline(0, 1, lty = 2)
		ord <- order(fit)
		lines(fit[ord], lower[ord])
		lines(fit[ord], upper[ord])
		list(lower = lower, upper = upper)
	}
	else list(lower = lower, upper = upper)
}