Learn About Multinomial Logit in Stata With Data From the Behavioral Risk Factor Surveillance System (2013)
Výsledky multinomického logistického regresního modelu jsou uvedeny na obrázku 3.
Výstup logistického regresního modelu vypadá následovně:
Řádek 1: tečka mlogit aktivní 1 žena do 30 let věk 65 let plus artritida
Řádek 2: Iterace 0: logaritmická pravděpodobnost je záporná 379353,27
Řádek 3: Iterace 1: logaritmická pravděpodobnost je záporná 357272,01
Řádek 4: Iterace 2: logaritmická pravděpodobnost je záporná 356455,81
Řádek 5: Iterace 3: logaritmická pravděpodobnost je záporná 356439,91
Řádek 6: Iterace 4: logaritmická pravděpodobnost je záporná 356439,9
Řádek 7: Počet pozorování se rovná 359,925
Řádek 8: LR chi2 (8) se rovná 45826,74
Řádek 9: Pravděpodobnost větší než chi2 se rovná 0,0000
Řádek 10: Logaritmická pravděpodobnost se rovná záporné 356439,9; Pseudo R2 se rovná 0,0604
Pod touto informací se objevuje tabulka, v níž jsou různé hodnoty pro „aktivní 1“ uvedeny takto:
Řádek 1: Žádné podtržítko nebo podtržítko pod podtržítkem střední, žena: Koeficient, záporný 0,0137624; standardní chyba, 0,0081021; z, záporný 1,70; P větší než absolutní hodnota z, 0,089; 95% interval spolehlivosti, záporný 0,0296422, 0,0021173.
Řádek 2: Žádné podtržítko nebo podtržítko pod podtržítkem střední, mladší 30 let: Koeficient, záporný 0,2132395; standardní chyba, 0,0136507; z, záporný 15,62; P větší než absolutní hodnota z, 0,000; 95procentní interval spolehlivosti, 0,7998839, 0,8377535.
Řádek 3: Žádné podtržítko nebo podtržítko pod podtržítkem střední, věk 65 let a více: Koeficient, 0,8188187; standardní chyba, 0,0096608; z, 84,76; P větší než absolutní hodnota z, 0,000; 95 procentní interval spolehlivosti, 0,7998839, 0,8377535.
Řádek 4: Žádné podtržítko nebo podtržítko pod podtržítkem střední, artritida: Koeficient, 0,4168419; standardní chyba, 0,008863; z, 47,03; P větší než absolutní hodnota z, 0,000; 95procentní interval spolehlivosti, 0,3994706, 0,4342131.
Řádek 5: Žádné podtržítko nebo podtržítko pod podtržítkem mírné, podtržítko cons: Koeficient, záporný 0,7812104; standardní chyba, 0,007175; z, záporný 108,88; P větší než absolutní hodnota z, 0,000; 95% interval spolehlivosti, záporný 0,7952731, 0,7671477.
Řádek 6: Středně silný, základní výsledek
Řádek 7: Silný, žena: Koeficient, záporný 0,0456556; standardní chyba, 0,0090931; z, záporný 5,02; P větší než absolutní hodnota z, 0,000; 95% interval spolehlivosti, záporný 0,0634778, 0,0278334.
Řádek 8: Vigorní, mladší 30 let: Koeficient, záporný 1,333241; standardní chyba, 0,0247722; z, záporný 53,82; P větší než absolutní hodnota z, 0,000; 95procentní interval spolehlivosti, záporný 1,381794, záporný 1,284689.
Řádek 9: Vitální, věk 65 let a více: Koeficient, 1,663682; standardní chyba, 0,0099416; z, 167,35; P větší než absolutní hodnota z, 0,000; 95procentní interval spolehlivosti, 1,644197, 1,683167.
Řádek 10: Vigorní, artritida: Koeficient, 0,0555993; standardní chyba, 0,0098662; z, 5,64; P větší než absolutní hodnota z, 0,000; 95procentní interval spolehlivosti, 0,036262, 0,0749366.
Řádek 11: Vigorózní, podtržené zápory: Koeficient, záporný 1,168047; standardní chyba, 0,0081195; z, záporný 143,86; P větší než absolutní hodnota z, 0,000; 95procentní interval spolehlivosti, záporný 1,183961, záporný 1,152133.
Multinomiální logistická regrese odhaduje modely pro všechny kategorie závislé proměnné kromě jedné. Ve výchozím nastavení Stata vynechává nejčastější výsledek a používá jej jako základ, od kterého se odvíjí porovnání ostatních kategorií.
Ve výsledcích je mnoho informací. V tomto příkladu zaměříme svou pozornost na jednotlivé odhady koeficientů spojujících nezávislé proměnné se závislou proměnnou a jejich odpovídající úroveň statistické významnosti. Vidíme, že každý odhad koeficientu se statisticky významně liší od nuly. To by nás vedlo k zamítnutí nulové hypotézy o koeficientu rovném nule pro všechny odhady.
Interpretace výsledků multinomického logitového modelu vyžaduje více než jen zkoumání směru a úrovně statistické významnosti samotných odhadů koeficientů.
Nejběžnějším způsobem interpretace výsledků multinomického logitového modelu je výpočet předpokládaných pravděpodobností na základě výsledků analýzy. Protože předpokládaná pravděpodobnost zařazení do některé z kategorií závislé proměnné je nelineární funkcí nezávislých proměnných, vyžaduje výpočet předpokládaných pravděpodobností nastavení každé nezávislé proměnné na určitou hodnotu. Představte si to jako vytvoření popisného profilu pro případ v souboru dat a výpočet předpokládané pravděpodobnosti, že se někdo s tímto profilem zapojí do tří různých úrovní zde popsané činnosti. Změnou vlastností profilu a opětovným výpočtem těchto pravděpodobností můžete vyhodnotit, jak různé nezávislé proměnné ovlivňují změny v předpovídaných pravděpodobnostech.
Na příkladu vypočítáme předpovídané pravděpodobnosti zařazení do jednotlivých úrovní namáhavosti činnosti pro respondenty s artritidou a bez ní. Pro každý z nich nastavíme hodnotu proměnné indikátoru ženy na její střední hodnotu a obě proměnné věku na nulu, což znamená, že tyto profily popisují respondenty ve věku 30 až 64 let.
Pro výpočet předpokládané pravděpodobnosti zařazení do jednotlivých úrovní aktivity pro respondenty s artritidou zadejte do příkazového okna programu Stata následující příkaz:
margins, at(mean(female) under30=(0) age65plus=(0) arthritis=(1))
Stisknutím klávesy Enter získáte předpokládanou pravděpodobnost zařazení do jednotlivých skupin úrovní aktivity.
Pro osoby bez artritidy je kód následující:
margins, at(mean(female) under30=(0) age65plus=(0) arthritis=(0))
Stiskněte Enter pro vytvoření predikované pravděpodobnosti pro zařazení do každé skupiny úrovně aktivity.
Chcete-li místo toho vypočítat tyto předpovězené pravděpodobnosti pomocí možností nabídky, po spuštění modelu multinomiální logistické regrese vyberte v nabídce programu Stata následující možnosti:
Statistika → Postestimace
V otevřeném dialogovém okně „Výběr postestimace“, zobrazeném na obrázku 4, rozbalte položku „Marginální analýza“. Pod položkou „Marginální efekty všech kovariát“ zvýrazněte kliknutím na „Vlastní a vysokorozměrné analýzy“.
Stisknutím tlačítka Spustit otevřete dialogové okno „Margins – Marginal means, predictive margins, and marginal effects“.
Dialogové okno „Postestimation Selector“ se skládá ze seznamu příkazů pod názvem „Postestimation commands“. Položky se zobrazují následovně:
- Marginální analýza (rozšířené zobrazení)
- Marginální prostředky a marginální efekty, základní analýzy
- Marginální průměry a analýza interakcí (sbalené zobrazení)
- Kontrasty marginálních průměrů (sbalené zobrazení)
- Marginální účinky kovariáty při hodnotách jiné kovariáty (rozšířené zobrazení)
- –
Průměrná hodnota populace (vzorek průměrného nadhodnocení)
- -.
Při výběrových průměrech
- –
- Marginální účinky všech kovariát (rozšířené zobrazení)
- –
Vlastní a high-dimenzionální analýzy (zvýrazněné)
- –
Profilové grafy po marginální analýze
- –
- Testy, kontrasty a porovnání odhadů parametrů (sbalené zobrazení)
- Specifikace, diagnostika, a analýza dobré shody (sbalené zobrazení)
- Předpovědi (sbalené zobrazení)
- Další sestavy (sbalené zobrazení)
- Správa výsledků odhadů (sbalené zobrazení)
V pravém horním rohu okna se zobrazí tlačítko „Spustit“ a v pravém dolním rohu tlačítko „Zrušit“.
V horní části dialogového okna „Margins – Marginal means, predictive margins, and marginal effects“ vyberte záložku označenou „At“. Na kartě „At“ by mělo být zvýrazněno tlačítko vedle položky „Všechny kovariáty při pozorovaných hodnotách ve vzorku“.
Klikněte na tlačítko „Vytvořit“ dole a otevře se nové dialogové okno „Specifikace 1“.
V sekci „Statistika“ v horní části dialogového okna „Specifikace 1“ lze nastavit příslušné proměnné na jejich střední hodnoty, v tomto případě ženy. Zaškrtněte „1“ kliknutím na políčko vlevo od něj a v textovém poli Statistika vpravo se objeví šipka s rozbalovací nabídkou. Vyberte možnost „Means“. V poli „Covariates“ (Průběžné proměnné) vyberte pomocí rozevírací nabídky ženu.
V části „Fixed values“ (Pevné hodnoty) můžete nastavit hodnoty pro zbývající proměnné. Klikněte na políčko vedle čísla 1. V rozevíracích možnostech „Covariate“ vyberte ze seznamu proměnných hodnotu under30. Do textového pole „Numlist“ zadejte ručně hodnotu proměnné na „0“. V poli 2 nastavte hodnotu age65plus na „0“. Nakonec v rámečku 3 nastavte hodnotu artritidy na „1“.
Obrázek 5 ukazuje, jak to vypadá v programu Stata.
Dialogové okno s názvem „Specifikace 1“ lze rozdělit na tři části nazvané jako „Statistika“, „Pevné hodnoty“ a „Výrazy“. Horní část s názvem „Statistika“ je dále rozdělena na dva panely, levý a pravý, které se zobrazují následovně.
Levý panel, Statistika:
- (zaškrtávací políčko) 1: text „Means,“ (rozbalovací tlačítko)
- (prázdné zaškrtávací políčko) 2: prázdný textový řádek, 50 (rozbalovací tlačítko), v režimu stupňů šedi
- (prázdné zaškrtávací políčko) 3: text „Means,“ (rozbalovací tlačítko): (rozbalovací tlačítko), v režimu stupňů šedi
- (prázdné zaškrtávací políčko) 4: prázdný textový řádek, 50 (rozbalovací tlačítko), v režimu stupňů šedi
Pravý panel, kovariáty:
- Text „žena“, rozbalovací tlačítko
- Prázdné textové pole, rozbalovací tlačítko
- Prázdné textové pole, rozbalovací tlačítko
- Prázdné textové pole, rozbalovací tlačítko
V poznámce dole je uvedeno: „Kovariáty mohou být také podtržítko všechny (všechny kovariáty) tečka podtržítko faktorové (všechny faktorové kovariáty) tečka a podtržítko spojité (všechny spojité kovariáty).“
Střední část nazvaná „Pevné hodnoty“ je dále rozdělena na dva panely, levý a pravý, které vypadají takto:
Levý panel, Kovariáty:
- (zaškrtávací políčko) 1: mladší 30 let (rozbalovací tlačítko)
- (zaškrtávací políčko) 2: věk 65 let a více (rozbalovací tlačítko)
- (zaškrtávací políčko) 3: věk 65 let a více (rozbalovací tlačítko)
- (zaškrtávací políčko) 4: věk 65 let a více (rozbalovací tlačítko) (rozbalovací tlačítko)
- (prázdné zaškrtávací políčko) 4: prázdný textový řádek (rozbalovací tlačítko)
Pravý panel, Numlist:
- 0
- 0
- 1 (vybraný textový řádek)
- Prázdný textový řádek
Dolní část s názvem „Výrazy“ je dále rozdělena na dva panely, levý a pravý, které se zobrazují následovně.
Levý panel, kovariát:
- (prázdné zaškrtávací políčko) 1: prázdný textový řádek (rozbalovací tlačítko)
- (prázdné zaškrtávací políčko) 2: prázdný textový řádek (rozbalovací tlačítko)
- (prázdné zaškrtávací políčko) 3: (rozbalovací tlačítko)
- (prázdné zaškrtávací políčko) 4: prázdný textový řádek (rozbalovací tlačítko)
Pravý panel, Expression:
- Prázdný textový řádek, s tlačítkem „Vytvořit“ v režimu stupňů šedi
- Prázdný textový řádek, s tlačítkem „Vytvořit“ v režimu stupňů šedi
- Prázdný textový řádek, s tlačítkem „Vytvořit“ v režimu stupňů šedi
- Prázdný textový řádek, s tlačítkem „Vytvořit“ v režimu stupňů šedi
V pravém dolním rohu dialogového okna se zobrazí dvě tlačítka, „OK“ a „Zrušit“, s tlačítkem OK ve vybraném režimu.
Klepnutím na OK se vrátíte do dialogového okna „margins – Marginal means, predictive margins, and marginal effects“. To by nyní mělo v okně „At specifications (modify settings from above)“ zobrazit „Specification 1“ a pod ním popis podrobností specifikace.
Znovu klikněte na tlačítko „Create“ a otevře se další dialogové okno „Specification 2“. Použijte možnosti jako dříve, tentokrát nastavte hodnotu artritidy na 0.
Znovu klikněte na tlačítko OK a vraťte se do dialogového okna „Margins – Marginal means, predictive margins, and marginal effects“. Nyní se v okně „At specifications (modify settings from above)“ zobrazí „Specification 1“ a „Specification 2“.
Obrázek 6 ukazuje, jak to vypadá ve Stata.
Dialogové okno má název „margins – Marginal means, predictive margins, and marginal effects“. Bezprostředně pod záhlavím se objeví sada 10 záložek jako: „Hlavní“, „Při“, „pokud nebo v nebo nad“, „V rámci“, „Kontrast“, „Párová srovnání“, „Váhy“, „SE“, „Pokročilé“ a „Hlášení“, přičemž „Při“ je ve zvoleném režimu. Obsah dialogového okna se dále zobrazuje následovně:
- (Vybrané přepínač) Všechny kovariáty při pozorovaných hodnotách ve vzorku
- (Prázdný přepínač) Všechny kovariáty při jejich středních hodnotách ve vzorku
- (Prázdný přepínač) Pro každou další faktorovou proměnnou, zacházet se všemi úrovněmi, jako by byly stejně pravděpodobné (ruší obě výše uvedená nastavení pro faktorové proměnné)
Dále se zobrazí seznamové pole s názvem „Při specifikacích (upravit nastavení z výše uvedeného)“ obsahující následující 2 položky:
- Specifikace 1
- Specifikace 2 (vybraná)
V pravé části pole je sada čtyř záložek jako: „Vytvořit“, „Upravit“, „Zakázat“ a „Povolit“, přičemž „Vytvořit“, „Upravit“ a „Zakázat“ jsou v aktivním režimu a „Povolit“ v režimu stupňů šedi. Text v dolní části zní: „u levého rodiče levý rodič průměrný pravý rodič žena do 30 let rovná se levý rodič 0 pravý rodič věk 65 let plus rovná se levý rodič 0 pravý rodič artritida rovná se levý rodič 0 pravý rodič pravý rodič.“
V pravém dolním rohu dialogového okna se zobrazí tři tlačítka: „OK“, „Storno“ a „Odeslat“, přičemž tlačítko OK je ve zvoleném režimu.
Znovu klikněte na tlačítko OK a zobrazí se výsledky, které jsou na obrázku 7.
Výsledky získané z předpokládané pravděpodobnosti jednotlivých skupin namáhavosti pohybové aktivity v posledních 30 dnech pro respondenty vypadají následovně:
Řádek 1. Na řádku 2 se zobrazuje pravděpodobnost, že respondenti budou vykonávat pohybovou aktivitu v posledních 30 dnech: Upravené předpovědi; počet pozorování se rovná 359 925
Řádek 2: Model VCE dvojtečka OIM
Řádek 3: 1 tečka podtržítko predikce dvojtečka Pr levý rodič aktivní 1 se rovná Žádný podtržítko nebo podtržítko Pod podtržítkem mírný pravý rodič, predikce levý rodič pr výsledek levý rodič 0 pravý rodič pravý rodič.
Line 4:
Řádek 5: 2 tečky podtržítko predikce dvojtečka Pr levá-rodič aktivní 1 rovná se rovná se mírná pravá-rodič, predikce levá-rodič pr výsledek levá-rodič 1 pravá-rodič pravá-rodič.
Řádek 5:
Řádek 6: 3 tečky podtržítko predikce dvojtečka Pr levá-parentéza aktivní 1 rovná se Silná pravá-parentéza, predikce levá-parentéza pr výsledek levá-parentéza 2 pravá-parentéza pravá-parentéza: 1 tečka podtržítka u dvojtečky žena se rovná 0,5679239 (průměr)
Řádek 7: do 30 let se rovná 0
Řádek 8: věk 65 let plus se rovná 0
Řádek 9: artritida 1 se rovná 1
Řádek 10: 2 tečky podtržítka u dvojtečky žena se rovná 0.5679239 (průměr)
Řádek 11: do 30 let se rovná 0
Řádek 12: věk 65 let plus se rovná 0
Řádek 13: aktivní 1 se rovná 0
Tabulka dále uvádí hodnoty „rozpětí“, „standardní chyba delta-metody“, „z“, „pravděpodobnost větší než z“ a 95procentní interval spolehlivosti“. Jednotlivé hodnoty predikované pravděpodobnosti jednotlivých skupin namáhavosti pohybové aktivity za posledních 30 dní pro respondenty jsou následující:
Řádek 1: Podskupina predikuje 1 1: Marže, 0,3429731; Delta-metoda Standardní chyba, 0,0016608; z, 206.51; Pravděpodobnost větší než z, 0,000; 95procentní interval spolehlivosti, 0,339718, 0,3462283.
Řádek 2: Underscore predict 1 underscore at 2: Margin, 0,258514; Delta-method Standard Error, 0,0010418; z, 248,14; Pravděpodobnost větší než z, 0.000; 95procentní interval spolehlivosti, 0,256472, 0,2605559.
Řádek 3: Podtržítko předpovídá 2 podtržítka na 1: rozpětí, 0,49762; standardní chyba Delta-metody, 0,0017772; z, 280,00; pravděpodobnost větší než z, 0,000; 95procentní interval spolehlivosti, 0.4941367, 0,5011032.
Řádek 4: Podtržítko předpovídá 2 podtržítka při 2: rozpětí, 0,5690545; standardní chyba delta-metody, 0,0011953; z, 476,06; pravděpodobnost větší než z, 0,000; 95 procentní interval spolehlivosti, 0,5667117, 0,57133973.
Řádek 5: podtržítko předpovídá 3 podtržítka na 1: rozpětí, 0,1594069; standardní chyba delta-metody, 0,00117; z, 136,25; pravděpodobnost větší než z, 0,000; 95procentní interval spolehlivosti, 0,1571138, 0,1617.
Řádek 6: podtržítko předpovídá 3 podtržítka u 2: rozpětí, 0,1724315; standardní chyba delta-metody, 0,0008886; z, 194,05; pravděpodobnost větší než z, 0,000; 95procentní interval spolehlivosti, 0,17068991, 0,1741731.
V levém sloupci tabulky dole jsou pod slovy „_predict#_at“ dva sloupce čísel. Čísla 1-3 vlevo jsou jednotlivé kategorie pro různé úrovně namáhavosti činnosti. Čísla 1 a 2 vpravo jsou dva různé profily, které jsme vytvořili pro to, zda respondent má či nemá artritidu. Podrobnosti o těchto profilech jsou uvedeny v informacích nad tabulkou.
Předpokládané pravděpodobnosti a jejich intervaly spolehlivosti jsou odhadnuty pomocí postestimulační simulace. Úplné pojednání o tomto procesu přesahuje rámec tohoto příkladu, ale stručně lze říci, že tento proces vypočítá 1 000 sad předpovídaných pravděpodobností pomocí simulace hodnot koeficientů modelu na základě jejich odhadovaných hodnot, rozptylů a kovariancí. Další informace naleznete v části „Využití statistických analýz:
Nejlepším způsobem, jak prozkoumat změny v předpovídaných pravděpodobnostech zahrnující nezávislou proměnnou s více než dvěma hodnotami, je vytvořit obrázek, jaký je znázorněn na obr. 8. Obrázek 8 znázorňuje predikovanou pravděpodobnost zařazení do jedné ze tří úrovní aktivity pro respondenty v každé z těchto tří věkových skupin při zachování ostatních dvou nezávislých proměnných v modelu na jejich příslušných středních hodnotách. Výška každého sloupce vystihuje velikost předpovídané pravděpodobnosti.
První sloupcový graf má název „do 30 let rovná se Ne, věk 65 plus rovná se Ne“. Svislá osa označuje pravděpodobnost v rozmezí od 0 do 0,6 a vodorovná osa označuje kategorizaci BMI v rozmezí od 0 do 2.
Odhadovaná predikovaná pravděpodobnost respondentů při různých klasifikacích BMI je následující: (0, 0,28); (1, 0,55; (2, 0,18).
Druhý sloupcový graf má název „do 30 let rovná se Ne, věk 65 plus rovná se Ano“. Odhadované předpokládané pravděpodobnosti respondentů při různých klasifikacích BMI jsou následující: (0, 0,30); (1, 0,25; (2, 0,41).
Třetí sloupcový graf má název „do 30 let rovná se Ano, věk 65 plus rovná se Ne“. Odhadované předpokládané pravděpodobnosti respondentů při různých klasifikacích BMI jsou následující: (0, 0,28); (1, 0,66); (2, 0,05).
Pro vytvoření tohoto grafu se vraťte do dialogového okna „Postestimation Selector“, které by mělo být stále otevřené. Ujistěte se, že je stále zvýrazněna položka „Vlastní a vysokorozměrné analýzy“.
Stisknutím tlačítka Spustit znovu otevřete dialogové okno „Margins – Marginal means, predictive margins, and marginal effects“.
Znovu vyberte kartu „At“ jako dříve. V poli „At specifikace (upravit nastavení z předchozího)“ se zobrazí obě předchozí specifikace. Označte každou z nich postupně kliknutím na ni a zakažte ji pomocí tlačítka „Zakázat“ vpravo. Kliknutím na tlačítko „Vytvořit“ otevřete nové dialogové okno „Specifikace“.
Stejně jako dříve jej použijeme k zadání hodnot proměnných pro zvolený profil. V části „Statistika“ v horní části dialogového okna „Specifikace 1“ lze nastavit příslušné proměnné (v tomto případě žena a artritida) na jejich střední hodnoty. Zaškrtněte „1“ kliknutím na políčko vlevo od něj a v textovém poli Statistika vpravo se zobrazí šipka s rozbalovací nabídkou. Vyberte možnost „Means“. V poli „Covariates“ (Průběžné proměnné) vyberte pomocí rozevíracích možností ženu. Totéž proveďte pro artritidu v řádku níže.
Dále použijte oddíl „Fixní hodnoty“ jako předtím a vyberte hodnoty pro věk do 30 let a věk65plus. Protože sledujeme změnu předpokládaných pravděpodobností ve třech věkových kategoriích, musíme vytvořit tři specifikace pro každou z věkových skupin s věkem nastaveným na následující hodnoty:
- Mladší než 30 let-do 30 let 1 age65plus 0
- 30-64 let-do 30 let 0 age65plus 0
- 65 let a starší než 30 let 0 age65plus 1
První z nich je znázorněna na obrázku 9. Na obrázku 9 je znázorněn věk a věk.
Jakmile se všechny tři nové specifikace zobrazí v dialogovém okně „Margins – Marginal means, predictive margins, and marginal effects“, klikněte na tlačítko OK, aby se vytvořily předpovězené pravděpodobnosti.
Chcete-li zadat příkaz přímo do příkazového okna programu Stata, kód je následující:
margins, at((mean) female (mean) arthritis under30=(1) age65plus=(0)) at((mean) female (mean) arthritis under30=(0) age65plus=(0)) at((mean) female (mean) arthritis under30=(0) age65plus=(1))
Stiskněte Enter pro vytvoření predikovaných pravděpodobností pro vybrané profily.
Dialogové okno s názvem „Specifikace 3“ lze rozdělit na tři části nazvané jako „Statistiky“, „Pevné hodnoty“ a „Výrazy“. Horní oddíl s názvem „Statistika“ je dále rozdělen na dva panely, levý a pravý, které se zobrazují následovně:
Levý panel, Statistika:
- (zaškrtnuté políčko) 1. Zaškrtněte políčko „Statistika“: Text „Means,“ (rozbalovací tlačítko)
- (zaškrtávací políčko) 2: Text „Means,“ (rozbalovací tlačítko)
- (prázdné zaškrtávací políčko) 3: prázdný textový řádek, 50 (rozbalovací tlačítko), v režimu stupňů šedi
- (prázdné zaškrtávací políčko) 4: Text „Means,“ (rozbalovací tlačítko): prázdný textový řádek, 50 (rozbalovací tlačítko), v režimu stupňů šedi
Pravý panel, kovariáty:
- Text „žena,“ rozbalovací tlačítko
- Text „artritida,“ rozbalovací tlačítko
- Prázdné textové pole, rozbalovací tlačítko
- Prázdné textové pole, rozevírací tlačítko
Poznámka dole zní: „Kovariáty mohou být také podtržítko všechny (všechny kovariáty) tečka podtržítko faktorové (všechny faktorové kovariáty) tečka a podtržítko spojité (všechny spojité kovariáty).“
Prostřední část nazvaná „Fixní hodnoty“ je dále rozdělena na dva panely, levý a pravý, které vypadají následovně:
Levý panel, kovariáty:
- (zaškrtnuté políčko) 1: do 30 let (rozbalovací tlačítko)
- (zaškrtávací políčko) 2: věk 65 let a více (rozbalovací tlačítko)
- (prázdné zaškrtávací políčko) 3: prázdný textový řádek (rozbalovací tlačítko)
- (prázdné zaškrtávací políčko) 4: (rozbalovací tlačítko)
Pravý panel, Numlist:
- 1
- 0 (vybraný textový řádek)
- Prázdný textový řádek
- Prázdný textový řádek
Dolní část s názvem „Výrazy“ je dále rozdělena na dva panely, levý a pravý, které se zobrazují takto:
Levý panel, Kovariata:
- (prázdné zaškrtávací pole) 1: prázdný textový řádek (rozbalovací tlačítko)
- (prázdné zaškrtávací pole) 2: prázdný textový řádek (rozbalovací tlačítko)
- (prázdné zaškrtávací pole) 3: prázdný textový řádek (rozbalovací tlačítko)
- (prázdné zaškrtávací pole) 4: prázdný textový řádek (rozbalovací tlačítko)
Pravý panel, Výraz:
- Prázdný textový řádek, s tlačítkem „Vytvořit“ v režimu stupňů šedi
- Prázdný textový řádek, s tlačítkem „Vytvořit“ v režimu stupňů šedi
- Prázdný textový řádek, s tlačítkem „Vytvořit“ v režimu stupňů šedi
- Prázdný textový řádek, s tlačítkem „Vytvořit“ v režimu stupňů šedi
V pravém dolním rohu dialogového okna se zobrazí dvě tlačítka, „OK“ a „Zrušit“, s tlačítkem OK ve vybraném režimu.
Pro vykreslení výsledků se vraťte do dialogového okna Postestimation Selector. Vyberte možnost „Profilové grafy po marginální analýze“, jak je znázorněno na obrázku 10.
Dialogové okno s názvem „Postestimation Selector“ se skládá ze seznamu příkazů pod názvem „Postestimation commands“. Položky se zobrazují následovně:
- Marginální analýza (rozšířené zobrazení)
- Marginální prostředky a marginální efekty, základní analýzy
- Marginální průměry a analýza interakcí (sbalené zobrazení)
- Kontrasty marginálních průměrů (sbalené zobrazení)
- Marginální účinky kovariáty při hodnotách jiné kovariáty (rozšířené zobrazení)
- –
Průměrná hodnota populace (vzorek průměrného nadhodnocení)
- -.
Při výběrových průměrech
- –
- Marginální účinky všech kovariát (rozšířené zobrazení)
- –
Vlastní a high-dimenzionální analýzy
- –
Grafy profilů po marginální analýze (zvýrazněné)
- –
- Testy, kontrasty a porovnání odhadů parametrů (sbalené zobrazení)
- Specifikace, diagnostika, a analýza dobré shody (sbalené zobrazení)
- Předpovědi (sbalené zobrazení)
- Další sestavy (sbalené zobrazení)
- Správa výsledků odhadů (sbalené zobrazení)
V pravém horním rohu okna se zobrazí tlačítko „Spustit“ a v pravém dolním rohu tlačítko „Zrušit“.
Stisknutím tlačítka Launch otevřete dialogové okno „marginsplot – Graph results from margins“. Na kartě „Main“ (Hlavní) ve třetím textovém poli označeném „Create subgraphs for groups defined by variables:“ (Vytvořit dílčí grafy pro skupiny definované proměnnými) vyberte z rozevíracích možností dvě proměnné věku, jak je znázorněno na obrázku 11.
Klikněte na tlačítko „Options“ (Možnosti) po straně textového pole. Tím se otevře nové dialogové okno „Podle popisku rozměru“. Zaškrtněte „Zahrnout název proměnné a znaménka rovnosti do všech popisků“, jak je znázorněno na obrázku 12.
Dialogové okno je nadepsáno „marginsplot – Graf výsledků z okrajů“. Hned pod titulkovým pruhem se objeví panel nabídek s 12 kartami jako: „Hlavní“, „Štítky“, „Graf“, „CI graf“, „Párový“, „Přidat grafy“, „Osa Y“, „Osa X“, „Nadpisy“, „Legenda“, „Celkově“ a „Podle možností“, přičemž ve vybraném režimu je karta „Hlavní“.
Další obsah dialogového okna se zobrazuje pod názvem „Rozměry“ takto:
- Proměnné, které definují osu x, s prázdným textovým pruhem a rozbalovací nabídkou a tlačítkem „Možnosti“ vedle ní.
- Vytvoření grafů pro skupiny definované proměnnými, s prázdným textovým pruhem a rozbalovací nabídkou a tlačítkem „Možnosti“ vedle ní.
- Vytvoření dílčích grafů pro skupiny definované proměnnými s textovým pruhem ve vybraném režimu, který zní: „ve věku (do 30 let) ve věku (nad 65 let)“. Textový pruh má také rozevírací nabídku a vedle ní tlačítko „Možnosti“.
- Vytvořte grafy pro skupiny definované proměnnými, s prázdným textovým pruhem a rozevírací nabídkou a vedle ní tlačítko „Možnosti“.
Pod ním se zobrazí následující kontrolní seznam:
- (prázdné zaškrtávací políčko) Prohodit osy x a y (vodorovný graf)
- (prázdné zaškrtávací políčko) Nevykreslovat intervaly spolehlivosti
- (prázdné zaškrtávací políčko) Název grafu, nebo stub v případě více grafů, s prázdným textovým pruhem pod ním
- (prázdné zaškrtávací políčko) Nahradit graf, pokud již existuje v paměti, v režimu stupňů šedi
V pravém dolním rohu dialogového okna se zobrazí tři tlačítka: „OK“, „Zrušit“ a „Odeslat“, přičemž tlačítko OK je ve zvoleném režimu.
Dialogové okno má název „Podle označení dimenze“. Další obsah dialogového okna se zobrazí následovně:
Řádek 1:
Další obsah dialogového okna se zobrazí následovně: Poznámka: Tato nastavení budou mít přednost před výchozím nastavením a nastavením karty Štítek
Řádek 2: Seznam uvozovkových řetězců pro označení jednotlivých úrovní dimenze
Řádek 3: (prázdný textový řádek ve vybraném režimu)
Řádek 4: Seznam vyjmenovaných označených, příklad 1 „jedna“ 2 „dvě“ elipsa
Řádek 5: (prázdný textový řádek) Označte skupiny jejich hodnotami, nikoli jejich popisky
Řádek 6: (prázdné zaškrtávací pole) Vynechat název proměnné a znaménka rovnosti ve všech popiscích, v režimu stupňů šedi
Řádek 7: (zaškrtnuté prázdné políčko) Zahrnout název proměnné a znaménka rovnosti do všech štítků
Řádek 8: (zaškrtnuté prázdné políčko) Oddělovač štítků při zadání více proměnných
Linka 9: (prázdný textový řádek) v režimu stupňů šedi
Řádek 10: (prázdné zaškrtávací políčko) Nepoužívat oddělovač mezi štítky
V pravém dolním rohu dialogového okna se zobrazí tři tlačítka: „OK“, „Zrušit“ a „Odeslat“, přičemž tlačítko OK je ve zvoleném režimu.
Na kartě „Plot“, zobrazené na obrázku 13, vyberte v horní části z rozbalovacích možností „Plot type“ možnost „Bar“.
Dialogové okno má název „marginsplot – Graph results from margins“. Hned pod titulkovým pruhem se objeví lišta nabídek s 12 kartami jako: „Hlavní“, „Popisky“, „Graf“, „CI graf“, „Párový“, „Přidat grafy“, „Osa Y“, „Osa X“, „Nadpisy“, „Legenda“, „Celkově“ a „Podle možností“, přičemž ve vybraném režimu je karta „Graf“.
Pod lištou nabídek se objeví rozbalovací nabídka s názvem „Typ grafu“, přičemž z nabídky je vybrána možnost „Sloupec“. Dále se zobrazí seznamové pole s názvem „Vybrat plot“ s následujícím seznamem položek:
- Všechny parcely (vybrané)
- Parcela 1
- Parcela 2
- Parcela 3
- Parcela 4
- Parcela 4 5
- Plot 6
- Plot 7
- Plot 8
Souběžně se seznamem se objeví záložka „Vlastnosti sloupce“. V dolní části je text: „Zde zobrazené parcely nemusí nutně odrážet počet skutečných parcel grafu.“
V pravém dolním rohu dialogového okna se objeví tři tlačítka: „OK“, „Storno“ a „Odeslat“, přičemž tlačítko OK je ve vybraném režimu.
Stata poskytuje další možnosti, jak graf lépe prezentovat. V tomto případě vybereme jen jednu další možnost. Vybereme kartu „Osa X“ a do textového pole „Název“ napíšeme „Úroveň aktivity“, čímž označíme osu x, jak je znázorněno na obrázku 14.
Klepnutím na tlačítko OK vytvoříme graf předpokládaných pravděpodobností zařazení do jednotlivých kategorií úrovně aktivity pro rozsah vybraných hodnot, jak je znázorněno na obrázku 8.
Pro zadání příkazu přímo do příkazového okna Stata je kód následující:
marginsplot, bydimension(at(under30) at(age65plus), nosimplelabels) recast(bar) recastci(rbar) xtitle(BMI category)
Stiskněte Enter pro vytvoření sloupcového grafu.
Dialogové okno má název „marginsplot – Graf výsledků z okrajů“. Bezprostředně pod titulkovým pruhem se objeví lišta nabídek s 12 kartami jako: „Hlavní“, „Popisky“, „Graf“, „CI graf“, „Párový“, „Přidat grafy“, „Osa Y“, „Osa X“, „Nadpisy“, „Legenda“, „Celkově“ a „Podle možností“, přičemž ve vybraném režimu je karta „Osa X“.
V textovém pruhu pod názvem „Nadpis“ se objeví položka „Úroveň aktivity“ a vedle ní karta „Vlastnosti“. Pod textovým pruhem se objevuje sada dalších pěti karet jako: „Vlastnosti hlavního zaškrtávacího políčka/nápisu“, „Vlastnosti vedlejšího zaškrtávacího políčka/nápisu“, „Vlastnosti měřítka osy“ a „Referenční čáry“. Následují dvě možnosti se zaškrtávacími políčky takto:
- (prázdné zaškrtávací políčko) Skrýt osu
- (prázdné zaškrtávací políčko) Umístit osu na opačnou stranu grafu
V pravém dolním rohu dialogového okna se objeví tři tlačítka: „OK“, „Storno“ a „Odeslat“, přičemž tlačítko OK je ve vybraném režimu.
Kompletní interpretace výsledků multinomického logitového modelu by představovala podobné tabulky nebo obrázky pro každou nezávislou proměnnou v modelu.
.