# leo los datos datos = read.table(file="datos.txt", header=TRUE, dec=",", sep=";") # permito usar dichos nombres directamente attach(datos) # expreso la relación funcional funcion = Consumo ~ Renta + Deuda + Hijos GUIME.LM <- function(funcion, datos) { # resultados de la regresión modelo = lm(funcion, data=datos) analisis = summary(modelo) # estudio de la normalidad residuos = modelo$resid desviacióntípica = sd(residuos) normalidad = ks.test(residuos, pnorm, 0, desviacióntípica) # heteroscedasticidad library(car) heteroscedasticidad = ncv.test(modelo) # autocorrelación library(lmtest) autocorrelacion = dwtest(modelo) # multicolinealidad multicolinealidad = vif(modelo) resultado = list(funcion, analisis, normalidad, heteroscedasticidad, autocorrelacion, multicolinealidad) names(resultado) = c("Modelo de regresión", "Estimación y validación del modelo", "Normalidad residuos", "Heteroscedasticidad residuos", "Autocorrelación residuos", "Multicolinealidad variables") # salida función sink("AnálisisR.txt") resultado } GUIME.LM(funcion, datos)