Lineare Regression mit dem Paket lsquares

load(lsquares)$

daten3:matrix([192,99],[181,76],[174,80],[191,98],[185,89],[177,72],
              [181,81],[178,80],[198,106],[169,62],[171,72],[171,75])$

Berechnung der Koeffizienten

lin:lsquares_estimates(daten3,[x,y],y=a*x+b,[a,b]);

Einsetzen der Koeffizienten in eine Funktion f(x):=a*x+b

lin_term:a*x+b,lin[1];

define(f3(x),lin_term)$

Umwandlung der Daten, die für die Berechnung mit lsquares als Matrix gegeben waren,
in eine Liste für den Plot.

daten3_list:makelist(daten3[i],i,1,matrix_size(daten3)[1])$

Plot der Daten und der Regressionsgerade

wxplot2d([f3(x),[discrete,daten3_list]],[x,150,200],[y,60,130],
  [style,[lines,1,6],[points,2,2,2]]);