Rastgeleleştirme yoğunluğu
100 tekrar, permütasyon mekanizmasını anlamana yardımcı olur. Ancak, oranlardaki sıfır hipotezi farklarının olası değerlerinin tüm aralığını gözlemlemek için 100 yeterli değildir.
Çıkarımın dört adımını hatırla. Bu adımlar, bu kurstaki ve gelecekteki istatistiksel çıkarım kurslarındaki tüm çıkarım egzersizlerinde kullanılacak aynı dört adımdır. Analiz sürecini hatırlamana yardımcı olması için fonksiyon adlarını kullan.
specifyyanıt ve açıklayıcı değişkenleri belirler.hypothesizesıfır hipotezini tanımlar.generateyeniden örneklemeleri, permütasyonları veya benzetimleri (simülasyonları) üretir.calculateözet istatistikleri hesaplar.
Bu egzersizde, oranlardaki sıfır hipotezi farklarının tam dağılımını görmek için işlemi 1000 kez tekrarlayacaksın.
Bu egzersiz, kursun bir parçasıdır
R ile Çıkarımın Temelleri
Egzersiz talimatları
dplyr, ggplot2, NHANES ve infer paketleri senin için yüklendi.
infersözdizimini kullanarakHomeOwndeğişkenini karıştırıp 1000 oran farkı üret.infersözdizimini hatırla:- İlgilendiğin ilişkiyi
HomeOwnileGenderarasındaspecifyet ve bu bağlamda başarı ev sahipliğidir:success = "Own". null = "independence"olduğu (yani cinsiyet ile ev sahipliğinin ilişkili olmadığı) sıfır hipotezinihypothesizeet.- 1000 permütasyon
generateet;repsdeğerini 1000 olarak ayarla. stat = "diff in props"istatistiğinic("male", "female")sırasıylacalculateet.
- İlgilendiğin ilişkiyi
- Farkların dağılımının yumuşatılmış görsel temsilini oluşturmak için yoğunluk grafiği kodunu çalıştır. Eğrinin şekli nasıl görünüyor?
Uygulamalı etkileşimli egzersiz
Bu egzersizi bu örnek kodu tamamlayarak deneyin.
# 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()