Densitatea de randomizare

Utilizarea a 100 de repetări îți permite să înțelegi mecanismul de permutare. Totuși, 100 de repetări nu sunt suficiente pentru a observa întreaga gamă de valori posibile pentru diferențele nule în proporții.

Amintește-ți cei patru pași ai inferenței. Acești pași sunt utilizați în toate exercițiile de inferență din acest curs și din cursurile viitoare de inferență statistică. Folosește numele funcțiilor pentru a reține mai ușor procesul de analiză.

specify va specifica variabilele de răspuns și explicative.
hypothesize va declara ipoteza nulă.
generate va genera reeșantioane, permutări sau simulări.
calculate va calcula statistici rezumative.

În acest exercițiu, vei repeta procesul de 1000 de ori pentru a înțelege distribuția completă a diferențelor nule în proporții.

Acest exercițiu face parte din cursul

Bazele inferenței în R

Vezi cursul

Instrucțiuni pentru exercițiu

Pachetele dplyr, ggplot2, NHANES și infer au fost deja încărcate pentru tine.

Generează 1000 de diferențe în proporții prin amestecarea variabilei HomeOwn folosind sintaxa infer. Amintește-ți sintaxa infer:
- specify că relația de interes este HomeOwn față de Gender și că un succes în acest context înseamnă deținerea unei locuințe, success = "Own".
- hypothesize că ipoteza nulă este adevărată, unde null = "independence" (adică genul și deținerea unei locuințe nu sunt corelate).
- generate 1000 de permutări; setează reps la 1000.
- calculate statistica stat = "diff in props" cu ordinea c("male", "female").
Rulează codul pentru graficul de densitate pentru a crea o reprezentare vizuală netezită a distribuției diferențelor. Ce formă are curba?

Exercițiu interactiv practic

Încearcă acest exercițiu completând acest cod de exemplu.

# Perform 1000 permutations
homeown_perm <- homes %>%
  # Specify HomeOwn vs. Gender, with `"Own" as success
  ___(___ ~ ___, success = "___") %>%
  # Use a null hypothesis of independence
  ___(___) %>% 
  # Generate 1000 repetitions (by permutation)
  ___(reps = ___, type = "permute") %>% 
  # Calculate the difference in proportions (male then female)
  ___(___, order = ___))

# Density plot of 1000 permuted differences in proportions
ggplot(homeown_perm, aes(x = stat)) + 
  geom_density()

Editează și rulează codul