Densité de randomisation

Utiliser 100 répétitions vous aide à comprendre le mécanisme de permutation. Toutefois, 100 ne suffit pas pour observer toute la plage des valeurs plausibles pour les différences nulles de proportions.

Rappelez-vous les quatre étapes de l'inférence. Ce sont les mêmes quatre étapes qui seront utilisées dans tous les exercices d'inférence de ce cours et des cours ultérieurs sur l'inférence statistique. Appuyez-vous sur les noms des fonctions pour vous rappeler le processus d'analyse.

specify sert à définir les variables réponse et explicative.
hypothesize sert à énoncer l'hypothèse nulle.
generate sert à produire des rééchantillonnages, des permutations ou des simulations.
calculate sert à calculer les statistiques sommaires.

Dans cet exercice, vous répéterez le processus 1000 fois pour avoir un aperçu de la distribution complète des différences nulles de proportions.

Cette activité fait partie du cours

Fondements de l'inférence en R

Voir le cours

Instructions de l’exercice

Les bibliothèques dplyr, ggplot2, NHANES et infer ont été chargées pour vous.

Générez 1000 différences de proportions en brassant la variable HomeOwn avec la syntaxe d'infer. Rappel de la syntaxe d'infer :
- specify que la relation d'intérêt est HomeOwn par rapport à Gender et qu'un succès ici correspond à être propriétaire, success = "Own".
- hypothesize que l'hypothèse nulle est vraie avec null = "independence" (ce qui signifie que le genre et le fait d'être propriétaire ne sont pas liés).
- generate 1000 permutations; définissez reps à 1000.
- calculate la statistique stat = "diff in props" avec l'ordre c("male", "female").
Exécutez le code du diagramme de densité pour créer une représentation lissée de la distribution des différences. Quelle forme a la courbe?

Exercice interactif pratique

Essayez cet exercice en complétant ce code d’exemple.

# 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()

Modifier et exécuter le code