Pokud jste jako já, tak nemáte rádi v rybách kosti a první věc, která vás napadne, je si zkrátka koupit nějakou tu rybu, která nese název “bez kostí”. Takže asi máme vyhráno, což? No, ne až tak úplně …
Testovaný produkt
I přesto, že u nás často kupujeme filety bez kostí, tak se občas v některé z nich nějaká ta kůstka najde. Jednou jsem měl takové období, kdy jsem měl pocit, že je kost snad v každé druhé z nich a když jsem si přečetl informace na obalu produktu, tak žádného % počtu filet s kostmi jsem se nedohledal. A tak jsem se rozhodl to zjistit.
V tomto výzkumu se podíváme na produkt “Mylord Premium – Aljašská treska bez kosti”, který jsem prováděl 1 rok a 4 měsíce. Za celou dobu jsem zkonzumoval 35 filet, ze kterých 5 mělo kost.
> summary(df)
kost
ano: 5
ne :30 Z toho plyne, že v průměru má kost 14,3% z nich:
> prop.table(df |> table())[1]
ano
0.1428571 Confidence interval
A jak moc jsme si jistí, že je to zrovna cca 14%? To zjistíme pomoci bootstrappingu na 5000 vzorcích:
BOOTSTRAP CONFIDENCE INTERVAL CALCULATIONS
Based on 5000 bootstrap replicates
CALL :
boot.ci(boot.out = b)
Intervals :
Level Normal Basic
95% ( 0.0300, 0.2598 ) ( 0.0286, 0.2571 )
Level Percentile BCa
95% ( 0.0286, 0.2571 ) ( 0.0286, 0.2571 )
Calculations and Intervals on Original Scale
Jelikož v případě % poměru nedostaneme normální rozdělení, použijeme údaj BCa, ze kterého nám vychází, že s 95% jistotou se skutečný poměr filet s kostmi nachází v rozmezí 2,9 – 25,7%.
Jak je vidno, tak kvůli malému počtu vzorků je nepřesnost docela vysoká, ale holt přecpávat se kvůli tomu filetami se rozhodně nechystám 🙂 . Tak tedy konzumaci rybích produktů zdar a hlavně co nejméně kostí!