Programa SAS para Outliers

 

data habit;

input Perfil $ IDH;

datalines;

Cafe 0.778

Cafe 0.731

Cafe 0.779

Cafe 0.761

Cafe 0.751

Cafe 0.729

Cafe 0.696

Cafe 0.682

Cafe 0.689

Cafe 0.706

Cafe 0.697

Cafe 0.661

Cafe 0.724

Cafe 0.678

Cafe 0.735

Cafe 0.78

Cafe 0.73

Cafe 0.761

Cafe 0.787

Cafe 0.714

Cafe 0.797

Cafe 0.778

Cafe 0.767

Cafe 0.551

Cafe 0.678

Cafe 0.56

Cafe 0.545

Cafe 0.718

Cafe 0.7

Cafe 0.714

Cafe 0.702

Cafe 0.743

Cafe 0.716

Cafe 0.713

Cafe 0.71

Cafe 0.687

Cafe 0.7

Cafe 0.697

Bovino 0.594

Bovino 0.7

Bovino 0.736

Bovino 0.708

Soja 0.744

Soja 0.618

Soja 0.754

Soja 0.754

Soja 0.755

Cana  0.772

Cana  0.712

Cana  0.74

Cana  0.788

Frango 0.782

Frango 0.719

Frango 0.671

Frango 0.768

Milho 0.744

Milho 0.758

Milho 0.734

Milho 0.757

 

;

proc print;

run;

/* input Perfil $ IDH; */

proc glm;

   class Perfil;

   model IDH = Perfil;

   output out=diagnosticos rstudent=res_stud;

     run;

proc print data=diagnosticos;

run;

data dados_limpos;

    set diagnosticos;

    if abs(res_stud) <= 2.57; /* Critério comum: remover valores > 2.57 desvios

                                Equivale a 99% de aceitação. descarto 1%. 0.5%

                                de cada lado */

run;

proc print data=dados_limpos;

run;