私の因子「フォーム」の複数のレベルで、複数の評価者の異人種間合意(Fleiss Kappa)を計算したいデータフレームがあります。コードを実行すると、次のエラーが発生します。
**Error in factor(ratings[i, ], levels = lev) : subscript out of bounds**
これがどのように見えるかのサンプルです:
form <- c("B01","B01","B01","B02","B02","B02","B03","B03","B03")
word<-c("arch","avenue","concrete","devotion","equations","moth","piling","relate","cover")
B01200<-c(1,2,3,NA,NA,NA,NA,NA,NA)
B01209<-c(3,2,1,NA,NA,NA,NA,NA,NA)
B01214<-c(1,3,2,NA,NA,NA,NA,NA,NA)
B01228<-c(NA,NA,NA,NA,NA,NA, 2,2,1)
B01235<-c(NA,NA,NA,NA,NA,NA ,1,2,3)
B01242<-c(NA,NA,NA,NA,NA,NA, 3,3,1)
B02215<-c(NA,NA,NA ,2,1,1,NA,NA,NA)
B02217<-c(NA,NA,NA,2,1,1,NA,NA,NA)
B02222<-c(NA,NA,NA, 3,1,3,NA,NA,NA)
x <- data.frame(form,word,B01200,B01209,B01214,B01228,B01235,B01242,B02215,B02217,B02222)
x$form <-factor(x$form)
by(x,x[,"form"], function(x) kappam.fleiss(x[,-c(1:2)], detail=TRUE))
str(x)を使用して基礎となる構造を見ると、「form」変数の問題が何であるかがわかりません。
表示される内容は次のとおりです。
'data.frame': 9 obs. of 11 variables:
$ form : Factor w/ 3 levels "B01","B02","B03": 1 1 1 2 2 2 3 3 3
$ word : Factor w/ 9 levels "architects","avenue",..: 1 2 3 5 6 7 8 9 4
$ B01200: num 1 2 3 NA NA NA NA NA NA
$ B01209: num 3 2 1 NA NA NA NA NA NA
$ B01214: num 1 3 2 NA NA NA NA NA NA
$ B01228: num NA NA NA NA NA NA 2 2 1
$ B01235: num NA NA NA NA NA NA 1 2 3
$ B01242: num NA NA NA NA NA NA 3 3 1
$ B02215: num NA NA NA 2 1 1 NA NA NA
$ B02217: num NA NA NA 2 1 1 NA NA NA
$ B02222: num NA NA NA 3 1 3 NA NA NA
私は何が間違っているのですか?
ありがとう!