hazard fuction 그래프를 그리고자 합니다.

[통계초짜]

data.csv

SPSS 로 통계를 겨우겨우 하고 있는 초짜 1인입니다.

현재 survival analysis를 하고 있는데요,,, 특정 네그룹간의 생존분석을 시행하였고, Cox regression을 시행하여 F/U 기간동안의hazard ratio를 구하였지만 누적시간의 흐름에 따른 생존양상의 추이를 보는 것이 중요하겠다고 여겨져서(정확히 이야기하자면 early recur인가 late recur인가,, 사건발생이 recurrence입니다) 여차여차 알아보다가 time에 따른 hazard function graph를 그리는 것이 좋겠다는 결론을 얻었습니다. 그리고 STATA 또는 R에서 이것이 가능한 것을 알았고요.

참고 논문에서 R program을 통해 kernel smooth 기법으로 hazard function graph를 그렸다고 되어 있던데...제가 구글을 통해 알아본 바로는 muhaz 기법을 통한 smooth 기법으로 hazard function graph를 그릴수 있다고 되어 있던데... 맞나요? R에서 hazard function graph를 그릴 수 있는 방법에 대해 알고 싶습니다. 기본적인 survival labrary안에서는 제공이 안되는 것같고... 다른 package를 깔아야 하는 것 같은데..

그리고 제가 아직 R을 안지 일주일밖에 안되고 오늘 책을 사본 초짜로서.... CRAN에 떠있는 muhaz package의 예제파일은 그룹내 전체 환자의 위험함수를 그리는 예제이던데... (아직 예제파일을 복사해서 저의 데이터파일화 시켜서 할 줄 밖에 몰라서) 저는 그룹내에 존재하는 네그룹간의 위험함수를 한 그래프에 그려서 그 추이를 보는 것이 저의 목표라서.... 그룹을 세팅하는 것과 각 그룹간 그래프에서 선의 모양과 선의 색을 지정하는 것을 어떻게 하는지.. 를 알아야 한답니다. 고수님들의 답변을 부탁드립니다.

저의 데이터 파일은 이렇습니다. cox regression 할때는 subtype_nstage가 LA_nnegative 인 그룹을 referrence로 한 상태에서 다른 그룹의 recur 및 survival 분석하였습니다.

	subtype	nstage	subtype_nstage	recur	dfs	survival	os
