Fitted Lines

xxxxxxxxxx

 

1

A,B=1+round(random(),2),round(.01+random(),2)

2

N=500; X=[.1*i*random() for i in range(N)]

3

Y=[A*x+B+10*random()*(-1)^randint(0,1) for x in X]

4

var('a,b,x'); F=a*x+b

5

ff=find_fit(list(zip(X,Y)),F,variables=[x]) 

6

a,b=float(str(ff[0])[5:]),float(str(ff[1])[5:])

7

F=[a*x+b for x in X]

8

pf=points(zip(X,F),color='#ff36ff',size=5,

9

          legend_label='fitted line')

10

pr=points(zip(X,Y),color='#3636ff',marker='*',

11

          size=9,legend_label='real data') 

12

(pf+pr).show(figsize=(6,5),gridlines=True,

13

             title=['[A == %s,B == %s]'%(A,B),ff])

Language: Sage Gap GP HTML Macaulay2 Maxima Octave Python R Singular

---

Messages