Tarea 1 - Hardy-Weinberg

---

Hardy Weinberg test

props <- c(189, 89, 9)
total <- sum(props)

Interpretación: Hay 189 individuos, homocigotos dominantes con el genotipo AA, 89 individuos, heterocigotos con el genotipo Aa y 9 individuos, homocigotos recesivos con el genotipo aa, el total de la muestra de la población es de 287 individuos

pA <- (2 * props[1] + props[2]) / (2 * total)
pa <- 1 - pA


print(pA)

[1] 0.8135889

print(pa)

[1] 0.1864111

interpretación: La frecuencia del alelo dominante A es 0.81 y la frecuencia del alelo recesivo a es de 0.18.

{+{total <- sum(props)} print(props) print(total)

expected <- c(
  pA^2 * total,
  2 * pA * pa * total,
  pa^2 * total
)

print(expected)

[1] 189.972997  87.054007   9.972997

Interpretación: Bajo el sistema de referencia de equilibrio Hardy-weinberg, los valores en los conteos esperados fueron los siguientes: AA=189.97, Aa=87.05 y aa=9.97.

chisq <- sum((props - expected)^2 / expected)
p.value <- pchisq(chisq, df = 1, lower.tail = FALSE)

print(chisq)

[1] 0.1434125

print(p.value)

[1] 0.7049118

Interpretación: El valorcalculado de chi-cuadrado fue 0.143, con un p ≈ 0.7050 > 05, lo que implica queno se descarta la hipótesis nula. Esto indica que las frecuencias genotípicas que se observaron coinciden con las que se esperaban según el modelo de Hardy-Weinberg, sugiriendo que la población se encuentra en equilibrio.

Tarea:Ejercicio de suma y frecuencia en R

sumar_elementos <-function(vector_numeros)
{suma_total <- sum(vector_numeros)
return(suma_total)}

vec_exercise <- c(15, 35, 10, 27, 40, 5)

resultado_suma <- sumar_elementos(vec_exercise)

print(paste("La suma de los elementos es:", resultado_suma))

[1] "La suma de los elementos es: 132"

vec_frecuencia <- function(vector_datos) {
  frecuencia_a <- sum(vector_datos == "A")
  return(frecuencia_a)}


mis_datos <- c("A", "A", "A","B")

frecuencia_resultado <- vec_frecuencia(mis_datos)

print(paste("La frecuencia de 'A' es:", frecuencia_resultado))

[1] "La frecuencia de 'A' es: 3"