[R-br] teste kappa por categorias

Éder Comunello comunello.eder em gmail.com
Quinta Março 6 14:43:29 BRT 2014


Prezado 'branco.sk8',

Sugiro que você envie um *dataset *básico para apoiar sua questão. Um
código mínimo mostrando o que você já conseguiu também é bem-vindo. Feito
isso, acredito que será mais fácil entender o que você precisa e conseguir
ajuda...

Talvez o código abaixo pode lhe servir como referência:

### <code r>
if (!require(fmsb)) {install.packages('fmsb', dep=T); require(fmsb)}
if (!require(vcd))  {install.packages('vcd', dep=T);  require(vcd)}

data <- data.frame(
          id=rep(c('MALE','FEMALE'),each=20),
          resp=c("N", "S", "N", "N", "S", "S", "S", "S", "S", "N", "S", "S",
                 "S", "S", "SR", "SR", "N", "SR", "S", "S", "S", "N", "N",
"S",
                 "S", "N", "S", "S", "S", "S", "S", "N", "S", "S", "S",
"N", "S",
                 "SR", "S", "SR"))
table(data$id)

MALE   <- data$resp[data$id=='MALE']
FEMALE <- data$resp[data$id=='FEMALE']
tab <- table(FEMALE,MALE); tab

tmp <- rbind(tab, FEMALE=margin.table(tab, 1))
tabtot <- cbind(tmp, MALE=margin.table(tmp, 1)); tabtot

Kappa(tab)
Kappa.test(tab)
### </code>





Éder Comunello <c <comunello.eder em gmail.com>omunello.eder em gmail.com>
Dourados, MS - [22 16.5'S, 54 49'W]
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140306/8209bc88/attachment.html>


Mais detalhes sobre a lista de discussão R-br