produce lowess smoothed functions
Lowess(X, Y, options)
L := Lowess(...)
L(x1, x2, ..., xN)
Matrix, Array, DataFrame, or listlist of data convertible to numeric. If XY is a Matrix, DataFrame, or Array, it must have at least two columns. If XY is a listlist, then each inner list is considered one row, with the kth of each inner list being the elements of the kth column, and the same requirement of at least two columns applies. Each row of XY is interpreted as one data point. If XY has m columns, then the first m−1 columns contain the values of the m−1 independent variables. The last column contains the corresponding values of the dependent variable.
Vector, list, Matrix, Array, DataSeries, or listlist of data convertible to numeric with the number of rows equal to the length of Y. The columns are interpreted as the values of the independent variables.
Vector, list, DataSeries, or Array of data convertible to numeric with length equal to the number of rows in X. The values in Y are the corresponding values of the dependent variable.
(optional) equation(s) of the form option=value where option is one of fitorder, bandwidth, or iters
x1, x2, ..., xN
evaluates L at (x1, x2, ..., xN), where N is equal to m−1 or the number of columns in X.
a matrix with number of columns equal to the number of columns in X or m−1. Returns a Vector where the ith element is L evaluated with the ith row of M as arguments
The degree of the polynomial used in each local regression. The default value is 1.
The proportion of the input data points used in each local regression. The default value depends on fitorder and the number of input data points.
The number of iterations when smoothing data of one independent variable. Each iteration makes the data smoother by eliminating outliers, thus making the computation more robust. This option has no effect when the data has more than one independent variable. The default is 2.
The Lowess command creates a function whose values represent the input data smoothed with the lowess algorithm.
Suppose the input data set XY is of m independent variables and has n data points, the lowess smoothed value at x≔x1,x2,…,xm is computed as follows.
Take the n⋅bandwidth points in XY that are closest to x.
Fit a polynomial P of m variables and degree fitorder to the points using weighted linear least squares, where the weight for a point w is computed by applying the tri-cube weight function to the distance between x and w.
Running one or more iterations, as specified by the iters option, will produce a set of weights to reduce the influence of outliers (that is, make the computation more robust). At each iteration, the weight of a point depends on the residual of the Lowess curve at that point in the previous iteration. These weights are combined with the weights given by the distance, as described previously.
L will return unevaluated if the arguments are non-convertible to numerics. But if the first and only argument is a Matrix with number of columns equal to the number of parameters of L, a Vector will be returned where the ith element is the L applied with the ith row of the Matrix as arguments.
Create a data sample and apply to it some error.
X ≔ Sample⁡Uniform⁡−2,2,200:
Y ≔ Sample⁡Uniform⁡−2,2,200:
Zerror ≔ Sample⁡Normal⁡0,0.1,200:
Z ≔ `~``*`⁡X,` $`,map⁡exp,−`~``^`⁡X,` $`,2−`~``^`⁡Y,` $`,2+Zerror:
XYZ ≔ Matrix⁡X,Y,Z,datatype=float8%T
Create the function whose graph is the smoothed surface.
L ≔ Lowess⁡XYZ,fitorder=2,bandwidth=0.3:
Plot the data sample, smoothed surface, and the region between the plane z=0.4 and the surface for −1.5≤x≤−0.5 and −1≤y≤1.
P ≔ ScatterPlot3D⁡XYZ:
Q ≔ plot3d⁡L,−2..2,−2..2,grid=25,25:
R ≔ plots:-shadebetween⁡L⁡x,y,0.4,x=−1.5..−0.5,y=−1..1,showboundary=false,negativeonly:
Find the volume of the shaded region.
For a two dimensional example we will create another data sample.
X ≔ Sample⁡Uniform⁡0,Pi,200
Yerror ≔ Sample⁡Normal⁡0,0.1,200
Y ≔ map⁡sin,X+Yerror
Create the function whose graph is the smoothed curve.
L ≔ CurveFitting:-Lowess⁡X,Y,fitorder=1,bandwidth=0.3:
Plot the data sample, smoothed curve, and the region between the x-axis and the curve for π8≤x≤3⁢π8.
P ≔ ScatterPlot⁡X,Y
Q ≔ plot⁡L⁡x,x=0..Pi
R ≔ plots:-shadebetween⁡L⁡x,0,x=Pi8..3⁢Pi8,showboundary=false,positiveonly
Find the area of the shaded region.
And find the maximum.
The Statistics[Lowess] command was introduced in Maple 2015.
For more information on Maple 2015 changes, see Updates in Maple 2015.
Download Help Document