LOESS Regression
The new release provides the following worksheet function to perform LOESS regression.
LOESS(Rx, Ry, Rx0, npts, deg): returns a column array with the fitted values of the LOESS regression defined by Rx, Ry, and npts for the x values in Rx0.
Here we assume that Rx is a column array or range containing x data values, Ry is a column array or range with the corresponding y values, Rx0 is a column array with any number of x values, and npts = the number of points in the span.
If Rx0 is omitted it defaults to the array values in Rx. deg = 1 (default) for weighted linear regression and deg = 2 for weighted quadratic regression. If npts is omitted (or set to zero) then npts is reset to (n+deg)/3 or the next larger integer where n = the # of elements in Rx (or Ry). Rx0 can also take a scalar x value, in which case the output is the fitted y value corresponding to this x value.
For more information, see LOESS Regression using Excel.
Runs Test with more than two categories
- RUNS_TEST(s, tails) = the p-value of the runs test for the sequence defined in the string s.
- RUNS_TEST(R1, tails) = the p-value of the runs test for the sequence defined by the column or row array or range R1.
- RUNS(s) = the number of runs in string s
- RUNS(R1) = the number of runs in the row or column array or range R1.
See Runs Test with more than Two Categories for more details.
Seasonal Kendall Test
The new release supports the Seasonal Kendall Test, which determines whether there is a monotonic trend for data that has seasonality. It is essentially the Mann-Kendall test for seasonal data. The following new function is now supported:
SK_TEST(R1, per, lab, tails, alpha): returns a column array with the values test statistic, s.e., z-stat, p-value, and trend.
The arguments for this function as well as the output are the same as for the MK_TEST function except that the per argument specifies the period of seasonality (e.g. per = 4 for quarterly, per = 12 for monthly, etc.).
See Seasonal Kendall Test for more details.
Queueing Models (analytic approach)
Exponential service rate
- MM1X(λ, μ, lab): M/M/1 queueing model
- MMs(λ, μ, s, lab): M/M/s queueing model
- MM1K(λ, μ, k, lab): M/M/1 model with a finite queue of size k
- MMsK(λ, μ, s, k, lab): M/M/s model with a finite queue of size k
- MM1N(λ, μ, nn, lab): M/M/1 model with a finite # of customers equal to nn
- MMsN(λ, μ, s, nn, lab): M/M/s model with a finite # of customers equal to nn
Non-exponential service rate
- MD1X(λ, μ, lab): M/D/1 queueing model; constant (deterministic) service rate
- MG1X(λ, μ, σ, lab): M/G/1 queueing model; general service rate; σ = standard deviation of the service time
- MGinf(λ, μ, lab): M/G/∞ queueing model; general service time, infinite # of servers
Priority queueing
- MM1P(R1, μ, lab): M/M/1 queueing model with non-preemptive priority queueing
- MMsP(R1, μ, s, lab): M/M/s queueing model with non-preemptive priority queueing
- MM1PS(R1, R2, lab): M/M/1 queueing model with non-preemptive priority queueing with variable service means
- MM1PP(R1, μ, lab): M/M/1 queueing model with preemptive priority queueing
For the following new function, the output consists of a row array with the Wq values for each priority class.
- MMsPP(R1, μ): M/M/s queueing model with preemptive priority queueing
Probability of n customers being served
- MM1Pn(λ, μ, n, cum): M/M/1 queueing model
- MMsPn(λ, μ, s, n, cum): M/M/s queueing model
- MM1KPn(λ, μ, k, n, cum): M/M/1 model with a finite queue of size k
- MMsKPn(λ, μ, s, k, n, cum): M/M/s model with a finite queue of size k
- MM1NPn(λ, μ, nn, n, cum): M/M/1 model with a finite # of customers equal to n
- MMsNPn(λ, μ, s, nn, n, cum): M/M/s model with a finite # of customers equal to n
- MD1Pn(λ, μ, n, cum): M/D/1 queueing model
- MGinfPn(λ, μ, n, cum): M/G/∞ queueing model
Probability of service or queueing times
Finally, the following new worksheet functions have been added. In the following, w = the time a customer remains in the system, and wq = the time a customer remains in the queue.
- MM1W(λ, μ, t) = P(w > t) for an M/M/1 queueing model
- MM1Wq(λ, μ, t) = P(wq > t) for an M/M/1 queueing model
- MMsW(λ, μ, s, t) = P(w > t) for an M/M/s queueing model
- MMsWq(λ, μ, s, t) = P(wq > t) for an M/M/s queueing model
- MD1W(λ, μ, t) = P(w > t) for an M/D/1 queueing model
- MD1Wq(λ, μ, t) = P(wq > t) for an M/D/1 queueing model
- MM1KW(λ, μ, k, t) = P(w > t) for an M/M/1 queueing model with a finite queue of size k
- MM1KWq(λ, μ, k, t) = P(wq > t) for an M/M/1 queueing model with a finite queue of size k
- MMsKWq(λ, μ, s, k, t) = P(wq > t) for an M/M/s queueing model with a finite queue of size k
See Queueing Theory for more details.
Queueing Models (simulation approach)
- QUEUE_SIM(R1, R2, lab); single server model for R1 and R2
- QUEUE_SIMX(R1, R2, s, lab): multi-server model for R1 and R2 where s = # of severs
Since these functions employ a simulation, they can be used with a wide variety of G/G/1 and G/G/s queueing systems. For example, if you want arrivals to follow a gamma distribution, you would place the formula =GAMMA.INV(RAND(), alpha, beta) in all the cells of R1 where alpha and beta are the parameters of the gamma distribution. The accuracy of these functions depends on the size of the R1 and R2 arrays.
POISSON_INV Enhancement
The existing POISSON_INV worksheet function has been revised to improve the processing speed.
Bug Fixes
- Eliminated errors related to TRUE/FALSE values in the output of some of the data analysis tools that were inadvertently introduced in Rel 8.5 (e.g. the Correlation data analysis tool). This includes the elimination of extraneous punctuation marks.