안녕하세요
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로 지정하는 것이군요. 감사합니다. !!!
집단변수가 2군이라서 col=1:3을 1:2로도 변경해봤는데, 같은 오류가 뜹니다... DIS 변수를 못 읽는거같은데... str(KMcurve1)로 자료를 보면 변수로 인식되어 있긴 합니다.