#Practico 2

#Ejercicio 3

attach(tabla_de_datos_practico_2)

Estudiantes=tabla_de_datos_practico_2

#Estudiantes1 = Estudiantes [order(Estudiantes$`Curso de nivelación`), ]

Estudiantes

#Seleccionar variables y cambiar nombres en una columna

Estudiantes$`Primera prueba`[`Primera prueba`=="Suficiente"]="Suficiente1"

Estudiantes$`Primera prueba`[`Primera prueba`=="Insuficiente"]="Insuficiente1"

Estudiantes ##llamar tabla original y ver los cambios

#Reordenar columnas

names(Estudiantes)#ver nombre de las columnas

Estudiantes1 = Estudiantes [ , c(3,2,4,1)] #comando para decir el nuevo orden de columnas

names(Estudiantes1)

#Para quitar una columna. Se descarga el paquete dplyr.

#Renombre el vector resultante y puse el comando select (datos, -la columna que se quiere sacar)

library(dplyr)

Estudiantes2=select(Estudiantes1, -Estudiante) 

Estudiantes2

#Diagrama de Sankey con ggplot

library(ggplot2)

library(highcharter)

library(htmlwidgets)

attach(Estudiantes2)

Tabla1= hchart(data_to_sankey(Estudiantes2), "sankey", name="Diagrama de Sankey")

Tabla1

#PRACTICO 3

#Ingresar datos

alfa=0.05

n=30

pestMont1=10/30 # 0.33

pestMont2=15/30 # 0.5

pestMont3=14/30 #0.46

Z=qnorm (1-alfa/2) #1.95

Z

#error estandar

EE1= sqrt(pestMont1*(1-pestMont1)/30) #0.0860663

EE1         

EE2= sqrt(pestMont2*(1-pestMont2)/30) #0.09128709

EE3= sqrt(pestMont3*(1-pestMont3)/30) #0.09108401

EE2

EE3

#Intervalos de confianza

I1=c(pestMont1-EE1*Z,pestMont1+EE1*Z)

I2=c(pestMont2-EE2*Z,pestMont2+EE2*Z)

I3=c(pestMont3-EE3*Z,pestMont3+EE3*Z)

#Hacer grafico de barras con las proporciones y agregar linea de IdC

library(ggplot2)

Tipos=c("Monte1", "Monte2", "Monte3") 

Presencia=c(pestMont1,pestMont2,pestMont3)

datoseje=data.frame(Tipos,Presencia) #tabla hecha con los vectores de arriba

IDCizq=c(0.16, 0.32,0.28) #Intervalos izquierdos de cada monte

IDCder=c(0.50, 0.67,0.64) #Intervalos derechos de cada monte

grafmonte<- ggplot(datoseje, aes(x=Tipos, y=Presencia)) +

geom_bar(position="dodge",fill=c("#7CAE00","#0CB702","#00BE67"), stat = "identity") +

geom_errorbar(aes(ymin=IDCizq, ymax=IDCder), width = .1, 

position=position_dodge(.7)) + ylab("Probabilidad de presencia") + 

ggtitle("Proporciones de los montes y sus intervalos de confianza") #geom_bar es la barra y geom_errorbar es la linea de los IdC

grafmonte #llamar grafica

#Binom test 

binom.test(10, 30, p = 0.05) #Monte 1

binom.test(15, 30, p = 0.05) #Monte 2

binom.test(14, 30, p = 0.05) #Monte 3