Apostila Introdução Ao Sas.pdf 1j6r1y

INTRODUÇÃO AO

SAS

Euclides Braga MALHEIROS Liliane Karla FIGUEIRA

FCAV/UNESP Campus de Jaboticabal

ii

ÍNDICE I. APRESENTAÇÃO........................................................................................................1 I.1. Módulos do SAS......................................................................................................1 I.2. Janelas do SAS ........................................................................................................1 I.3. Barra de ferramentas .............................................................................................4 I.4. Help do SAS.............................................................................................................5 II. PROGRAMAÇÃO SAS ...............................................................................................6 II.1. Alguns comandos do SAS .....................................................................................6 OPTIONS ...................................................................................................................6 TITLE .........................................................................................................................7 FOOTNOTE ...............................................................................................................7 COMENTÁRIOS .......................................................................................................7 DATA .........................................................................................................................7 LIBNAME ..................................................................................................................8 FILENAME ................................................................................................................8 RUN............................................................................................................................8 COMANDOS ARITMÉTICOS..................................................................................8 Constantes..............................................................................................................8 Variáveis ................................................................................................................8 Operadores .............................................................................................................8 Aritméticos .......................................................................................................9 Lógicos..............................................................................................................9 De comparação..................................................................................................9 Funções ..................................................................................................................9 II.2. Algumas estruturas usadas na programação....................................................10 IF ... THEN ... ELSE ................................................................................................10 DO ... END ..............................................................................................................10 II.3. Algumas formas para criar um SAS-DATA-SET ............................................11 Quando os dados são colocados no próprio programa .............................................11 Quando os dados estão gravados como um arquivo ASCII......................................11 Comando INPUT ...............................................................................................11 Listado.............................................................................................................11 Colunado .........................................................................................................11 Formatado .......................................................................................................11

iii

II.4. Procedimentos do SAS ........................................................................................12 Estrutura básica dos procedimentos .........................................................................12 Comandos que podem ser usados nos procedimentos..............................................13 VAR.....................................................................................................................13 BY........................................................................................................................13 CLASS .................................................................................................................13 ID .........................................................................................................................13 TITLE ..................................................................................................................13 FOOTNOTE ........................................................................................................13 WHERE ...............................................................................................................13 OUTPUT OUT ....................................................................................................13 III. Tutorial do SAS .........................................................................................................14 IV. Programas Exemplos.................................................................................................14 V. INFORMAÇÕES (Procedimentos e comandos) separadas por aula .....................15 V.1. Informações necessárias para a aula 1.............................................................15 Comandos .................................................................................................................15 DROP...................................................................................................................15 KEEP ...................................................................................................................15 LABEL.................................................................................................................15 Procedimentos ..........................................................................................................15 PROC PRINT ......................................................................................................15 PROC SORT........................................................................................................15 PROC MEANS ....................................................................................................16 Como criar um arquivo ASCII a partir do Excel ......................................................16 Lista de Exercícios – Aula 1...................................................................................17 V.2. Informações necessárias para a aula 2.............................................................19 Comandos .................................................................................................................19 SET ......................................................................................................................19 MERGE ...............................................................................................................19 INFILE .................................................................................................................20 Lista de Exercícios – Aula 2...................................................................................21 V.3. Informações necessárias para a aula 3.............................................................22 Procedimentos ..........................................................................................................22 FREQ ...................................................................................................................22 TTEST .................................................................................................................22 Lista de Exercícios – Aula 3...................................................................................23 V.4. Informações necessárias para a aula 4.............................................................24

iv

Procedimentos ..........................................................................................................24 UNIVARIATE.....................................................................................................24 INSIGHT .............................................................................................................24 Lista de Exercícios – Aula 4...................................................................................25 V.5. Informações necessárias para a aula 5.............................................................26 Gráficos de dispersão ...............................................................................................26 PROC PLOT ........................................................................................................26 PROC GPLOT .....................................................................................................27 Outros gráficos .........................................................................................................27 Histogramas (PROC GCHART)..........................................................................27 Gráficos de 3 dimensões (PROC G3D) ...............................................................28 Lista de Exercícios – Aula 5...................................................................................29 V.6. Informações necessárias para a aula 6.............................................................30 Regressão Linear .....................................................................................................30 PROC REG..........................................................................................................30 Lista de Exercícios – Aula 6...................................................................................31 V.7. Informações necessárias para as aulas 7 e 8....................................................32 PROC ANOVA .......................................................................................................32 CLASS .................................................................................................................32 MODEL ...............................................................................................................32 MEANS ...............................................................................................................33 TEST H=... E=... ...............................................................................................33 PROC GLM .............................................................................................................33 LSMEANS...........................................................................................................33 CONTRAST ........................................................................................................34 RANDOM............................................................................................................34 Lista de Exercícios – Aula 7...................................................................................35 V.8. Informações necessárias para a aula 8 .............................................................37 Considerações estatísticas .......................................................................................37 Esperanças dos quadrados médios .......................................................................37 Regras práticas para obtenção das esperanças dos quadrados médios E(QM) dados balanceados................................................................................................38 Procedimentos para estimativas de componentes da variância ................................39 PPROC VARCOMP ............................................................................................39 Lista de Exercícios – Aula 8...................................................................................40

v

V.9. Informações necessárias para a aula 9 .............................................................41 Procedimentos para estimar componentes da variância ..........................................41 PROC NPAR1WAY............................................................................................41 Lista de Exercícios – Aula 9...................................................................................42 V.10. Informações necessárias para a aula 10 .........................................................43 Procedimento para Regressão Não Linear ...............................................................43 PROC NLIN ........................................................................................................43 Lista de Exercícios – Aula 10.................................................................................44

VI. RESOLUÇÃO DAS LISTAS DE EXERCÍCIOS...................................................45 VI.1. Resolução da Lista de Exercícios – Aula 1.......................................................45 VI.2. Resolução da Lista de Exercícios – Aula 2.......................................................47 VI.3. Resolução da Lista de Exercícios – Aula 3.......................................................49 VI.4. Resolução da Lista de Exercícios – Aula 4.......................................................51 VI.5. Resolução da Lista de Exercícios – Aula 5.......................................................52 VI.6. Resolução da Lista de Exercícios – Aula 6.......................................................53 VI.7. Resolução da Lista de Exercícios – Aula 7.......................................................55 VI.8. Resolução da Lista de Exercícios – Aula 8.......................................................58 VI.9. Resolução da Lista de Exercícios – Aula 9.......................................................59 VI.10. Resolução da Lista de Exercícios – Aula 10...................................................60

INTRODUÇÃO AO SAS Euclides Braga MALHEIROS1 Liliane Karla FIGUEIRA2

I. Apresentação O SAS é um sistema composto por vários módulos, com aplicações diversas. Além do módulo para análises estatísticas, o SAS/STAT, vários outros módulos fazem parte do sistema.

I.1. Módulos do SAS ♦ ♦ ♦ ♦ ♦ ♦ ♦ ♦ ♦ ♦ ♦ ♦ ♦ ♦ ♦

SAS/BASE - Procedimentos básicos SAS/ACCESS - Banco de Dados SAS/AF - Gerador de Aplicativos SAS/CALC - Análise e modelagem financeira/Planilha tridimensional SAS/CONNECT - Conectividade do SAS com diferentes ambientes SAS/ETS – Análises econômicas e de séries temporais SAS/FSP - Processador de textos SAS/GIS - Geoprocessamento SAS/GRAPH - Gráfico SAS/IML - Linguagem matricial interativa SAS/INSIGHT – Módulo interativo do tipo “point and click” SAS/OR - Pesquisa Operacional SAS/QC - Controle de qualidade SAS/STAT - Estatística SAS/SHARE – Sistema multiusuário

I.2. Janelas do SAS O SAS tem a característica de trabalhar simultaneamente com 3 janelas, com finalidades específicas, que são:

1

Departamento de Ciências Exatas - FCAV/UNESP - Campus de Jaboticabal. Aluna do Curso de Pós-Graduação em Agronomia, área de Concentração em Entomologia Agrícola, em nível de doutorado - FCAV/UNESP - Campus de Jaboticabal.

2

2

♦ EDITOR DE PROGRAMAS “PGM” → Janela onde o usuário edita os programas. Quando o programa é processado (RUN), o SAS gera o conteúdo das janelas LOG e OUTPUT. ♦ LOG → Janela onde é apresentada uma depuração do programa, ou seja, informações sobre o processamento como: tempo gasto; listagem de erros, se houver; etc. (gerada pelo SAS). ♦ OUTPUT → Janela para a saída dos resultados (gerada pelo SAS). A tela de abertura do SAS, quando as janelas são apresentadas em colunas, é mostrada na Figura 1.

Figura 1. Tela de abertura do SAS.

A opção WINDOW, na barra de menus, permite minimizar as janelas ou apresentá-las de diferentes formas, ou seja: em cascata (Cascade), em forma de ladrilhos em colunas (Tile) ou reorganiza no formato padrão do SAS (Resize). Permite também ar cada uma das janelas. Ver a Figura 2.

Figura 2. Opção WINDOW da Barra de Menus.

3

os às janelas, bem como outras operações, podem ser feitos utilizando-se as Teclas de Atalhos, como: CTRL + E → Limpa a Janela. F1 → HELP. F4 → Recarrega o último programa rodado (Janela PGM). F9 → Janela com as definições das funções do teclado (KEYS). Pode ser usado para redefinir essas funções. F12 → Limpa a janela (não é padrão do SAS, precisa ser definido na janela KEYS). Essas e outras Teclas de Atalho podem ser vistas na Janela KEYS. Ver a Figura 3:

Figura 3. Teclas de Atalho do SAS.

Toda vez que o Programa existente na janela PGM é rodado (“RUN”), são geradas linhas nas janelas LOG e OUTPUT, conforme apresentado na Figura 4.

LOG PGM

RUN

OUTPUT Figura 4. Esquema usado pelo SAS para gerar as Janelas LOG e OUTPUT.

4

É interessante destacar que: 1. Ao rodar um programa, a) caso as janelas LOG e OUTPUT não estejam limpas, os resultados serão acumulados nos conteúdos dessas janelas. b) A janela PGM é limpa automaticamente. Assim sendo, toda vez que um programa for rodado e precisar ser editado, para a correção de algum erro ou para a inclusão de outros comandos, deve-se: ar a janela OUTPUT e apagá-la .............................................(F7 / F12), ar a janela LOG e apagá-la .....................................................(F6 / F12), ar a janela PGM e recarregar o último programa rodado ..........(F5 / F4). 2. O conteúdo de cada janela é salvo individualmente. O SAS usa as extensões: .SAS para os programas (janela PGM), .LOG para o conteúdo da janela LOG, .LST para o conteúdo do OUTPUT .

I.3. Barra de ferramentas A maioria das ferramentas disponíveis na Barra de Ferramentas do SAS (Figura 5) são conhecidas pelos usuários do Microsoft Word.

HELP SAS/ASSISTENTE PASTAS PROMPT DO DOS DESFAZER COLAR COPIAR RECORTAR IMPRIMIR ARQUIVO NOVO ABRIR ARQUIVO RODAR (RUN) COMANDOS DE OPERAÇÃO Figura 5. Barra de Ferramentas do SAS.

5

Na ferramenta “Comandos de Operação” são permitidos comandos de operações nas janelas como: ♦ CLEAR: Limpa a janela , ♦ TOP: Vai para o início da janela, ♦ BOTTOM: Vai para o fim da janela, ♦ LEFT: Rola a janela à esquerda, ♦ RIGHT: Rola a janela à direita, ♦ RUN: Roda o programa (só na janela PGM), ♦ BYE ou ENDSAS: Sai do SAS, ♦ Nomes de Janelas, para á-las, como: LOG, OUT, PGM, KEYS, HELP (a janela KEYS é a que contem as teclas de atalho), ♦ etc. À medida que vão sendo usados comandos na ferramenta “Comandos de Operação”, o SAS vai criando uma lista que, ao repetir um comando já usado basta selecioná-lo. A lista é obtida clicando na seta à direita da ferramenta. Veja o exemplo na Figura 6.