1	LA	negative	LA_nnegative	0	73	0	73
2	LA	negative	LA_nnegative	0	134	0	134
3	LA	negative	LA_nnegative	0	57	0	57
4	LA	negative	LA_nnegative	0	98	0	98
5	LA	negative	LA_nnegative	1	12	0	130
6	LA	negative	LA_nnegative	0	82	0	82
7	LA	negative	LA_nnegative	0	125	0	125
8	LA	negative	LA_nnegative	0	126	0	126
9	LA	negative	LA_nnegative	0	86	0	86
10	LA	negative	LA_nnegative	0	123	0	123
11	LA	negative	LA_nnegative	0	116	0	116
12	LA	negative	LA_nnegative	0	117	0	117
13	LA	negative	LA_nnegative	0	119	0	119
14	LA	negative	LA_nnegative	0	121	0	121
15	LA	negative	LA_nnegative	0	104	0	104
16	LA	negative	LA_nnegative	0	126	0	126
17	LA	negative	LA_nnegative	0	126	0	126
18	LA	negative	LA_nnegative	0	123	0	123
19	LA	negative	LA_nnegative	0	82	0	82
20	LA	negative	LA_nnegative	0	69	0	69
21	LA	negative	LA_nnegative	0	123	0	123
22	LA	negative	LA_nnegative	0	124	0	124
23	LA	negative	LA_nnegative	0	110	0	110
24	LA	negative	LA_nnegative	0	123	0	123
25	LA	negative	LA_nnegative	0	114	0	114
26	LA	negative	LA_nnegative	0	121	0	121
27	LA	negative	LA_nnegative	0	74	0	74
28	LA	negative	LA_nnegative	0	114	0	114
29	LA	negative	LA_nnegative	0	88	0	88
30	LA	negative	LA_nnegative	0	121	0	121
31	LA	negative	LA_nnegative	0	12	0	12
32	LA	negative	LA_nnegative	0	82	0	82
33	LA	negative	LA_nnegative	0	110	0	110
34	LA	negative	LA_nnegative	0	71	0	71
35	LA	negative	LA_nnegative	0	84	0	84
36	LA	negative	LA_nnegative	1	7	0	117
37	LA	negative	LA_nnegative	0	75	0	75
38	LA	negative	LA_nnegative	0	86	0	86
39	LA	negative	LA_nnegative	0	116	0	116
40	LA	negative	LA_nnegative	0	92	0	92
41	LA	negative	LA_nnegative	0	111	0	111
42	LA	negative	LA_nnegative	0	114	0	114
43	LA	negative	LA_nnegative	0	109	0	109
44	LA	negative	LA_nnegative	0	115	0	115
45	LA	negative	LA_nnegative	0	74	0	74
46	LA	negative	LA_nnegative	1	51	0	70
47	LA	negative	LA_nnegative	0	68	0	68
48	LA	negative	LA_nnegative	0	111	0	111
49	LA	negative	LA_nnegative	0	111	0	111
50	LA	negative	LA_nnegative	1	27	1	77
51	LA	negative	LA_nnegative	0	98	0	98
52	LA	negative	LA_nnegative	0	110	0	110
53	LA	negative	LA_nnegative	0	109	0	109
54	LA	negative	LA_nnegative	0	104	0	104
55	LA	negative	LA_nnegative	0	14	0	14
56	LA	negative	LA_nnegative	0	93	0	93
57	LA	negative	LA_nnegative	0	96	0	96
58	LA	negative	LA_nnegative	0	44	0	44
59	LA	negative	LA_nnegative	0	94	0	94
60	LA	negative	LA_nnegative	0	104	0	104
61	LA	negative	LA_nnegative	0	93	0	93
62	LA	negative	LA_nnegative	0	61	0	61
63	LA	negative	LA_nnegative	0	98	0	98
64	LA	negative	LA_nnegative	0	99	0	99
65	LA	negative	LA_nnegative	0	113	0	113
66	LA	negative	LA_nnegative	0	95	0	95
67	LA	negative	LA_nnegative	0	97	0	97
68	LA	negative	LA_nnegative	0	32	0	32
69	LA	negative	LA_nnegative	0	95	0	95
70	LA	negative	LA_nnegative	0	94	0	94
71	LA	negative	LA_nnegative	0	74	0	74
72	LA	negative	LA_nnegative	0	93	0	93
73	LA	negative	LA_nnegative	0	73	0	73
74	LA	negative	LA_nnegative	0	95	0	95
75	LA	negative	LA_nnegative	0	85	0	85
76	LA	negative	LA_nnegative	0	80	0	80
77	LA	negative	LA_nnegative	0	71	0	71
78	LA	negative	LA_nnegative	0	87	0	87
79	LA	negative	LA_nnegative	0	85	0	85
80	LA	negative	LA_nnegative	0	64	1	64
81	LA	negative	LA_nnegative	0	80	0	80
82	LA	negative	LA_nnegative	0	65	0	65
83	LA	negative	LA_nnegative	0	72	0	72
84	LA	negative	LA_nnegative	0	81	0	81
85	LA	negative	LA_nnegative	0	36	0	36
86	LA	negative	LA_nnegative	0	71	0	71
87	LA	negative	LA_nnegative	0	72	0	72
88	LA	negative	LA_nnegative	0	74	0	74
89	LA	negative	LA_nnegative	0	77	0	77
90	LA	negative	LA_nnegative	0	64	0	64
91	LA	negative	LA_nnegative	0	72	0	72
92	LA	negative	LA_nnegative	0	67	0	67
93	LA	negative	LA_nnegative	0	61	0	61
94	LA	negative	LA_nnegative	0	76	0	76
95	LA	negative	LA_nnegative	0	67	0	67
96	LA	negative	LA_nnegative	0	70	0	70
97	LA	negative	LA_nnegative	1	6	0	55
98	LA	negative	LA_nnegative	0	60	0	60
99	LA	negative	LA_nnegative	0	71	0	71
100	LA	negative	LA_nnegative	0	73	0	73
101	LA	negative	LA_nnegative	0	60	0	60
102	LA	negative	LA_nnegative	0	74	0	74
103	LA	negative	LA_nnegative	0	73	0	73
104	LA	negative	LA_nnegative	0	73	0	73
105	LA	negative	LA_nnegative	0	69	0	69
106	LA	negative	LA_nnegative	0	72	0	72
107	LA	negative	LA_nnegative	0	69	0	69
108	LA	negative	LA_nnegative	0	70	0	70
109	LA	negative	LA_nnegative	0	59	0	59
110	LA	negative	LA_nnegative	0	72	0	72
111	LA	negative	LA_nnegative	0	19	0	19
112	LA	negative	LA_nnegative	0	58	0	58
113	LA	negative	LA_nnegative	0	67	0	67
114	LA	negative	LA_nnegative	0	68	0	68
115	LA	negative	LA_nnegative	0	66	0	66
116	LA	negative	LA_nnegative	0	67	0	67
117	LA	negative	LA_nnegative	0	41	0	41
118	LA	negative	LA_nnegative	1	53	0	66
119	LA	negative	LA_nnegative	0	64	0	64
120	LA	negative	LA_nnegative	0	65	0	65
121	LA	negative	LA_nnegative	0	65	0	65
122	LA	negative	LA_nnegative	0	59	0	59
123	LA	negative	LA_nnegative	0	65	0	65
124	LA	negative	LA_nnegative	0	64	0	64
125	LA	negative	LA_nnegative	0	61	0	61
126	LA	negative	LA_nnegative	0	59	0	59
127	LA	negative	LA_nnegative	0	62	0	62
128	LA	negative	LA_nnegative	0	59	0	59
129	LA	negative	LA_nnegative	0	62	0	62
130	LA	negative	LA_nnegative	0	62	0	62
131	LA	negative	LA_nnegative	0	56	0	56
132	LA	negative	LA_nnegative	1	15	0	36
133	LA	negative	LA_nnegative	0	60	0	60
134	LA	negative	LA_nnegative	0	59	0	59
135	LA	negative	LA_nnegative	0	59	0	59
136	LA	negative	LA_nnegative	0	29	1	29
137	LA	negative	LA_nnegative	0	57	0	57
138	LA	negative	LA_nnegative	0	52	0	52
139	LA	negative	LA_nnegative	0	57	0	57
140	LA	negative	LA_nnegative	0	57	0	57
141	LA	negative	LA_nnegative	0	56	0	56
142	LA	negative	LA_nnegative	0	55	0	55
143	LA	negative	LA_nnegative	0	53	0	53
144	LA	negative	LA_nnegative	0	53	0	53
145	LA	negative	LA_nnegative	0	53	0	53
146	LA	negative	LA_nnegative	0	52	0	52
147	LA	positive	LA_npositive	0	130	0	130
148	LA	positive	LA_npositive	0	109	0	109
149	LA	positive	LA_npositive	0	101	0	101
150	LA	positive	LA_npositive	0	128	0	128
151	LA	positive	LA_npositive	0	127	0	127
152	LA	positive	LA_npositive	0	109	0	109
153	LA	positive	LA_npositive	0	125	0	125
154	LA	positive	LA_npositive	1	82	0	124
155	LA	positive	LA_npositive	0	107	0	107
156	LA	positive	LA_npositive	0	116	0	116
157	LA	positive	LA_npositive	0	116	0	116
158	LA	positive	LA_npositive	0	114	0	114
159	LA	positive	LA_npositive	0	100	0	100
160	LA	positive	LA_npositive	0	111	0	111
161	LA	positive	LA_npositive	0	101	0	101
162	LA	positive	LA_npositive	0	107	0	107
163	LA	positive	LA_npositive	0	56	0	56
164	LA	positive	LA_npositive	0	72	0	72
165	LA	positive	LA_npositive	0	63	0	63
166	LA	positive	LA_npositive	0	103	0	103
167	LA	positive	LA_npositive	0	95	0	95
168	LA	positive	LA_npositive	0	96	0	96
169	LA	positive	LA_npositive	1	26	0	27
170	LA	positive	LA_npositive	1	45	1	100
171	LA	positive	LA_npositive	0	89	0	89
172	LA	positive	LA_npositive	0	70	0	70
173	LA	positive	LA_npositive	1	25	1	33
174	LA	positive	LA_npositive	0	84	0	84
175	LA	positive	LA_npositive	0	90	0	90
176	LA	positive	LA_npositive	0	90	0	90
177	LA	positive	LA_npositive	0	72	0	72
178	LA	positive	LA_npositive	0	62	0	62
179	LA	positive	LA_npositive	0	74	0	74
180	LA	positive	LA_npositive	0	76	0	76
181	LA	positive	LA_npositive	0	83	0	83
182	LA	positive	LA_npositive	0	81	0	81
183	LA	positive	LA_npositive	0	60	0	60
184	LA	positive	LA_npositive	1	6	0	72
185	LA	positive	LA_npositive	0	69	0	69
186	LA	positive	LA_npositive	0	78	0	78
187	LA	positive	LA_npositive	0	74	1	74
188	LA	positive	LA_npositive	0	66	0	66
189	LA	positive	LA_npositive	0	71	0	71
190	LA	positive	LA_npositive	0	59	0	59
191	LA	positive	LA_npositive	0	67	0	67
192	LA	positive	LA_npositive	0	66	0	66
193	LA	positive	LA_npositive	0	67	0	67
194	LA	positive	LA_npositive	0	71	0	71
195	LA	positive	LA_npositive	0	67	0	67
196	LA	positive	LA_npositive	0	65	0	65
197	LA	positive	LA_npositive	0	67	0	67
198	LA	positive	LA_npositive	0	62	0	62
199	LA	positive	LA_npositive	0	59	0	59
200	LA	positive	LA_npositive	0	59	0	59
201	LA	positive	LA_npositive	0	56	0	56
202	LA	positive	LA_npositive	1	33	0	61
203	LA	positive	LA_npositive	0	5	0	5
204	LA	positive	LA_npositive	0	59	0	59
205	LA	positive	LA_npositive	0	58	0	58
206	LA	positive	LA_npositive	0	52	0	52
207	LA	positive	LA_npositive	0	56	0	56
208	TNBC	negative	TNBC_nnegative	0	32	0	32
209	TNBC	negative	TNBC_nnegative	1	18	1	42
210	TNBC	negative	TNBC_nnegative	0	95	0	95
211	TNBC	negative	TNBC_nnegative	1	51	1	60
212	TNBC	negative	TNBC_nnegative	0	130	0	130
213	TNBC	negative	TNBC_nnegative	1	8	0	36
214	TNBC	negative	TNBC_nnegative	0	27	0	27
215	TNBC	negative	TNBC_nnegative	0	100	0	100
216	TNBC	negative	TNBC_nnegative	0	94	0	94
217	TNBC	negative	TNBC_nnegative	1	50	0	53
218	TNBC	negative	TNBC_nnegative	0	92	0	92
219	TNBC	negative	TNBC_nnegative	0	76	0	76
220	TNBC	negative	TNBC_nnegative	0	116	0	116
221	TNBC	negative	TNBC_nnegative	0	122	0	122
222	TNBC	negative	TNBC_nnegative	0	102	0	102
223	TNBC	negative	TNBC_nnegative	0	113	0	113
224	TNBC	negative	TNBC_nnegative	0	100	0	100
225	TNBC	negative	TNBC_nnegative	0	108	0	108
226	TNBC	negative	TNBC_nnegative	1	24	0	38
227	TNBC	negative	TNBC_nnegative	1	25	1	37
228	TNBC	negative	TNBC_nnegative	0	109	0	109
229	TNBC	negative	TNBC_nnegative	1	39	1	75
230	TNBC	negative	TNBC_nnegative	0	116	0	116
231	TNBC	negative	TNBC_nnegative	0	30	0	30
232	TNBC	negative	TNBC_nnegative	0	116	0	116
233	TNBC	negative	TNBC_nnegative	0	102	0	102
234	TNBC	negative	TNBC_nnegative	0	78	0	78
235	TNBC	negative	TNBC_nnegative	0	105	0	105
236	TNBC	negative	TNBC_nnegative	0	86	0	86
237	TNBC	negative	TNBC_nnegative	0	56	0	56
238	TNBC	negative	TNBC_nnegative	0	108	0	108
239	TNBC	negative	TNBC_nnegative	0	63	0	63
240	TNBC	negative	TNBC_nnegative	0	108	0	108
241	TNBC	negative	TNBC_nnegative	0	54	0	54
242	TNBC	negative	TNBC_nnegative	0	108	0	108
243	TNBC	negative	TNBC_nnegative	0	70	0	70
244	TNBC	negative	TNBC_nnegative	0	73	0	73
245	TNBC	negative	TNBC_nnegative	0	68	0	68
246	TNBC	negative	TNBC_nnegative	0	105	0	105
247	TNBC	negative	TNBC_nnegative	0	61	0	61
248	TNBC	negative	TNBC_nnegative	0	78	0	78
249	TNBC	negative	TNBC_nnegative	1	36	0	62
250	TNBC	negative	TNBC_nnegative	0	88	0	88
251	TNBC	negative	TNBC_nnegative	0	60	0	60
252	TNBC	negative	TNBC_nnegative	0	73	0	73
253	TNBC	negative	TNBC_nnegative	0	94	0	94
254	TNBC	negative	TNBC_nnegative	0	92	0	92
255	TNBC	negative	TNBC_nnegative	0	86	0	86
256	TNBC	negative	TNBC_nnegative	0	80	0	80
257	TNBC	negative	TNBC_nnegative	0	23	0	23
258	TNBC	negative	TNBC_nnegative	0	38	0	38
259	TNBC	negative	TNBC_nnegative	0	28	0	28
260	TNBC	negative	TNBC_nnegative	0	84	0	84
261	TNBC	negative	TNBC_nnegative	0	79	0	79
262	TNBC	negative	TNBC_nnegative	0	83	0	83
263	TNBC	negative	TNBC_nnegative	0	82	0	82
264	TNBC	negative	TNBC_nnegative	0	80	0	80
265	TNBC	negative	TNBC_nnegative	0	74	0	74
266	TNBC	negative	TNBC_nnegative	0	70	0	70
267	TNBC	negative	TNBC_nnegative	0	73	0	73
268	TNBC	negative	TNBC_nnegative	0	62	0	62
269	TNBC	negative	TNBC_nnegative	0	72	0	72
270	TNBC	negative	TNBC_nnegative	0	69	0	69
271	TNBC	negative	TNBC_nnegative	1	31	0	73
272	TNBC	negative	TNBC_nnegative	0	58	0	58
273	TNBC	negative	TNBC_nnegative	0	56	0	56
274	TNBC	negative	TNBC_nnegative	0	49	0	49
275	TNBC	negative	TNBC_nnegative	0	63	0	63
276	TNBC	negative	TNBC_nnegative	0	49	0	49
277	TNBC	negative	TNBC_nnegative	0	60	0	60
278	TNBC	negative	TNBC_nnegative	0	60	0	60
279	TNBC	negative	TNBC_nnegative	1	17	0	63
280	TNBC	negative	TNBC_nnegative	0	62	0	62
281	TNBC	negative	TNBC_nnegative	0	57	0	57
282	TNBC	negative	TNBC_nnegative	0	60	0	60
283	TNBC	negative	TNBC_nnegative	0	61	0	61
284	TNBC	negative	TNBC_nnegative	0	61	0	61
285	TNBC	negative	TNBC_nnegative	0	60	0	60
286	TNBC	negative	TNBC_nnegative	0	24	0	24
287	TNBC	negative	TNBC_nnegative	0	29	0	29
288	TNBC	negative	TNBC_nnegative	0	54	0	54
289	TNBC	negative	TNBC_nnegative	0	52	0	52
290	TNBC	negative	TNBC_nnegative	0	52	0	52
291	TNBC	negative	TNBC_nnegative	0	42	0	42
292	TNBC	positive	TNBC_npositive	1	29	0	134
293	TNBC	positive	TNBC_npositive	0	75	0	75
294	TNBC	positive	TNBC_npositive	0	99	0	99
295	TNBC	positive	TNBC_npositive	0	86	0	86
296	TNBC	positive	TNBC_npositive	0	102	0	102
297	TNBC	positive	TNBC_npositive	1	19	1	23
298	TNBC	positive	TNBC_npositive	0	64	0	64
299	TNBC	positive	TNBC_npositive	1	43	0	109
300	TNBC	positive	TNBC_npositive	0	109	0	109
301	TNBC	positive	TNBC_npositive	1	5	1	17
302	TNBC	positive	TNBC_npositive	0	12	1	12
303	TNBC	positive	TNBC_npositive	0	89	0	89
304	TNBC	positive	TNBC_npositive	1	32	0	47
305	TNBC	positive	TNBC_npositive	1	83	0	98
306	TNBC	positive	TNBC_npositive	1	26	0	35
307	TNBC	positive	TNBC_npositive	0	81	0	81
308	TNBC	positive	TNBC_npositive	0	79	0	79
309	TNBC	positive	TNBC_npositive	0	57	0	57
310	TNBC	positive	TNBC_npositive	0	3	1	3
311	TNBC	positive	TNBC_npositive	0	59	0	59
312	TNBC	positive	TNBC_npositive	0	59	0	59
313	TNBC	positive	TNBC_npositive	1	12	0	29
314	TNBC	positive	TNBC_npositive	1	30	0	55
315	TNBC	positive	TNBC_npositive	0	53	0	53
316	TNBC	positive	TNBC_npositive	0	53	0	53

[cardiomoon]

아님 엑셀화일을 올려주셔도 되구요.

[통계초짜]

데이터 파일을 올렸습니다. 본 글에 데이터 붙인건 지워지지가 않네요^^

[cardiomoon]

이런 그래프를 원하시나요?