안녕하세요
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 |
---|---|---|---|
1234 | 요즘 아주 잘 쓰고 있습니다. [1] | 김승욱 | 2015.08.17 |
1233 | 자료 (xlsx or CSV) [2] | medi79 | 2015.08.31 |
1232 | 혹시 이런 질문 드려도 되는지 모르겠습니다. [5] | 고든 | 2015.09.07 |
1231 | RM ANOVA [2] | flocky | 2015.09.09 |
1230 | ANCOVA에 대해서 여쭤보고 싶습니다. [5] | 고든 | 2015.09.13 |
1229 | sensitivity, specificity 계산가능여부 여쭤봅니다. [2] | 곽미리 | 2015.09.16 |
1228 | 표만들기에서 질문드립니다. [3] | wiswis | 2015.09.22 |
1227 | Standard error [2] | R까기 | 2015.09.24 |
1226 | 결측값 처리와 관련하여 [3] | 화모기 | 2015.10.07 |
1225 | 표만들기_2차그룹변수로 했을 때 결과물 이상 [1] | 화모기 | 2015.10.10 |
1224 | 로지스틱 회귀분석 반응변수 수 제한 [3] | 화모기 | 2015.10.11 |
1223 | 향후 5년간 핵심 화두가 Cognitive Business가 될 것 같아 먼저 여기에 글을 올립니다. [2] | Ben | 2015.10.14 |
1222 | 파일업로드에 대해 문의드립니다. [3] | flapper | 2015.10.15 |
1221 | 스피어만 상관계수(Spearman correlation coefficient) [1] | ㅎㅎㅎ | 2015.10.15 |
1220 | PSM에 관해 문의 드립니다. [3] | sodium75 | 2015.10.19 |
1219 | ROC 분석과 관련하여 질문합니다. [1] | 송민지 | 2015.10.21 |
1218 | 파일 업로드 자체가 안되네요.. [3] | 주민호 | 2015.11.13 |
1217 | 사소한 건의사항이 있습니다. [1] | EPfunny | 2015.11.15 |
1216 | 나무분석 [3] | 조윤성 | 2015.12.02 |
1215 | 'x' and 'units' must have length > 0 에러 메세지 [4] | 빙하 | 2015.12.03 |
집단변수가 2군이라서 col=1:3을 1:2로도 변경해봤는데, 같은 오류가 뜹니다... DIS 변수를 못 읽는거같은데... str(KMcurve1)로 자료를 보면 변수로 인식되어 있긴 합니다.