Figura 6. Lista de comandos SAS.

I.4. Help do SAS O SAS possui um Auxílio completo e interativo. Para utilizá-lo, basta ar a janela HELP pela tecla F1, ou pela Barra de Ferramentas. Os procedimentos do SAS são encontrados em: SAS 6.11: HELP / INDEX. SAS 6.12: HELP / Extended help / SAS System help: Main Menu / Limited Index.

6

II. Programação SAS Os programas SAS são escritos na janela PGM, em linhas de 256 caracteres. Os comandos de programação terminam com ponto e vírgula (;), podendo ter tantos comandos por linha quantos se desejar, desde que respeitado o tamanho da linha. Os programas SAS possuem duas fases: a Fase 1, que constitui na criação do arquivo SAS (SAS-DATA-SET) e a Fase 2, que inclui os Procedimentos do SAS (PROC) para trabalhar com os dados (impressão, representações gráficas, análises, etc). Veja o esquema da Figura 7. Janelas como

Fase 1

SAS-DATA-SET

Fase 2

OUTPUT, GRAPHICS, etc Ou

ARQUIVOS DE DADOS Figura 7. Fases de um programa SAS.

O SAS-DATA-SET é o arquivo de dados a ser usado pelo SAS. Ele é formado por linhas e colunas, sendo nas colunas as variáveis e nas linhas os registros. Exemplo: Nome Maria José Rosa João

Endereço Rua 8, 567 Rua Tiradentes, 123 Av. Marajás, 750 Rua Poty, 1250

Cidade Itu Matão Jaboticabal Maringá

Estado SP SP SP PR

Peso 59.50 85.00 . 79.00

Observação perdida Ponto para separar as casas decimais

II.1. Alguns comandos do SAS OPTIONS: Define, dentre muitas opções, parâmetros para configurar a página do OUTPUT. Ex: OPTIONS LS=80 PS=66 NOCENTER NODATE; LS = Line Size e PS = Page Size

7

Obs: a) Os parâmetros não definidos assumem os valores “Default”. b) Os valores dos parâmetros definidos com o comando OPTIONS, só são válidos para o programa em questão. c) Para mudar o “Default” use: GLOBAL / OPTIONS / GLOBAL OPTIONS. TITLE: Define a n-ésima (1 ≤ n ≤ 10) linha do cabeçalho para o OUTPUT. sintaxe: TITLE[n] < título >; Ex: TITLE1 Análise do Experimento de Fertilidade; TITLE3 Coleta do Ano 1998; Obs: a) Se for omitido o “n” o 1o TITLE define a 1ª linha do Cabeçalho, o 2o a 2ª linha, e assim por diante. b) Se o cabeçalho não for definido, o SAS usa o definido como “Default”. c) Para mudar o “Default” use: GLOBAL / OPTIONS / TITLE. FOOTNOTE: Define a n-ésima (1 ≤ n ≤ 10) linha do rodapé para o OUTPUT. Uso equivalente ao uso do TITLE. COMENTÁRIOS: Usados na documentação dos programas. sintaxe: * ; ou: /* */ Obs: a) Os comentários não são executados. b) Os comentários podem ser colocados em qualquer parte do programa. c) Pode-se usar tantos comentários quantos se desejar. DATA: Atribui nome ao SAS_DATA_SET. sintaxe: DATA [< nome do SAS_DATA_SET >]; Obs: a) O nome do SAS_DATA_SET pode ser: simples – quando apenas o nome do arquivo é especificado (Ex: NNNN). Neste caso será criado um arquivo temporário (até fechar o SAS), de nome NNNN, na pasta SASWORK, também temporária. composto – quando o nome tem duas partes, separadas por “ponto” (Ex: CCCC.NNNN). Neste caso será criado um arquivo permanente, de nome NNNN, no caminho nomeado por CCCC. b) O nome do SDS deve ter de 1 a 8 caracteres alfanuméricos (alfabéticos ou numéricos) e o primeiro deve ser alfabético. c) Se o nome for omitido, o SAS usa DATA1 para o 1o, DATA2 para o 2o, e assim por diante.

8

LIBNAME: Atribui nome a um caminho. Ex: LIBNAME CAM “C:\SAS”; Com este comando, o caminho C:\SAS é nomeado por CAM. FILENAME: Atribui nome a um arquivo. Ex: FILENAME ARQ “C:\SAS\DADOS.ASC”; Com este comando, o arquivo C:\SAS\DADOS.ASC é nomeado por ARQ. RUN: Determina ao SAS para Rodar (Executar) o programa até este ponto. COMANDOS ARITMÉTICOS: Usados para criar novas variáveis. Sintaxe: =<EXPRESSÃO>; Ex: PORCG=PGORD/PLEITE*100; As expressões podem conter: Constantes, Variáveis, Operadores e Funções. Constantes: Podem ser: Numéricas, Caracter, Data e Hora. Numéricas: Sistema decimal no modo padrão (Ex: -125.45) ou em notação científica (Ex: -1.456E-03); sistema octal (Ex: 0B123X). Caracter: Até 200 caracteres entre aspas (Ex: “Femea”). Data: Data seguida por D, e entre aspas (Ex: “25/10/98D”). Hora: Hora seguida por T, e entre aspas (Ex: “10:15:45T”). Para maiores detalhes consultar o HELP, usando a seqüência: SAS 6.11: HELP / SAS LANGUAGE / SAS EXPRESSIONS / SAS CONSTANTS SAS 6.12: HELP / Extended help / SAS System help: Main Menu / SAS Language / SAS Expressions / SAS Constants. Variáveis: As variáveis podem ser numéricas ou caracter. Os nomes das variáveis devem ter de 1 a 8 caracteres alfanuméricos (alfabéticos ou numéricos) e o primeiro deve ser alfabético. Operadores: Os mais usuais são os operadores aritméticos, lógicos e de comparação. Operadores aritméticos:

9

Operador + * / **

Operação Adição Subtração ou Negação Multiplicação Divisão Potenciação

Exemplos Usual SAS A+B A+B A-B ou -B A-B ou -B A.B A*B A:B A/B AB A**B

Operadores lógicos: Operador & ou AND | ou OR ^ ou NOT

Operação Negação

Exemplos SAS A&B ou A AND B A|B ou A OR B ^A ou NOT A

Usual AeB AouB Não A

Operadores de comparação: Operador = ou EQ > ou GT < ou LT >= ou GE <= ou LE ^= ou NE

Operação Igual a Maior que Menor que Maior ou igual que Menor ou igual que Diferente de

Usual A=B A>B A
Exemplos SAS A=B ou A EQ B A>B ou A GT B A =B ou A GE B A<=B ou A LE B A^=B ou A NE B

Para maiores detalhes consultar o HELP, usando a seqüência: 6.11: HELP / Extended help / SAS Language / SAS Expressions / SAS Operator. 6.12: HELP / Extended help / SAS System help: Main Menu / SAS Language / SAS Expressions / SAS Operator. Funções: As funções SAS são escritas na forma: Função(lista de argumentos separados por vírgulas) Existe no SAS um número muito grande de funções disponíveis, como exemplos: SQRT(X) – calcula a raiz quadrada de X, MDY(MES,DIA,ANO) – resulta uma data correspondente ao Dia, Mês e Ano especificados como argumentos. Obs: Para armazenar datas o SAS armazena um número que corresponde ao total de dias entre a data e o 1º dia do século.

10

Para maiores detalhes consultar o HELP, usando a seqüência: 6.11: HELP / Extended help / SAS Language / SAS Functions. 6.12: HELP / Extended help / SAS System help: Main Menu / SAS Language / SAS Functions.

II.2. Algumas estruturas usadas na programação IF ... THEN ... ELSE: Usado para testar se uma condição é verdadeira ou falsa. Sintaxe: IF THEN ; [ELSE < comando>;] Ex: IF MES<7 THEN SEMESTRE=”PRIMEIRO”; ELSE SEMESTRE=”SEGUNDO”; Obs: Se o ELSE for omitido, o SAS executa o próximo comando na seqüência.

DO ... END: Usado para controle de “loops”, ou seja, para executar a mesma seqüência de comandos um número finito de vezes. Sintaxe: DO x=a TO b [BY c]; Comando 1; Comando 2; .... END; Ex: DATA EXDO; DO X=5 TO 30 BY 5; X2=X**2; X3=X**3; OUTPUT; END; Com esse “looping”, teremos no SAS-DATA-SET

11

Obs . 1 2 3 4 5 6

X 5 10 15 20 25 30

X2 25 100 225 400 625 900

X3 125 1.000 3.375 8.000 15.625 27.000

II.3 Algumas formas para criar o SAS-DATA-SET Quando os dados são colocados no próprio programa. Seqüência: DATA <sds>; INPUT <lista de variáveis>; CARDS; ;

Define o nome do SAS-DATA-SET. Informa as variáveis a serem lidas. Informa que iniciarão os dados. Informa que terminaram os dados.

Quando os dados estão gravados num arquivo ASCII (American Standard Code – Interchange Information). Seqüência: DATA <sds>; INFILE “<arquivo ASC>”; INPUT <lista de variáveis>;

Define o nome do SAS-DATA-SET. Arquivo ASC a ser importado. Lista de variáveis importadas.

Comando INPUT O comando INPUT relaciona os nomes das variáveis a serem criadas no SAS-DATA-SET. Ele pode ser de uma das formas: listado, colunado ou formatado. Listado: quando se especifica no INPUT apenas a lista de variáveis (nomes separados por pelo menos um espaço em branco). Usado quando os dados a serem lidos estão separados por, pelo menos, um espaço em branco. Ex: INPUT NOME$ PESO ALTURA; Colunado: quando se especifica a variável e a coluna inicial e final onde se encontram os dados a serem lidos. Ex: INPUT NOME$ 1-20 PESO 20-25 ALTURA 25-30; Formatado: quando se especifica a variável e formato que o dado a ser lido deve ser considerado.

12

Ex: INPUT NOME$ 1-20 DATANAS DDMMYY6. PESONAS 8.2 ♦ DDMMYY6. especifica que a data a ser lida contém 2 caracteres para o dia, 2 para o mês e 2 para o ano, num total de 6. ♦ 8.2 especifica que o peso a ser lido tem 8 casas, sendo 2 decimais. Obs.: Quando se usa o INPUT, o número de valores por linha de dados deve coincidir com o número de variáveis do INPUT, a menos que a lista de variáveis seja seguida por @@. Exemplo: DATA EXEMPLO; INPUT X Y; CARDS; 10 20 12 22 15 35 ;

DATA EXEMPLO; INPUT X Y@@; CARDS; 10 20 12 22 15 35 ;

Para maiores detalhes consultar o HELP, usando a seqüência: 6.11: HELP / Extended Help / SAS Language / SAS Data Step / Syntax / Input. 6.12: HELP / Extended help / SAS System help: Main Menu / SAS Language / SAS DATA step / Syntax / INPUT.

II.4. Procedimentos do SAS Existem vários Procedimentos, dentro de cada módulo, disponíveis a usuários de uma grande diversidade de áreas de aplicações. Estrutura básica dos procedimentos A 1ª linha de cada procedimento tem a seguinte estrutura:

PROC NOME < opções >;

Lista de Opções. Nome do Procedimento. Prefixo de todo Procedimento (PROCEDURE).

Os procedimentos podem incluir outras linhas para especificações especiais. A estrutura geral de um procedimento é:

13

PROC NOME < opções >; ...

