UnaMuestra1 <- function()
{
		cat("\n")
		op = c("Intervalo para la media de una población normal.", "Intervalo para la varianza de una población normal.", "Contraste para la media de una población normal.", "Contraste para la varianza de una población normal.")
		submenu = menu(op, title="¿Qué deseas hacer?")
	
		if (submenu == 1)
		{
			cat("\n")
			cat("Introduce la muestra (para terminar pulsa dos veces enter): \n")
			muestra = scan(quiet=T) 
			media = mean(muestra)
			varianza = var(muestra)
			tamano = length(muestra)
			radio = sqrt(varianza/tamano)
			cat("\n")
			cat("Introduce el nivel de confianza (para terminar pulsa dos veces enter): \n")
			confianza = scan(quiet=T)
			confianza = 1-confianza
			alfa = 1 - (confianza/2)
			t = qt(alfa, tamano-1)
			extremoinferior = media - t * radio
			extremosuperior = media + t * radio
			cat("\n")
			cat("Intervalo de confianza al nivel", 1-confianza,"para la media de una población normal con varianza desconocida: \n")
			cat("[",extremoinferior,", ",extremosuperior,"] \n")
			cat("\n")
			resultado = list(muestra, extremoinferior, extremosuperior, media, varianza, tamano, 1-confianza, t)
			names(resultado) = c("Muestra", "Extremo inferior del intervalo", "Extremo superior del intervalo", "Media muestral", "Cuasivarianza muestral", "Tamaño de la muestra", "Nivel de confianza", "Estadístico teórico")
		}

		if (submenu == 2)
		{
			cat("\n")
			cat("Introduce la muestra (para terminar pulsa dos veces enter): \n")
			muestra = scan(quiet=T) 
			media = mean(muestra)
			varianza = var(muestra)
			tamano = length(muestra)
			radio = (tamano-1)*varianza
			cat("\n")
			cat("Introduce el nivel de confianza (para terminar pulsa dos veces enter): \n")
			confianza = scan(quiet=T)
			confianza = 1-confianza
			alfa = 1 - (confianza/2)
			chi1 = qchisq(alfa, tamano-1)
			chi2 = qchisq(1-alfa, tamano-1)
			extremoinferior = radio/chi1
			extremosuperior = radio/chi2
			cat("\n")
			cat("Intervalo de confianza al nivel", 1-confianza,"para la varianza de una población normal con media desconocida: \n")
			cat("[",extremoinferior,", ",extremosuperior,"] \n")
			cat("\n")
			resultado = list(muestra, extremoinferior, extremosuperior, media, varianza, tamano, 1-confianza, chi1, chi2)
			names(resultado) = c("Muestra", "Extremo inferior del intervalo", "Extremo superior del intervalo", "Media muestral", "Cuasivarianza muestral", "Tamaño de la muestra", "Nivel de confianza", "Estadístico teórico", "Estadístico teórico")
		}

		if (submenu == 3)
		{
			cat("\n")
			cat("Introduce la muestra (para terminar pulsa dos veces enter): \n")
			muestra = scan(quiet=T) 
			media = mean(muestra)
			varianza = var(muestra)
			tamano = length(muestra)
			radio = sqrt(varianza/tamano)
			cat("\n")
			cat("Introduce el nivel de significación (para terminar pulsa dos veces enter): \n")
			confianza = scan(quiet=T)
			cat("\n")
			cat("Introduce el valor de comparación (para terminar pulsa dos veces enter): \n")
			mu = scan(quiet=T)
			texp = (media-mu)/radio
			cat("\n")
			opc = c("Igual.", "Menor.", "Mayor.")
			sub = menu(opc, title="¿De qué tipo es la hipótesis nula?")
			if (sub == 1)
			{
				alfa = 1 - (confianza/2)
				t = qt(alfa, tamano-1)
				if (abs(texp) > t)
				{
					cat("\n")
					cat("Se rechaza la hipótesis nula de que la media poblacional sea igual a",mu)
					cat("\n")
					cat("\n")
				}
				if (abs(texp) < t)
				{
					cat("\n")
					cat("No se rechaza la hipótesis nula de que la media poblacional sea igual a",mu)
					cat("\n")
					cat("\n")
				}
			}
			if (sub == 2)
			{
				alfa = 1 - confianza
				t = qt(alfa, tamano-1)
				if (texp > t)
				{
					cat("\n")
					cat("Se rechaza la hipótesis nula de que la media poblacional sea menor a",mu)
					cat("\n")
					cat("\n")
				}
				if (texp < t)
				{
					cat("\n")
					cat("No se rechaza la hipótesis nula de que la media poblacional sea menor a",mu)
					cat("\n")
					cat("\n")
				}
			}
			if (sub == 3)
			{
				alfa = confianza
				t = qt(alfa, tamano-1)
				if (texp < t)
				{
					cat("\n")
					cat("Se rechaza la hipótesis nula de que la media poblacional sea mayor a",mu)
					cat("\n")
					cat("\n")
				}
				if (texp > t)
				{
					cat("\n")
					cat("No se rechaza la hipótesis nula de que la media poblacional sea mayor a",mu)
					cat("\n")
					cat("\n")
				}
			}
			resultado = list(muestra, texp, t, media, varianza, tamano, confianza)
			names(resultado) = c("Muestra", "Estadístico muestral", "Estadístico teórico", "Media muestral", "Cuasivarianza muestral", "Tamaño de la muestra", "Nivel de confianza")
		}

		if (submenu == 4)
		{
			cat("\n")
			cat("Introduce la muestra (para terminar pulsa dos veces enter): \n")
			muestra = scan(quiet=T) 
			media = mean(muestra)
			varianza = var(muestra)
			tamano = length(muestra)
			cat("\n")
			cat("Introduce el nivel de significación (para terminar pulsa dos veces enter): \n")
			confianza = scan(quiet=T)
			cat("\n")
			cat("Introduce el valor de comparación (para terminar pulsa dos veces enter): \n")
			var = scan(quiet=T)
			chiexp = ((tamano-1)*varianza)/var
			cat("\n")
			opc = c("Igual.", "Menor.", "Mayor.")
			sub = menu(opc, title="¿De qué tipo es la hipótesis nula?")
			if (sub == 1)
			{
				alfa = 1 - (confianza/2)
				chi1 = qchisq(alfa, tamano-1)
				chi2 = qchisq(1-alfa, tamano-1)
				if ((chiexp < chi2) || (chiexp > chi1))
				{
					cat("\n")
					cat("Se rechaza la hipótesis nula de que la varianza poblacional sea igual a",var)
					cat("\n")
					cat("\n")
				}
				if ((chiexp > chi2) & (chiexp < chi1))
				{
					cat("\n")
					cat("No se rechaza la hipótesis nula de que la varianza poblacional sea igual a",var)
					cat("\n")
					cat("\n")
				}
			resultado = list(muestra, chiexp, chi1, chi2, media, varianza, tamano, confianza)
			names(resultado) = c("Muestra", "Estadístico muestral", "Estadístico teórico", "Estadístico teórico", "Media muestral", "Cuasivarianza muestral", "Tamaño de la muestra", "Nivel de confianza")
			}
			if (sub == 2)
			{
				alfa = 1 - confianza
				chi = qchisq(alfa, tamano-1)
				if (chiexp > chi)
				{
					cat("\n")
					cat("Se rechaza la hipótesis nula de que la varianza poblacional sea menor a",var)
					cat("\n")
					cat("\n")
				}
				if (chiexp < chi)
				{
					cat("\n")
					cat("No se rechaza la hipótesis nula de que la varianza poblacional sea menor a",var)
					cat("\n")
					cat("\n")
				}
			resultado = list(muestra, chiexp, chi, media, varianza, tamano, confianza)
			names(resultado) = c("Muestra", "Estadístico muestral", "Estadístico teórico", "Media muestral", "Cuasivarianza muestral", "Tamaño de la muestra", "Nivel de confianza")
			}
			if (sub == 3)
			{
				alfa = confianza
				chi = qchisq(alfa, tamano-1)
				if (chiexp < chi)
				{
					cat("\n")
					cat("Se rechaza la hipótesis nula de que la varianza poblacional sea mayor a",var)
					cat("\n")
					cat("\n")
				}
				if (chiexp > chi)
				{
					cat("\n")
					cat("No se rechaza la hipótesis nula de que la varianza poblacional sea mayor a",var)
					cat("\n")
					cat("\n")
				}
			resultado = list(muestra, chiexp, chi, media, varianza, tamano, confianza)
			names(resultado) = c("Muestra", "Estadístico muestral", "Estadístico teórico", "Media muestral", "Cuasivarianza muestral", "Tamaño de la muestra", "Nivel de confianza")
			}
		}
		resultado
}


