Jitn kvality technologickch proces 12 Statistick analza dat

Jištění kvality technologických procesů 12. Statistická analýza dat Richard Březina

Statistická analýza dat o Základní pojmy o Průzkumová analýza dat o Vyhodnocování trendů

Statistická analýza dat þ Základní pojmy o Statistické šetření o Charakteristiky polohy a variability o Grafické metody o Průzkumová analýza dat o Vyhodnocování trendů

Statistická analýza dat o Základní pojmy þ Průzkumová analýza dat o Výpočet popisných a intervalových charakteristik o Grafické metody o Vyhodnocování trendů

Statistická analýza dat o Základní pojmy o Průzkumová analýza dat þ Vyhodnocování o trendů Tvorba a použití regulačních diagramů

Základní pojmy l Statistické zkoumání l Ve statistickém zkoumání se sledují hromadné jevy, tj. jevy vyskytující se u velkého množství prvků. Z hlediska účelu zkoumání hromadných jevů se předpokládá definování vymezené množiny objektů, prvků zkoumání neboli statistického souboru (soubor podniků, soubor obyvatelstva, soubor událostí apod. ). Počet objektů, prvků, jednotek statistického souboru se nazývá rozsah souboru. l Základní soubor l Statistické soubory, které jsou předmětem statistického zkoumání, jsou označovány jako základní soubor. l Výběrový soubor l Z důvodu ekonomické náročnosti, praktické složitosti a velkých rozsahů se v praxi zpravidla nepracuje s celým rozsahem statistického souboru, ale jen se vzorkem statistických jednotek neboli s výběrovým souborem. Výběrové soubory se používají ke zkoumání zejména proto, že zobecnění provedené z dat výběrového souboru je považováno pro daný účel zkoumání za dostatečně přesné a z hlediska poznání za reprezentativní a že zkoumání celého statistického souboru by bylo nákladné, časově zdlouhavé nebo z jiných praktických ohledů neuskutečnitelné. Zkoumané vlastnosti statistického souboru sleduje statistika prostřednictvím měřitelných vlastností statistických jednotek. tzv. statistických znaků. Pracuje-li se u každé statistické jednotky jen s jedním znakem (s jednou proměnnou), jedná se o jednorozměrný soubor. Zjišťuje-li se současně více znaků a jejich vzájemné vztahy, jde o dvou-, tří-, resp. obecně vícerozměrný soubor. l

Základní pojmy l Kvalitativní a kvantitativní znaky l Základním tříděním statistických znaků je rozlišování znaků číselných (kvantitativních, numerických) a znaků slovních (kvalitativních, alfabetických, kategoriálních). Číselné (kvantitativní, numerické) statistické znaky vyjadřují sledované vlastnosti statistických jednotek čísly. Číselné znaky se dělí na spojité (kontinuální), které mohou teoreticky nabývat libovolných reálných číselných hodnot v určitém intervalu (průtok vody, hmotnost výrobku, výška, peněžní obrat apod. ) a znaky nespojité (diskrétní), které mohou v oboru reálných čísel nabývat pouze určitých číselných hodnot l l Cíle statistického zjišťování l Cílem statistického šetření je získávání statistických dat o hodnotách statistických znaků u prvků, které tvoří statistický soubor. Každé statistické zjišťování má určitý, předem vymezený a konkrétní cíl. Z definice cílů šetření vyplývá, l l jak vymezíme statistický soubor jaké proměnné statistické znaky budeme zjišťovat. l Zpracování rutinních dat l Pro účely této přednášky se omezíme se na zpracování jednorozměrného souboru dat, u něhož sledujeme jednu vlastnost, měřením jednoho spojitého znaku

Základní pojmy l Postup průzkumové analýzy dat l Při statistickém vyhodnocení rutinních dat se předpokládá, že l l l l jde o nezávislé veličiny stejně rozdělené veličiny pocházející z normálního rozdělení Test nezávislosti prvků výběru l Závislost měření je obvykle způsobena nestabilitou měřicího zařízení, nekonstantností podmínek, zanedbáním faktorů, které významně ovlivňují výsledek měření Test homogenity výběru l K nehomogenitě naměřených dat dochází všude tam, kde se vyskytuje výrazná nestejnoměrnost měřených vlastností vzorků, mění se náhle podmínky experimentu a data obsahují vybočující měření Test normality dat Určení výběrového rozdělení a jeho konstrukce l Stanovení parametrů polohy a rozptýlení Průzkumová analýza dat pomocí grafických metod

Základní pojmy l Rozdělení pravděpodobnosti l Normální rozdělení Rozdělení je charakterizováno parametrem polohy a rozptýlení funkcí hustoty pravd. l Střední hodnotou m l Rozptylem s 2 l l Studentovo rozdělení l l Pro výběry pocházející z normálního rozdělení platí, že náhodná veličina t má Studentovo rozdělení s(n-1) stupni volnosti l C 2 - rozdělení l Pro výběry pocházející z normálního rozdělení platí, náhodná veličina C 2 má C 2 – rozdělení s (n-1) stupni volnosti

