DosMuestras1 <- function()
{
		cat("\n")
		op = c("Intervalo para la diferencia de medias de dos poblaciones normales con varianzas desconocidas e iguales.", "Intervalo para el cociente de varianzas de dos poblaciones normales.", "Contraste para la diferencia de medias de dos poblaciones normales con varianzas desconocidas e iguales.", "Contraste para el cociente de varianzas de dos poblaciones normales.")
		submenu = menu(op, title="¿Qué deseas hacer?")
	
		if (submenu == 1)
		{
			cat("\n")
			cat("Introduce la primera muestra (para terminar pulsa dos veces enter): \n")
			muestra1 = scan(quiet=T) 
			media1 = mean(muestra1)
			varianza1 = var(muestra1)
			tamano1 = length(muestra1)
			cat("\n")
			cat("Introduce la segunda muestra (para terminar pulsa dos veces enter): \n")
			muestra2 = scan(quiet=T) 
			media2 = mean(muestra2)
			varianza2 = var(muestra2)
			tamano2 = length(muestra2)
			radio1 = sqrt(((tamano1-1)*varianza1+(tamano2-1)*varianza2)/(tamano1+tamano2-2))
			radio2 = sqrt((tamano1+tamano2)/(tamano1*tamano2))
			radio = radio1*radio2
			media = media1-media2
			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, tamano1+tamano2-2)
			extremoinferior = media - t * radio
			extremosuperior = media + t * radio
			cat("\n")
			cat("Intervalo de confianza al nivel", 1-confianza,"para la diferencia de medias de dos poblaciones normales con varianzas desconocidas e iguales: \n")
			cat("[",extremoinferior,", ",extremosuperior,"] \n")
			cat("\n")
			resultado = list(muestra1, muestra2, extremoinferior, extremosuperior, media1, varianza1, media2, varianza2, tamano1, tamano2, 1-confianza, t)
			names(resultado) = c("Primera muestra", "Segunda muestra", "Extremo inferior del intervalo", "Extremo superior del intervalo", "Primera media muestral", "Primera cuasivarianza muestral", "Segunda media muestral", "Segunda cuasivarianza muestral", "Tamaño primera muestra", "Tamaño segunda muestra", "Nivel de confianza", "Estadístico teórico")
		}

		if (submenu == 2)
		{
			cat("\n")
			cat("Introduce la primera muestra (para terminar pulsa dos veces enter): \n")
			muestra1 = scan(quiet=T) 
			media1 = mean(muestra1)
			varianza1 = var(muestra1)
			tamano1 = length(muestra1)
			cat("\n")
			cat("Introduce la segunda muestra (para terminar pulsa dos veces enter): \n")
			muestra2 = scan(quiet=T) 
			media2 = mean(muestra2)
			varianza2 = var(muestra2)
			tamano2 = length(muestra2)
			radio = varianza2/varianza1
			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)
			F2 = qf(alfa, tamano1-1, tamano2-1)
			F1 = qf(1-alfa, tamano1-1, tamano2-1)
			extremoinferior = radio*F1
			extremosuperior = radio*F2
			cat("\n")
			cat("Intervalo de confianza al nivel", 1-confianza,"para el cociente de varianzas de dos poblaciones normales con medias desconocidas: \n")
			cat("[",extremoinferior,", ",extremosuperior,"] \n")
			cat("\n")
			resultado = list(muestra1, muestra2, extremoinferior, extremosuperior, media1, varianza1, media2, varianza2, tamano1, tamano2, 1-confianza, F1, F2)
			names(resultado) = c("Primera muestra", "Segunda muestra", "Extremo inferior del intervalo", "Extremo superior del intervalo", "Primera media muestral", "Primera cuasivarianza muestral", "Segunda media muestral", "Segunda cuasivarianza muestral", "Tamaño primera muestra", "Tamaño segunda muestra", "Nivel de confianza", "Estadístico teórico", "Estadístico teórico")
		}

		if (submenu == 3)
		{
			cat("\n")
			cat("Introduce la primera muestra (para terminar pulsa dos veces enter): \n")
			muestra1 = scan(quiet=T) 
			media1 = mean(muestra1)
			varianza1 = var(muestra1)
			tamano1 = length(muestra1)
			cat("\n")
			cat("Introduce la segunda muestra (para terminar pulsa dos veces enter): \n")
			muestra2 = scan(quiet=T) 
			media2 = mean(muestra2)
			varianza2 = var(muestra2)
			tamano2 = length(muestra2)
			radio1 = sqrt(((tamano1-1)*varianza1+(tamano2-1)*varianza2)/(tamano1+tamano2-2))
			radio2 = sqrt((tamano1+tamano2)/(tamano1*tamano2))
			radio = radio1*radio2
			media = media1-media2
			cat("\n")
			cat("Introduce el nivel de significación (para terminar pulsa dos veces enter): \n")
			confianza = scan(quiet=T)
			texp = media/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, tamano1+tamano2-2)
				if (abs(texp) > t)
				{
					cat("\n")
					cat("Se rechaza la hipótesis nula de igualdad de medias poblacionales")
					cat("\n")
					cat("\n")
				}
				if (abs(texp) < t)
				{
					cat("\n")
					cat("No se rechaza la hipótesis nula de igualdad de medias poblacionales")
					cat("\n")
					cat("\n")
				}
			}
			if (sub == 2)
			{
				alfa = 1 - confianza
				t = qt(alfa, tamano1+tamano2-2)
				if (texp > t)
				{
					cat("\n")
					cat("Se rechaza la hipótesis nula de que la media de la primera población sea menor que la de la segunda")
					cat("\n")
					cat("\n")
				}
				if (texp < t)
				{
					cat("\n")
					cat("No se rechaza la hipótesis nula de que la media de la primera población sea menor que la de la segunda")
					cat("\n")
					cat("\n")
				}
			}
			if (sub == 3)
			{
				alfa = confianza
				t = qt(alfa, tamano1+tamano2-2)
				if (texp < t)
				{
					cat("\n")
					cat("Se rechaza la hipótesis nula de que la media de la primera población sea mayor que la de la segunda")
					cat("\n")
					cat("\n")
				}
				if (texp > t)
				{
					cat("\n")
					cat("No se rechaza la hipótesis nula de que la media de la primera población sea mayor que la de la segunda")
					cat("\n")
					cat("\n")
				}
			}
			resultado = list(muestra1, muestra2, texp, t, media1, varianza1, media2, varianza2, tamano1, tamano2, confianza)
			names(resultado) = c("Primera muestra", "Segunda muestra", "Estadístico muestral", "Estadístico teórico", "Primera media muestral", "Primera cuasivarianza muestral", "Segunda media muestral", "Segunda cuasivarianza muestral", "Tamaño primera muestra", "Tamaño segunda muestra", "Nivel de confianza")
		}

		if (submenu == 4)
		{
			cat("\n")
			cat("Introduce la primera muestra (para terminar pulsa dos veces enter): \n")
			muestra1 = scan(quiet=T) 
			media1 = mean(muestra1)
			varianza1 = var(muestra1)
			tamano1 = length(muestra1)
			cat("\n")
			cat("Introduce la segunda muestra (para terminar pulsa dos veces enter): \n")
			muestra2 = scan(quiet=T) 
			media2 = mean(muestra2)
			varianza2 = var(muestra2)
			tamano2 = length(muestra2)
			cat("\n")
			cat("Introduce el nivel de significación (para terminar pulsa dos veces enter): \n")
			confianza = scan(quiet=T)
			Fexp = varianza1/varianza2
			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)
				F2 = qf(alfa, tamano1-1, tamano2-1)
				F1 = qf(1-alfa, tamano1-1, tamano2-1)
				if ((Fexp < F1) || (Fexp > F2))
				{
					cat("\n")
					cat("Se rechaza la hipótesis nula de igualdad de varianzas poblacionales")
					cat("\n")
					cat("\n")
				}
				if ((Fexp > F1) & (Fexp < F2))
				{
					cat("\n")
					cat("No se rechaza la hipótesis nula de igualdad de varianzas poblacionales")
					cat("\n")
					cat("\n")
				}
			resultado = list(muestra1, muestra2, Fexp, F1, F2, media1, varianza1, media2, varianza2, tamano1, tamano2, confianza)
			names(resultado) = c("Primera muestra", "Segunda muestra", "Estadístico muestral", "Estadístico teórico", "Estadístico teórico", "Primera media muestral", "Primera cuasivarianza muestral", "Segunda media muestral", "Segunda cuasivarianza muestral", "Tamaño primera muestra", "Tamaño segunda muestra", "Nivel de confianza")
			}
			if (sub == 2)
			{
				alfa = 1 - confianza
				F = qf(alfa, tamano1-1, tamano2-1)
				if (Fexp > F)
				{
					cat("\n")
					cat("Se rechaza la hipótesis nula de que la varianza de la primera población sea menor que la de la segunda")
					cat("\n")
					cat("\n")
				}
				if (Fexp < F)
				{
					cat("\n")
					cat("No se rechaza la hipótesis nula de que la varianza de la primera población sea menor que la de la segunda")
					cat("\n")
					cat("\n")
				}
			resultado = list(muestra1, muestra2, Fexp, F, media1, varianza1, media2, varianza2, tamano1, tamano2, confianza)
			names(resultado) = c("Primera muestra", "Segunda muestra", "Estadístico muestral", "Estadístico teórico", "Primera media muestral", "Primera cuasivarianza muestral", "Segunda media muestral", "Segunda cuasivarianza muestral", "Tamaño primera muestra", "Tamaño segunda muestra", "Nivel de confianza")
			}
			if (sub == 3)
			{
				alfa = confianza
				F = qf(alfa, tamano1-1, tamano2-1)
				if (Fexp < F)
				{
					cat("\n")
					cat("Se rechaza la hipótesis nula de que la varianza de la primera población sea mayor que la de la segunda")
					cat("\n")
					cat("\n")
				}
				if (Fexp > F)
				{
					cat("\n")
					cat("No se rechaza la hipótesis nula de que la varianza de la primera población sea mayor que la de la segunda")
					cat("\n")
					cat("\n")
				}
			resultado = list(muestra1, muestra2, Fexp, F, media1, varianza1, media2, varianza2, tamano1, tamano2, confianza)
			names(resultado) = c("Primera muestra", "Segunda muestra", "Estadístico muestral", "Estadístico teórico", "Primera media muestral", "Primera cuasivarianza muestral", "Segunda media muestral", "Segunda cuasivarianza muestral", "Tamaño primera muestra", "Tamaño segunda muestra", "Nivel de confianza")
			}
		}
		resultado
}


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

