해색주의 조용한 공간

a <- c(8, 2, 0, 18, 1, 1, 28, 1, 1

     , 10, 2, 0, 20, 3, 2, 32, 1, 0

     , 12, 3, 0, 22, 2, 1, 34, 1, 1

     , 14, 3, 0, 24, 1, 0, 38, 3, 2

     , 16, 3, 0, 26, 1, 1)

b <- matrix(a, ncol = 3, byrow=T)

head(b)

c <- data.frame(b)

colnames(c) = c('li', 'n', 'y')

c <- transform(c,

               b = n -y,

               t = y/n,

               li_sq = li ^ 2)

head(c)

x <- c$li

outcome <- cbind(c$y, c$b)

model.c <- glm(outcome ~ c$li, family = binomial(link = "logit"))

summary(model.c)

model.c1 <- glm(outcome ~ c$li + c$li_sq, family = binomial(link = "logit"))

summary(model.c1)

 원래는 다양한 테스트를 해보려고 했는데, 아무래도 R은 많이 사용하지 않으니 익숙하지가 않다. 시간 내서 대량 데이터 핸들링을 해봐야 하는데, 그런 것은 아무래도 데이터를 많이 돌려봐야 할듯 하네.

medicine = c(rep('A', 500), rep('B', 200), rep('C', 300))

test <- data.frame(medicine)

test

help(transform)

test <- transform(test,

                  cure = ifelse(medicine == "A"

                               , 1

                               , ifelse(medicine == "B", -1, 0)))

test

toread <- "id sex age inc r1 r2 r3 
1  F  35 17  7 2 2 
17  M  50 14  5 5 3 
33  F  45  6  7 2 7 
49  M  24 14  7 5 7 
65  F  52  9  4 7 7 
81  M  44 11  7 7 7 
2   F  34 17  6 5 3 
18  M  40 14  7 5 2 
34  F  47  6  6 5 6 
50  M  35 17  5 7 5" 

survey <- read.table(textConnection(toread), header = TRUE) 
closeAllConnections() 
survey 

http://r.789695.n4.nabble.com/SAS-quot-datalines-quot-or-quot-cards-quot-statement-equivalent-in-R-td949489.html

손차장님이 말씀하신 것이 이것인가 시간날 때 회사에서 해봐야겠다.

http://support.sas.com/documentation/cdl/en/lrdict/64316/HTML/default/viewer.htm#a000161869.htm