rm(list=ls())

# Ejercicio 1
#Operaciones en R
2+3
15-7
4*6 #*multiplicacion
13/5
13%/%5 #Cociente de la division entera
13%%5  #Resto de la division entera
1/3+1/5
sqrt(25) 
sqrt(26)
sin(pi)
sin(3.14)

#Errores
log(0)
2/0
3*
3/*4
sqrt(-4)

#Limpieza de la consola de comandos Ctrl+L


# Ejercicio 2  
#Variables y vectores
a = 2
b = 3
c = a + b
a = b * c
b = (c - a)^2
c = a * b
a;b;c

#Ejercicio 3
edades=c(22,21,18,17,21,17,19,25)
edades2=c(17,20,21)
edades=c(edades,edades2)


# Ejercicio 4
#setear directorio de trabajo en Session
setwd("C:/Curso ESTADISTICA/ClasesR")
scan("../datos/Tut02-Edades.csv")

errorDeLectura = scan(file = "EsteFicheroNoExiste.csv")

scan("../datos/Tut02-Edades.csv",sep=",")


# Ejercicio 5
muestra = c(29, 17, 63, 31, 55, 9, 92, 61, 10, 16, 63, 6, 61, 59, 66, 41, 68, 6, 99, 21, 87, 68, 52, 83, 66, 98, 45, 50,
            24, 100, 83, 37, 44, 4, 97, 67, 56, 74, 75, 71, 55, 22, 86, 22, 93, 65, 38, 84, 54, 83, 100, 71, 99, 19,
            63, 11, 11, 62, 91, 20, 79, 42, 59, 95, 70, 74, 8, 25, 45, 58, 57, 75, 81, 34, 70, 68, 39, 12, 14, 21)
#Para guardar este vector en un fichero usamos la funcion cat

cat(muestra,file="muestra.csv")
#Mirarlo

cat(muestra,file="muestra.csv",fill=1)
#Mirarlo y ver la diferencia



# Ejercicio 6


vector1 = c(8, 5, 19, 9, 17, 2, 28, 18, 3, 4, 19, 1)
# Para cargar el vector2 asegurate de seleccionar el directorio de trabajo.

vector2 = scan(file="../datos/Tut02-vector2.csv",sep=";") 
vector1 + vector2
vector1*vector2
vector1+5
2*vector1
vector1^2
log(vector1)
logVector1=log(vector1)
1:15
-7:13
1.5:7.5
1.5:7.4


# Ejercicio 7

(cuadrados = (1:10)^2 )

cuadrados2 = (2:11)^2 

cuadrados2 - cuadrados
length(cuadrados)
sum(cuadrados)

# Ejercicio 8

sum((1:10)^2)
# No confundir con 
sum(1:10)^2


# Ejercicio 9
vectorEdades = scan(file="../datos/Tut02-Edades.csv")
( n = length(vectorEdades) )
( mediaEdades = sum(vectorEdades) / n )
( varEdades = sum( ( vectorEdades - mediaEdades )^2) / n )
( desvTipicaEdades = sqrt(varEdades))
( cuasiVarEdades = sum( ( vectorEdades - mediaEdades )^2) / (n - 1) )
( cuasiDesvTipicaEdades = sqrt(cuasiVarEdades))
sum( vectorEdades - mediaEdades )^2



# Ejercicio 10
rm(list=ls())
var3 = scan("../datos/Tut02-var3.csv")

# Ejercicio 11 

(minimo = min(var3) ) 
(maximo = max(var3) ) 
(rango = range(var3) )
( tablaFrecAbs = table(var3) )
length(tablaFrecAbs)
sum(tablaFrecAbs)
min(tablaFrecAbs)
max(tablaFrecAbs)
range(tablaFrecAbs)

#Ejercicio 12

class(var3) 
class(tablaFrecAbs) 
( n = length( var3 ) ) 
( tablaFrecRel = tablaFrecAbs / n )
sum(tablaFrecRel)

# Ejercicio 13

( tablaFrecAcu = cumsum(tablaFrecAbs)  )
( tablaFrecRelAcu = cumsum(tablaFrecRel) ) 
( tablaFrecRelAcu = (tablaFrecAcu / n) ) 

# Ejercicio 14

barplot(tablaFrecAbs, col = heat.colors(15)) 

boxplot(var3) 

boxplot(var3, horizontal=TRUE) 

boxplot(var3, outline = FALSE) 

boxplot(var3, horizontal = TRUE, outline = FALSE)

box(lwd=2)
# y ver que el resultado del segundo no es lo que se esperaba.

# Ejercicio 15

boxplot(var3, horizontal = TRUE, outline = FALSE,lwd=2)
box(lwd=2)
boxplot(var3, outline = FALSE, horizontal = TRUE)


# Ejercicio 16

vectorEdades = scan(file="../datos/Tut02-Edades.csv")

(mediaEdades = mean(vectorEdades))

( mediaVar3 = mean(var3)) # Primera forma
( mediaVar3 = sum(var3)/length(var3)) # Segunda forma

# Ejercicio 17

( varMuestral = var(var3))
( desvTipMuestral = sd(var3))
sqrt(varMuestral)
n=length(vectorEdades)
( varPobl= ( (n-1) / n ) * varMuestral )

( desvTipPobl = sqrt( varPobl ) )

( mediana = median(var3))

# Ejercicio 18
summary(var3)
class(summary(var3))

( rangoIntCuart= IQR(var3))

( percentiles = quantile(var3, c(0.05, 0.15, 0.58, 0.75)) )



# Ejercicio 22

sample(1:50)
sample(1:6, size=1)
dado = sample(1:6, size=1000, replace = TRUE) 
table(dado)
class(dado)


# Ejercicio 23

(edades=c(22,21,18,19,17,21,18,20,17,18,17,22,20,19,18,19,18,22,20,10))
(edadesAlAzar=sample(edades,size=7,replace=FALSE))
(edadesDistintas=unique(edades))
(edadesDistintasOrdenadas=sort(unique(edades)))
(edadesDistintasOrdenadas=sort(unique(edades),decreasing=TRUE))


( Clave = sample(c(LETTERS, letters, 0:9), size=15, replace=TRUE) )
paste(Clave, collapse="")

(deTresEnTres=seq(from=1,to=100,by=3))
(nodos=seq(from=15,to=27,length.out=20))

valores=c(1,2,3,4,5)
rep(valores,times=3)

# Ejercicio 24
set.seed(2014)
dado = sample(1:6, size=1000, replace = TRUE)
table(dado)
dado[1]
dado[length(dado)]
tail(dado, 1)
tail(dado, 10)
head(dado)

# Ejercicio 25
masoIgualQue2 = dado[dado>=2]
head(masoIgualQue2)
pares = dado[dado%%2==0]
head(pares)


# Ejercicio 26
(7 < 5) | (2 < 4)
(7<5)&(2<4)
3==5
dado[dado==3]

# Ejercicio 27

set.seed(2014)
dado2 = sample(1:6, 1500, replace = TRUE)
sum(dado2 == 4)