Link - Frau Prof. Dr. rer. nat. Kira Klenke - Hochschule Hannover

SAS-Programme
TABLES (f1_10beh_w_a f1_10beh_w_h -- f1_10beh_w_j f1_11verhalt_a --
f1_11verhalt_c) * f1_3legelinie/nocum NOROW nopercent;
%FormLab;
RUN;
TITLE 'Futter pro 1000 Durchschnittsh. u. Wo. nach Legelinie';
PROC TABULATE DATA=ADS;
CLASS f1_3legelinie;
VAR futter1000wokg;
TABLE futter1000wokg*(n='Anzahl Durchgänge' Mean='arithm. Mittel'
median='Median' std='Standardabw.'cv='Variationskoeffizient'
min='Minimum' p10='10%-Quantil' p90='90%-Quantil'
max='Maximum'),f1_3legelinie/rts=20 CONDENSE;
%FormLab;
RUN;
ODS RTF CLOSE;
B.4 Ausreißertest
/*************************************************************************
Programmdatei: Ausreissertest.sas
Autor: Jan Schäl
Austeißertest nach Grubbs/Beck
**************************************************************************/
*Logarithmierung der Variablen verluste1000wo;
DATA neu;
SET ADS_A;
LG_Verl=LOG10(verluste1000wo);
RUN;
PROC SORT DATA=neu;
BY f1_4haltung;
RUN;
*Berechnung der Residuen für die Variable LG_Verl;
PROC GLM DATA=neu noprint;
CLASS f1_4haltung;
MODEL LG_Verl = f1_4haltung;
OUTPUT OUT=Resi R=Residuen;
RUN;
QUIT;
*Berechnung der erforderlichen Kenngrößen für den Grubbs-Test;
PROC MEANS DATA=Resi NOPRINT;
BY f1_4haltung;
VAR Residuen;
OUTPUT OUT=Grubbs MEAN=xq STD=s;
RUN;
PROC SORT DATA=Grubbs;
BY f1_4haltung;
RUN;
PROC SORT DATA=Resi;
BY f1_4haltung;
RUN;
*Berechnung des Kritischen Werts T getrennt nach Haltungssystemen;
134