Uma das opções é: DATA = < Arquivo SAS >; Especifica o arquivo SAS que será usado. Se for omitido, o SAS usa o último criado

VAR < lista de variáveis >; BY < lista de variáveis >; CLASS < lista de variáveis>; ID < lista de variáveis >; TITLE < “título” >; FOOTNOTE < “rodapé”>; WHERE < condição >;

Opcionais e, no geral, independem da ordem.

OUTPUT OUT= <SDS> = ;

...

RUN; Comandos que podem ser usados nos procedimentos VAR < lista de variáveis >; Especifica as variáveis para as quais será executado o procedimento. BY < lista de variáveis >; Especifica que o procedimento deve ser executado para cada valor da(s) variável(is) especificada(s). Este comando exige que o SAS-DATA-SET esteja classificado por essa(s) variavel(is). CLASS < lista de variáveis >; Especifica que o procedimento deve ser executado para cada valor da(s) variável(is) especificada(s). ID < lista de variáveis >; Especifica a variável de identificação dos registros. Se não for especificado, o SAS usa o número do Registro. TITLE < “título” >; Especifica o título para todas as páginas geradas por este procedimento. FOOTNOTE < “rodapé >; Especifica o rodapé para todas as páginas geradas por este procedimento. WHERE < condição >; Especifica que o procedimento deve ser executado somente para os registros para os quais a condição for satisfeita. OUTPUT OUT= <SDS> = ;

Coloca no SAS-DATA-SET “SDS”, além das variáveis do SAS-DATA-SET em uso, variáveis criadas pelo SAS “VS” (cujos nomes são criados pelo SAS e devem ser especificados corretamente), com os nomes das variáveis do usuário “VU”.

14

III. Tutorial do SAS O SAS possui um tutorial bastante interessante para os principiantes. São algumas aulas que permitem aos usuários uma familiarização com a terminologia, alguns conceitos básicos e a estrutura do Sistema. O Tutorial é um módulo à parte, que precisa ser selecionado durante a instalação do SAS. Para entrar no Tutorial siga os os: 1) Selecione a opção de menu: Help / Online training. Com esta opção, será apresentada uma lista dos cursos disponíveis. 2) Selecione o curso e a opção “Start Course”. Com esta opção, será apresentada uma lista de lições do curso selecionado, e as possíveis seções (Instructions - que corresponde às instruções da lição, Summary - que corresponde a um resumo da lição e Quiz - que é um teste de avaliação da lição). 3) Selecione a lição e a seção, e siga em frente.

IV. Programas Exemplos O SAS possui uma lista de programas prontos (programas exemplos) que contribuem em muito para o aprendizado. Para ar esta lista de programas, siga os os: 1) 2) 3) 4) 5) 6) 7)

Selecione a opção de menu: Help / Sample Programs. Selecione a opção: Conteúdo / SAS Sample Library. Selecione o módulo do qual deseja o exemplo. Selecione o programa. Marque o programa como um bloco e copie (jogue para a área de transferência). e a janela “Program Editor”, e nela cole o conteúdo da área de transferência. Rode o Programa (Run).

15

V. INFORMAÇÕES (Procedimentos e Comandos) separadas por aula V.1. INFORMAÇÕES NECESSÁRIAS PARA A AULA 1 COMANDOS DROP: especifica a lista de variáveis a serem eliminadas do SAS-DATA-SET. As variáveis não incluídas na lista serão mantidas. Sintaxe: DROP <lista de variáveis>; KEEP: especifica a lista de variáveis a serem mantidas do SAS-DATA-SET. As variáveis não incluídas na lista serão eliminadas. Sintaxe: KEEP <lista de variáveis>; LABEL: atribui label’s a variáveis. Sintaxe: LABEL =”label1” =”label2” ... ; Ex.: LABEL X=”DOSES DE N” Y=”PRODUÇÃO”; PROCEDIMENTOS

PROC CONTENTS Finalidade: Apresenta uma descrição do SAS-DATA-SET. Sintaxe: PROC CONTENTS < opções >;

PROC PRINT Finalidade: Imprimir o SAS-DATA-SET no OUTPUT. Sintaxe: PROC PRINT < opções >; Algumas opções são: NOOBS – Não imprime o número da observação (Registro). DOUBLE – Imprime com espaçamento duplo. ROUND – Imprime os números reais arredondados para 2 casas decimais. LABEL – Imprime tendo como cabeçalho das colunas os label’s no lugar dos nomes das variáveis.

PROC SORT Finalidade: Classificar o SAS-DATA-SET.

16

Sintaxe: PROC SORT < opções >; BY [DESCENDING] Var1 [DESCENDING] Var2 ...; Obs:. As variáveis precedidas por DESCENDING são classificadas em ordem decrescente, caso contrário, em ordem crescente. Algumas opções: NODUP – Elimina os registros para os quais o(s) valor(es) da(s) chave(s) de classificação forem iguais.

PROC MEANS Finalidade: Calcular estatísticas descritivas de variáveis numéricas. Sintaxe: PROC MEANS < opções > < lista de estatísticas >; Algumas opções: ALPHA=valor – Especifica o nível de significância para o intervalo de confianca. Algumas estatísticas: N – Número de observações. NMISS – Número de observações perdidas. MEAN - Média aritmética. VAR – Variância. STD - Desvio padrão. STDERR - Erro padrão da média. MIN - Valor mínimo. MAX - Valor máximo. CLM - Intervalo de confiança para a média. RANGE - Amplitude (Máximo-Mínimo). T - Teste t para a hipótese H0:média=0 vs H1:média≠0. PRT - Nível Mínimo de significância para o teste t. COMO CRIAR UM ARQUIVO ASCII, A PARTIR DO EXCEL a) b) c) d) e)

Digite o arquivo, uma variável em cada coluna, sem títulos das colunas. Preencha as caselas vazias com vírgulas “ ,”. Salve o arquivo especificando o tipo: Texto (OS/2 ou MS-DOS)(*.txt). Feche o EXCEL. Entre no EDIT do DOS. INICIAR / PROGRAMAS / PROMPT MS-DOS / EDIT: Abrir o Arquivo. Trocar as vírgulas por pontos. Salvar o arquivo. Sair do Edit. f) Sair do DOS (EXIT).

17

LISTA DE EXERCÍCIOS - AULA 1 1. Fazer um programa no SAS para: a) Colocar comentários informando: o número da lista, do exercício, a data e o nome do programador. b) Definir as opções para o OUTPUT: tamanho da linha=78, tamanho da página=64. c) Definir o Cabeçalho: CURSO SAS - ANÁLISE DE DADOS. d) Definir o Rodapé: FCAV/UNESP. d) Criar um arquivo SAS (SAS-DATA-SET) temporário, a partir dos dados da Tabela a seguir, usando a opção CARDS e digitando todas as variáveis. TRAT. T1 T2 T3

1 35 40 39

2 19 35 27

REPETIÇÕES 3 31 46 20

4 15 41 29

5 30 33 45

f) Imprimir o arquivo. g) Mostrar uma descrição do conteúdo do arquivo. 2. Fazer um programa SAS para: a) Colocar comentários informando: o número da lista, do exercício e a data. b) Definir as opções para o OUTPUT: tamanho da linha=78, tamanho da página=64. c) Criar um arquivo SAS, a partir do arquivo L1E2.DAT, que é um arquivo ASCII, e está codificado como segue: COLUNAS CAMPO 1-4 vaca 5-8 paiv 9-12 maev 13-14 diap 15-16 mesp 17-18 anop 19-20 orparto 21 sexo 22-23 anod 24-25 mesd 26-27 diad 28-33 pesonasc 34-39 pesodesm Obs: Sexo = M ou F e Pesos com duas casas decimais d) Imprimir o arquivo.

18

e) Criar 4 novas variáveis: Idade à desmama, Ganho de peso total até a desmama, Ganho de peso diário até a desmama, Período do Nascimento (ÁGUAS out-mar, SECA abr-set). f) Eliminar do arquivo SAS as variáveis utilizadas apenas como valores intermediários. g) Classificar o arquivo pela variável sexo. h) Classificar o arquivo pelas variáveis: sexo e período do nascimento. 3. Fazer um programa no SAS para: a) Colocar comentários informando: o número da lista, do exercício, a data e o nome do programador. b) Definir as opções para o OUTPUT: tamanho da linha=78, tamanho da página=64. c) Criar um arquivo SAS, permanente, na pasta C:\SAS, a partir do arquivo L1E3.DAT, no drive A, que é um arquivo ASCII, e contém 5 variáveis: Raças, Origem, Repetições, Peso ao nascer (PN) e Ganho de peso em 30 dias (GP30). d) Estabelecer “label’s” associados às variáveis. e) Imprimir o arquivo. f) Salvar o Programa. 4. Fazer um programa no SAS para, usando o arquivo SAS criado anteriormente, a) Calcular a média, variância e desvio padrão para as variáveis PN e GP30. b) Calcular as médias para PN e GP30, por raça. c) Calcular as médias para PN e GP30, por raça e origem. 5. Fazer um programa no SAS para criar um Arquivo SAS, a partir dos dados do Exercício 1, considerando-se que a segunda repetição do tratamento 1 tenha sido perdida, e digitando apenas os valores observados (Usar comandos DO ... END). 6. Fazer um programa no SAS para criar um Arquivo SAS, a partir dos dados do Exercício 1, digitando os dados no EXCEL.

19

V.2. INFORMAÇÕES NECESSÁRIAS PARA A AULA 2 Comandos SET: Carrega os registros (linhas) no SAS-DAT-SET. Sintaxe: SET < lista de Arquivo SAS >; Esse comando é normalmente utilizado para reabrir um SAS-DATA-SET ou para concatenar linhas. Ex. 1. Reabrir o Sas-Data-Set de nome SDS DATA SDS; SET SDS; ... Ex. 2. Criar um Sas-Data-Set de nome SDT contendo as linhas do Sas-Data-Set SDS1 seguida das linhas do SDS2 DATA SDT; SET SDS1 SDS2; ... Obs.: Neste caso SDS1 e SDS2 devem conter as mesmas variáveis. MERGE: Carrega variáveis (colunas) no SAS-DAT-SET. Sintaxe: MERGE < lista de Arquivo SAS >; BY < variável >; Ex. Sejam os arquivos SAS: SDS1 Nome Antônio Maria Pedro Rosa

Sexo M F M F

SDS2 Função Aux. Escritório Secretária Servente Cozinheira

Com os comandos: PROC SORT DATA=SDS1; BY NOME; RUN; PROC SORT DATA=SDS2; BY NOME; RUN; DATA SDS; MERGE SDS1 SDS2; BY NOME;

Nome Antônio Maria Rosa

Salário 600.00 800.00 300.00

20

Resulta o arquivo SAS : SDS Nome Antônio Maria Pedro Rosa

Sexo M F M F

Função Aux. Escritório Secretária Servente Cozinheira

Salário 600.00 800.00 . 300.00

Obs.: Os arquivos SDS1 e SDS2 devem estar classificados pela variável NOME.

INFILE “< arquivo ASC >” < opções >; Uma opção é a LRECL, que define o Tamanho Lógico do Registro – “LOGICAL RECORD LENGTH”. Se o tamanho não for definido, o SAS usa o tamanho 256.

21

