#---------------------------------------------------------------------------------------------

# Instalacion de las librerias necesarias: SOLO es necesario hacerlo UNA vez

cosa = winDialog("yesno", "¿Quieres instalar la libreria fdth?")
if (cosa=="YES"){
  local({r <- getOption("repos"); r["CRAN"] <- "http://cran.r-project.org/"; options(repos=r)})
  utils:::install.packages("fdth",dependencies=TRUE)
}

cosa = winDialog("yesno", "¿Quieres instalar la libreria gmodels?")
if (cosa=="YES"){
  local({r <- getOption("repos"); r["CRAN"] <- "http://cran.r-project.org/"; options(repos=r)})
  utils:::install.packages("gmodels",dependencies=TRUE)
}

#---------------------------------------------------------------------------------------------

# Cargas las librerias
library(fdth)
library(gmodels)

# Lees los  Datos  
nombre = file.choose()
cosa <- read.csv2(nombre, header=TRUE, sep=";" ,quote="\"",
dec=".", fill=TRUE, comment.char="")

attach(cosa)

# Tablas de Frecuencias
d = fdt(Renta) 
d

# Histograma
X11()
hist(Politica, main="Histograma de tendencias politicas", 
breaks=3, col="red")

# Diagrama de tartas
table(Define)
X11()
pie(table(Define), col=rainbow(4))

# Diagrama de cajas
X11()
boxplot(Renta~Zona, col=c("violet","yellow"))

# Tablas de frecuencias bidimensionales
ftable(xtabs(~Politica+Cree))
CrossTable(xtabs(~Politica+Cree), format="SPSS")

# Medidas de resumen, cuartiles y medidas de dispersion
summary(Politica)
cat("La moda es: ",which.max(table(Politica))[[1]],"\n")

# Aqui no hay relacion
# Correlacion
cor(Amor,Sexo, use="complete.obs")

# Regresion
rela = lm(Amor~Sexo)
summary(rela)
X11()
plot(Sexo,Amor)
abline(coef=coef(rela))

# Aqui si hay relacion
# Correlacion
cor(Toros,Politica, use="complete.obs")

# Regresion
relb = lm(Toros~Politica)
summary(relb)
X11()
plot(Politica,Toros)
abline(coef=coef(relb))

X11()
plot(Politica,Toros)
abline(coef=coef(relb))
segments(Politica,fitted(relb),Politica,Toros,col="blue")

# Intervalo de confianza para la media y test de hipotesis
t.test(Renta, mu=50)