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.brasil@ipec.fiocruz.br
email: emmanuel.brasil@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]]