DosMuestras2 <- function()
{
		cat("\n")
		op = c("Intervalo para la diferencia de medias de dos poblaciones normales con varianzas conocidas.", "Intervalo para la diferencia de medias de dos poblaciones normales con varianzas desconocidas e iguales.", "Intervalo para el cociente de varianzas de dos poblaciones normales con medias desconocidas.", "Contraste para la diferencia de medias de dos poblaciones normales con varianzas conocidas.", "Contraste para la diferencia de medias de dos poblaciones normales con varianzas desconocidas e iguales.", "Contraste para el cociente de varianzas de dos poblaciones normales con medias desconocidas.")
		submenu = menu(op, title="¿Qué deseas hacer?")
	
		if (submenu == 1)
		{
			cat("\n")
			cat("Introduce la media de la primera muestra (para terminar pulsa dos veces enter): \n")
			media1 = scan(quiet=T)
			cat("\n")
			cat("Introduce el tamaño de la primera muestra (para terminar pulsa dos veces enter): \n")
			tamano1 = scan(quiet=T)
			cat("\n")
			cat("Introduce la varianza de la primera población (para terminar pulsa dos veces enter): \n")
			varianza1 = scan(quiet=T)
			cat("\n")
			cat("Introduce la media de la segunda muestra (para terminar pulsa dos veces enter): \n")
			media2 = scan(quiet=T)
			cat("\n")
			cat("Introduce el tamaño de la segunda muestra (para terminar pulsa dos veces enter): \n")
			tamano2 = scan(quiet=T)
			cat("\n")
			cat("Introduce la varianza de la segunda población (para terminar pulsa dos veces enter): \n")
			varianza2 = scan(quiet=T)
			radio = sqrt((varianza1/tamano1)+(varianza2/tamano2))
			media = media1-media2
			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 diferencia de medias de dos poblaciones normales con varianzas conocidas: \n")
			cat("[",extremoinferior,", ",extremosuperior,"] \n")
			cat("\n")
			resultado = list(extremoinferior, extremosuperior, media1, media2, 1-confianza, z)
			names(resultado) = c("Extremo inferior del intervalo", "Extremo superior del intervalo", "Media primera muestra", "Media segunda muestra", "Nivel de confianza", "Estadístico teórico")
		}

		if (submenu == 2)
		{
			cat("\n")
			cat("Introduce la media de la primera muestra (para terminar pulsa dos veces enter): \n")
			media1 = scan(quiet=T)
			cat("\n")
			cat("Introduce la cuasivarianza de la primera muestra (para terminar pulsa dos veces enter): \n")
			varianza1 = scan(quiet=T)
			cat("\n")
			cat("Introduce el tamaño de la primera muestra (para terminar pulsa dos veces enter): \n")
			tamano1 = scan(quiet=T)
			cat("\n")
			cat("Introduce la media de la segunda muestra (para terminar pulsa dos veces enter): \n")
			media2 = scan(quiet=T)
			cat("\n")
			cat("Introduce la cuasivarianza de la segunda muestra (para terminar pulsa dos veces enter): \n")
			varianza2 = scan(quiet=T)
			cat("\n")
			cat("Introduce el tamaño de la segunda muestra (para terminar pulsa dos veces enter): \n")
			tamano2 = scan(quiet=T)
			radio1 = sqrt(((tamano1-1)*varianza1+(tamano2-1)*varianza2)/(tamano1+tamano2-2))
			radio2 = sqrt((tamano1+tamano2)/(tamano1*tamano2))
			radio = radio1*radio2
			media = media1-media2
			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, tamano1+tamano2-2)
			extremoinferior = media - t * radio
			extremosuperior = media + t * radio
			cat("\n")
			cat("Intervalo de confianza al nivel", 1-confianza,"para la diferencia de medias de dos poblaciones normales con varianza desconocidas e iguales: \n")
			cat("[",extremoinferior,", ",extremosuperior,"] \n")
			cat("\n")
			resultado = list(extremoinferior, extremosuperior, media1, varianza1, media2, varianza2, 1-confianza, t)
			names(resultado) = c("Extremo inferior del intervalo", "Extremo superior del intervalo", "Media primera muestra", "Cuasivarianza primera muestra", "Media segunda muestra", "Cuasivarianza segunda muestra", "Nivel de confianza", "Estadístico teórico")
		}


		if (submenu == 3)
		{
			cat("\n")
			cat("Introduce la cuasivarianza de la primera muestra (para terminar pulsa dos veces enter): \n")
			varianza1 = scan(quiet=T) 
			cat("\n")
			cat("Introduce el tamaño de la primera muestra (para terminar pulsa dos veces enter): \n")
			tamano1 = scan(quiet=T)
			cat("\n")
			cat("Introduce la cuasivarianza de la segunda muestra (para terminar pulsa dos veces enter): \n")
			varianza2 = scan(quiet=T) 
			cat("\n")
			cat("Introduce el tamaño de la segunda muestra (para terminar pulsa dos veces enter): \n")
			tamano2 = scan(quiet=T) 
			radio = varianza2/varianza1
			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)
			F2 = qf(alfa, tamano1-1, tamano2-1)
			F1 = qf(1-alfa, tamano1-1, tamano2-1)
			extremoinferior = radio*F1
			extremosuperior = radio*F2
			cat("\n")
			cat("Intervalo de confianza al nivel", 1-confianza,"para el cociente de varianzas de dos poblaciones normales con medias desconocidas: \n")
			cat("[",extremoinferior,", ",extremosuperior,"] \n")
			cat("\n")
			resultado = list(extremoinferior, extremosuperior, varianza1, varianza2, 1-confianza, F1, F2)
			names(resultado) = c("Extremo inferior del intervalo", "Extremo superior del intervalo", "Primera cuasivarianza muestral", "Segunda cuasivarianza muestral", "Nivel de confianza", "Estadístico teórico", "Estadístico teórico")
		}

		if (submenu == 4)
		{
			cat("\n")
			cat("Introduce la media de la primera muestra (para terminar pulsa dos veces enter): \n")
			media1 = scan(quiet=T)
			cat("\n")
			cat("Introduce el tamaño de la primera muestra (para terminar pulsa dos veces enter): \n")
			tamano1 = scan(quiet=T)
			cat("\n")
			cat("Introduce la varianza de la primera población (para terminar pulsa dos veces enter): \n")
			varianza1 = scan(quiet=T)
			cat("\n")
			cat("Introduce la media de la segunda muestra (para terminar pulsa dos veces enter): \n")
			media2 = scan(quiet=T)
			cat("\n")
			cat("Introduce el tamaño de la segunda muestra (para terminar pulsa dos veces enter): \n")
			tamano2 = scan(quiet=T)
			cat("\n")
			cat("Introduce la varianza de la segunda población (para terminar pulsa dos veces enter): \n")
			varianza2 = scan(quiet=T)
			radio = sqrt((varianza1/tamano1)+(varianza2/tamano2))
			media = media1-media2
			cat("\n")
			cat("Introduce el nivel de significación (para terminar pulsa dos veces enter): \n")
			confianza = scan(quiet=T)
			zexp = media/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 igualdad de medias poblacionales")
					cat("\n")
					cat("\n")
				}
				if (abs(zexp) < z)
				{
					cat("\n")
					cat("No se rechaza la hipótesis nula de igualdad de medias poblacionales")
					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 primera media poblacional sea menor a la segunda media poblacional")
					cat("\n")
					cat("\n")
				}
				if (zexp < z)
				{
					cat("\n")
					cat("No se rechaza la hipótesis nula de que la primera media poblacional sea menor a la segunda media poblacional")
					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 primera media poblacional sea mayor a la segunda media poblacional")
					cat("\n")
					cat("\n")
				}
				if (zexp > z)
				{
					cat("\n")
					cat("No se rechaza la hipótesis nula de que la primera media poblacional sea mayor a la segunda media poblacional")
					cat("\n")
					cat("\n")
				}
			}
			resultado = list(zexp, z, media1, varianza1, media2, varianza2, confianza)
			names(resultado) = c("Estadístico muestral", "Estadístico teórico", "Primera media muestral", "Primera cuasivarianza muestral", "Segunda media muestral", "Segunda cuasivarianza muestral", "Nivel de confianza")
		}


		if (submenu == 5)
		{
			cat("\n")
			cat("Introduce la media de la primera muestra (para terminar pulsa dos veces enter): \n")
			media1 = scan(quiet=T)
			cat("\n")
			cat("Introduce la cuasivarianza de la primera muestra (para terminar pulsa dos veces enter): \n")
			varianza1 = scan(quiet=T)
			cat("\n")
			cat("Introduce el tamaño de la primera muestra (para terminar pulsa dos veces enter): \n")
			tamano1 = scan(quiet=T)
			cat("\n")
			cat("Introduce la media de la segunda muestra (para terminar pulsa dos veces enter): \n")
			media2 = scan(quiet=T)
			cat("\n")
			cat("Introduce la cuasivarianza de la segunda muestra (para terminar pulsa dos veces enter): \n")
			varianza2 = scan(quiet=T)
			cat("\n")
			cat("Introduce el tamaño de la segunda muestra (para terminar pulsa dos veces enter): \n")
			tamano2 = scan(quiet=T)
			radio1 = sqrt(((tamano1-1)*varianza1+(tamano2-1)*varianza2)/(tamano1+tamano2-2))
			radio2 = sqrt((tamano1+tamano2)/(tamano1*tamano2))
			radio = radio1*radio2
			media = media1-media2
			cat("\n")
			cat("Introduce el nivel de significación (para terminar pulsa dos veces enter): \n")
			confianza = scan(quiet=T)
			texp = media/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, tamano1+tamano2-2)
				if (abs(texp) > t)
				{
					cat("\n")
					cat("Se rechaza la hipótesis nula de igualdad de medias poblacionales")
					cat("\n")
					cat("\n")
				}
				if (abs(texp) < t)
				{
					cat("\n")
					cat("No se rechaza la hipótesis nula de igualdad de medias poblacionales")
					cat("\n")
					cat("\n")
				}
			}
			if (sub == 2)
			{
				alfa = 1 - confianza
				t = qt(alfa, tamano1+tamano2-2)
				if (texp > t)
				{
					cat("\n")
					cat("Se rechaza la hipótesis nula de que la primera media poblacional sea menor a la segunda media poblacional")
					cat("\n")
					cat("\n")
				}
				if (texp < t)
				{
					cat("\n")
					cat("No se rechaza la hipótesis nula de que la primera media poblacional sea menor a la segunda media poblacional")
					cat("\n")
					cat("\n")
				}
			}
			if (sub == 3)
			{
				alfa = confianza
				t = qt(alfa, tamano1+tamano2-2)
				if (texp < t)
				{
					cat("\n")
					cat("Se rechaza la hipótesis nula de que la primera media poblacional sea mayor a la segunda media poblacional")
					cat("\n")
					cat("\n")
				}
				if (texp > t)
				{
					cat("\n")
					cat("No se rechaza la hipótesis nula de que la primera media poblacional sea mayor a la segunda media poblacional")
					cat("\n")
					cat("\n")
				}
			}
			resultado = list(texp, t, media1, varianza1, media2, varianza2, confianza)
			names(resultado) = c("Estadístico muestral", "Estadístico teórico", "Primera media muestral", "Primera cuasivarianza muestral", "Segunda media muestral", "Segunda cuasivarianza muestral", "Nivel de confianza")
		}

		if (submenu == 6)
		{
			cat("\n")
			cat("Introduce la cuasivarianza de la primera muestra (para terminar pulsa dos veces enter): \n")
			varianza1 = scan(quiet=T) 
			cat("\n")
			cat("Introduce el tamaño de la primera muestra (para terminar pulsa dos veces enter): \n")
			tamano1 = scan(quiet=T)
			cat("\n")
			cat("Introduce la cuasivarianza de la segunda muestra (para terminar pulsa dos veces enter): \n")
			varianza2 = scan(quiet=T) 
			cat("\n")
			cat("Introduce el tamaño de la segunda muestra (para terminar pulsa dos veces enter): \n")
			tamano2 = scan(quiet=T) 
			cat("\n")
			cat("Introduce el nivel de significación (para terminar pulsa dos veces enter): \n")
			confianza = scan(quiet=T)
			Fexp = varianza1/varianza2
			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)
				F1 = qf(alfa, tamano1-1, tamano2-1)
				F2 = qf(1-alfa, tamano1-1, tamano2-1)
				if ((Fexp < F2) || (Fexp > F1))
				{
					cat("\n")
					cat("Se rechaza la hipótesis nula de igualdad de varianzas poblacionales")
					cat("\n")
					cat("\n")
				}
				if ((Fexp > F2) & (Fexp < F1))
				{
					cat("\n")
					cat("No se rechaza la hipótesis nula de igualdad de varianzas poblacionales")
					cat("\n")
					cat("\n")
				}
			resultado = list(Fexp, F1, F2, varianza1, varianza2, confianza)
			names(resultado) = c("Estadístico muestral", "Estadístico teórico", "Estadístico teórico", "Primera cuasivarianza muestral", "Segunda cuasivarianza muestral", "Nivel de confianza")
			}
			if (sub == 2)
			{
				alfa = 1 - confianza
				F = qf(alfa, tamano1-1, tamano2-1)
				if (Fexp > F)
				{
					cat("\n")
					cat("Se rechaza la hipótesis nula de que la varianza de la primera población sea menor que la de la segunda")
					cat("\n")
					cat("\n")
				}
				if (Fexp < F)
				{
					cat("\n")
					cat("No se rechaza la hipótesis nula de que la varianza de la primera población sea menor que la de la segunda")
					cat("\n")
					cat("\n")
				}
			resultado = list(Fexp, F, varianza1, varianza2, confianza)
			names(resultado) = c("Estadístico muestral", "Estadístico teórico", "Primera cuasivarianza muestral", "Segunda cuasivarianza muestral", "Nivel de confianza")
			}
			if (sub == 3)
			{
				alfa = confianza
				F = qf(alfa, tamano1-1, tamano2-1)
				if (Fexp < F)
				{
					cat("\n")
					cat("Se rechaza la hipótesis nula de que la varianza de la primera población sea mayor que la de la segunda")
					cat("\n")
					cat("\n")
				}
				if (Fexp > F)
				{
					cat("\n")
					cat("No se rechaza la hipótesis nula de que la varianza de la primera población sea mayor que la de la segunda")
					cat("\n")
					cat("\n")
				}
			resultado = list(Fexp, F, varianza1, varianza2, confianza)
			names(resultado) = c("Estadístico muestral", "Estadístico teórico", "Primera cuasivarianza muestral", "Segunda cuasivarianza muestral", "Nivel de confianza")
			}
		}

		resultado
}


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


DosMuestras <- function()
{

	cat("\n")
	opciones = c("A partir de las muestras.", "A partir de las características poblacionales y muestrales.")
	menuprincipal = menu(opciones, title="Inferencia estadística para distribuciones procedentes de dos poblaciones normales")

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

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

	resultado
}