Luk Arbuckle
2012-Feb-01 18:47 UTC
[R] Function to compute multi-response, multi-rater kappa?
I'm looking for a function in R that extends kappa to multiple raters when there is more than one response per subject. For example, say a group of doctors have to assign diseases to patients. Each patient will be assigned one to many diseases, and the number of doctors assigning diseases to any one patient will be two to many. Here's an extremely simple example of the type of data I might have (two patients, three doctors, five diagnoses, and doctor 'x' assigns multiple diagnoses per patient): pat<-c('a','a','a','b','b','b','b') doc<-c('x','x','y','x','x','y','z') dx<-c('1','2','2','3','4','4','5') df<-data.frame(pat=pat,doc=doc,dx=dx) df I found a paper that can address this [1], although I can't find it referenced on cran. What I have found on cran doesn't address the multi-response case. I don't want to reinvent the wheel if this has already been implemented in R, and I'm still new to the language so I'm not a fast coder. If anyone can help it would be greatly appreciated. Cheers! Luk Arbuckle [1] Kraemer, Helena. “Extension of the Kappa Coefficient.” Biometrics 36, no. 2 (1980). http://dx.doi.org/10.2307/2529972. [[alternative HTML version deleted]]
Luk Arbuckle
2012-Feb-01 20:13 UTC
[R] Function to compute multi-response, multi-rater kappa?
I'm very sorry for double posting! My r-help setting Receive your own posts to the list? is set to Yes, and Mail delivery is Enabled. Yet I did not get a copy of my post (this message is a reply from my sent mail). I only learned of the double posting when I found it copied in an r-help archive. Again, my apologies. Luk Arbuckle On Wed, Feb 1, 2012 at 13:47, Luk Arbuckle wrote:> I'm looking for a function in R that extends kappa to multiple raters when > there is more than one response per subject. For example, say a group of > doctors have to assign diseases to patients. Each patient will be assigned > one to many diseases, and the number of doctors assigning diseases to any > one patient will be two to many. > > Here's an extremely simple example of the type of data I might have (two > patients, three doctors, five diagnoses, and doctor 'x' assigns multiple > diagnoses per patient): > > pat<-c('a','a','a','b','b','b','b') > doc<-c('x','x','y','x','x','y','z') > dx<-c('1','2','2','3','4','4','5') > df<-data.frame(pat=pat,doc=doc,dx=dx) > df > > I found a paper that can address this [1], although I can't find it > referenced on cran. What I have found on cran doesn't address the > multi-response case. I don't want to reinvent the wheel if this has > already been implemented in R, and I'm still new to the language so I'm not > a fast coder. If anyone can help it would be greatly appreciated. > > Cheers! > > Luk Arbuckle > > [1] Kraemer, Helena. “Extension of the Kappa Coefficient.” Biometrics 36, > no. 2 (1980). http://dx.doi.org/10.2307/2529972. >[[alternative HTML version deleted]]