안녕하세요
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 |
---|---|---|---|
1233 | 선형회귀분석 시 보정변수 사용 관련 질문드립니다. [2] | 전한울 | 2023.09.18 |
1232 | 그래프 관련 문의드립니다. | nanak | 2023.09.14 |
1231 | 영수증 메일로 받을수 있을까요? [1] | 딸기공주 | 2023.09.13 |
1230 | mycsv 관련 문의드립니다. [1] | 도삼 | 2023.09.12 |
1229 | matched group 을 대상으로 stratified cox | saga | 2023.09.06 |
1228 | IPTW 시 SMD 계산 관련 여쭙습니다. | 장한나 | 2023.09.04 |
1227 | 회귀분석 관련 문의드립니다. [2] | 전한울 | 2023.09.04 |
1226 | PSM matching 실행시 disconnected 됩니다... [1] | 맹구모션 | 2023.08.25 |
1225 | 메타분석_출간오류분석 관련 질문드립니다. [2] | 블루레몬 | 2023.08.21 |
1224 | plot 문의 드립니다. | 1122 | 2023.08.20 |
1223 | semMediation 설치문의 [2] | 마술샘 | 2023.08.11 |
1222 | 로지스틱 회귀분석 문의드립니다. [2] | JKJ | 2023.08.09 |
1221 | 지출 문의 드립니다. [1] | 블루레몬 | 2023.08.08 |
1220 | autoReg 실행 중 에러 [1] | pkjsdf123 | 2023.08.08 |
1219 | 영수증 처리 [1] | britnepak | 2023.07.28 |
1218 | multipleROC 커브 글을 따라가다 궁금한 것이 생겨서 질문드립니다. [2] | rgppwls | 2023.07.20 |
1217 | disconnected from server 오류 [2] | JKJ | 2023.07.13 |
1216 | mytable에서 변수 순서와 ztable적용시 변수가 숫자로 변형 [1] | 둥니맨 | 2023.07.12 |
1215 | box plot 시 X 축 label이 다 나오지 않아 질문드립니다. [2] | 난둥 | 2023.07.03 |
1214 | 안녕하십니까 교수님 PSM 결과에 대해 궁금한 것이 있어 여쭤봅니다 [2] | 김도훈 | 2023.06.29 |
집단변수가 2군이라서 col=1:3을 1:2로도 변경해봤는데, 같은 오류가 뜹니다... DIS 변수를 못 읽는거같은데... str(KMcurve1)로 자료를 보면 변수로 인식되어 있긴 합니다.