R friends,
Im trying to construct a table with results of agreement analysis of two
observers. Most of the variables I have are categorical answers and Im using
the epicalc::kap function. Initially I was intending to use a function
estimate the ICCC (Intraclass correlation coefficient) for continuous
variables, but the in epi packages I didnt find any function for it.
Instead, I found epi.ccc that estimate the Lin concordance correlation
agreement, which I was not aware of. I did look for some bibliography and it
seems that this statistic is intended to estimate reliability as ICCC and
give very similar results.
However, this function is giving NAs instead of the desired estimate. Below
are the columns of my data and the function output. Perhaps someone would
suggest something to workaround, a better alternative to make it work or
identify something Im doing wrong.
> clin[,c('Q6.A1','Q6.A2')]
Q6.A1 Q6.A2
1 53 53
2 52 52
3 19 19
4 20 NA
5 53 53
6 32 32
7 9 10
8 49 49
9 64 64
10 77 NA
11 14 14
12 52 52
13 12 13
14 54 54
15 79 79
16 43 43
17 38 39
18 45 46
19 41 41
20 64 64
21 53 54
22 66 NA
23 56 56
24 58 58
25 47 48
26 62 63
27 42 42
28 50 50
29 14 14
30 80 81
31 60 60
32 87 88
33 46 47
34 62 62
35 30 30
36 24 24
37 49 49
38 27 27
39 22 22
40 25 25
41 54 54
42 68 69
43 36 37
44 73 73
45 56 55
46 36 36
47 30 30
48 59 58
49 53 54
50 30 30
51 42 42
52 38 0
53 45 46
54 29 29
55 55 55
56 26 27
57 44 44
58 48 48
59 52 52
60 20 20
61 50 51
62 41 41
63 44 44
64 48 48
65 40 40
66 45 45
67 65 65
68 60 60
69 38 37
70 0 29
71 51 52
72 49 49
73 46 46
74 19 19
75 42 42
76 31 31
77 24 25
78 54 54
79 38 38
80 45 46
81 38 38
82 69 69
83 51 51
84 79 80
85 59 60
86 47 47
87 65 65
88 40 40
89 28 28
90 33 33
91 29 29
92 47 47
93 58 NA
94 41 41
95 46 46
96 56 56
97 54 54
98 63 63
99 62 62
100 52 52
101 NA 51
102 24 25
103 54 54
104 70 70
105 46 46
106 29 30
107 28 28
108 50 50
109 52 52
110 61 61
111 51 51
112 52 52
113 53 53
114 54 54
115 36 36
116 61 61
117 14 15
118 52 52
119 54 54
120 56 0
121 68 69
122 58 64
123 49 49
124 29 29
125 45 45
126 63 63
127 67 67
128 43 44
129 62 62
130 60 61
131 18 18
132 38 37
133 44 43
134 27 27
135 73 73
136 63 60
137 27 27
138 59 59
139 31 31
140 47 48
141 44 54
142 45 45
143 55 55
144 54 54
145 49 49
146 67 67
147 46 46
148 59 59
149 57 57
150 34 34
151 24 24
152 54 54
153 83 82
154 59 59
155 83 84
156 55 55
157 42 42
158 54 54
159 36 36
160 56 56
161 72 72
162 20 54
163 65 65
164 62 62
165 27 27
166 58 58
167 54 54
168 67 67
169 65 65
170 43 43
171 70 70
172 68 68
173 62 62
174 39 40
175 20 20
176 55 56
177 0 65
178 63 64
179 40 40
180 29 0
181 57 57
182 24 25
183 58 58
184 46 946
185 74 NA
186 53 53
187 67 67
188 54 54
189 61 61
190 28 28
191 46 46
192 20 20
193 44 44
194 45 45
195 44 44
196 35 35
197 42 42
198 43 43
199 36 37
200 69 69
201 69 69
202 44 44
203 49 NA
204 49 NA
205 49 NA
206 49 NA
207 19 19
208 68 68
209 28 28
210 49 NA
211 69 69
212 74 74
213 59 61
214 74 74
215 15 15
216 26 26
217 36 36
218 39 39
219 54 54
220 42 42
221 62 62
222 17 17
223 52 52
224 64 64
225 65 65
226 35 35
227 69 69
228 32 32
229 42 42
230 34 35
231 42 43
232 55 55
233 17 17
234 46 46
235 34 34
236 45 45
237 56 56
238 28 27
239 34 34
240 47 47
241 42 39
242 52 51
243 45 45
244 37 37
245 54 54
246 26 26
247 68 69
248 66 66
249 49 49
250 58 58
251 53 54
252 26 26
253 75 75
254 37 37
255 80 80
256 54 54
257 55 55
258 38 38
259 70 70
260 48 49
261 43 43
262 88 88
263 31 32
264 61 61
265 57 57
266 36 36
267 48 49
268 17 18
269 53 54
270 75 75
271 60 60
272 69 69
273 49 49
274 60 60
275 69 69
276 50 50
277 67 67
278 66 66
279 50 51
280 0 58
281 27 27
282 61 61
283 56 56
284 40 40
285 43 44
286 45 44
287 35 35
288 46 46
289 35 35
290 64 64
291 50 50
292 23 23
293 63 63
294 38 38
295 42 42
296 63 64
297 65 65
298 77 77
299 58 51
300 19 NA
301 63 63
302 40 40
303 65 65
304 58 58
305 54 54
306 52 53
307 81 81
308 61 61
309 60 60
310 24 24
311 43 43
312 46 46
313 62 62
314 54 54
315 68 68
316 64 65
317 39 39
318 53 52
319 59 58
> epi.ccc(clin$Q6.A1,clin$Q6.A2)
$rho.c
est lower upper
1 NA NA NA
$s.shift
[1] NA
$l.shift
[1] NA
$C.b
[1] NA
$blalt
mean delta
1 53.0 0
2 52.0 0
3 19.0 0
4 NA NA
5 53.0 0
6 32.0 0
7 9.5 -1
8 49.0 0
9 64.0 0
10 NA NA
11 14.0 0
12 52.0 0
13 12.5 -1
14 54.0 0
15 79.0 0
16 43.0 0
17 38.5 -1
18 45.5 -1
19 41.0 0
20 64.0 0
21 53.5 -1
22 NA NA
23 56.0 0
24 58.0 0
25 47.5 -1
26 62.5 -1
27 42.0 0
28 50.0 0
29 14.0 0
30 80.5 -1
31 60.0 0
32 87.5 -1
33 46.5 -1
34 62.0 0
35 30.0 0
36 24.0 0
37 49.0 0
38 27.0 0
39 22.0 0
40 25.0 0
41 54.0 0
42 68.5 -1
43 36.5 -1
44 73.0 0
45 55.5 1
46 36.0 0
47 30.0 0
48 58.5 1
49 53.5 -1
50 30.0 0
51 42.0 0
52 19.0 38
53 45.5 -1
54 29.0 0
55 55.0 0
56 26.5 -1
57 44.0 0
58 48.0 0
59 52.0 0
60 20.0 0
61 50.5 -1
62 41.0 0
63 44.0 0
64 48.0 0
65 40.0 0
66 45.0 0
67 65.0 0
68 60.0 0
69 37.5 1
70 14.5 -29
71 51.5 -1
72 49.0 0
73 46.0 0
74 19.0 0
75 42.0 0
76 31.0 0
77 24.5 -1
78 54.0 0
79 38.0 0
80 45.5 -1
81 38.0 0
82 69.0 0
83 51.0 0
84 79.5 -1
85 59.5 -1
86 47.0 0
87 65.0 0
88 40.0 0
89 28.0 0
90 33.0 0
91 29.0 0
92 47.0 0
93 NA NA
94 41.0 0
95 46.0 0
96 56.0 0
97 54.0 0
98 63.0 0
99 62.0 0
100 52.0 0
101 NA NA
102 24.5 -1
103 54.0 0
104 70.0 0
105 46.0 0
106 29.5 -1
107 28.0 0
108 50.0 0
109 52.0 0
110 61.0 0
111 51.0 0
112 52.0 0
113 53.0 0
114 54.0 0
115 36.0 0
116 61.0 0
117 14.5 -1
118 52.0 0
119 54.0 0
120 28.0 56
121 68.5 -1
122 61.0 -6
123 49.0 0
124 29.0 0
125 45.0 0
126 63.0 0
127 67.0 0
128 43.5 -1
129 62.0 0
130 60.5 -1
131 18.0 0
132 37.5 1
133 43.5 1
134 27.0 0
135 73.0 0
136 61.5 3
137 27.0 0
138 59.0 0
139 31.0 0
140 47.5 -1
141 49.0 -10
142 45.0 0
143 55.0 0
144 54.0 0
145 49.0 0
146 67.0 0
147 46.0 0
148 59.0 0
149 57.0 0
150 34.0 0
151 24.0 0
152 54.0 0
153 82.5 1
154 59.0 0
155 83.5 -1
156 55.0 0
157 42.0 0
158 54.0 0
159 36.0 0
160 56.0 0
161 72.0 0
162 37.0 -34
163 65.0 0
164 62.0 0
165 27.0 0
166 58.0 0
167 54.0 0
168 67.0 0
169 65.0 0
170 43.0 0
171 70.0 0
172 68.0 0
173 62.0 0
174 39.5 -1
175 20.0 0
176 55.5 -1
177 32.5 -65
178 63.5 -1
179 40.0 0
180 14.5 29
181 57.0 0
182 24.5 -1
183 58.0 0
184 496.0 -900
185 NA NA
186 53.0 0
187 67.0 0
188 54.0 0
189 61.0 0
190 28.0 0
191 46.0 0
192 20.0 0
193 44.0 0
194 45.0 0
195 44.0 0
196 35.0 0
197 42.0 0
198 43.0 0
199 36.5 -1
200 69.0 0
201 69.0 0
202 44.0 0
203 NA NA
204 NA NA
205 NA NA
206 NA NA
207 19.0 0
208 68.0 0
209 28.0 0
210 NA NA
211 69.0 0
212 74.0 0
213 60.0 -2
214 74.0 0
215 15.0 0
216 26.0 0
217 36.0 0
218 39.0 0
219 54.0 0
220 42.0 0
221 62.0 0
222 17.0 0
223 52.0 0
224 64.0 0
225 65.0 0
226 35.0 0
227 69.0 0
228 32.0 0
229 42.0 0
230 34.5 -1
231 42.5 -1
232 55.0 0
233 17.0 0
234 46.0 0
235 34.0 0
236 45.0 0
237 56.0 0
238 27.5 1
239 34.0 0
240 47.0 0
241 40.5 3
242 51.5 1
243 45.0 0
244 37.0 0
245 54.0 0
246 26.0 0
247 68.5 -1
248 66.0 0
249 49.0 0
250 58.0 0
251 53.5 -1
252 26.0 0
253 75.0 0
254 37.0 0
255 80.0 0
256 54.0 0
257 55.0 0
258 38.0 0
259 70.0 0
260 48.5 -1
261 43.0 0
262 88.0 0
263 31.5 -1
264 61.0 0
265 57.0 0
266 36.0 0
267 48.5 -1
268 17.5 -1
269 53.5 -1
270 75.0 0
271 60.0 0
272 69.0 0
273 49.0 0
274 60.0 0
275 69.0 0
276 50.0 0
277 67.0 0
278 66.0 0
279 50.5 -1
280 29.0 -58
281 27.0 0
282 61.0 0
283 56.0 0
284 40.0 0
285 43.5 -1
286 44.5 1
287 35.0 0
288 46.0 0
289 35.0 0
290 64.0 0
291 50.0 0
292 23.0 0
293 63.0 0
294 38.0 0
295 42.0 0
296 63.5 -1
297 65.0 0
298 77.0 0
299 54.5 7
300 NA NA
301 63.0 0
302 40.0 0
303 65.0 0
304 58.0 0
305 54.0 0
306 52.5 -1
307 81.0 0
308 61.0 0
309 60.0 0
310 24.0 0
311 43.0 0
312 46.0 0
313 62.0 0
314 54.0 0
315 68.0 0
316 64.5 -1
317 39.0 0
318 52.5 1
319 58.5 1
Abraço forte e que a força esteja com você,
Dr. Pedro Emmanuel A. A. do Brasil
Instituto de Pesquisa Clínica Evandro Chagas
Fundação Oswaldo Cruz
Rio de Janeiro - Brasil
Av. Brasil 4365
Tel 55 21 3865-9648
email: pedro.bra...@ipec.fiocruz.br
email: emmanuel.bra...@gmail.com
---Apoio aos softwares livres
www.zotero.org - gerenciamento de referências bibliográficas.
www.broffice.org ou www.openoffice.org - textos, planilhas ou apresentações.
www.epidata.dk - entrada de dados.
www.r-project.org - análise de dados.
www.ubuntu.com/ - sistema operacional
[[alternative HTML version deleted]]
_______________________________________________
R-sig-Epi@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-sig-epi
