안녕하세요
R 통계와 그래프를 보면서 KM curve를 그리는데 자꾸 오류가 뜨는데 뭐가 문제인지 잘 모르겠습니다.
fit1=survfit(Surv(fuduration,dyslipidemia==1)~DIS,data=KMcurve1)
plot(fit1,col=1:3,lty=1:3)
legend("topright",legend=levels(DIS),col = 1:3,lty = 1:3)
위와 같이 명령어 넣었습니다. 이벤트까지의 시간이 "fuduration"이고, 이벤트는 고지혈증 발생 이라 "dyslipidemia" 입니다.
집단 변수를 "DIS"로 잡고, 두군간의 고지혈증 발생 KM curve를 그리고 싶어 위와 같이 돌렸는데
> legend("topright",legend=levels(DIS),col = 1:3,lty = 1:3)
Error in levels(DIS) : 객체 'DIS'를 찾을 수 없습니다
이렇게 뜨면서 오류가 납니다.
집단 변수 빼고 전체 그룹 KM curve는 잘 그려집니다..
어떤게 잘못되었을까요.. ?? 질문드립니다.
Comment 10
-
스날푸
2020.04.15 16:46
-
cardiomoon
2020.04.15 20:13
KMcurve1$DIS 로 해보시겠어요?
-
스날푸
2020.04.16 09:27
fit4=survfit(Surv(fuduration,dyslipidemia==1)~DIS,data=KMcurve2$DIS)plot(fit4,col=1:3,lty=1:3)legend("topright",legend=levels(DIS),col(1:3),lty = 1:3)(fuduration이 number로 안 되어 있어서 숫자로 변환하여 KMcurve2로 새로운 데이터 넣고 돌려봤습니다.)Error in levels(DIS) : 객체 'DIS'를 찾을 수 없습니다로 계속뜹니다.변수를 잘못 넣은걸까요??str(KMcurve2) 해보면$ DIS : Factor w/ 2 levels "biologics","cDMARDs": 1 1 1 1 1 1 1 1 1 1 ...$ fuduration : num 48 49 50 48 24 12 48 44 47 49 ...로 나오긴 합니다 ㅠㅠ;; -
cardiomoon
2020.04.16 09:48
fit4=survfit(Surv(fuduration,dyslipidemia==1)~DIS,data=KMcurve2)plot(fit4,col=1:3,lty=1:3)legend("topright",legend=levels(KMcurve2$DIS),col(1:3),lty = 1:3)이렇게 한번 해보셔요 -
스날푸
2020.04.16 11:09
교수님... 위에 해준신데로 해봤는데
> legend("topright",legend=levels(KMcurve5$DIS),col = 1:3,lty = 1:3)
Error in legend("topright", legend = levels(KMcurve5$DIS), col = 1:3, :
'legend'의 길이가 0입니다
로 뜹니다.. 변수 이름도 DIS -> dis로 바꿔서 해보고 했는데... 계속 에러가 뜹니다
-
cardiomoon
2020.04.16 12:52
데이터파일을 올려봐주시겠어요?
-
스날푸
2020.04.16 15:47
이 파일입니다.
event를 dyslipidemia로 잡고, 그룹은 DIS로 잡았습니다.
fuduration 이 event까지의 시간 입니다.
감사합니다 !!
-
cardiomoon
2020.04.16 23:05
KMcurve1=read.csv("~/Downloads/KMcurve_20200416.csv")KMcurve1$DIS=factor(KMcurve1$DIS)fit1=survfit(Surv(fuduration,dyslipidemia==1)~DIS,data=KMcurve1)plot(fit1,col=1:2,lty=1:2)legend("topright",legend=levels(KMcurve1$DIS),col = 1:2,lty = 1:2) -
cardiomoon
2020.04.16 23:06
-
스날푸
2020.04.17 09:04
감사합니다.
factor로 지정하는 것이군요. 감사합니다. !!!
No. | Subject | Author | Date |
---|---|---|---|
914 | R을 이용한 조건부과정분석 [1] | MJ89 | 2020.08.18 |
913 | 로그인 에러가 자주뜹니다. [1] | 피터팬79 | 2021.01.11 |
912 | Plot inset 만들기 | cardiomoon | 2015.06.09 |
911 | 생존분석 중에 문의 드립니다. [1] | hahachristi | 2017.05.03 |
910 | R에서 그룹별 사분위수 계산 [2] | 버블 | 2021.05.04 |
909 | Sigmoid curve (회귀분석) 질문드립니다. [1] | chanheebox | 2017.09.19 |
908 | 메타분석 질문 [2] | 남일성 | 2017.10.06 |
907 | WGCNA 분석 [1] | zkbi | 2020.04.23 |
906 | Univariate and Multiple Logistic Regression: predictor의 선택 | 신희순 | 2018.10.04 |
» | Survival 패키지를 이용한 KM curve 작성시 오류 문의 [10] | 스날푸 | 2020.04.15 |
904 | 교수님 안녕하세요. multipleROC package에 관해 질문드리고 싶습니다. [4] | 포카리 | 2021.03.09 |
903 | effect size 문의 | 보라도리 | 2019.09.28 |
902 | 안녕하세요 r 구조방정식 질문이 있습니다 [6] | 마술샘 | 2016.02.23 |
901 | standard error 에 관한 질문입니다. [6] | R까기 | 2016.01.28 |
900 | 생존분석에서 연속형변수의 cutpoint에 대해 질문드립니다. [2] | zepiros15 | 2018.12.17 |
899 | 생존 분석시 문제 [3] | selteen | 2016.06.09 |
898 | 결측치 제거에 관한 문의입니다. [1] | 방랑자 | 2015.03.26 |
897 | 8개의 군에서 Cox regression 에러 [5] | 정암 | 2017.05.21 |
896 | R 3.2.3 is released!!! [1] | fermat39 | 2015.12.15 |
895 | scientific reports 투고 관련 문의 드립니다. [2] | 날아라랭지 | 2018.09.16 |
집단변수가 2군이라서 col=1:3을 1:2로도 변경해봤는데, 같은 오류가 뜹니다... DIS 변수를 못 읽는거같은데... str(KMcurve1)로 자료를 보면 변수로 인식되어 있긴 합니다.