안녕하세요 교수님. 멋진 패키지 개발해 주셔서 감사합니다.
덕분에 논문 작업 및 통계에 큰 도움을 받고 있습니다.
autoReg(fit, uni=TRUE)를 시행하면, univariable 결과에서 p values가 0.2 미만인 변수들만 선택해서, multi를 돌리는데요.
혹시 multi로 가게 되는 변수를 p-value threshhold로 선택하는 것이 아닌, mannual로 선택할 수 있는 방법도 있을 까요?
autoReg(fit, uni=TRUE) %>% myft() 를 시행했을때 보여지는 표가 바로 논문에 사용할 수 있을 정도로 좋은데,
해당 표에서 multivariable에 들어가는 변수들을 매뉴얼로 선택할 수 있다면 더할 나위 없을 것 같은데, 방법이 있는지 여쭈어봅니다.
항상 감사드리고, 이제 연말이 다가오는데 남은 한해 행복한 일들만 가득하시길 빕니다.
P.S. autoreg의 gaze는 mytable과 거의 동일해 보이는데, gaze 함수를 따로 개발하신 이유가 있으실까요? mytable 대비 어떠한 장점이 있는 것인지 궁금합니다.
Comment 6
-
cardiomoon
2022.11.16 12:01
-
김찬식
2022.11.21 12:05
와우....엄청나군요. 감사합니다!!!!
-
cardiomoon
2022.11.16 12:02
gaze 함수를 따로 만든 이유는 mytable() 함수를 autoReg() 함수에서 쓸수 있도록 새로 만든 것입니다.^^
myft()함수를 쓸때 warning 이 나오는 분은 github에 있는 autoReg 0.3.1을 사용하시기 바랍니다.
-
김찬식
2022.11.21 12:07
그렇군요! 그렇다면 mytable 함수와 gaze 함수 내부에서 돌아가는 로직은 동일하다고 보면 될까요?
동일하다면 gaze 함수가 myft()와 함께 쓸 수 있어 좋아보여 향후 통계를 돌릴때 mytable을 gaze로 대체하여 사용할까 합니다.
-
cardiomoon
2022.11.21 15:53
거의 동일하다고 보시면 됩니다.
-
김찬식
2022.11.21 16:19
답변 감사합니다. 연말이 다가오는데, 행복한 일들만 가득하시길 빕니다^^
1) 먼저 univariate 모형만 만드신 후
2) 원하시는 설명변수를 골라 multivariate 모형을 만드시고
3) addFitSummary() 함수로 합치시면 됩니다.
예를 들어 다음을 실행해보셔요
require(autoReg)
require(survival)
require(dplyr)
data(cancer,package="survival")
fit=coxph(Surv(time,status)~rx+age+sex+nodes+obstruct+perfor,data=colon)
ft1<-autoReg(fit,uni=TRUE,multi=FALSE)
fit2=coxph(Surv(time,status)~rx+age+sex,data=colon)
addFitSummary(ft1,fit2,statsname="HR (multivariable)") %>% myft()