LISTA DE EXERCÍCIOS - AULA 2 1. Fazer um programa no SAS para: a) Colocar comentários informando: o número da lista, do exercício e a data. b) Definir as opções para o OUTPUT: tamanho da linha=78, tamanho da página=64. e) Criar um arquivo SAS (SAS-DATA-SET), a partir do Arquivo L2E1.DAT, que é um arquivo ASCII e contém as variáveis: A B C D E Y1-Y32. Nota: registro muito grande - usar a opção “lrecl=nº” no INFILE. 2. Fazer um programa no SAS para: a) Colocar comentários informando: o número da lista, do exercício e a data. b) Definir as opções para o OUTPUT: tamanho da linha=78, tamanho da página=64. e) Criar um arquivo SAS, a partir do arquivo L2E2.XLS, que é um arquivo Excel com 9 variáveis e muitas informações perdidas. 3. Fazer um programa no SAS para: a) Colocar comentários informando: o número da lista, do exercício e a data. b) Definir as opções para o OUTPUT: tamanho da linha=78, tamanho da página=64. c) ler o arquivo L2E3A.DAT; que é um arquivo ASCII, e contem as variáveis: ANO (92-94), MES, Y. d) ler o arquivo L2E3B.DAT; que é um arquivo ASCII, e contem as variáveis: ANO (9596), MES, Y. e) criar um arquivo contendo os registros de L2E3A mais os de L2E3B. L2E3A L2E3B f) calcular a média, variância e desvio padrão de Y para os 5 anos (92 a 96). 4. Fazer um programa no SAS para: a) Colocar comentários informando: o número da lista, do exercício e a data. b) Definir as opções para o OUTPUT: tamanho da linha=78, tamanho da página=64. c) ler o arquivo L2E4A.DAT; que é um arquivo ASCII, e contém as variáveis: Tratam., Repet., Y1, Y2 e Y3. d) ler o arquivo L2E4B.DAT; que é um arquivo ASCII, e contém as variáveis: Tratam., Repet., Y4, Y5. e) criar um arquivo contendo as colunas de L2E4A.DAT mais as de L2E4B.DAT L2E4A L2E4B f) calcular o coeficiente de correlação de Pearson entre as variáveis Y1 a Y5. 5. Fazer um programa no SAS para: a) Colocar comentários informando: o número da lista, do exercício e a data. b) Definir as opções para o OUTPUT: tamanho da linha=78, tamanho da página=64. c) ler o arquivo L2E5.DAT; que é um arquivo ASCII, e contém as variáveis NOME SEXO IDADE PESO e ALTURA. d) Criar 2 arquivos SAS, um para os registros em que SEXO=M e outro com os registros em que SEXO=F.

22

V.3. INFORMAÇÕES NECESSÁRIAS PARA A AULA 3 Procedimentos PROC FREQ Finalidade: Criar tabelas de freqüência. A Tabela conterá a freqüência, as porcentagens em relação ao total geral, as porcentagens em relação ao total das linhas e as porcentagens em relação ao total de colunas. Sintaxe: PROC FREQ < opções1 >; TABLE < variável A > [* < variável B >] / < opções2 >; Algumas das são: DATA= <SDS> - especifica O Sas-Data-Set a ser usado PAGE – especifica para imprimir, no OUTPUT, uma tabela por página. Algumas das são: NOFREQ – Não apresenta as freqüências NOPERCENT – Não apresenta as porcentagens em relação ao total geral NOROW – Não apresenta as porcentagens em relação ao total de linhas NOCOL – Não apresenta as porcentagens em relação ao total de linhas EXPECTED – Apresenta os valores esperados CHISQ – Realiza o teste Qui-Quadrado para a independência das variáveis EXACT - Realiza o teste Exato de Fischer para a independência das variáveis

PROC TTEST Finalidade: Teste t de Student para comparação das médias de 2 grupos. Sintaxe: PROC TTEST < opções >; CLASS < variável >;

Variável que classifica os grupos

Uma das : DATA= <SDS> - especifica O Sas-Data-Set a ser usado

23

LISTA DE EXERCÍCIOS - AULA 3 1. Fazer um programa no SAS para: a) Ler os dados L3E1.DAT (Diâmetro do folículo 1 a 3; Qualidade do folículo 1 a 3 e qualidade do oócito 1 a 2). Observe que os dados não estão separados por espaços. b) Fazer o Teste Qui-quadrado. 2. Fazer um programa no SAS para fazer o teste Qui-quadrado e o teste Exato de Fisher para a tabela de contingência: RESULTADO DO ACASALAMENTO fecundo infecundo 52 128 77 83

RAÇA Charolesa Zebu

3. Digitar os Dados a seguir no Excel e criar um arquivo ASC (L3E3.TXT) Y1 Y2

72,5 72

69 75 69,1 75

70,8 71,2 73 70 67,1 71 72 73 70,9 71 72,9 69,9 37 71,1 71,9 73

68 69 69,9 70,1 71 74 68,2 68,9 70 70 71,1 74

4. Fazer um programa no SAS para: a) Colocar comentários informando: o número da lista e do exercício, a data e o nome do programador. b) Definir as opções para o OUTPUT: Tamanho da linha=78, tamanho da página=64. c) Criar um arquivo SAS, a partir do arquivo L3E3.TXT. d) Imprimir o arquivo. e) Obter o intervalo de confiança para a média de Y1 (Y1 ). f) Testar a Hipótese: H0: Y1 = 70 . 5. Fazer um programa no SAS para: a) Colocar comentários informando: o número da lista e do exercício e a data. b) Definir as opções para o OUTPUT: Tamanho da linha=78, tamanho da página=64. c) Criar um arquivo SAS a partir dos dados: T1 T2

760 758

755 748

758 757

761 753

755 755

d) Imprimir o arquivo. e) Aplicar o test t para testar a hipótese: H0: T1 = T2 . e) Substituir as 2 primeiras repetições de T2 por 708 e 798, respectivamente, e ver o que acontece.

72 72

24

V.4. INFORMAÇÕES NECESSÁRIAS PARA A AULA 4 Procedimentos PROC UNIVARIATE Finalidade: Apresenta as estatísticas descritivas necessárias para um diagnóstico da distribuição dos dados. Sintaxe: PROC UNIVARIATE < opções >; Algumas opções: NORMAL – Faz o teste de normalidade. Teste de Shapiro-Wilk (W) se o número de observações for <2000, caso contrário o teste de Kolmogorov (D). PLOT – Gera um gráfico de ramos e folhas que permite interpretar certas características da distribuição e um gráfico contendo os valores esperados, caso a distribuição for normal (+), e os valores observados (*). PROC INSIGHT Finalidade: Entra no modo interativo do SAS. 1. Sintaxe via programa: PROC INSIGHT DATA=<SDS>; DIST Var ; RUN; 2. Sintaxe via menu: Seleciona as opções: Globals / Analyze / Interactive data analysis a) Seleção do arquivo: Selecione o caminho. Campo LIBRARY. Selecione o arquivo. Campo DATA SET. Selecione Open b) Selecione na planilha a variável que vai estudar Um clique no cabeçalho da Coluna c)

Selecione as opções: Analyze / Distribution(Y)

Até este ponto, em qualquer uma das opções chegarão no mesmo local Aqui pode-se pedir, entre muitas coisas: Teste de Normalidade: Curves / Test for Distribution Gráfico da normal: Curves / Parametric density Tabela de Freqüências: Tables / Frequence table Intervalo de confiança para a média: Tables / C. L. for means

25

LISTA DE EXERCÍCIOS - AULA 4 1. Fazer um programa no SAS para: a) Ler os dados L3E3.TXT (variáveis Y1 Y2). b) Fazer teste de normalidade dessas variáveis (Usando o PROC UNIVARIATE e usando o SAS INSIGHT). 2. Fazer um programa no SAS para: a) Gerar 1000 dados com distribuição Normal com média 0 e Variância 1. c) Fazer uma análise exploratória da distribuição (Usando o PROC UNIVARIATE e usando o SAS INSIGHT).

26

V.5. INFORMAÇÕES NECESSÁRIAS PARA A AULA 5 Os gráficos no SAS podem ser feitos no OUTPUT ou numa janela apropriada GRAPHICS. Na janela GRAPHICS os gráficos são apresentados individualmente e com maior precisão que no OUTPUT.

Gráficos de dispersão Os procedimentos utilizados para gráficos de dispersão são: PROC PLOT – para gráficos no OUTPUT e PROC GPLOT – para gráficos na janela GRAPHICS. PROC PLOT Finalidade: Editar gráficos de dispersão no OUTPUT. Sintaxe: PROC PLOT < opções1 >; PLOT < Vertical >* < Horizontal > [=”caracter”] / < opções2 >; Algumas das opções1 são: caracter a ser usado NOLEGEND – Não apresenta as legendas do gráfico na representação dos HPERCENT – Especifica a porcentagem da página, pontos na horizontal, deve ser usada pelo gráfico VPERCENT – Especifica a porcentagem da página, na vertical, deve ser usada pelo gráfico Algumas das opções2 são: VPOS=Nº – especifica o número de linhas que serão usadas HPOS=Nº – especifica o número de colunas que serão usadas HAXIS =

TO BY - especifica a escala do eixo horizontal (X), sendo “li” o limite inferior, “ls” o limite superior e “inc” o incremento VAXIS =

TO BY - especifica a escala do eixo vertical (Y), sendo “li” o limite inferior, “ls” o limite superior e “inc” o incremento HREF = n1 n2 ... – coloca linhas de referências verticais, nas posições n1, n2, etc, do eixo horizontal VREF = n1 n2 ... – coloca linhas de referências horizontais, nas posições n1, n2, etc, do eixo vertical OVERLAY – usado para sobrepor 2 ou mais gráficos VZERO – Faz com que a origem do eixo vertical seja no ponto zero HZERO – Faz com que a origem do eixo horizontal seja no ponto zero

27

PROC GPLOT Finalidade: Editar gráficos de dispersão na janela GRAPHICS. Sintaxe: PROC GPLOT < opções1 >; PLOT < Vertical >* < Horizontal > / < opções2 >; SYMBOL[n] ; 1 – primeira sequência de pontos 2 – Segunda etc...... Uma das opções1 é: UNIFORM – Especifica que todos os gráficos feitos por este procedimento devem ter a mesma escala Algumas das opções2 são: Além das opções: VPOS; HPOS; HAXIS; VAXIS; HREF; VREF; OVERLAY; HZERO e VZERO, descritas no PROC PLOT, pode-se usar: GRID – Coloca linhas de grade (horizontais e verticais) FRAME – Coloca um contorno na área do gráfico LEGEND – Inclui a legenda no gráfico Algumas das opções3 são: COLOR = – Especifica a cor dos símbolos INPERPOL = <none | | R | spline | ... > - Especifica a forma de ligação dos pontos, que pode ser, respectivamente, nenhum, linha segmentada, equação de regressão, curva suavizada, etc. VALUE = - Especifica a forma do símbolo, que pode ser, respectivamente, +, , etc. WIDTH = < 1, 2, 3, ... > - Especifica a espessura da linha HEIGHT = < 1, 2, 3, ... > - Especifica o tamanho dos símbolos

Outros gráficos Os procedimentos utilizados são: PROC CHART ou GCHART– para gráficos no output ou na janela GRAPHICS, respectivamente. VBAR ou HBAR ou PIE, etc – define o tipo do gráfico, que pode ser: Barras verticais, Barras Horizontais, tipo Pizza, etc. Histograma (PROC GCHART) Exemplo 1 – Fazer um histograma, com 10 classes de freqüências, por sexo. DATA EX1; INPUT SEX$ PESO @@; CARDS; M 18.3 F 22 M 19 F 26 M 20 F 15.5 M 25.5 F 12 M 18 F 21.7 M 13.5 M 16 M 24

28

M 22 F 23.2 M 20 F 21 M 31 F 28 F 20 M 13.5 F 23 M 19 M 18 M 16 F 17 ; RUN; GOPTIONS DEVICE=WIN; PATTERN V=SOLID COLOR=BLUE; PROC GCHART; VBAR PESO/GROUP=SEX LEVELS=10 ; RUN;

