Fernando, agradeço a ajuda, mas infelizmente já tentei com table, xtabs, ftable, e clogit e não consegui o que preciso. Os resultados não são contagens de pares. Tenho certeza de que estou fazendo algo errado.

Devido à minha dificuldade, recorri à lista.

Não postei as tentativas erradas que fiz, para deixar o script mais limpo.

Abraço,

Maurício Cardeal

UFBA

help(xtabs)

2017-04-07 12:59 GMT-05:00 Mauricio Cardeal via R-br <r-br@listas.c3sl.ufpr.br>:

Boa tarde!

Por favor, como faço para construir uma tabela como a do modelo do final do script,

onde o conteúdo das 4 caselas são contagens do número de pares que combinam

melhora e não melhora de 2 tratamentos aplicados simultaneamente no mesmo

individuo (perna direita e perna esquerda), por exemplo ?

# CMR:

par <- c(1,1,2,2,3,3,4,4,5,5,6,6,7,7)  # refere-se ao número do individuo
trata <-  c(1,0,1,0,1,0,1,0,1,0,1,0,1,0) # tratamento 1 (A) ou 0 (B)
melhora <-  c(1,0,1,0,1,0,0,1,0,1,0,0,1,1) # resultado do tratamento 1 (melhora) 0 (não melhora)
tratac <- factor(trata)
levels(tratac)<-list("A"=1,"B"=0 )
melhorac <- factor(melhora)
levels(melhorac)<-list("+"=1,"-"=0)

data.frame(par,tratac,melhorac)

   par tratac melhorac
1    1      A        +
2    1      B        -
3    2      A        +
4    2      B        -
5    3      A        +
6    3      B        -
7    4      A        -
8    4      B        +
9    5      A        -
10   5      B        +
11   6      A        -
12   6      B        -
13   7      A        +
14   7      B        +

# COMO FAÇO PARA GERAR ESSA TABELA A PARTIR DOS DADOS ACIMA?

#                      tratamento B
#                      melhora melhora
#                         +       -
#tratamento A melhora +   1       3
#             melhora -   2       1

Obrigado,

Maurício Cardeal

UFBA


_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e fornea cdigo mnimo reproduzvel.