로지스틱 분포

R 중급 회귀

Richie Cotton

Data Evangelist at DataCamp

가우시안 확률밀도함수(PDF)

gaussian_distn <- tibble(
  x = seq(-4, 4, 0.05),
  gauss_pdf_x = dnorm(x)
)

ggplot(gaussian_distn, aes(x, gauss_pdf_x)) +
  geom_line()

선형-가우시안-pdf.png

가우시안 누적분포함수(CDF)

gaussian_distn <- tibble(
  x = seq(-4, 4, 0.05),
  gauss_pdf_x = dnorm(x),
  gauss_cdf_x = pnorm(x)
)

ggplot(gaussian_distn, aes(x, gauss_cdf_x)) +
  geom_line()

선형-가우시안-cdf.png

가우시안 역 CDF

gaussian_distn_inv <- tibble(
  p = seq(0.001, 0.999, 0.001),
  gauss_inv_cdf_p = qnorm(p)
)

ggplot(gaussian_distn_inv, aes(p, gauss_inv_cdf_p)) +
  geom_line()

선형-가우시안-icdf.png

분포 함수 이름

곡선	접두사	정규분포	로지스틱	암기법
PDF	d	`dnorm()`	`dlogis()`	"d"는 미분. CDF를 미분하면 PDF
CDF	p	`pnorm()`	`plogis()`	"p"는 q의 반대 → 역 CDF의 역함수
역 CDF	q	`qnorm()`	`qlogis()`	"q"는 분위수

glm()의 family 인수

lm(response ~ explanatory, data = dataset)

glm(response ~ explanatory, data = dataset, family = gaussian)

glm(response ~ explanatory, data = dataset, family = binomial)

¹ https://campus.datacamp.com/courses/introduction-to-regression-in-r/simple-logistic-regression?ex=1

gaussian()

str(gaussian())

List of 11
 $ family    : chr "gaussian"
 $ link      : chr "identity"
 $ linkfun   :function (mu)  
 $ linkinv   :function (eta)  
 $ variance  :function (mu)  
 $ dev.resids:function (y, mu, wt)  
 $ aic       :function (y, n, mu, wt, dev)  
 $ mu.eta    :function (eta)  
 $ initialize:  expression({  n <- rep.int(1, nobs)  if (is.null(etastart) && is.null(start) &&
     is.null(mustart) &&  ((family$link| __truncated__
 $ validmu   :function (mu)  
 $ valideta  :function (eta)  
 - attr(*, "class")= chr "family"

linkfun과 linkinv

링크 함수(link function)는 반응변수를 변환합니다.

gaussian()$linkfun

function (mu) 
mu

gaussian()$linkinv

function (eta) 
eta

로지스틱 PDF

logistic_distn <- tibble(
  x = seq(-6, 6, 0.05),
  logistic_pdf_x = dlogis(x)
)

ggplot(logistic_distn, aes(x, logistic_pdf_x)) +
  geom_line()

선형-로지스틱-pdf.png

로지스틱 분포

로지스틱 분포의 CDF는 로지스틱 함수라고도 합니다.
$\text{cdf}(x) = \frac{1}{(1 + exp(-x))}$
로지스틱 분포의 역 CDF는 로짓(logit) 함수라고도 합니다.
$\text{inverse\_cdf}(p) = log(\frac{p}{(1 - p)})$

연습해 봅시다!

R 중급 회귀