Gráficos de 3 dimensões (PROC G3D) O procedimento utilizado é: PROC G3D – para gráficos na janela GRAPHICS Exemplo 2. Fazer o gráfico da equação Z=X2 + Y2 , -5<X<5 E –5 <5. DATA EX2; DO X=-5 TO 5 BY 0.2; DO Y=-5 TO 5 BY 0.2; Z=X**2+Y**2; OUTPUT; END; END; PROC G3D; PLOT X*Y=Z/GRID CAXIS=RED CBOTTOM=BLUE; RUN;

29

LISTA DE EXERCÍCIOS - AULA 5 1. Fazer um programa SAS para ler os dados da Tabela abaixo e representar graficamente o perfil dos 3 tratamentos ao longo do tempo. Usar as opções que julgar serem interessantes. Semanas Tratamento 1 2 3 4 5 6 7 8 9 10 TR1 10 11 11 11.5 12 11 12.3 12.5 12.3 12.5 TR2 12 12.1 12.5 12.3 12.4 13 13.1 12.9 13.1 13.4 TR3 13.1 13 12.9 12.5 11.7 12.5 11.5 11 10.5 10 2. Fazer um programa SAS para gerar 100 dados com distribuição normal e fazer um histograma (gráfico de barras verticais). Usar as opções que você julgar serem interessantes. 3. Fazer um programa SAS para fazer o gráfico da equação: Z = 1 / (X2 + Y2 + 0.5), 10<X<10 e –10 <10. Usar as opções que você julgar serem interessantes.

-

30

V.6. INFORMAÇÕES NECESSÁRIAS PARA A AULA 6 Regressão Linear O procedimento utilizado para regressão linear e o PROC REG:

PROC REG Finalidade:

Ajuste da equação de regressão linear.

Sintaxe: PROC REG < opções1 >; MODEL < Var. Dep > = < Lista Var. Indep. > / < opções2 >; RUN; A < Lista Var. Indep. > é tal que: Se conter apenas 1 variável – Reg Linear simples Se conter 2 ou mais variáveis – Reg. linear Múltipla Algumas das são: P – apresenta os valores estimados (preditos) CLI – apresenta os intervalos de confiança dos valores preditos CLM – apresenta os intervalos de confiança das médias SELECTION= BACKWARD (B) | FORWARD (F) | STEPWISE | SLSTAY= n1 | SLENTRY=n2 As opções BACKWARD (B) | FORWARD (F) | STEPWISE são usadas para definir o método para seleção do modelo de regressão As opções SLSTAY= n1 | SLENTRY=n2 são usadas para definir os Níveis de significância para ficar ou entrar, respectivamente, uma variável no modelo

31

LISTA DE EXERCÍCIOS - AULA 6 1. REGRESSÃO LINEAR SIMPLES. Fazer um programa no SAS para: a) Ler os dados no SAS usando a opção CARDS: Valores Valores de Y de X R1 R2 R3 1 8.2 6.6 9.8 2 19.7 15.7 16.0 3 28.6 25.0 31.9 4 30.8 37.8 40.2 5 40.3 42.9 32.6 b) Representar (no OUTPUT) os pontos graficamente. c) Fazer a análise de regressão linear y=f(x)=a+bx. d) Fazer a análise de regressão linear com opções para obter: valores estimados (preditos) de Y ( y i ), Intervalos de confiança para y i , e representar graficamente (na janela GRAPHICS) pontos observados e os estimados com respectivos intervalos. 2. REGRESSÃO LINEAR MÚLTIPLA. Fazer um programa no SAS para: a) Ler o arquivo L6E2.DAT, que é um arquivo ASC, e contém as variáveis: X1 = N, X2 = P, X3 = K, X4 = CA, X5 = MG, X6 = S e Y = % terra retida na peneira 18. b) Fazer a análise de regressão linear múltipla: Y=f(X1,X2, X3,X4, X5,X6). Obter os valores preditos com respectivos intervalos de confiança. c) Testar se o intercepto difere de 50 e se o coeficiente de N difere de 1. d) Representar graficamente os valores observados vs estimados. e) Fazer a análise de regressão linear múltipla e selecionar o modelo usando os métodos: BACKWARD, FORWARD, STEPWISE e R-MÁXIMO. Definir o Nível de Significância para a seleção dos modelos. 3. REGRESSÃO POLINOMIAL. Fazer um programa no SAS para: a) Ler os dados no SAS usando a opção CARDS: quantidade Tempo para Gelar (Y) de pó (X) Rep.1 Rep.2 2.50 7.39 7.30 2.55 7.00 7.03 2.60 6.90 6.95 2.65 6.85 6.80 2.70 6.70 6.30 2.75 6.33 6.20 2.80 5.97 5.90 2.85 5.90 5.82 2.90 5.80 5.80 2.95 6.15 6.00 3.00 6.30 6.15 b) Representar graficamente Y=f(X). c) Fazer a análise de regressão polinomial, até 3º grau; de Y em função de X.

32

V.7. INFORMAÇÕES NECESSÁRIAS PARA AS AULAS 7 E 8 As análises de variância no SAS podem ser realizadas pelos procedimentos: PROC ANOVA e PROC GLM. O PROC ANOVA (ANalyze Of VAriance) é específico para dados balanceados e o PROC GLM (General Linear Model) pode ser utilizado em qualquer caso (dados balanceados ou não).

PROC ANOVA A sintaxe do PROC ANOVA é:

PROC ANOVA ; CLASS ; MODEL <lista de var, dep,>=efeitos / ; MEANS <efeitos> / ; TEST H= <lista de efeitos> E= <efeito>; RUN; Os comandos em negrito são obrigatórios, Nos Comandos:

CLASS ; corresponde às variáveis de classificação que serão utilizadas na análise

MODEL <lista de var, dep,>=<efeitos> / ; <lista de var, dep,> - corresponde às variáveis a serem analisadas <efeitos> - corresponde aos efeitos envolvidos no modelo, eqüivale às fontes

de variação do esquema da análise da variância, exceto o resíduo - Efeitos cruzados (interações) são representados incluindo-se um * entre as variáveis da interação, para a interação de A com B, usa-se: A*B - Efeitos aninhados são representados colocando-se a variável aninhada entre parênteses, para representar A aninhado em B usa-se A(B) Obs.: Efeitos incluídos no modelo, que não constam da lista de variáveis de classificação são consideradas como Covariáveis. - Algumas da opções são: - SS1 | SS2 | SS3 | SS4 - Para escolher o tipo das Somas de Quadrados a ser usado

33

MEANS <efeitos> / ;

Apresenta as médias para os níveis dos efeitos especificados, sendo: <efeitos> - lista de efeitos para os quais deseja as médias - Algumas da opções são: - Teste para Comparações de médias (TUKEY | DUNCAN | T | BON | DUNNETT | SCHEFFE | etc.) - Teste para Homogeneidade de Variância (HOVTEST = BARTLETT | BF (Brown Forsythe) LEVENE | OBRIEN )

TEST H= <lista de efeitos> E= <efeito>; H= <lista de efeitos> - especifica os efeitos que devem ser testados (teste F), usando como termo do erro o efeito especificado em E= <efeito>

PROC GLM A sintaxe do PROC GLM é:

PROC GLM ; CLASS ; MODEL <lista de var, dep,>=<efeitos> / ; MEANS <efeitos> / ; LSMEANS efeitos / ; CONTRAST ' label'<efeito> / ; RANDOM <efeitos> / ; TEST H= <lista de efeitos> E= <efeito>; RUN; Os comandos em Negrito são Obrigatórios, Apenas os comandos não descritos no PROC ANOVA (em azul) serão descritos aqui, uma vez que para os outros a descrição seria igual,

LSMEANS efeitos / ;

Apresenta as médias ajustadas para os níveis dos efeitos especificados, sendo: <efeitos> - lista de efeitos para os quais deseja as médias - Algumas da opções são: - PDIFF | TDIFF – requer a apresentação do valor da probabilidade “p-value” ou do valor da estatística T, respectivamente, para as comparações de médias duas a duas,

34

- ADJUST = (TUKEY | BON | SCHEFFE | etc) – requer o ajuste do “pvalue” para o teste especificado

CONTRAST “label” <efeito> / ; Faz o teste F para o contraste especificado: label – qualquer comentário (até 20 caracteres) que caracterize o contraste <efeito> - efeito ao qual corresponde o contraste coeficientes do contraste – coeficientes do contraste, um para cada nível do efeito, separados por pelo menos um espaço - Algumas da opções são: - E – requer a impressão dos coeficientes do contraste - E=efeito – especifica efeito que será usado como termo do erro, para o teste F RANDOM <efeitos> / ;

Especifica os efeitos aleatórios do modelo, Sempre que este comando for incluído, o SAS apresenta as Esperanças dos Quadrados Médios do Modelo. <efeitos> - lista de efeitos aleatórios - Algumas da opções são: - Q – requer que todas as formas quadráticas dos efeitos fixos sejam apresentadas no OUTPUT - TEST – requer que os testes F para a análise da variância sejam feitos não com o resíduo geral, e sim com os resíduos apropriados

35

LISTA DE EXERCÍCIOS - AULA 7 1.A. Fazer um programa no SAS para: a) Criar um arquivo SAS, a partir dos dados a seguir (Arquivo ASC- L7E1.TXT): REPETIÇÕES 1 2 3 4 5 1 2370,0 1687,0 2592,0 2283,0 2910,0 2 1282,0 1527,0 871,0 1025,0 825,0 3 562,0 321,0 636,0 317,0 485,0 4 173,0 127,0 132,0 150,0 129,0 5 193,0 71,0 82,0 62,0 96,0 Fonte: BANZATTO, D.A., KRONKA, S.N. Experimentação Agrícola. Jaboticabal: Funep, 1992. TRATAM.

6 3020,0 920,0 842,0 227,0 44,0

b) Verificar as condições de Normalidade do Erro c) Verificar se há homogeneidade da variância d) Com a finalidade de estudar se a heterocedasticidade da variância é regular, verificar se a regressão linear Log(Var)=a + b Log(Med) é significativa. 1.B. Com os dados do exercício anterior, fazer um programa no SAS para: a) Transformar os dados em Z=Ln(Y) - logarítmo neperiano b) Verificar as condições de Normalidade do Erro dos dados transformados c) Verificar se há homogeneidade da variância d) Fazer a Análise da Variância (DIC) e) Comparar as médias pelo teste de Duncan 2.A. Fazer um programa no SAS para: a) Criar um arquivo SAS, a partir dos dados a seguir (Arquivo ASC- L7E2.TXT): REPETIÇÕES 1 2 3 4 1 142,36 144,78 145,19 138,88 2 139,28 137,77 144,44 130,61 3 140,73 134,06 136,07 144,11 4 150,88 135,83 136,97 136,36 5 153,49 165,02 151,75 150,22 Fonte: BANZATTO, D.A., KRONKA, S.N. Experimentação Agrícola. Jaboticabal:Funep, 1992. TRATAM,

b) Fazer a Análise da Variância (DBC) c) Comparar as médias pelo teste de Tukey d) Fazer um desdobramento dos graus de liberdade, por contraste, considerando que T1, T2 e T3 possuem uma característica não observada em T4 e T5 2.B. Fazer um programa no SAS para a Análise de Variância dos dados do exercício 2.A., supondo que a observação (Trat. 5 e Bloco 2) tenha sido perdida.

36

3)

Fazer um programa no SAS para: a) Criar um arquivo SAS, a partir dos dados a seguir (Arquivo ASC – L7E3TXT)

