File:Linear least squares example2.png

Summary

Description
English: Illustration of linear least squares.
Date
Source Own work
Author Oleg Alexandrov
Other versions
File:Linear least squares example2.svg is a vector version of this file. It should be used in place of this PNG file when not inferior.Category:Vector version available

File:Linear least squares example2.png → File:Linear least squares example2.svg

For more information, see Help:SVG.

In other languages
Alemannisch  العربية  беларуская (тарашкевіца)  български  বাংলা  català  нохчийн  čeština  dansk  Deutsch  Ελληνικά  English  British English  Esperanto  español  eesti  euskara  فارسی  suomi  français  Frysk  galego  Alemannisch  עברית  हिन्दी  hrvatski  magyar  հայերեն  Bahasa Indonesia  Ido  italiano  日本語  ქართული  한국어  lietuvių  македонски  മലയാളം  Bahasa Melayu  မြန်မာဘာသာ  norsk bokmål  Plattdüütsch  Nederlands  norsk nynorsk  norsk  occitan  polski  prūsiskan  português  português do Brasil  română  русский  sicilianu  Scots  slovenčina  slovenščina  српски / srpski  svenska  தமிழ்  ไทย  Türkçe  татарча / tatarça  українська  vèneto  Tiếng Việt  中文  中文(中国大陆)  中文(简体)  中文(繁體)  中文(马来西亚)  中文(新加坡)  中文(臺灣)  +/−
New SVG image

PNG development
InfoField

Licensing

Public domain This work has been released into the public domain by its author, Oleg Alexandrov. This applies worldwide.
In some countries this may not be legally possible; if so:
Oleg Alexandrov grants anyone the right to use this work for any purpose, without any conditions, unless such conditions are required by law.
Category:Self-published work#Linear%20least%20squares%20example2.pngCategory:PD-self#Linear%20least%20squares%20example2.png

Source code (MATLAB)

% Illustration of linear least squares.
function main()
 
   % KSmrq's colors
   red    = [0.867 0.06 0.14];
   blue   = [0, 129, 205]/256;
   green  = [0 100 0]/256;
   green=[0 139 0]/256;
   
   yellow = [254, 194,   0]/256;
   white = 0.99*[1, 1, 1];
   gray = 0.1*white;
 
   Xr = [1, 2, 3, 4];
   Yr = round(1+1.5*Xr+4*rand(1, 4))
   %Yr = [9  13  15 14];
   Yr = [6 5 7 10];
   
   A = min(Xr)-1; B = max(Xr)+1;
   C = min(Yr)-1; D = max(Yr)+1;
   N = 100;
   small = 0.2;
   X = linspace(A, B, N);
 
 
   % least squares fitting
   Yr = Yr';
   Xr=Xr';
   Mat = [(0*Xr+1) Xr];
   V=Mat'*Yr;
   V
   Mat
   V=(Mat'*Mat)\V;
   pe = V(1); qe = V(2);

   V
   
   % Set up the figure
   lw = 5; % linewidth
   lw2 = lw/2;
   fs = 22; % font size
   figure(1); clf; hold on;
   set(gca, 'fontsize', fs);
   set(gca, 'linewidth', lw2)
   hold on; grid on;
   axis equal;
 

   grid off;

   for i=1:length(Xr)
      plot([Xr(i), Xr(i)], [Yr(i) pe+qe*Xr(i)], 'color', green, 'linewidth', 0.95*lw);
   end

   myrad = 0.12;
   for i=1:length(Xr)
      ball(Xr(i), Yr(i), myrad, red);
   end

   plot(X, pe+qe*X, 'b', 'linewidth', lw);
   
   set(gca, 'GridLineStyle', '-', 'xcolor', gray);
   set(gca, 'GridLineStyle', '-', 'ycolor', gray);
   set(gca, 'XTick', A:1:B);
   set(gca, 'YTick', C:1:D);

   
   Yr - (pe+qe*Xr)
   
   axis equal;
   axis([A B C D]);
 
   saveas(gcf, 'Linear_least_squares_example.eps', 'psc2'); % save as eps
 
 
function ball(x, y, r, color)
   Theta=0:0.1:2*pi;
   X=r*cos(Theta)+x;
   Y=r*sin(Theta)+y;
   H=fill(X, Y, color);
   set(H, 'EdgeColor', 'none')
Category:Linear regression Category:Files by User:Oleg Alexandrov from en.wikipedia Category:Images with MATLAB source code
Category:Files by User:Oleg Alexandrov from en.wikipedia Category:Images with MATLAB source code Category:Linear regression Category:PD-self Category:PNG created with MATLAB Category:Self-published work Category:Vector version available