Regressionsanalyse

Die Daten liegen als ASCII-Datei in zwei Spalten vor. Die erste Spalte enthält die x-Werte und die zweite Spalte die Funktionswerte.

restart:
with(stats):with(statplots):with(plots):
with(share):readshare(fit):
readlib(readdata):
data:=readdata(`datei`,zahlenformat,2):
`(* Angabe der Fitfunktion z.B. a0*1+a1*x+a2*x^2+a3*x^3 in Form einer Liste *)`:
lst:=[1,x,x^2,x^3]:
n:=nops(data):
dat:=array(data):
datx:=[seq(dat[i,1],i=1..n)]:
daty:=[seq(dat[i,2],i=1..n)]:
p1:=scatter2d(datx,daty):
f:=fit(datx,daty,lst,x);
p2:=plot(f,x=0..dat[n,1]):
display([p1,p2]);