獨斷論

# 두 개의 데이터가 있다고 가정 data1

27 Survival analysis 27.1 Overview Survival analysis focuses on describing for a given individual or group of individuals, a defined point of event called the failure (occurrence of a disease, cure from a disease, death, relapse after response to treatment…) that occurs after a period of time called failure time (or follow-up time in cohort/population-based studies) during which individuals are ..

### ### Tutorial of survival analysis ### ## subject time event ## 1 1 3 0 ## 2 2 5 1 ## 3 3 7 1 ## 4 4 2 1 ## 5 5 18 0 ## 6 6 16 1 ## 7 7 2 1 ## 8 8 9 1 ## 9 9 16 1 ## 10 10 5 0 ## ## where: ## ## subject is the individual’s identifier ## time is the time to event (in years)2 ## event is the event status (0 = censored, 1 = event happened) # library(survival) ## Example 1 # Data input survDat = ..

100명을 무작위로 나이와 심장병이 있는지 조사하였을때 데이터가 아래와 같다고 가정하자. 데이터파일: 위 데이터파일을 내려받고 아래 R script를 실행 chd_df = read.csv("chdage.csv", header = TRUE) plot(CHD ~ Age, data = chd_df) 대략 20대에서 9명은 심장병이 없고 1명만 심장병을 가지고 있는 반면, 60대에서는 심장병이 없는 사람은 2명이고 심장병을 가지고 있는 사람은 8명이다. 달리말하면 20대에서 심장병을 갖고 있을 확률은 0.1이고 60대가 심장병을 가질 확률은 0.8이다. 40대 초반이 심장병을 가질 확률은 대략 0.5이다. 이러한 확률은 아래와 같은 함수로 나타낼수가 있다(https://en.wikipedia.org/wiki/L..

There are a number of extractor functions (functions that extract or derive specific information from a model) available including: residuals() Extracts the residuals from the model rstandard() Extracts the standardized residuals from the model rstudent() Extracts the studentized residuals from the model fitted() Extracts the predicted (expected) response values (on the link scale) at the observ..

9M1 풀이 9장 본문의 R 스크립트 중 sigma 부분만 바꾸면 된다. 아래 m91과 m911을 비교 library(rethinking) rm(list=ls()) data(rugged) d = rugged d$log_gdp = log(d$rgdppc_2000) dd = d[complete.cases(d$rgdppc_2000), ] dd$log_gdp_std = dd$log_gdp / mean(dd$log_gdp) dd$rugged_std = dd$rugged / max(dd$rugged) dd$cid = ifelse(dd$cont_africa == 1, 1, 2) dat_slim = list( log_gdp_std = dd$log_gdp_std, rugged_std = dd$rugged_std, cid..

R의 lm( )과 python의 statsmodel을 이용하여 단순회귀분석을 비교해보도록 한다. R의 lm( ) 독립변수 x와 종속변수 y를 console에서 직접입력한다. 총 50개의 sample size이다. x = c(2.655,3.721,5.729,9.082,2.017,8.984,9.447,6.608,6.291,0.618, 2.06,1.766,6.87,3.841,7.698,4.977,7.176,9.919,3.8,7.774, 9.347,2.121,6.517,1.256,2.672,3.861,0.134,3.824,8.697,3.403, 4.821,5.996,4.935,1.862,8.274,6.685,7.942,1.079,7.237,4.113, 8.209,6.471,7.829,5.53,5.297,7.89..

Introduction to Bayesian inference with JAGS 강의 정리 3 channel9.msdn.com/Events/useR-international-R-User-conferences/useR-International-R-User-2017-Conference/Introduction-to-Bayesian-inference-with-JAGS-II Introduction to Bayesian inference with JAGS II Introduction to Bayesian inference with JAGS II channel9.msdn.com 병원수술 완치 예제 n