Re: [問題] GLM回歸請教

[jasonfghx (工人)]

※ 引述《clansoda (小笨)》之銘言：
idx <- c(1, 3, 5) #假設你要1、3、5這三個variables

my_formula <-

as.formula(paste("Classkd ~ ", paste(colnames(training.data)[idx],
           collapse = "+")))

剩下來的就是把你的logistic function裡面的formula改成上面這個就好了

--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 1.165.133.189
※ 文章網址: https://www.ptt.cc/bbs/R_Language/M.1522983688.A.D9C.html


謝謝這位大大回應

但因為我的變數中的有空格

用你的語法會出問題

所以我改成這樣

as.formula(paste("Classkd ~ ", paste("'",new[c(1,3,4)] ,collapse = "+","'")))

加了一些'就正常了

謝謝

--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 180.217.254.88
※ 文章網址: https://www.ptt.cc/bbs/R_Language/M.1523001997.A.4A2.html
→ clansoda: 對的，做法只是一個概念，還是要依照資料做變化
再請教一下

my_formula <-as.formula(paste("Classkd ~ ", paste("'",new[idx] ,collapse =
"+","'")))

data.glm <- glm(my_formula,data.training,family="binomial")

我這樣做他有警示
Error in terms.formula(formula, data = data) :
  invalid model formula in ExtractVars

想請教我可以怎麼改進?
※ 編輯: jasonfghx (180.217.254.88), 04/06/2018 16:40:25

data.glm <- glm(formula=Classckd~.,data=data.training,family="binomial")
我用這樣是全部都讀

Age                          -2.441e-09  1.827e+03       0        1
`Blood Pressure`              4.457e-09  2.665e+03       0        1
`Specific Gravity`            5.912e-06  7.377e+06       0        1
Albumin                      -7.880e-07  4.209e+04       0        1
Sugar                        -1.612e-07  8.159e+04       0        1
`Red Blood Cellsabnormal`     7.015e-07  7.672e+04       0        1
`Red Blood Cellsnormal`              NA         NA      NA       NA
`Pus Cellabnormal`           -3.664e-07  8.803e+04       0        1
`Pus Cellnormal`                     NA         NA      NA       NA
`Pus Cell clumpsnotpresent`  -6.867e-07  1.370e+05       0        1
`Pus Cell clumpspresent`             NA         NA      NA       NA
Bacterianotpresent            1.807e-06  2.058e+05       0        1
Bacteriapresent                      NA         NA      NA       NA
`Blood Glucose Random`        2.400e-10  7.645e+02       0        1
`Blood Urea`                  7.269e-09  1.332e+03       0        1
`Serum Creatinine`           -1.812e-07  2.287e+04       0        1

但我現在希望說我不要全部都去回歸

只是就語法上有的有` 有的沒有

就讓我很困擾
※ 編輯: jasonfghx (180.217.254.88), 04/06/2018 16:49:45
→ clansoda: 不知道你的new是什麼東西
→ clansoda: 因為我是用colnames來抓的
> my_formula <-as.formula(paste("Classckd ~ ",
paste("`",colnames(data.training)[idx],collapse ="+","`")))
> data.glm <- glm(my_formula,data.training,family="binomial")
Error in eval(predvars, data, env) : object ' Age ' not found

我跟你用一樣來討論好了  謝謝!!
※ 編輯: jasonfghx (180.217.254.88), 04/06/2018 17:23:41
→ jasonfghx: 我覺得又是" Age " 中的space造成錯誤的
※ 編輯: jasonfghx (180.217.254.88), 04/06/2018 17:26:26
→ clansoda: 那你要不要把你的colnames修正一下呢?
→ clansoda: colnames(training.data) <- trimws(colnames(train.dat
→ clansoda: a))
→ clansoda: 因為手上也沒有你的data，有時候作業系統也會有影響
→ clansoda: 如果還是not work，就拿大絕招出來
→ clansoda: training.data %>% select(c(1, 3, 5))
→ clansoda: 直接這樣做配上.來選你的model也可以，只是醜一點而已
→ jasonfghx: 我覺得這題就乾脆用刪掉法比較快