Základní pojmy l Výběrové charakteristiky polohy l Výběrový průměr l l l Je maximálně věrohodným odhadem střední hodnoty m Je velmi citlivý na odlehlé hodnoty Kvantily l Kvantily jsou přirozenou mírou polohy. Daný výběr se nejdříve seřadí od nejmenší hodnoty po největší a poté určíme pro daný p% kvantil pořadové číslo jednotky np , pro které platí: l Medián l l l padesátiprocentní kvantil, který charakterizuje z hlediska velikosti prostřední hodnotu souboru. Ve statistickém souboru odděluje polovinu hodnot menších od poloviny hodnot větších. Kvartily l robustní 25% a 75% kvantil Modus l robustní Modus je hodnota , která se vyskytuje nejčastěji. robustní

Základní pojmy l Výběrové charakteristiky rozptýlení l Výběrový rozptyl a směrodatná odchylka l je nevychýleným odhadem rozptylu s 2 pro normální rozdělení. Je definován jako aritmetický průměr ze čtverců odchylek jednotlivých hodnot od průměru, tj. , s – směrodatná odchylka l Varianční rozpětí l l Variační rozpětí R je rychlou, jednoduchou, ale jen orientační charakteristikou variability založenou na informaci o maximální a minimální hodnotě v souboru, tj. Interkvartilové rozpětí l robustní Kvantilový odhad směrodatné odchylky robustní

Základní pojmy l Charakteristika tvaru rozdělení l Šikmost l l slouží k zjištění míry symetrie (souměrnosti) rozdělení Jednoduchou charakteristikou šikmosti je Pearsonův koeficient α. K hodnocení stupně šikmosti využívá vztah mezi velikostí aritmetického průměru a mediánu daný ve tvaru Jiný přístup k měření šikmosti je založen na aplikaci tzv. momentových charakteristik. Při práci s daty uspořádanými do rozdělení četností je vhodná tzv. momentová míra šikmosti Špičatost l Míry špičatosti jsou založeny na srovnání stupně nahuštěnosti hodnot prostřední velikosti se stupněm nahuštěnosti ostatních

Základní pojmy l Intervalové odhady statistických parametrů l l l Koeficient spolehlivosti, hladina významnosti Ze statistického hlediska mají bodové odhady malý význam, protože neříkají nic o tom, kde leží skutečné hodnoty parametrů. Více informací poskytuje intervalový odhad, který určuje interval, v němž se bude se zadanou pravděpodobností ( 1 - a )nacházet skutečná hodnota daného parametru. Pravděpodobnost nazýváme koeficient spolehlivosti ( 1 - a ), volíme nejčastěji 95% nebo 99% a parametr a se nazývá hladina významnosti. l Intervalový odhad střední hodnoty l Intervalový odhad rozptylu l Intervalový odhad mediánu

Základní pojmy Grafické metody Diagram rozptýlení a krabicový graf

Základní pojmy Grafické metody Histogram Q-Q graf Odhad hustoty pravděpodobnosti Graf rozptýlení s kvantily

Průzkumová analýza dat Příklad: Homogenita tabletoviny l Obsahová stejnoměrnost l l Naměřené hodnoty l l Obsah účinné látky v tabletovině Výsledky analýzy obsahu účinné látky ve vzorcích odebraných z homogenizačního kontejneru při prospektivní validaci na 3 validačních šaržích Limity pro obsahovou stejnoměrnost v tabletovině : l l l Interval 93 -107 % pro jednotlivé naměřené hodnoty Interval 18, 6 – 21, 4 mg/tbl RSD ≤ 3, 2%

Průzkumová analýza dat Naměřené hodnoty Označení vzorku Obsah účinné látky v tabletovině [mg/tbl. ] Š 1 Š 2 Š 3 A 19, 70 19, 50 19, 60 B 19, 70 19, 51 19, 69 C 19, 60 19, 84 19, 47 D 19, 60 19, 66 18, 94 E 19, 60 19, 70 19, 52 F 19, 90 19, 47 19, 52 G 19, 80 19, 54 19, 52 H 19, 90 20, 10 19, 67 CH 19, 40 19, 69 19, 59 I 19, 80 19, 48 19, 45

Průzkumová analýza dat Test normality : Název sloupce : šarže 1 šarže 2 šarže 3 Normalita : Přijata Zamítnuta Vypočtený : 0, 573 4, 012 7, 690 Teoretický : 5, 991 Pravděpodobnost : 0, 751 0, 134 0, 021 Vybočující body : Název sloupce : šarže 1 šarže 2 šarže 3 Homogenita : Přijata Zamítnuta Počet vybočujících bodů : 0, 000 1, 000 Spodní mez : 19, 222 19, 083 19, 185 Horní mez : 20, 178 20, 087 19, 855 Autokorelace : Řád autokorelace : 4 Název sloupce : šarže 1 šarže 2 šarže 3 Řád autokorelace 1 Závěr : Nevýznamný

