quarta-feira, 22 de maio de 2019
MANOVA
MANOVA
data imc_dat;
input cat $ imc corr kcal;
cards;
AT 20.2 60.7 3200
AT 21.3 54.8 3100
AT 19.3 49.6 2800
AT 21.1 52.3 3300
SEM 22.4 14.9 2600
SEM 21.9 17.8 2700
SEM 23.8 18.6 3200
SEM 24.1 15.1 3300
SE 27.3 2.5 2700
SE 23.4 4.3 2300
SE 25.2 2.3 2600
SE 26.4 2.6 3200
PR 26.2 4.1 2600
PR 24.2 2.1 2700
PR 25.4 1.9 2650
;
proc print;
run;
proc glm;
class cat;
model imc corr kcal = cat;
contrast " Atl e Semiat Vs Seden e Prof" cat 1 -1 -1 1;
contrast " Professor Vs Sedentario" cat 0 1 -1 0;
manova h=_all_ / printe printh;
contrast " Atl e Semiat Vs Seden e Prof" 1 -1 -1 1;
contrast " Professor Vs Sedentario" cat 0 1 -1 0;
run;
Correlação Residual MANOVA:
Outro Exemplo de Contrastes, Ordem Alfabética:
- Análise Multivariada Exemplos:
var imc corr kcal; ==> Kruskal Wallis
Arquivo Word para Download:
data imc_dat;
input cat $ imc corr kcal;
cards;
AT 20.2 60.7 3200
AT 21.3 54.8 3100
AT 19.3 49.6 2800
AT 21.1 52.3 3300
SEM 22.4 14.9 2600
SEM 21.9 17.8 2700
SEM 23.8 18.6 3200
SEM 24.1 15.1 3300
SE 27.3 2.5 2700
SE 23.4 4.3 2300
SE 25.2 2.3 2600
SE 26.4 2.6 3200
PR 26.2 4.1 2600
PR 24.2 2.1 2700
PR 25.4 1.9 2650
;
proc print;
run;
proc glm;
class cat;
model imc corr kcal = cat;
contrast " Atl e Semiat Vs Seden e Prof" cat 1 -1 -1 1;
contrast " Professor Vs Sedentario" cat 0 1 -1 0;
manova h=_all_ / printe printh;
contrast " Atl e Semiat Vs Seden e Prof" 1 -1 -1 1;
contrast " Professor Vs Sedentario" cat 0 1 -1 0;
run;
/*
contrast " Atl e Semiat Vs Seden e Prof" cat 1 -1 -1 1;
Pega a ordem alfabetica dos niveis dos fatores para os coeficientes
dos contrastes
Se colocar o contraste depois do comando MANOVA faz os contrastes univariados
tambem
contrast " Atl e Semiat Vs Seden e Prof" cat 1 -1 -1 1;
Pega a ordem alfabetica dos niveis dos fatores para os coeficientes
dos contrastes
Se colocar o contraste depois do comando MANOVA faz os contrastes univariados
tambem
*/
Correlação Residual MANOVA:
Partial Correlation Coefficients from the Error SSCP Matrix / Prob > |r|
|
|||
DF = 11
|
imc
|
corr
|
kcal
|
imc
|
1.000000
|
0.017729
0.9564
|
0.697770
0.0116
|
corr
|
0.017729
0.9564
|
1.000000
|
0.132496
0.6814
|
kcal
|
0.697770
0.0116
|
0.132496
0.6814
|
Outro Exemplo de Contrastes, Ordem Alfabética:
- Análise Multivariada Exemplos:
Comparar resultados ANOVA com Kruskal-Wallis.
data imc_dat;
input cat $ imc corr kcal;
cards;
AT 20.2 60.7 3200
AT 21.3 54.8 3100
AT 19.3 49.6 2800
AT 21.1 52.3 3300
SEM 22.4 14.9 2600
SEM 21.9 17.8 2700
SEM 23.8 18.6 3200
SEM 24.1 15.1 3300
SE 27.3 2.5 2700
SE 23.4 4.3 2300
SE 25.2 2.3 2600
SE 26.4 2.6 3200
PR 26.2 4.1 2600
PR 24.2 2.1 2700
PR 25.4 1.9 2650
;
proc print;
run;
proc npar1way data=imc_dat wilcoxon dscf;
class cat;
var imc corr kcal = cat;
run;
var imc corr kcal = cat; ==> ANOVA, GLM, MANOVA
Tem um erro nesse programa, onde esta? Compare com os slides apresentados antes do programa SAS (Dois slides de Kruskal Wallis com trechos de programa SAS).
var imc corr kcal; ==> Kruskal Wallis
Resultados de Kruskal Wallis
Resultados ANOVA dados brutos, transformados e da Estat. Robusta. Tese doutorado do Gabriel.
quarta-feira, 15 de maio de 2019
Cluster Analysis Exemplo Gabriela Trivellato
data propriedades;
input cat $ DMERCALI DMERCINSU MELHORAGUA CIRCCURTOS ORGAGRICUL ORGSOCIEDADE;
/*P1 = produtor agroecológico
P2 = produtor orgânico
P3 = grande produtor agroindústria produção convencional
P4 = pequeno produtor agricultura familiar produção convencional
DMERCALI = dependência do mercado para aquisição de alimentos
DMERCINSU = dependência do mercado para aquisição de insumos
MELHORAGUA = papel dos agricultores para melhoria de serviços ecossitêmicos - água
CIRCCURTOS = circuitos curtos de comercialização
ORGAGRICUL = produção de alimentos orgânicos para os produtores
ORGSOCIEDADE = produção de alimentos orgânicos para a sociedade
FOI CONSIDERADA ESCALA DE LIKERT PARA AS RESPOSTAS DE 1 A 5
5 = IDENTIFICO MUITO
4 = IDENTIFICO
3 = IDENTIFICO POUCO
2 = IDENTIFICO POUQUÍSSIMO
1 = NÃO IDENTIFICO */
cards;
P1 2 1 4 5 5 5
P2 4 4 4 2 1 5
P3 5 5 1 1 1 1
P4 4 4 4 3 1 1
;
proc cluster data=propriedades outtree = arvore method = average;
var DMERCALI DMERCINSU MELHORAGUA CIRCCURTOS ORGAGRICUL ORGSOCIEDADE;
id cat;
run;
PROC TREE DATA = arvore;
RUN;
input cat $ DMERCALI DMERCINSU MELHORAGUA CIRCCURTOS ORGAGRICUL ORGSOCIEDADE;
/*P1 = produtor agroecológico
P2 = produtor orgânico
P3 = grande produtor agroindústria produção convencional
P4 = pequeno produtor agricultura familiar produção convencional
DMERCALI = dependência do mercado para aquisição de alimentos
DMERCINSU = dependência do mercado para aquisição de insumos
MELHORAGUA = papel dos agricultores para melhoria de serviços ecossitêmicos - água
CIRCCURTOS = circuitos curtos de comercialização
ORGAGRICUL = produção de alimentos orgânicos para os produtores
ORGSOCIEDADE = produção de alimentos orgânicos para a sociedade
FOI CONSIDERADA ESCALA DE LIKERT PARA AS RESPOSTAS DE 1 A 5
5 = IDENTIFICO MUITO
4 = IDENTIFICO
3 = IDENTIFICO POUCO
2 = IDENTIFICO POUQUÍSSIMO
1 = NÃO IDENTIFICO */
cards;
P1 2 1 4 5 5 5
P2 4 4 4 2 1 5
P3 5 5 1 1 1 1
P4 4 4 4 3 1 1
;
proc cluster data=propriedades outtree = arvore method = average;
var DMERCALI DMERCINSU MELHORAGUA CIRCCURTOS ORGAGRICUL ORGSOCIEDADE;
id cat;
run;
PROC TREE DATA = arvore;
RUN;
Horários de Consulta:
Horários de
Consulta:
Skype:
hc_gabriel
Segundas: 18-19
Quintas: 18-19
Segundas: 18-19
Quintas: 18-19
Sábados:
14:15- 15:15
Domingos: 21:15-23:15Podemos marcar outros horários se necessário, conversamos por Whatsapp: o meu é 019-988-627-438.
quarta-feira, 8 de maio de 2019
Aula 15/5 - Cluster Cerveja Envelhecida
Aula 15/5 - Cluster An. Cerveja Envelhecida
Outros Cases
MANOVA se der tempo
Outros Cases
MANOVA se der tempo
Assinar:
Postagens (Atom)