TRATAMENTOS TRAT1 TRAT2 TRAT3 TRAT4 TRAT5 Y X Y X Y X Y X Y X 1 74 9 58 7 118 9 41 6 95 8 2 51 9 67 8 48 9 38 9 41 8 3 95 8 40 5 49 9 77 8 39 9 4 62 9 58 8 64 9 92 9 114 9 5 60 9 29 6 67 8 57 7 35 6 6 47 9 64 8 51 8 77 7 49 8 7 14 6 55 9 15 8 59 8 39 9 8 19 8 47 8 29 9 32 8 100 9 Fonte: PIMENTELGOMES, F. Estatística Experimental. Piracicaba:ESALQ/USP, 12. Ed., 1987. BLOCOS

b) Fazer análise de Covariância 4.A. Fazer um programa SAS para: a) Criar um arquivo SAS com os dados a seguir (Arquivo ASC – L7E4.TXT): REPETIÇÕES 1 2 3 4 1-R1E1 26,2 26,0 25,0 25,4 2-R1E2 24,8 24,6 26,7 25,2 3-R2E1 25,7 26,3 25,1 26,4 4-R2E2 19,6 21,1 19,0 18,6 5-R3E1 22,8 19,4 18,8 19,2 6-R3E2 19,8 21,4 22,8 21,3 Fonte: BANZATTO, D.A., KRONKA, S.N. Experimentação Agrícola. Jaboticabal:Funep, 1992, TRATAM.

b) Fazer a Análise da variância considerando um delineamento DBC e um esquema fatorial 3x2. c) Fazer o desdobramento dos graus de liberdade da interação. 4.B. Fazer um programa no SAS para a Análise de Variância dos dados do exercício 4.A., supondo que B esta aninhado em A. 5) Fazer um programa no SAS para: a) Criar um arquivo SAS, a partir dos dados (Arquivo ASC – L7E5.TXT):: TRATAM, A1B1 A1B2 A1B3 A2B1 A2B2 A2B3

BLOCOS 1 42,9 53,8 49,5 53,3 57,6 59,8

2 41,6 58,5 53,8 69,6 69,6 65,8

3 28,9 43,9 40,7 45,4 42,4 41,4

b) Fazer a análise de variância considerando um Delineamento em Parcelas Subdivididas

4 30,8 46,3 39,4 35,1 51,9 45,4

37

V.8. INFORMAÇÕES NECESSÁRIAS PARA A AULA 8 Considerações estatísticas Esperanças dos quadrados médios A esperança do quadrado médio é utilizada para verificar as hipóteses que estão sendo testadas e os denominadores apropriados para os testes F da análise da variância e para obter as estimativas dos componentes da variância. Num Delineamento Interiamente Casualizado, com nt tratamentos e nr repetições, as esperanças dos quadrados médios são dadas por: F. V. Trat, Resíduo Total

G.L.

S.Q.

nt – 1 nt(nr – 1) nt nr – 1

SQ(Trat,) SQ(Res,) SQ(Total)

Onde φt =

E(QM) Ef. Aleatório Ef. Fixo 2 2 σ + nt σt σ2+ nt φt 2 σ2 σ

F QM(Trat)/QM(Res)

1 nt (µ − µ i ) 2 nt − 1 i =1

hipóteses que estão sendo testadas Observe que o Teste F esta testando a hipótese: a) Efeito aleatório: H0: σt2=0 vs H1: σt2 ≠ 0, b) Efeito fixoo: H0: µ1= µ2 = µ3 = ... = µnt = µ vs H1: µi ≠ µj, i ≠ j.

estimativas dos componentes da variância Uma forma de obter as estimativas dos componentes da variância é igualar os Quadrados Médios observados às suas respectivas esperanças, assim:

σˆ 2 = QM (Re s) σˆ 2 + nt σˆ t2 = QM (Trat ) logo:

σˆ 2 = QM (Re s ) σˆ t2 =

1 [QM (Trat ) − QM (Re s )] nt

38

Regras práticas para obtenção das esperanças dos quadrados médios E(QM) dados balanceados 1. Determine as Fontes de Variação que comporão a Análise de Variância. Os efeitos principais e todas as possíveis interações. Se a letra ocorre em ambos os lados, a interação não existe. Se a letra se repete do lado direito, coloque a letra apenas uma vez. Ex: Fatores Interação ----------------------------------------AeB AB A e C:B AC:B A:B e C:B AC:BB = AC:B A:B e B:DE AB:BDE (não existe) ----------------------------------------2. A análise de Variância terá uma linha para cada Causa de Variação, inclusive para o resíduo. 3. O nº de graus de liberdade de uma linha genérica é obtida pelo produto dos (ni-1) associados aos fatores à esquerda dos dois pontos (:) com os nj aos à direita. Ex: Para AB:XY ------> (na-1)(nb-1)nxny graus de liberdade. 4. Denote os Componentes da Variância por s2 com respectivos índices Ex: Para A:B use s2a:b. 5. Os coeficientes dos Componentes da Variância é o produto de todos os ni cujos índices não aparecem em s2. 6. Cada E(QM) conterá aqueles s2 (com respectivos coeficientes) que possuem como índices, no mínimo, todas as letras presentes naquela linha. 7. Para os modelos mistos, sem restrições nos parâmetros, tratá-los como aleatórios, substituindo os s2 dos efeitos fixos pelas funções quadráticas correspondentes. 8. Para modelos mistos, com restrições nos parâmetros, eliminar da E(QM) aqueles Componentes da Variância que têm entre suas letras, qualquer letra de efeito fixo, exceto as da linha.

39

Procedimentos para estimar componentes da variância PROC VARCOMP Finalidade: Estimar componentes da variância Sintaxe: PROC VARCOMP ;

CLASS ; MODEL <lista de var, dep,>=<efeitos> / ; RUN;

CLASS e MODEL – como definidas no PROC ANOVA e PROC GLM - Algumas da opções são: - SAS-data-set – especifica o SAS-DATA-SET - EPSILON = n – para especificar o valor de ε (epsilon) utilizado no teste de convergência para os métodos ML e REML, Default: 1*E-8 - MAXITER = n - para especificar o número máximo de iterações para os métodos ML e REML, Default: 50 - METHOD = TYPE1|MIVQUE0|ML|REML - para especificar o método de estimação dos componentes da variância, Default: MIVQUE0 - Uma opção é: - FIXED = n - para especificar que os n primeiros efeitos são fixos e, automaticamente, os demais são considerados aleatórios, Default: Todos os efeitos aleatórios.

40

LISTA DE EXERCÍCIOS - AULA 8 1. Fazer um programa em SAS para : a) Ler o arquivo L8E1,TXT, que é um arquivo ASC, O arquivo contém as variáveis A, B, Rep, Y. FATORES REPETIÇÕES 1 2 3 A B 1 1 86,8 90,20 76,80 1 2 94,00 88,60 83,00 1 3 44,00 83,80 56,50 2 1 72,20 52,40 55,00 2 2 72,40 120,20 90,80 2 3 104,60 92,00 101,00 3 1 68,40 91,00 55,20 3 2 78,80 49,00 56,00 34,00 57,20 32,40 3 3 4 1 54,00 24,40 35,00 33,00 33,60 34,80 4 2 4 3 33,20 32,00 33,60 b) Fazer análise de variância considerando que B está aninhado em A, sendo A fixo e B(A) aleatório. Comparar as médias para o efeito fixo. Estimar os componentes da variância para o efeito aleatório. 2. Fazer a análise da variância de um experimento contendo os fatores: A, B e C(A), todos aleatórios. a) Obter as EQM e montar a análise. b) Pedir ao SAS as EQM, e então montar a análise. c) Pedir ao SAS as EQM e a montagem da análise. Dados: A B C ----- Y ----1 1 1 5 6 9 6 7 1 1 2 6 8 9 8 7 1 1 3 4 5 8 6 6 1 2 1 8 6 9 7 5 1 2 2 7 8 9 6 8 1 2 3 10 8 7 9 6 2 1 1 5 7 6 8 6 2 1 2 9 7 6 9 7 2 1 3 6 8 7 9 5 2 2 1 6 5 6 7 8 2 2 2 7 8 7 9 6 2 2 3 6 5 7 9 8

41

V.9. INFORMAÇÕES NECESSÁRIAS PARA A AULA 9 Procedimentos para estimar componentes da variância PROC NPAR1WAY Finalidade: Análise não paramétrica para experimentos com 1 fator, Sintaxe: PROC NPAR1WAY ;

CLASS ; VAR ; RUN;

CLASS e MODEL – como definidas no PROC ANOVA e PROC GLM - Algumas das opções são:

- SAS-data-set – especifica o SAS-DATA-SET - ANOVA – requer uma análise de variância padrão, além da análise não paramétrica - VW – requer que os escores de Van der Waerden sejam analisados - WILCOXON – requer uma análise dos números de ordem dos dados, ou escores de Wilcoxon, para 2 níveis, eqüivale ao teste da soma das ordens de Wilcoxon, para qualquer número de níveis, é o teste de Kruskal-Wallis

42

LISTA DE EXERCÍCIOS - AULA 9 1. Fazer um programa no SAS para: a) Criar um arquivo SAS, a partir dos dados a seguir: TRATAM, 1 2 3

REPETIÇÕES 1 2212,8 2195,2 1770,3

2 2025,3 2031,5 1800,0

3 1989,0 1876,5 1852,8

4 2232,8 1750,3 1769,0

5 2027,8 1060,3

Fonte: CAMPOS, H. Estatística Experimental não paramétrica. Piracicaba:ESALQ/USP. 3. Ed., 1979.

b) Verificar as condições de Homogeneidade da Variância (Teste de Bartlett). c) Fazer o teste não paramétrico de Wilcoxon para comparação dos tratamentos.

43

V.10. INFORMAÇÕES NECESSÁRIAS PARA A AULA 10 Procedimentos para regressão não linear PROC NLIN (Non LINear regression) Finalidade: Sintaxe:

-

-

Estimativa de parâmetros de equações de regressão não-lineares.

PROC NLIN ; PARAMETERS <param>= ... ; MODEL = expressão; BOUNDS expressões ... ; DER<.param> = expression;

- Algumas da opções são:

METHOD= GAUSS | MARQUARDT | NEWTON | GRADIENT | DUD – Especifica o método a ser usado. . Se o METHOD não for especificado e as derivadas “DER” forem, o SAS usa o método : GAUSS . Se o METHOD=DUD for selecionado, as derivadas “DER” não precisam ser especificadas. OUTTEST=SAS-data-set – cria um SAS-DATA-SET com os parâmetros prodizidos em cada iteração MAXITER = n - para especificar o número máximo de iterações, Default: 50

PARAMETERS param,= - identifica os parâmetros a serem estimados e atribui valores iniciais aos mesmos MODEL = expressão – identifica a relação funcional entre a

variável dependente e a independente

BOUNDS expressões – define limites para as estimativas dos parâmetros DER<.param> = expressão – define as derivadas da equação de

regressão, em relação a cada um dos parâmetros

44

LISTA DE EXERCÍCIOS - AULA 10 1. Fazer um programa no SAS para: a) Os dados a seguir referem-se à digestibilidade da matéria seca de capim coast-cross, em função do tempo de Incubação. Criar um arquivo SAS, a partir desses dados (Arquivo ASC – L19E1). Tempo de Incub. 0 3 6 12 24 48 72

1 14,85 17,39 22,32 29,21 40,73 47,34 50,41

2 14,80 18,11 21,04 30,20 40,96 51,34 54,62

Repetições 3 4 15,07 14,27 18,50 16,03 23,43 19,53 29,78 32,52 41,36 43,88 51,61 53,35 55,72 56,21

5 20,46 22,83 27,94 36,27 46,00 57,07 59,93

6 20,07 21,75 23,98 31,78 39,13 49,49 55,12

b) Representar graficamente Y=f(X), onde Y=Degradabilidade e X=Tempo de Incubação. c) Estimar os parâmetros do modelo de ORSKOV & McDONALD (1979): Deg(t)=a+b(1-e-ct). d) Representar graficamente os valores observados e o modelo estimado.