Průzkumová analýza dat Klasické parametry : Název sloupce : šarže 1 šarže 2 šarže 3 Průměr : 19, 7 19, 649 19, 497 Spodní mez : 19, 588 19, 506 19, 346 Horní mez : 19, 812 19, 792 19, 648 Rozptyl : 0, 024 0, 040 0, 045 Směr. odchylka : 0, 156 0, 200 0, 211 -0, 368 1, 195 -1, 974 Šikmost Odchylka od 0 : Nevýznamná Špičatost : Odchylka od 3 : Nevýznamná 2, 438 Nevýznamná Významná 3, 549 Nevýznamná 6, 134 Významná Polosuma 19, 650 19, 785 19, 315 Modus : 19, 700 19, 520 19, 558 Robustní parametry : Název sloupce : šarže 1 šarže 2 šarže 3 Medián : 19, 700 19, 600 19, 520 IS spodní : 19, 527 19, 392 19, 393 IS horní : 19, 873 19, 808 19, 647 Mediánová směr. odchylka : 0, 077 0, 092 0, 056 Mediánový rozptyl : 0, 006 0, 008 0, 003

Průzkumová analýza dat Grafické metody s 1 s 2 s 3

Průzkumová analýza dat

Vyhodnocování trendů l Vyhodnocování trendů zahrnuje metody a postupy sledování a vyhodnocování vybraných veličin procesu. Všechny tyto činnosti mají za cíl odhalit v získaných datech signály, že existují podmínky, které mohou způsobit odchylku od ustáleného stavu procesu a pomoci tak identifikovat nepříznivé vlivy, které mohou mít negativní dopad na produkt. Prakticky to znamená, že předpokládáme nějaké statistické vlastnosti procesu (obecně dat) a vyhodnocování trendů slouží (přinejmenším) k testování, zda se tyto statistické vlastnosti nemění. l Postupy: l Sledování středních hodnot a rozptylů l Grafické zobrazení dat l Regulační diagramy l

Vyhodnocování trendů Regulační diagramy l Diagram x-průměr l Použití tohoto grafu se doporučuje pokud je možné vytvořit tzv. podskupiny (dvě nebo více hodnot pro jeden časový okamžik). l Diagram x-individual l Používá se tam, kde z nějakého důvodu není možné vytvořit podskupiny, do grafu se vynáší přímo naměřené hodnoty l Diagram c l Je vhodný pro sledování počtu jednotlivých vad na definovaném množství výrobku, který má Poissonovo rozdělení

Vyhodnocování trendů Pravidla regulačních diagramů l l l l 1. Jedna hodnota mimo kontrolní meze, 2. Devět hodnot na téže straně centrální linie, 3. Šest hodnot za sebou roste nebo klesá, 4. Čtrnáct hodnot má střídavé znaménko svých diferencí, 5. Dvě ze tří hodnot jsou od základní linie dále než 2 s, 6. Čtyři z pěti hodnot jsou na téže straně základní linie dále než 1 s, 7. Patnáct hodnot je uvnitř intervalu ±s od základní linie, 8. Osm hodnot po sobě je mimo interval ±s od základní linie.

Vyhodnocování trendů Příklad: Sledování výtěžnosti l Výtěžek l l Množství přijaté na sklad z výrobní zakázky Naměřené hodnoty l Množství přijaté na sklad z výrobních zakázek 9 sledovaných přípravků za období 2003 -2004 byly získány ze systému SAP

Vyhodnocení trendů l Sledování středních hodnot a rozptylů l Z tabulky je patrné, že Přípravek 5 vykazuje za rok 2004 pokles výtěžnosti a změnu rozptylu. Proto provedeme analýzu pomocí regulačního diagramu X-individual. Nejdříve vytvoříme regulační diagram pro data za rok 2003.

Vyhodnocení trendů l Regulační diagram l Z grafu je patrné, že došlo zhruba po 100 šaržích ke změně ve výrobním postupu, proto je potřeba pro výpočet parametrů ustáleného stavu procesu vzít data pro šarže s pořadovým číslem větším než 100.

Vyhodnocení trendů l Při tvorbě regulačního diagramu byly vypočteny následující hodnoty parametrů: l ZL = 280, 318 (výběrový průměr) l l s = 2, 675 (výběrová směrodatná odchylka) LCL = 271, 937 UCL = 288, 700 Tyto vypočtené parametry byly použity pro regulační diagram z dat naměřených v roce 2004

Vyhodnocení trendů l l l Červeně označené body indikují porušení jednoho nebo více z 8 pravidel regulačních diagramů. Kromě bodů mimo meze (kde se jednalo o výrobu na výjimku), se také jedná o šarže kolem pořadového čísla 300, kde je pro více než dvacet bodů porušeno pravidlo č. 7 Patnáct hodnot je uvnitř intervalu ±s od základní linie Tento jev byl odůvodněn výrobou pod dohledem technologa z důvodů testování nového dodavatele suroviny.