|
|
|
|
|
|
|
|
Die Daten können als Zahlenwerte in einer Spaltendatei in ASCII-code abgelegt sein.
restart:
readlib(readdata):
with(stats): with(transform): with(describe):
with(statplots):with(plots):
data:=readdata(`dateiname`,Zahlenformat,1);
`(* Sortieren der Daten *)`:
d1:=sort(data);
`(* Zahl der Daten *)`:
n:=count(data);
`(* Kumulative Summe *)`:
cs:=seq(sum(d1[i],i=1..m),m=1..n);
`(* Berechnung der Häufigkeit Intervallgröße 1 *)`;
dl:=op(d1):
ma:=max(dl):
mi:=min(dl):
d2:=tallyinto(d1,[seq(i..(i+1),i=mi..(ma+1))]);
`(* Häufigkeitstabelle *)`;
d3:=frequency(d2);
histdata:=scaleweight[1./count(d2)](d2):
histogram(histdata);
`(* Klasseneinteilung mit den Klassengrenzen k1,k2,k3 usw. *)`:
k1:= :
k2:= :
k3:= :
ma:=max(op(d1)):
tallyinto(d1,[k1..k2,k2..k3,k3..(ma+1)]);
frequency(");
`(* Klasseneinteilung mit äquidistanten Klassengrenzen *)`;
schritt:= ; `(* Breite einer Klasse *)`:
start:=min(op(d1)) ; `(* Start der Klasseneinteilung *)`:
fin:=round(max(op(d1))/schritt) ; `(* Ende der Klasseneinteilung *)`:
st:=round(start/schritt);
fn:= fin-st:
h1:=tallyinto(d1,[seq((i*schritt+start)..((i+1)*schritt+start),i=0.. fn)]) ;
frequency(h1);
histgram2:=scaleweight[1./count(h1)](h1):
`(* Gesamtfläche des Histogramms ist auf eins normiert *)`;
histogram(histgram2);
