Program:

#import the required header files

from graphics import *

import math as m

#definition of "graphWin" method

def graphWin(Title):

    #Create a graphics window

    win = GraphWin(Title, 400, 400)

    #set the coords

    win.setCoords(-10, -10, 10, 10)

    #tell user to select multiple locations on the screen to designate points

    message = Text(Point(-3, 8), "Click to delineate points on the graph.")

    #draw the point

    message.draw(win)

    #call the "Line" method

    axisX = Line(Point(-10,0), Point(10,0))

    #call the "draw" method

    axisX.draw(win)

    #call the "Line" method

    axisY = Line(Point(0,10), Point(0,-10))

    #call the "draw" method

    axisY.draw(win)

    #call the "Rectangle" method

    r = Rectangle(Point(-9, -9), Point(-7,-8))

    #call the "draw" method

    r.draw(win)

    #get the centre position

    rCenter = r.getCenter()

    #set the text

    stopMouse = Text(rCenter, "Done")

    #call the "draw" method

    stopMouse.draw(win)

    #accept input from user until <done> button is pressed

    click = Point(0,0)

    #set the list

    allPoints = []

    #check the condition

    while True:

        #get the mouse action

        click = win.getMouse()

        #check the condition

        if ((-9 <= click.getX() <= -7) and (-9 <= click.getY() <= -8)):

            #break the loop

            break

        #otherwise

        else:

            #store user points in an appended list

            allPoints.append(click)

            #call the "draw" method

            click.draw(win)

    #return the values

    return allPoints, win

#definition of "average" method

def average(allPoints):

    #set the values for the variables

    sumX = 0

    sumY = 0

    count = 0

    sumXiYi = 0

    sumSqXi = 0

    sumSqYi = 0

    #iterate "i" until it reaches "allPoints" value

    for i in allPoints:

        #get the "x" and "y" value

        x = i.getX()

        y = i.getY()

        #calculate the "sumX" value

        sumX = sumX + x

        #calculate the "sumY" value

        sumY = sumY + y

        #increment the value

        count = count + 1

        #calculate the "xy" value

        xy = x * y

        #calculate the "sumXiYi" value

        sumXiYi = sumXiYi + xy

        #calculate the "SqX" value

        SqX = x * x

        #calculate the "sumSqXi" value

        sumSqXi = sumSqXi + SqX

        #calculate the "SqY" value

        SqY = y * y

        #calculate the "sumSqYi" value

        sumSqYi = sumSqYi + SqY

    #calculate the "a" and "b" value

    a = ((sumY * sumSqXi) - (sumX * sumXiYi)) / (count * (sumSqXi) - sumX ** 2)

    b = ((count * sumXiYi) - (sumX * sumY)) / (count * (sumSqXi) - sumX ** 2)

    #return the values

    return a, b

#definition of "main" method

def main():

    #get the text for the window

    allPoints, win = graphWin("Regression Line")

    #call the method

    a, b = average(allPoints)

    #set the values

    x1 = -10

    x2 = 10

    #call the "Line" method

    regressLine = Line(Point(x1, (a + b * x1)), Point(x2, (a + b * x2)))

    #call the "draw" method

    regressLine.draw(win)

    #get the mouse action

    win.getMouse()

#call the "main" method

main()