wkappa | R Documentation |
Computes a weighted Kappa for agreement in the case of 2 raters. The diagnosis (the object of the rating) may have k possible ordered values.
wkappa(r,weights="squared")
r |
n*2 matrix or dataframe, n subjects and 2 raters |
weights |
weights="squared" to obtain squared weights. If not, absolute weights are computed (see details) |
Diagnoses have to be coded by numbers (ordered naturally). For weigths="squared", weights are related to squared differences between rows and columns indices (in this situation wkappa is close to an icc). For weights!="squared", weights are related to absolute values of differences between rows and columns indices. The function is supposed to deal with the case where the two raters have not exactly the same scope of rating. Missing value are omitted.
A list with :
$table |
the 2*k table of raw data (first rater in rows, second rater in columns) |
$weights |
"squared" or "absolute" |
$kappa |
Weighted Kappa |
Bruno Falissard
Cohen, J. Weighted kappa: nominal scale agreement with provision for scaled disagreement or partial credit. Psychological Bulletin 70 (1968): 213-220.
data(expsy) wkappa(expsy[,c(11,13)]) # weighted kappa (squared weights) #to obtain a 95%confidence interval: #library(boot) #wkappa.boot <- function(data,x) {wkappa(data[x,])[[3]]} #res <- boot(expsy[,c(11,13)],wkappa.boot,1000) #quantile(res$t,c(0.025,0.975)) # two-sided bootstrapped confidence interval of weighted kappa #boot.ci(res,type="bca") # adjusted bootstrap percentile (BCa) confidence interval (better)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.