#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

UnaMuestra2 <- function()
{
		cat("\n")
		op = c("Intervalo para la media de una población normal con varianza conocida.", "Intervalo para la media de una población normal con varianza desconocida.", "Intervalo para la varianza de una población normal con media desconocida.", "Contraste para la media de una población normal con varianza conocida.", "Contraste para la media de una población normal con varianza desconocida.", "Contraste para la varianza de una población normal con media desconocida.")
		submenu = menu(op, title="¿Qué deseas hacer?")
	
		if (submenu == 1)
		{
			cat("\n")
			cat("Introduce la media de la muestra (para terminar pulsa dos veces enter): \n")
			media = scan(quiet=T)
			cat("\n")
			cat("Introduce el tamaño de la muestra (para terminar pulsa dos veces enter): \n")
			tamano = scan(quiet=T)
			cat("\n")
			cat("Introduce la varianza de la población (para terminar pulsa dos veces enter): \n")
			varianza = scan(quiet=T)
			radio = sqrt(varianza/tamano)
			cat("\n")
			cat("Introduce el nivel de confianza (para terminar pulsa dos veces enter): \n")
			confianza = scan(quiet=T)
			confianza = 1-confianza
			alfa = 1 - (confianza/2)
			z = qnorm(alfa)
			extremoinferior = media - z * radio
			extremosuperior = media + z * radio
			cat("\n")
			cat("Intervalo de confianza al nivel", 1-confianza,"para la media de una población normal con varianza conocida: \n")
			cat("[",extremoinferior,", ",extremosuperior,"] \n")
			cat("\n")
			resultado = list(extremoinferior, extremosuperior, media, 1-confianza, z)
			names(resultado) = c("Extremo inferior del intervalo", "Extremo superior del intervalo", "Media muestral", "Nivel de confianza", "Estadístico teórico")
		}

		if (submenu == 2)
		{
			cat("\n")
			cat("Introduce la media de la muestra (para terminar pulsa dos veces enter): \n")
			media = scan(quiet=T)
			cat("\n")
			cat("Introduce la cuasivarianza de la muestra (para terminar pulsa dos veces enter): \n")
			varianza = scan(quiet=T)
			cat("\n")
			cat("Introduce el tamaño de la muestra (para terminar pulsa dos veces enter): \n")
			tamano = scan(quiet=T)
			radio = sqrt(varianza/tamano)
			cat("\n")
			cat("Introduce el nivel de confianza (para terminar pulsa dos veces enter): \n")
			confianza = scan(quiet=T)
			confianza = 1-confianza
			alfa = 1 - (confianza/2)
			t = qt(alfa, tamano-1)
			extremoinferior = media - t * radio
			extremosuperior = media + t * radio
			cat("\n")
			cat("Intervalo de confianza al nivel", 1-confianza,"para la media de una población normal con varianza desconocida: \n")
			cat("[",extremoinferior,", ",extremosuperior,"] \n")
			cat("\n")
			resultado = list(extremoinferior, extremosuperior, media, varianza, 1-confianza, t)
			names(resultado) = c("Extremo inferior del intervalo", "Extremo superior del intervalo", "Media muestral", "Cuasivarianza muestral", "Nivel de confianza", "Estadístico teórico")
		}


		if (submenu == 3)
		{
			cat("\n")
			cat("Introduce la cuasivarianza de la muestra (para terminar pulsa dos veces enter): \n")
			varianza = scan(quiet=T) 
			cat("\n")
			cat("Introduce el tamaño de la muestra (para terminar pulsa dos veces enter): \n")
			tamano = scan(quiet=T) 
			radio = (tamano-1)*varianza
			cat("\n")
			cat("Introduce el nivel de confianza (para terminar pulsa dos veces enter): \n")
			confianza = scan(quiet=T)
			confianza = 1-confianza
			alfa = 1 - (confianza/2)
			chi1 = qchisq(alfa, tamano-1)
			chi2 = qchisq(1-alfa, tamano-1)
			extremoinferior = radio/chi1
			extremosuperior = radio/chi2
			cat("\n")
			cat("Intervalo de confianza al nivel", 1-confianza,"para la varianza de una población normal con media desconocida: \n")
			cat("[",extremoinferior,", ",extremosuperior,"] \n")
			cat("\n")
			resultado = list(extremoinferior, extremosuperior, varianza, 1-confianza, chi1, chi2)
			names(resultado) = c("Extremo inferior del intervalo", "Extremo superior del intervalo", "Cuasivarianza muestral", "Nivel de confianza", "Estadístico teórico", "Estadístico teórico")
		}

		if (submenu == 4)
		{
			cat("\n")
			cat("Introduce la media de la muestra (para terminar pulsa dos veces enter): \n")
			media = scan(quiet=T) 
			cat("\n")
			cat("Introduce el tamaño de la muestra (para terminar pulsa dos veces enter): \n")
			tamano = scan(quiet=T)
			cat("\n")
			cat("Introduce la varianza de la población (para terminar pulsa dos veces enter): \n")
			varianza = scan(quiet=T)
			radio = sqrt(varianza/tamano)
			cat("\n")
			cat("Introduce el nivel de significación (para terminar pulsa dos veces enter): \n")
			confianza = scan(quiet=T)
			cat("\n")
			cat("Introduce el valor de comparación (para terminar pulsa dos veces enter): \n")
			mu = scan(quiet=T)
			zexp = (media-mu)/radio
			cat("\n")
			opc = c("Igual.", "Menor.", "Mayor.")
			sub = menu(opc, title="¿De qué tipo es la hipótesis nula?")
			if (sub == 1)
			{
				alfa = 1 - (confianza/2)
				z = qnorm(alfa)
				if (abs(zexp) > z)
				{
					cat("\n")
					cat("Se rechaza la hipótesis nula de que la media poblacional sea igual a",mu)
					cat("\n")
					cat("\n")
				}
				if (abs(zexp) < z)
				{
					cat("\n")
					cat("No se rechaza la hipótesis nula de que la media poblacional sea igual a",mu)
					cat("\n")
					cat("\n")
				}
			}
			if (sub == 2)
			{
				alfa = 1 - confianza
				z = qnorm(alfa)
				if (zexp > z)
				{
					cat("\n")
					cat("Se rechaza la hipótesis nula de que la media poblacional sea menor a",mu)
					cat("\n")
					cat("\n")
				}
				if (zexp < z)
				{
					cat("\n")
					cat("No se rechaza la hipótesis nula de que la media poblacional sea menor a",mu)
					cat("\n")
					cat("\n")
				}
			}
			if (sub == 3)
			{
				alfa = confianza
				z = qnorm(alfa)
				if (zexp < z)
				{
					cat("\n")
					cat("Se rechaza la hipótesis nula de que la media poblacional sea mayor a",mu)
					cat("\n")
					cat("\n")
				}
				if (zexp > z)
				{
					cat("\n")
					cat("No se rechaza la hipótesis nula de que la media poblacional sea mayor a",mu)
					cat("\n")
					cat("\n")
				}
			}
			resultado = list(zexp, z, media, varianza, confianza)
			names(resultado) = c("Estadístico muestral", "Estadístico teórico", "Media muestral", "Cuasivarianza muestral", "Nivel de confianza")
		}


		if (submenu == 5)
		{
			cat("\n")
			cat("Introduce la media de la muestra (para terminar pulsa dos veces enter): \n")
			media = scan(quiet=T) 
			cat("\n")
			cat("Introduce la cuasivarianza de la muestra (para terminar pulsa dos veces enter): \n")
			varianza = scan(quiet=T)
			cat("\n")
			cat("Introduce el tamaño de la muestra (para terminar pulsa dos veces enter): \n")
			tamano = scan(quiet=T)
			radio = sqrt(varianza/tamano)
			cat("\n")
			cat("Introduce el nivel de significación (para terminar pulsa dos veces enter): \n")
			confianza = scan(quiet=T)
			cat("\n")
			cat("Introduce el valor de comparación (para terminar pulsa dos veces enter): \n")
			mu = scan(quiet=T)
			texp = (media-mu)/radio
			cat("\n")
			opc = c("Igual.", "Menor.", "Mayor.")
			sub = menu(opc, title="¿De qué tipo es la hipótesis nula?")
			if (sub == 1)
			{
				alfa = 1 - (confianza/2)
				t = qt(alfa, tamano-1)
				if (abs(texp) > t)
				{
					cat("\n")
					cat("Se rechaza la hipótesis nula de que la media poblacional sea igual a",mu)
					cat("\n")
					cat("\n")
				}
				if (abs(texp) < t)
				{
					cat("\n")
					cat("No se rechaza la hipótesis nula de que la media poblacional sea igual a",mu)
					cat("\n")
					cat("\n")
				}
			}
			if (sub == 2)
			{
				alfa = 1 - confianza
				t = qt(alfa, tamano-1)
				if (texp > t)
				{
					cat("\n")
					cat("Se rechaza la hipótesis nula de que la media poblacional sea menor a",mu)
					cat("\n")
					cat("\n")
				}
				if (texp < t)
				{
					cat("\n")
					cat("No se rechaza la hipótesis nula de que la media poblacional sea menor a",mu)
					cat("\n")
					cat("\n")
				}
			}
			if (sub == 3)
			{
				alfa = confianza
				t = qt(alfa, tamano-1)
				if (texp < t)
				{
					cat("\n")
					cat("Se rechaza la hipótesis nula de que la media poblacional sea mayor a",mu)
					cat("\n")
					cat("\n")
				}
				if (texp > t)
				{
					cat("\n")
					cat("No se rechaza la hipótesis nula de que la media poblacional sea mayor a",mu)
					cat("\n")
					cat("\n")
				}
			}
			resultado = list(texp, t, media, varianza, confianza)
			names(resultado) = c("Estadístico muestral", "Estadístico teórico", "Media muestral", "Cuasivarianza muestral", "Nivel de confianza")
		}

		if (submenu == 6)
		{
			cat("\n")
			cat("Introduce la cuasivarianza de la muestra (para terminar pulsa dos veces enter): \n")
			varianza = scan(quiet=T) 
			cat("\n")
			cat("Introduce el tamaño de la muestra (para terminar pulsa dos veces enter): \n")
			tamano = scan(quiet=T)
			cat("\n")
			cat("Introduce el nivel de significación (para terminar pulsa dos veces enter): \n")
			confianza = scan(quiet=T)
			cat("\n")
			cat("Introduce el valor de comparación (para terminar pulsa dos veces enter): \n")
			var = scan(quiet=T)
			chiexp = ((tamano-1)*varianza)/var
			cat("\n")
			opc = c("Igual.", "Menor.", "Mayor.")
			sub = menu(opc, title="¿De qué tipo es la hipótesis nula?")
			if (sub == 1)
			{
				alfa = 1 - (confianza/2)
				chi1 = qchisq(alfa, tamano-1)
				chi2 = qchisq(1-alfa, tamano-1)
				if ((chiexp < chi2) || (chiexp > chi1))
				{
					cat("\n")
					cat("Se rechaza la hipótesis nula de que la varianza poblacional sea igual a",var)
					cat("\n")
					cat("\n")
				}
				if ((chiexp > chi2) & (chiexp < chi1))
				{
					cat("\n")
					cat("No se rechaza la hipótesis nula de que la varianza poblacional sea igual a",var)
					cat("\n")
					cat("\n")
				}
			resultado = list(chiexp, chi1, chi2, varianza, confianza)
			names(resultado) = c("Estadístico muestral", "Estadístico teórico", "Estadístico teórico", "Cuasivarianza muestral", "Nivel de confianza")
			}
			if (sub == 2)
			{
				alfa = 1 - confianza
				chi = qchisq(alfa, tamano-1)
				if (chiexp > chi)
				{
					cat("\n")
					cat("Se rechaza la hipótesis nula de que la varianza poblacional sea menor a",var)
					cat("\n")
					cat("\n")
				}
				if (chiexp < chi)
				{
					cat("\n")
					cat("No se rechaza la hipótesis nula de que la varianza poblacional sea menor a",var)
					cat("\n")
					cat("\n")
				}
			resultado = list(chiexp, chi, varianza, confianza)
			names(resultado) = c("Estadístico muestral", "Estadístico teórico", "Cuasivarianza muestral", "Nivel de confianza")
			}
			if (sub == 3)
			{
				alfa = confianza
				chi = qchisq(alfa, tamano-1)
				if (chiexp < chi)
				{
					cat("\n")
					cat("Se rechaza la hipótesis nula de que la varianza poblacional sea mayor a",var)
					cat("\n")
					cat("\n")
				}
				if (chiexp > chi)
				{
					cat("\n")
					cat("No se rechaza la hipótesis nula de que la varianza poblacional sea mayor a",var)
					cat("\n")
					cat("\n")
				}
			resultado = list(chiexp, chi, varianza, confianza)
			names(resultado) = c("Estadístico muestral", "Estadístico teórico", "Cuasivarianza muestral", "Nivel de confianza")
			}
		}

		resultado
}


#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


UnaMuestra <- function()
{

	cat("\n")
	opciones = c("A partir de la muestra.", "A partir de las características poblacionales y muestrales.")
	menuprincipal = menu(opciones, title="Inferencia estadística para una distribución procedente de una población normal")

	if (menuprincipal == 1)
	{
		resultado = UnaMuestra1()
	}

	if (menuprincipal == 2)
	{
		resultado = UnaMuestra2()
	}

	resultado
}