\mnb150<nb wrap="0" txtwidth="75" autowidth="0" pprint="1" typeset="1" xface="Times New Roman" xsize="360" xbold="0" xitalic="0" xcolor="#0000FF" slidents="1" sizemult="71" minsize="160" ifont="0 280 255 0 49 Lucida Sans Typewriter" ofont="0 280 16711680 0 49 Lucida Sans Typewriter" tfont="0 280 0 0 34 Trebuchet MS"></nb>ÿ{\rtf1\ansi\deff0\deftab720{\fonttbl{\f0\fswiss MS Sans Serif;}{\f1\froman\fcharset2 Symbol;}{\f2\fmodern\fprq1 Lucida Sans Typewriter;}{\f3\fswiss\fprq2 Trebuchet MS;}{\f4\fswiss\fprq2 System;}}
{\colortbl\red0\green0\blue0;\red255\green0\blue0;\red0\green0\blue255;}
\deflang1033\pard\ri4\plain\f3\fs72\cf0 Chapter  6 Data visualization and 
\par elementary statistics\plain\f3\fs36\cf0 
\par 
\par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f2\fs36\cf1 {\pntext\f1\'b7\tab}READPATH:="D:\\\\My Documents\\\\mupad":
\par \pard\li600\ri1\fi-300\plain\f2\fs36\cf1 read("statistics.mu"):
\par read("frequency.mu"):
\par \pard\ri4\plain\f3\fs36\cf0 
\par 6.1 Data formats
\par 
\par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f2\fs36\cf1 {\pntext\f1\'b7\tab}data1 := [
\par \pard\li600\ri1\fi-300\plain\f2\fs36\cf1 ["Name","fat","sodium","carbo","sugar","rating"], 
\par ["Bran", 1, 130, 5, 6, 68], 
\par ["All_bran", 1, 260, 7, 5, 34], 
\par ["Cheerios", 2, 290, 17, 1, 51],
\par ["Cocoa puffs", 1, 180, 12, 13, 23], 
\par ["Corn flakes", 0, 290, 21, 2, 36],
\par ["Corn pops", 0, 90, 12, 12, 36],
\par ["Crispix", 0, 220, 21, 3, 47],
\par ["Frosted flakes", 0, 200, 14, 11, 31],
\par ["Grains pecan", 3, 75, 13, 4, 46],
\par ["Muesli", 3, 95, 16, 11, 37],
\par ["Muesli II", 3, 150, 16, 11, 34],
\par ["Oatmeal", 2, 170, 13.5, 10, 30],
\par ["Wheaties", 1, 200, 17, 3, 52]]
\par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f2\fs36\cf1 {\pntext\f1\'b7\tab}data2 := [130, 260, 290, 180, 290, 90, 220, 200, 75, 95, 150, 170, 200]
\par {\pntext\f1\'b7\tab}data2[4]
\par {\pntext\f1\'b7\tab}data1[4]
\par {\pntext\f1\'b7\tab}data1[4][3]
\par {\pntext\f1\'b7\tab}[data1[i][4] $ i=1..nops(data1)]
\par {\pntext\f1\'b7\tab}A := table(
\par \pard\li600\ri1\fi-300\plain\f2\fs36\cf1    (i=data1[i]) $ i=2..5
\par )
\par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f2\fs36\cf1 {\pntext\f1\'b7\tab}A[2]
\par {\pntext\f1\'b7\tab}A[M] := ([abcd,1234]):
\par {\pntext\f1\'b7\tab}A
\par {\pntext\f1\'b7\tab}nops(A), nops(data1), nops(data2)
\par {\pntext\f1\'b7\tab}data3 := import::readdata("Book1.csv",",")
\par \pard\ri4\plain\f3\fs36\cf0 
\par >>>>>>>>>> Time for Exercise 6.1   
\par 
\par 6.2 Plotting tabulated data
\par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f2\fs36\cf1 {\pntext\f1\'b7\tab}fat := [data3[i][2] $ i=2..nops(data3)]
\par {\pntext\f1\'b7\tab}sodium := [data3[i][3] $ i=2..nops(data3)]
\par {\pntext\f1\'b7\tab}carbo := [data3[i][4] $ i=2..nops(data3)]
\par {\pntext\f1\'b7\tab}sugar := [data3[i][4] $ i=2..nops(data3)]
\par {\pntext\f1\'b7\tab}rating := [data3[i][5] $ i=2..nops(data3)]
\par {\pntext\f1\'b7\tab}sodiumBar := plot::Bars2d(sodium):
\par \pard\li600\ri1\fi-300\plain\f2\fs36\cf1 plot(sodiumBar)
\par \pard\ri4\plain\f3\fs36\cf0     
\par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f2\fs36\cf1 {\pntext\f1\'b7\tab}carboBar := plot::Bars2d(carbo,
\par \pard\li600\ri1\fi-300\plain\f2\fs36\cf1    DrawMode=Horizontal,
\par    BarStyle=LinesPoints
\par ):
\par plot(carboBar)
\par \pard\ri4\plain\f3\fs36\cf0     
\par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f2\fs36\cf1 {\pntext\f1\'b7\tab}sugarPoints := [[i,sugar[i]] $ i=1..nops(sugar)]
\par {\pntext\f1\'b7\tab}sugarPlot := plot::Polygon2d(sugarPoints,
\par \pard\li600\ri1\fi-300\plain\f2\fs36\cf1    GridVisible=TRUE,
\par    YSubgridVisible=TRUE
\par ):
\par plot(sugarPlot)
\par \pard\ri4\plain\f3\fs36\cf0     
\par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f2\fs36\cf1 {\pntext\f1\'b7\tab}spices := [2, 2, 2, 2, 3, 5, 3]
\par {\pntext\f1\'b7\tab}labels := ["ginger","cumin","coriander","paprika",
\par \pard\li600\ri1\fi-300\plain\f2\fs36\cf1           "chili", "tomato", "garlic"]
\par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f2\fs36\cf1 {\pntext\f1\'b7\tab}spicesPlot := plot::Piechart2d(spices, 
\par \pard\li600\ri1\fi-300\plain\f2\fs36\cf1    Titles=labels
\par ):
\par plot(spicesPlot)
\par \pard\ri4\plain\f3\fs36\cf0     
\par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f2\fs36\cf1 {\pntext\f1\'b7\tab}data4 := import::readdata("accidents.csv",","):
\par \pard\ri4\plain\f3\fs36\cf0 
\par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f2\fs36\cf1 {\pntext\f1\'b7\tab}data4 := [data4[i] $ i=2..nops(data4)]
\par {\pntext\f1\'b7\tab}plot(pareto(data4))
\par {\pntext\f1\'b7\tab}allInOne := plot::Bars2d([fat, carbo, sugar, rating]):
\par \pard\li600\ri1\fi-300\plain\f2\fs36\cf1 plot(allInOne)
\par \pard\ri4\plain\f3\fs36\cf0     
\par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f2\fs36\cf1 {\pntext\f1\'b7\tab}allInOne := plot::Bars3d(
\par \pard\li600\ri1\fi-300\plain\f2\fs36\cf1    [fat, carbo, sugar, rating]):
\par plot(allInOne)
\par \pard\ri4\plain\f3\fs36\cf0     
\par >>>>>>>>>> Time for Exercise 6.2 
\par 
\par 6.3 Frequency distributions
\par 
\par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f2\fs36\cf1 {\pntext\f1\'b7\tab}ageList := [52, 19, 36, 9, 7, 3, 75, 63, 52, 39, 47, 80, 
\par \pard\li600\ri1\fi-300\plain\f2\fs36\cf1    59, 73, 81, 37, 52, 63, 69, 7, 47, 50, 48, 18, 39, 71, 
\par    51, 5, 43, 30, 43, 40, 28, 29, 25, 23, 31, 45, 50, 33, 
\par    49, 37, 39, 40, 31, 50, 23, 33, 29, 32, 50, 27, 39, 31, 
\par    21, 41, 41, 33, 42, 34, 24, 33, 20, 38, 26, 23, 20, 33, 
\par    23, 29, 48, 20, 41, 31, 47]
\par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f2\fs36\cf1 {\pntext\f1\'b7\tab}min(ageList)
\par {\pntext\f1\'b7\tab}max(ageList)
\par {\pntext\f1\'b7\tab}nops(ageList)
\par {\pntext\f1\'b7\tab}sort(ageList)
\par {\pntext\f1\'b7\tab}plot(dotPlot(ageList, 0, 85))
\par \pard\ri4\plain\f3\fs36\cf0     
\par Technical comment:
\par =================   
\par In this chapter, while talking about classes, we will symbolize and define 
\par them using boundaries (not upper and lower class limit), and we will require 
\par that our observations do not coincide with class boundaries. In order to be 
\par consistent throughout the book, we will declare classes using only one 
\par method---the range operator a..b. for example, 0..5 for the class with the 
\par lower boundary 0 and the upper boundary 5.
\par 
\par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f2\fs36\cf1 {\pntext\f1\'b7\tab}clasSize := 5:
\par \pard\li600\ri1\fi-300\plain\f2\fs36\cf1 startFrom := 0:
\par endOn := max(ageList) div clasSize:
\par ageClasses := [
\par    startFrom+clasSize*i..startFrom+clasSize*(i+1) $       
\par    i=0..endOn
\par ]
\par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f2\fs36\cf1 {\pntext\f1\'b7\tab}clasSize := 10:
\par \pard\li600\ri1\fi-300\plain\f2\fs36\cf1 startFrom:= 0-0.5:
\par endOn := max(ageList) div clasSize:
\par ageClasses10 := [  
\par    startFrom+clasSize*i..startFrom+clasSize*(i+1) $  
\par    i=0..endOn
\par ]
\par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f2\fs36\cf1 {\pntext\f1\'b7\tab}clasSize  := 5:
\par \pard\li600\ri1\fi-300\plain\f2\fs36\cf1 startFrom := 0-0.5:
\par endOn := max(ageList) div clasSize:
\par ageClasses5 := [  
\par    startFrom+clasSize*i..startFrom+clasSize*(i+1) $  
\par   i=0..endOn
\par ]
\par \pard\ri4\plain\f3\fs36\cf0 
\par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f2\fs36\cf1 {\pntext\f1\'b7\tab}stats::frequency(ageList, ageClasses10)
\par {\pntext\f1\'b7\tab}stats::frequency(ageList, ageClasses5)
\par {\pntext\f1\'b7\tab}ageHist10 := plot::Histogram2d(
\par \pard\li600\ri1\fi-300\plain\f2\fs36\cf1     ageList, Cells=ageClasses10
\par ):
\par \pard\ri4\plain\f3\fs36\cf0 
\par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f2\fs36\cf1 {\pntext\f1\'b7\tab}plot(ageHist10)
\par {\pntext\f1\'b7\tab}ageHist5 := plot::Histogram2d(
\par \pard\li600\ri1\fi-300\plain\f2\fs36\cf1   ageList, Cells=ageClasses5
\par ):
\par plot(ageHist5)
\par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f2\fs36\cf1 {\pntext\f1\'b7\tab}totalFreq := nops(ageList):
\par {\pntext\f1\'b7\tab}ageDensityHist10 := plot::Histogram2d(ageList,
\par \pard\li600\ri1\fi-300\plain\f2\fs36\cf1    Cells=ageClasses10,
\par    Area=totalFreq
\par ):
\par plot(ageDensityHist10)
\par \pard\ri4\plain\f3\fs36\cf0     
\par 
\par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f2\fs36\cf1 {\pntext\f1\'b7\tab}ageRelativeHist10 := plot::Histogram2d(ageList,
\par \pard\li600\ri1\fi-300\plain\f2\fs36\cf1    Cells=ageClasses10,
\par    Area=1
\par ):
\par plot(ageRelativeHist10)
\par \pard\ri4\plain\f3\fs36\cf0     
\par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f2\fs36\cf1 {\pntext\f1\'b7\tab}ageList2 := [41, 19, 72, 21, 11, 55, 7, 3, 25, 35, 23, 3, 19, 40, 57, 5, 29, 
\par \pard\li600\ri1\fi-300\plain\f2\fs36\cf1       56, 24, 47, 38, 42, 20, 26, 20, 23, 45, 35, 22, 33, 27, 20, 33, 46, 29, 
\par       33, 48, 40, 31, 42, 38, 40, 37, 25, 50, 23, 41, 30, 22, 45, 42, 39, 48, 
\par      26, 31, 28, 47, 24, 38, 27]
\par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f2\fs36\cf1 {\pntext\f1\'b7\tab}plot(dotPlotDouble(ageList, ageList2, 0, 85))
\par \pard\ri4\plain\f3\fs36\cf0     
\par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f2\fs36\cf1 {\pntext\f1\'b7\tab}ageRelativeHist10 := plot::Histogram2d(ageList,
\par \pard\li600\ri1\fi-300\plain\f2\fs36\cf1    Cells=ageClasses10,
\par    Area=1,
\par    ViewingBoxYMax=0.035
\par ):
\par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f2\fs36\cf1 {\pntext\f1\'b7\tab}ageRelativeHist10b := plot::Histogram2d(ageList2,
\par \pard\li600\ri1\fi-300\plain\f2\fs36\cf1    Cells=ageClasses10,
\par    Area=1,
\par    ViewingBoxYMax=0.035
\par ):
\par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f2\fs36\cf1 {\pntext\f1\'b7\tab}S1 := plot::Scene2d(ageRelativeHist10):
\par \pard\li600\ri1\fi-300\plain\f2\fs36\cf1 S2 := plot::Scene2d(ageRelativeHist10b):
\par plot(S1, S2)
\par \pard\ri4\plain\f3\fs36\cf0     
\par >>>>>>>>>> Time for Exercise 6.3  
\par 
\par 6.4 Measures of central tendency
\par 
\par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f2\fs36\cf1 {\pntext\f1\'b7\tab}reset():  // after reset load libraries again
\par {\pntext\f1\'b7\tab}READPATH:="D:\\\\My Documents\\\\mupad":
\par \pard\li600\ri1\fi-300\plain\f2\fs36\cf1 read("statistics.mu"):
\par read("frequency.mu"):
\par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f2\fs36\cf1 {\pntext\f1\'b7\tab}ageList := [52, 19, 36, 9, 7, 3, 75, 63, 52, 39, 47, 80, 59, 73, 81, 37, 52, 63, 69, 7, 
\par \pard\li600\ri1\fi-300\plain\f2\fs36\cf1     47, 50, 48, 18, 39, 71, 51, 5, 43, 30, 43, 40, 28, 29, 25, 23, 31, 45, 50, 33, 49, 
\par    37, 39, 40, 31, 50, 23, 33, 29, 32, 50, 27, 39, 31, 21, 41, 41, 33, 42, 34, 24, 33, 
\par    20, 38, 26, 23, 20, 33, 23, 29, 48, 20, 41, 31, 47]:
\par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f2\fs36\cf1 {\pntext\f1\'b7\tab}ageList2 := [41, 19, 72, 21, 11, 55, 7, 3, 25, 35, 23, 3, 19, 40, 57, 5, 29, 
\par \pard\li600\ri1\fi-300\plain\f2\fs36\cf1       56, 24, 47, 38, 42, 20, 26, 20, 23, 45, 35, 22, 33, 27, 20, 33, 46, 29, 
\par       33, 48, 40, 31, 42, 38, 40, 37, 25, 50, 23, 41, 30, 22, 45, 42, 39, 48, 
\par      26, 31, 28, 47, 24, 38, 27]
\par 
\par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f2\fs36\cf1 {\pntext\f1\'b7\tab}nops(ageList)
\par {\pntext\f1\'b7\tab}S := sum(ageList[i], i=1..75)
\par {\pntext\f1\'b7\tab}S := sum(ageList[i], i=1..nops(ageList))
\par {\pntext\f1\'b7\tab}meanAge := sum(
\par \pard\li600\ri1\fi-300\plain\f2\fs36\cf1    ageList[i],i=1..nops(ageList)
\par )/nops(ageList)
\par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f2\fs36\cf1 {\pntext\f1\'b7\tab}data := [46.4, 22.2, 66.1, 68.8, 5.7, 94.2, 11.0, 52.3, 74.1, 80.7, 49.5, 
\par \pard\li600\ri1\fi-300\plain\f2\fs36\cf1 84.2, 26.3, 72.8, 66.0, 48.3, 62.7, 19.1, 34.6, 42.9, 42.7, 17.0, 30.8, 
\par 41.8, 75.6, 61.8, 94.5, 96.4, 17.8, 86.7]:
\par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f2\fs36\cf1 {\pntext\f1\'b7\tab}factors:=[18, 12, 5, 6, 14, 4, 5, 19, 15, 1, 9, 1, 0, 2, 6, 13, 
\par \pard\li600\ri1\fi-300\plain\f2\fs36\cf1 14, 11, 16, 18, 4, 6, 7, 20, 9, 19, 18, 8, 20, 7]:
\par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f2\fs36\cf1 {\pntext\f1\'b7\tab}sum(
\par \pard\li600\ri1\fi-300\plain\f2\fs36\cf1   data[i]*factors[i], i=1..nops(data)
\par )/sum(factors[i], i=1..nops(data))
\par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f2\fs36\cf1 {\pntext\f1\'b7\tab}ageList := sort(ageList)
\par \pard\li600\ri1\fi-300\plain\f2\fs36\cf1 
\par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f2\fs36\cf1 {\pntext\f1\'b7\tab}nops(ageList)
\par {\pntext\f1\'b7\tab}median := ageList[floor(75/2)+1]
\par {\pntext\f1\'b7\tab}datatemp:=[1,2,3,4,5,6]:
\par \pard\li600\ri1\fi-300\plain\f2\fs36\cf1 median := (datatemp[nops(datatemp)/2]+ 
\par            datatemp[nops(datatemp)/2+1])/2
\par \pard\ri4\plain\f3\fs36\cf0 
\par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f2\fs36\cf1 {\pntext\f1\'b7\tab}stats::modal(ageList)
\par \pard\ri4\plain\f3\fs36\cf0 
\par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f2\fs36\cf1 {\pntext\f1\'b7\tab}quartiles := quantile(ageList, 4, 1),
\par \pard\li600\ri1\fi-300\plain\f2\fs36\cf1 quantile(ageList, 4, 2),
\par quantile(ageList, 4, 3)
\par \pard\ri4\plain\f3\fs36\cf0 
\par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f2\fs36\cf1 {\pntext\f1\'b7\tab}deciles := quantile(ageList, 10, i) $ i=1..9
\par {\pntext\f1\'b7\tab}B := plot::Boxplot(ageList, ageList2, 
\par \pard\li600\ri1\fi-300\plain\f2\fs36\cf1    DrawMode=Horizontal
\par ):
\par plot(B, ViewingBoxXMin=0)
\par \pard\ri4\plain\f3\fs36\cf0     
\par >>>>>>>>>> Time for Exercise 6.4   
\par 
\par 6.5 Standard deviation
\par 
\par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f2\fs36\cf1 {\pntext\f1\'b7\tab}ageList := [3, 5, 7, 7, 9, 18, 19, 20, 20, 20, 21, 23, 23, 23, 23, 24, 25, 26, 
\par \pard\li600\ri1\fi-300\plain\f2\fs36\cf1 27, 28, 29, 29, 29, 30, 31, 31, 31, 31, 32, 33, 33, 33, 33, 33, 34, 36, 37, 
\par 37, 38, 39, 39, 39, 39, 40, 40, 41, 41, 41, 42, 43, 43, 45, 47, 47, 47, 48, 
\par 48, 49, 50, 50, 50, 50, 51, 52, 52, 52, 59, 63, 63, 69, 71, 73, 75, 80, 81]
\par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f2\fs36\cf1 {\pntext\f1\'b7\tab} float(stats::meandev(ageList))
\par {\pntext\f1\'b7\tab}float(stats::stdev(ageList, Population))
\par \pard\ri4\plain\f3\fs36\cf0 
\par >>>>>>>>>> Time for Exercise 6.5   
\par 
\par 6.6 Correlation and regression
\par 
\par 
\par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f2\fs36\cf1 {\pntext\f1\'b7\tab}temperature := [27.4, 27.5, 27.3, 26.8, 26.2, 25.5, 24.9, 25.0, 
\par \pard\li600\ri1\fi-300\plain\f2\fs36\cf1    25.5, 26.1, 26.8, 27.1]:
\par  rainfall := [274.8, 236.0, 316.8, 395.5, 389.4, 420.4, 507.1, 
\par     535.4, 436.5, 407.9, 319.7, 339.9]:
\par  pressure := [1007.0, 1006.9, 1007.4, 1008.5, 1009.7, 1010.8, 
\par     1011.1, 1011.3, 1010.9, 1010.2, 1008.8, 1007.5]:
\par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f2\fs36\cf1 {\pntext\f1\'b7\tab}A := plot::Scatterplot(temperature, rainfall):
\par \pard\li600\ri1\fi-300\plain\f2\fs36\cf1 plot(A)
\par \pard\ri4\plain\f3\fs36\cf0     
\par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f2\fs36\cf1 {\pntext\f1\'b7\tab}A := plot::Scatterplot(pressure, temperature):
\par \pard\li600\ri1\fi-300\plain\f2\fs36\cf1 plot(A)
\par \pard\ri4\plain\f3\fs36\cf0     
\par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f2\fs36\cf1 {\pntext\f1\'b7\tab}A := plot::Scatterplot(pressure, rainfall):
\par \pard\li600\ri1\fi-300\plain\f2\fs36\cf1 plot(A)
\par \pard\ri4\plain\f3\fs36\cf0 
\par \pard\li300\ri5\fi-300{\*\pn\pnlvlblt\pnf1\pnindent300{\pntxtb\'b7}}\plain\f2\fs36\cf1 {\pntext\f1\'b7\tab}stats::linReg(temperature, rainfall)
\par {\pntext\f1\'b7\tab}stats::correlation(temperature, rainfall)
\par {\pntext\f1\'b7\tab}stats::correlation(pressure, temperature)
\par {\pntext\f1\'b7\tab}stats::correlation(pressure,rainfall)
\par \pard\ri4\plain\f3\fs36\cf0 >>>>>>>>>> Time for Exercise 6.6  
\par }