45

VI. Resolução das listas de exercícios VI.1. RESOLUÇÃO DA LISTA DE EXERCÍCIOS - AULA 1 EXERCÍCIO 1 /* LISTA 1 – EXERCÍCIO 1 */ /* DD/MM/AA - FULANO DE TAL */ OPTIONS LS=78 PS=64; TITLE “CURSO SAS – ANÁLISE DE DADOS”; FOOTNOTE “FCAV/UNESP”; DATA L1E1; INPUT TR RP Y@@; CARDS; 1 1 35 1 2 19 1 3 31 1 4 15 1 5 30 2 1 40 2 2 35 2 3 46 2 4 41 2 5 33 3 1 39 3 2 27 3 3 20 3 4 29 3 5 45 ; PROC PRINT; PROC CONTENTS; RUN;

EXERCÍCIO 2 /* LISTA 1 - EXERCÍCIO 2 */ OPTIONS PS=64 LS=78; DATA L1E2; INFILE "A:\A1E2.DAT"; INPUT VACA 1-4 PAIV 5-8 MAEV 9-12 DIAP 13-14 MESP 15-16 ANOP 17-18 ORP 19-20 SEXO $ 21 ANOD 22-23 MESD 24-25 DIAD 26-27 PESON 28-33 .2 PESOD 34-39 .2; *PROC PRINT; *RUN; ID=MDY(MESD,DIAD,ANOD)-MDY(MESP,DIAP,ANOP); GPTD=PESOD-PESON; GPDD=GPTD/ID; IF 4<=MESP<=9 THEN PERN=”SECA”; ELSE PERN=”AGUAS”; DROP DIAP MESP ANOP DIAD MESD ANOD PESON PESOD; PROC PRINT; RUN; PROC SORT; BY SEXO; PROC PRINT; RUN; PROC SORT; BY SEXO PERN; PROC PRINT; RUN;

46

EXERCÍCIO 3 /* LISTA 1 - EXERCÍCIO 3 – DD/MM/AA – FULANO DE TAL */ OPTIONS PS=64 LS=78; LIBNAME PASTA “C:\SAS”; DATA PASTA.L1E3; INFILE "A:\A1E3.DAT"; INPUT RC OR RP PN GP30; LABEL RC="RAÇA" OR="ORIGEM" RP="REPETIÇÕOES" PN="PESO AO NASCER" GP30="GANHO DE PESO AOS 30 DIAS"; PROC PRINT; RUN;

EXERCÍCIO 4 /* LISTA 1 - EXERCÍCIO 4 */ OPTIONS PS=64 LS=78; LIBNAME PASTA “C:\SAS”; PROC MEANS DATA=PASTA.L1E3 MEAN VAR STD; VAR PN GP30; RUN; PROC SORT; BY RC; RUN; PROC MEANS MEAN; VAR PN GP30; BY RC; PROC MEANS MEAN; VAR PN GP30; CLASS RC OR; RUN;

EXERCÍCIO 5 /* LISTA 1 - EXERCÍCIO 5 */ OPTIONS PS=64 LS=78; DATA L1E5; DO TRAT=1TO 3; DO REP=1 TO 5; INPUT Y@@; OUTPUT; END; END; CARDS; 35 19 31 15 30 40 35 46 41 33 39 27 20 29 45 ; PROC PRINT; RUN;

EXERCÍCIO 6 SIGA OS OS ESPECIFICADOS NO ITEM: COMO CRIAR UM ARQUIVO ASC, A PARTIR DO EXCEL, DESCRITO ANTERIORMENTE

47

VI.2. RESOLUÇÃO DA LISTA DE EXERCÍCIOS - AULA 2 EXERCÍCIO 1 /* LISTA 2 - EX 1 */ OPTIONS LS = 78 PS = 64; DATA L2E1; INFILE "A:\L2E1.DAT" LRECL = 300; INPUT A B C D E Y1-Y32; PROC PRINT; RUN;

EXERCÍCIO 2 TITLE "LISTA 2 EX 2"; OPTIONS LS = 78 PS = 64; DATA L2E2; INFILE "A:\L2E2.TXT"; INPUT AL$ GRAD GRUPO MES COL HGU1 HGU2 HGAR1 HGAR2; HGU1T = LOG(HGU1 + 1); IF GRAD = . THEN DELETE; *PROC SORT; *BY DESCENDING HGU1; *RUN; IF HGU1 >= 15 THEN DELETE; PROC UNIVARIATE PLOT NORMAL; VAR HGU1T; RUN;

EXERCÍCIO 3 /*LISTA 2 EXERCICIO 3 */ OPTION LS = 78 PS = 64; DATA L2E3A; INFILE "A:\L2E3A.DAT"; INPUT ANO MES Y; DATA L2E3B; INFILE "A:\L2E3B.DAT"; INPUT ANO MES Y; DATA L2E3; SET L2E3A L2E3B; RUN; PROC PRINT DATA = L2E3A; PROC PRINT DATA = L2E3B; PROC PRINT DATA = L2E3; RUN; PROC MEANS MEAN VAR STD; VAR Y; RUN;

48

EXERCÍCIO 4 /*LISTA 2 EXERCICIO 4 */ OPTION LS = 78 PS = 64; DATA L2E4A; INFILE "A:\L2E4A.DAT"; INPUT TR RP Y1 Y2 Y3; DATA L2E4B; INFILE "A:\L2E4B.DAT"; INPUT TR RP Y4 Y5; DATA L2E4; MERGE L2E4A L2E4B; RUN; PROC PRINT DATA = L2E4A; PROC PRINT DATA = L2E4B; PROC PRINT DATA = L2E4; RUN; PROC CORR; VAR Y1 Y2; WITH Y3-Y5; RUN;

EXERCÍCIO 5 * LISTA 2 EXERCICIO 5; OPTIONS LS=78 PS=64; DATA L2E5; INFILE "A:L2E5.DAT"; INPUT NOME$ SEXO$ ID PS ALT; PROC PRINT; RUN; DATA MASC FEM; SET L2E5; IF SEXO="M" THEN OUTPUT MASC; ELSE OUTPUT FEM; RUN; PROC PRINT DATA=FEM; PROC PRINT DATA=MASC; RUN;

49

VI.3. RESOLUÇÃO DA LISTA DE EXERCÍCIOS - AULA 3 EXERCÍCIO 1 /* LISTA 3 EXERCICIO 1 */ OPTIONS LS=78 PS = 64; DATA L3E1; INFILE "A:\L3E1.DAT"; INPUT DF 1 QF 2 QO 3; *PROC PRINT; RUN; PROC FREQ ; *TABLE DF*QF/CHISQ; *TABLE DF*QO/CHISQ; *TABLE QF*QO/CHISQ; TABLE DF*QF*QO/CHISQ; RUN;

EXERCÍCIO 2 /* LISTA 3 EXERCICIO 2 */ OPTIONS LS=78 PS=64; DATA L3E2; DO I=1 TO 52; RC = "CH"; FC="F"; OUTPUT; END; DO I=1 TO 128; RC = "CH"; FC="I"; OUTPUT; END; DO I=1 TO 77; RC = "ZB"; FC="F"; OUTPUT; END; DO I=1 TO 83; RC = "ZB"; FC="I"; OUTPUT; END; DROP I; PROC PRINT NOOBS; RUN; PROC FREQ; TABLE RC*FC/CHISQ; RUN;

EXERCÍCIO 3 SIGA OS OS ESPECIFICADOS NO ITEM: COMO CRIAR UM ARQUIVO ASC, A PARTIR DO EXCEL, DESCRITO ANTERIORMENTE

EXERCÍCIO 4

50 /* LISTA 3 EXERCICIO 4 */ OPTIONS LS=78 PS=64; DATA L3E4; INFILE "A:L3E3.TXT"; INPUT Y1 Y2; Z = Y1-70; IF Y2=37 THEN Y2= .; *PROC PRINT; RUN; PROC UNIVARIATE NORMAL PLOT; RUN; PROC MEANS ALPHA=0.10 MEAL=CLM; VAR Y1; RUN; PROC MEANS T PRT; VAR Z; RUN;

EXERCÍCIO 5 /* LISTA 3 EXERCICIO 5 */ OPTIONS LS=78 PS=64; DATA L3E5; INPUT TR$ Y; CARDS; T1 760 T1 755 T1 758 T1 761 T1 755 T2 758 T2 748 T2 757 T2 753 T2 755 ; PROC PRINT; RUN; PROC TTEST; CLASS TR; RUN;

51

VI.4. RESOLUÇÃO DA LISTA DE EXERCÍCIOS - AULA 4 EXERCÍCIO 1 SIGA OS OS ESPECIFICADOS NO ITEM: SAS - INSIGHT

EXERCÍCIO 2 DATA L4E2; DO I=1 TO 1000; Y = NORMAL (0); OUTPUT; END; PROC PRINT; RUN;

52

VI.5. RESOLUÇÃO DA LISTA DE EXERCÍCIOS - AULA 5 EXERCÍCIO 1 /* LISTA 5 EXERCICIO 1 */ OPTION LS=78 PS=64; DATA LGE1; INPUT S T1 T2 T3 @@; CARDS; 1 10 12 13.1 2 11 12.1 13 3 11 12.5 12.9 4 11.5 12.3 12.5 5 12 12.4 11.7 6 11 13 12.5 7 12.3 13.1 11.5 8 12.5 12.9 11 9 12.3 13.1 10.5 10 12.5 13.4 10 ; PROC PRINT; RUN; PROC GPLOT; PLOT T1*S T2*S T3*S/OVERLAY LEGEND ; SYMBOL1 COLOR=RED INTERPOL=J VALUE=STAR HEIGHT=1 WIDTH=1; SYMBOL2 COLOR=BLUE INTERPOL=J VALUE=DOT HEIGHT=1 WIDTH=1; SYMBOL3 COLOR=GREEN INTERPOL=J VALUE=PLUS HEIGHT=1 WIDTH=1; RUN;

EXERCÍCIO 2 /* LISTA G EXERCICIO 2 */ OPTION LS=78 PS=64; DATA LGE2; DO I=1 TO 100 BY 1; Y=NORMAL(0); OUTPUT; END; PROC PRINT; RUN;

53

VI.6. RESOLUÇÃO DA LISTA DE EXERCÍCIOS - AULA 6 EXERCÍCIO 1 /* LISTA 6- EXERCÍCIO 1 */ OPTIONS LS=78 PS=64 NODATE; DATA L4E1; DO X=1 TO 5; DO RP=1 TO 3; INPUT Y@@; OUTPUT; END; END; CARDS; 8.2 6.6 9.8 19.7 15.7 16.0 28.6 25.0 31.9 30.8 37.8 40.2 40.3 42.9 32.6 ; PROC PRINT; RUN; PROC PLOT NOLEGEND; PLOT Y*X="*"/VPOS=20 HPOS=60; TITLE "GRÁFICO DE DISPERSÃO Y=F(X)"; RUN; PROC REG; MODEL Y=X/P CLI; /* CLI = INTERVALO DE CONFIANCA DO VALOR PREDITO*/ OUTPUT OUT=NSDS P=YE L95=LI U95=LS; /*P = PREDITO L95 = LIMITE INFERIOR U95 = LIMITE SUPERIOR */ RUN; PROC PRINT; RUN; PROC PLOT NOLEGEND; PLOT YE*X="*" LI*X="-" LS*X="-" Y*X="."/VPOS=20 HPOS=60 OVERLAY VAXIS=-10 TO 60 BY 10 HAXIS = 0 TO 5 BY 1 VREF=10 30; RUN;

