Exemplo:
Os dados do estudo Kawasaki encontram-se num conjunto de dados SAS com 167 observações (uma para cada criança) e três variáveis, um número de identificação, braço de tratamento (GG ou ASA), e uma variável indicadora de qualquer anomalia CA na visita 3 ou visita 4.
O ORDER= Opção
A declaração PROC FREQ tem uma opção que define a ordem em que os valores aparecem nas frequências e nas etiquetas cruzadas geradas por PROC FREQ.
O padrão é ORDER=INTERNAL, o que significa que os dados são ordenados (por ordem alfabética ou numérica) pelos valores não-formatados dos dados. Por exemplo, a variável ARM no exemplo acima assume um valor de ‘ASA’ ou ‘GG,’ e assim, por defeito, os valores ASA aparecerão antes dos valores GG na saída PROC FREQ.
A opção ORDER=FORMATTED ordenará os dados por valores formatados (ascendentes) das variáveis. Os impactos de outras opções ORDER= são dados no final deste módulo.
Formatando o resultado de modo a que o evento esteja na primeira coluna
Usando o formato abaixo, uma vez que “E” vem antes de “N” alfabeticamente, “Event” estará na coluna 1 e “No event” na coluna 2. No entanto, a ASA estará na linha 1, uma vez que a ASA está formatada “0-Aspirina” e a GG está formatada “1-Gamma Globulina”.
formatoproc;
valor $armf “ASA”=”0-Aspirina” “GG”=”1-Gama-Globulina”;
value eventf 0=’Nenhum evento’ 1=’Evento’;
run;
proc freq data=d.kawa; order=formatted;
format arm $armf. anyv34 eventf.;
tabela braço*anyv34;
correr;
Outras opções
Podemos manter a inclusão de uma declaração de formato em cada proc, mas em vez disso, vamos formatá-las num passo de dados.
data one;set d.kawa;
format arm $armf. anyv34 eventf.;
Existem várias opções que podem ser incluídas após a / na declaração da TABELA.
- As opções norow, nocol,e nopercent limitam o número de entradas na tabela.
- A opção de medidas estima o odds ratio e o risco relativo com os intervalos de confiança que os acompanham.
- A opção chisq solicita o teste do qui-quadrado. A opção esperada pede que as frequências celulares esperadas sejam incluídas nas células. Um aviso exibido na saída se mais de 20% das células tiverem contagens esperadas de menos de 5.
Suprimir a Coluna e Percentagem Global
proc freq data=uma ordem=formatada;
braço de tabela*anyv34 / nocol nopercent;
run;
Incluindo Frequências Esperadas
proc freq data=uma encomenda=formatada;
tabela braço*anyv34 / esperado ;
correr;
Requerendo o Teste do Quadrado do Chi
proc freq data=uma ordem=formatada;
tabela braço*anyv34 / chisq ;
run;
A tabela 2 x 2 é produzida como acima, mais a seguinte saída.
A linha realçada contém a estatística do quadrado chi e a sua associada p-value
br>
Nota: Se > 20% das frequências celulares são <5, SAS imprimirá um aviso, e não deve utilizar o teste do qui-quadrado. Em vez disso, utilize o Teste Exacto de Dois Lados de Fisher (impresso por defeito quando a tabela é 2 x 2).
0 comentários