python - Least square fit a 2D line -


i realize use numpy find line so:

import numpy np import matplotlib.pyplot plt = np.array([1,2,3,4,6,7]) b = np.array([5,4,3,2,-2,-1]) k,m = np.polyfit(a,b,1) plt.scatter(a,b) plt.plot([0,10],[m,10*k+m]) plt.show()

enter image description here

but i'd use raw python code instead. math rusty, if can done in few lines of code i'd appreciate help!

if looking simple linear regression based on minimizing quadratic error, pure python implementation pretty straightforward (check equations α , β on link above):

def linear_fit(x, y):     """for set of points `(xi, yi)`, return linear polynomial `f(x) = k*x + m`     minimizes sum of quadratic errors.     """     meanx = sum(x) / len(x)     meany = sum(y) / len(y)     k = sum((xi-meanx)*(yi-meany) xi,yi in zip(x,y)) / sum((xi-meanx)**2 xi in x)     m = meany - k*meanx     return k, m

for sample input:

>>> x = [1,2,3,4,6,7] >>> y = [5,4,3,2,-2,-1] >>> linear_fit(x, y) (-1.1614906832298135, 6.285714285714285)

Comments

Post a Comment

Popular posts from this blog

python - Operations inside variables -

this question has answer here: evaluating mathematical expression in string 12 answers i want make python script finds maximum possible output using 2 operations. problem says there syntax error @ 2nd last line, in b in there. how can fix it? x = 0 = int(input("1st number:")) c = int(input("2nd number:" )) e = int(input("3rd number:" )) in range(4): if == 0: b = "+" elif == 1: b = "-" elif == 2: b = "/" else: b = "*" j in range(4): if j == 0: d = "+" elif j == 1: d = "-" elif j == 2: d = "/" else: d = "*" k = b c d e print(k) if you're after variables e on single line, couldn...
Read more

Generic Map Parameter java -

i having lot of methods of type: public static list<eduusers> getdetails(class c, map<string, integer> params) { list<eduusers> details = null; session session = hibernateutil.getsessionfactory().opensession(); transaction tx = null; try { tx = session.begintransaction(); criteria cr = session.createcriteria(c); (map.entry<string, integer> entry : params.entryset()) { cr.add(restrictions.eq(entry.getkey(), entry.getvalue())); } details = cr.list(); tx.commit(); } catch (exception asd) { log.debug(asd.getmessage()); if (tx != null) { tx.commit(); } } { session.close(); } return details; } i trying have generic method them , have written far: public static <t> list<t> getdetails(class<t> c, class<t> b) { l...
Read more

arrays - What causes a java.lang.ArrayIndexOutOfBoundsException and how do I prevent it? -

what arrayindexoutofboundsexception mean , how rid of it? here code sample triggers exception: string[] name = {"tom", "dick", "harry"}; for(int = 0; i<=name.length; i++) { system.out.print(name[i] +'\n'); } your first port of call should documentation explains reasonably clearly: thrown indicate array has been accessed illegal index. index either negative or greater or equal size of array. so example: int[] array = new int[5]; int boom = array[10]; // throws exception as how avoid it... um, don't that. careful array indexes. one problem people run thinking arrays 1-indexed, e.g. int[] array = new int[5]; // ... populate array here ... (int index = 1; index <= array.length; index++) { system.out.println(array[index]); } that miss out first element (index 0) , throw exception when index 5. valid indexes here 0-4 inclusive. correct, idiomatic for statement here be: for (int index = 0; index <...
Read more