EXERCÍCIO 2 /* LISTA 6- EXERCÍCIO 2 */ OPTIONS LS=78 PS=64 ; DATA L4E2; INFILE "A:L4E2.DAT"; INPUT X1-X6 Y; LABEL X1= "N" X2="P" X3="K" X4="CA" X5="MG" X6="S" Y="%TR"; PROC REG; MODEL Y=X1-X6/SELECTION=FORWARD SLSTAY=0.15; /*PEGA O MODELO E VAI AJUSTANDO*/ RUN; /* LISTA 6- EXERCÍCIO 2 */ OPTIONS LS=78 PS=64 ; DATA L4E2;

54 INFILE "A:L4E2.DAT"; INPUT X1-X6 Y; LABEL X1= "N" X2="P" X3="K" X4="CA" X5="MG" X6="S" Y="%TR"; PROC REG; *MODEL Y=X1-X6/SELECTION=FORWARD SLENTRY=0.15; /*PEGA O MODELO E VAI AJUSTANDO*/ *MODEL Y=X1-X6/SELECTION=BACKWARD SLSTAY=0.15; MODEL Y=X1-X6/SELECTION=STEPWISE SLENTRY=0.15 SLSTAY=0.15; OUTPUT OUT=SDSN P=YEST; RUN; PROC PRINT; RUN; PROC GPLOT; PLOT YEST*Y; RUN;

EXERCÍCIO 3 /* LISTA 6 - EXERCICIO 3 */ OPTIONS LS = 78 PS = 64; DATA L4E3; INPUT X Y @@; X2=X**2; X3=X**3; CARDS; 2.5 7.39 2.5 7.30 2.55 7.0 2.55 7.03 2.6 6.9 2.6 6.95 2.65 6.85 2.65 6.8 2.7 6.7 2.7 6.3 2.75 6.33 2.75 6.2 2.8 5.97 2.8 5.9 2.85 5.9 2.85 5.82 2.9 5.8 2.9 5.8 2.95 6.15 2.95 6 3 6.3 3 6.15 ; PROC PRINT; RUN; PROC PLOT; PLOT Y*X /VPOS=20 HPOS=60; RUN; PROC REG; MODEL Y=X X2 X3 /SELECTION = B SLSTAY=0.10; RUN;

55

VI.7. RESOLUÇÃO DA LISTA DE EXERCÍCIOS - AULA 7 EXERCÍCIO 1.A. /* "LISTA 7 = EXERCÍCIO 1" */ OPTIONS LS=78 PS=64; DATA L7E1A; INFILE "A:L7E1.TXT"; INPUT TR RP Y; PROC PRINT; RUN; PROC GLM; CLASS TR RP; MODEL YT = TR; OUTPUT OUT=NSDS R = ERRO; MEANS TR/HOVTEST = BARTLETT; RUN; PROCC UNIVARIATE NORMAL PLOT; VAR ERRO; RUN; PROC SORT DATA=L7E1; BY TR; RUN; PROC MEANS MEAN VAR NOPRINT; VAR Y; OUTPUT OUT=THV MEAN=M VAR=V; BY TR; RUN; DATA THV; SET THV; DP=SQRT(V); LM=LOG10(M); LV=LOG10(V); KEEP TR M V DP LM LV; PROC PRINT; RUN; PROC REG; MODEL LV=LM; RUN;

EXERCÍCIO 1.B. /* "LISTA 7 = EXERCÍCIO 2" */ OPTIONS LS=78 PS=64; DATA L7E1B; INFILE "A:L5E1.TXT"; INPUT TR RP Y; YT=LOG(Y); PROC PRINT; RUN; PROC GLM; CLASS TR RP; MODEL YT = TR; OUTPUT OUT=NSDS R = ERRO; MEANS TR/HOVTEST = BARTLETT DUNCAN; RUN; PROCC UNIVARIATE NORMAL PLOT;

56 VAR ERRO; RUN;

EXERCÍCIO 2.A. /* LISTAS - EXERCÍCIO 2.A*/ OPTIONS LS=78 PS=64; DATA L7E2A; INFILE "A:\L7E2.TXT"; INPUT TR BL Y; PROC PRINT; RUN; PROC GLM; CLASS TR BL; MODEL Y = BL TR/SS3; CONTRAST "(T1 + T2 + T3) VS (T1 + T5)" TR 2 2 2 -3 -3; CONTRAST "T4 VS T5 " TR 0 0 0 1 -1; CONTRAST "T1 VS (T2 + T3) " TR 2 -1 -1 0 0 ; CONTRAST "T2 VS T3 " TR 0 1 -1 0 0 ; MEANS TR/TUKEY; RUN;

EXERCÍCIO 2.B /* LISTAS - EXERCÍCIO 2.B */ OPTIONS LS=78 PS=64; DATA L7E2B; INFILE "A:\L7E2.TXT"; INPUT TR BL Y; IF TR=5 AND BL=2 THEN Y=.; PROC PRINT; RUN; PROC GLM; CLASS TR BL; MODEL Y = BL TR/SS3; LSMEANS TR/PDIFF ADJUST = TUKEY; RUN;

EXERCÍCIO 3 /* LISTA 7 - EXERCICIO 4 */ OPTIONS LS=78 PS=64; DATA L7E3; INFILE "A:\L7E3.TXT"; INPUT TR BL Y X; PROC PRINT; RUN; PROC GLM; CLASS TR BL; MODEL Y = BL TR X/SS3; /* COMO X NAO ESTA NA CLASSE O SAS INTERPRETA COMO UMA COVARIAVEL */ RUN; MEANS TR; LSMEANS TR/PDIFF ADJUST=TUKEY; RUN; EXERCÍCIO 4.A.

57 /* LISTA 7 - EXERCICIO 4 */ OPTIONS LS=78 PS=64; DATA L7E4A; INFILE "A:\L7E4.TXT"; INPUT A B BL Y; PROC PRINT; RUN; PROC GLM; CLASS A B BL; MODEL Y = BL A B A*B/SS3; RUN; /* DESDOBRAMENTO DA INTERAÇAO (3x2) */ DATA L7E4A; SET L7E4A; IF A=1 THEN BdA1=B; ELSE BdA1=0; IF A=2 THEN BdA2=B; ELSE BdA2=0; IF A=3 THEN BdA3=B; ELSE BdA3=0; IF B=1 THEN AdB1=A; ELSE AdB1=0; IF B=2 THEN AdB2=A; ELSE AdB2=0; PROC PRINT; RUN; TITLE "DESDOBRAMENTO A d. B"; PROC GLM; CLASS BL B AdB1-AdB2; MODEL Y=BL B AdB1(B) AdB2(B)/SS3; LSMEANS AdB1(B) AdB2(B)/PDIFF ADJUST=TUKEY; RUN; TITLE "DESDOBRAMENTO B d. A"; PROC GLM; CLASS BL A BdA1-BdA3; MODEL Y= BL A BdA1(A) BdA2(A) BdA3(A)/SS3; LSMEANS BdA1(A) BdA2(A) BdA3(A) /PDIFF ADJUST=TUKEY; RUN; EXERCÍCIO 4B – SUPONDO B ANINHADO EM A /* LISTA 7 - EXERCICIO 5 - SUPONDO QUE B ESTA ANINHADO EM A*/ OPTIONS LS=78 PS=64; DATA L7E4B; INFILE "A:\L7E4.TXT"; INPUT A B BL Y; PROC PRINT; RUN; PROC GLM; CLASS A B BL; MODEL Y = BL A B(A)/SS3; MEANS A/TUKEY; RUN; /*ESTUDO DO EFEITO DE BdA POR CONTRAST */ DATA L7E5; SET L7E5; TR=(A-1)*2+B; /* FORMULA GERAL: (A-1)*NB +A */ *PROC PRINT; RUN; PROC GLM; CLASS BL TR;

58 MODEL Y=BL TR/SS3; MEANS TR; CONTRAST "BdA1" TR 1 -1 0 0 0 0; CONTRAST "BdA2" TR 0 0 1 -1 0 0; CONTRAST "BdA3" TR 0 0 0 0 1 -1; RUN;

EXERCÍCIO 5 /* LISTA 7 - EXERCICIO 5 */ OPTIONS LS=78 PS=64; DATA L7E5; INFILE "A:\L7E5.TXT"; INPUT A B BL Y; PROC PRINT; RUN; PROC GLM; CLASS A B BL; MODEL Y= BL A BL*A B A*B/SS3; TEST H=BL A E=BL*A; MEANS A/TUKEY E=BL*A; MEANS B/TUKEY; RUN;

59

VI.8. RESOLUÇÃO DA LISTA DE EXERCÍCIOS - AULA 8 EXERCÍCIO 1 /* LISTA 8 - EXERCICIO 1 */ OPTIONS LS=78 PS=64; DATA L8E1; INFILE "A:\L8E1.TXT"; INPUT A B RP Y; IF A=. THEN DELETE; PROC PRINT; RUN; PROC GLM; CLASS A B RP; MODEL Y=A B(A)/SS3; MEANS A/TUKEY E=B(A); RANDOM A B(A)/TEST; RUN;

EXERCÍCIO 2 /* LISTA 8 - EXERCICIO 2 */ OPTIONS LS=78 PS=64; DATA L8E2; DO A=1 TO 2; DO B=1 TO 2; DO C = 1 TO 3; DO RP = 1 TO 5; INPUT Y@@; OUTPUT; END; END; END; END; CARDS; 569676898745866 8 6 9 7 5 7 8 9 6 8 10 8 7 9 6 576869769768 795 656787879665798 ; PROC PRINT; RUN; PROC GLM; CLASS A B C RP; MODEL Y=A B C(A) A*B B*C(A)/SS3; RANDOM A B C(A) B B*C(A)/TEST; RUN;

60

VI.9. RESOLUÇÃO DA LISTA DE EXERCÍCIOS - AULA 9 EXERCÍCIO 1 /* LISTA 9 - EXERCICIO 1 */ OPTIONS LS=78 PS=64; DATA L9E1; INFILE "A:L9E1.TXT"; INPUT TR RP Y; PROC PRINT; RUN; /* PROC GLM; CLASS TR RP; MODEL Y=TR/SS3; MEANS TR/HOVTEST=BARTLETT; RUN; */ PROC NPAR1WAY ANOVA WILCOXON; CLASS TR; VAR Y; RUN;

61

VI.10. RESOLUÇÃO DA LISTA DE EXERCÍCIOS - AULA 10 EXERCÍCIO 1 /* LISTA 10 - EXERCICIO 1 */ OPTIONS LS=78 PS=64; DATA L10E1; INFILE "A:L10E1.TXT"; INPUT T RP Y; PROC PRINT; RUN; PROC PLOT; PLOT Y*T="*"/HPOS=60 VPOS=20; RUN; /* OBSERBE QUE: LIM DEG (T TENTE A INFINITO) = A+B DEG(0)=A OBSERVANDO O GRAFICO, A+B ~ 60 E A ~ 15 */ TITLE "COM AS DERIVADAS"; PROC NLIN; PARAMETERS A=15 B=45 C=0.15; U=EXP(-C*T); MODEL Y=A+B*(1-U); DER.A=1; DER.B=1-U; DER.C=B*T*U; RUN; TITLE "SEM AS DERIVADAS"; PROC NLIN METHOD=DUD; PARAMETERS A=15 B=45 C=0.15; U=EXP(-C*T); MODEL Y=A+B*(1-U); OUTPUT OUT=VPAR PARMS=A B C; RUN; DATA VPAR; SET VPAR; YE=A+B*(1-EXP(-C*T)); PROC PRINT; RUN; PROC GPLOT; PLOT Y*T YE*T/OVERLAY LEGEND; SYMBOL1 COLOR=BLUE INTERPOL=NONE VALUE=PLUS ; SYMBOL2 COLOR=RED INTERPOL=SPLINE VALUE=DOT ; TITLE "VALORES OBSERVADOS E ESTIMADOS"; RUN;