Viacnásobná lineárna regresia v programovom prostredí R
 
 
Step: Jednoduchá logistická regresia...
  Jednoduchá logistická regresia
Potom príkaz bude vyzerať nasledovne:
summary(glm(y~x,family=binomial, data=nazov_objektu))
Vysvetlivky k príkazom

y - závislá binárna premenná
x - nezávislá premenná
nazov_objektu - názov objektu/ databázy s údajmi
Príklad
Výskumník zisťoval, či má vplyv skóre prestíže zamestnania írskych otcov (fathocc) na ukončenie vzdelania ich detí s certifikátom (lvcert=1) alebo bez certifikátu (lvcert=0). Keďže premenná vyjadrujúca ukončenie vzdelania dieťaťa s alebo bez certifikátu je binárna, výskumník využil logistickú regresiu na bližšiu charakteristiku vzťahu medzi prestížou zamestnania otca a vzdelaním jeho dieťaťa. Výstupom príkazu
summary(glm(formula=lvcert~fathocc,
family=binomial,data=irished))

je regresný koeficient pre skóre prestíže zamestnania otca a p-hodnota.

H0: Prestíž zamestnania otca nemá štatisticky významný vplyv na skutočnosť, či dosiahne dieťa vzdelanie s certifikátom.
HA: Prestíž zamestnania otca má štatisticky významný vplyv na skutočnosť, či dosiahne dieťa vzdelanie s certifikátom.

Z výstupu logistickej regresie usúdil, že vplyv prestíže zamestnania otcov na ukončenie vzdelania ich detí je štatisticky významný (na hladine významnosti 95% p< 0,05).

Regresný koeficient kvôli lepšej interpretovateľnosti prepočítal na OR.
 
 
Step: Prepočet regresného koeficienta na OR pri jednoduchej logis...
  Prepočet regresného koeficienta na OR pri jednoduchej logistickej regresii, výpočet 95% CI
Potom príkaz bude vyzerať nasledovne:
# zistenie regresných koeficientov

sum.coef<-summary(glm (formula= y~x,family=binomial,data= nazov_objektu))$coef

# prepočet regresných koeficientov na OR

est<-exp(sum.coef[,1])

# výpočet horného konfidenčného intervalu

upper.ci<-exp(sum.coef[,1]+1.96*sum.coef[,2])

# výpočet dolného konfidenčného intervalu

lower.ci<-exp(sum.coef[,1]-1.96*sum.coef[,2])

# sumár (regresný koeficient, štandardná chyba, hodnota testovacej štatistiky z, p-hodnota, OR, dolný CI, horný CI)

cbind(sum.coef,est,lower.ci,upper.ci)

Vysvetlivky k príkazom

y - závislá binárna premenná
x - nezávislá premenná
Príklad
S využitím nasledujúcich príkazov výskumník prepočítal regresný koeficient pre skóre prestíže zamestnania otca na OR a vypočítal 95% konfidenčný interval pre OR.
sum.coef<-summary(glm (formula=lvcert~ fathocc,family=binomial,data=irished))$coef
est<-exp(sum.coef[,1])
upper.ci<-exp(sum.coef[,1]+1.96*sum.coef[,2])
lower.ci<-exp(sum.coef[,1]-1.96*sum.coef[,2])
cbind(sum.coef,est,lower.ci,upper.ci)

Keď sa zvýši skóre prestíže zamestnania otca o jednu jednotku, zvýši sa šanca dieťaťa dosiahnuť vzdelanie s certifikátom 1,035-násobne (95%CI: 1,023-1,047).