>
|
n := 20:
m := 5:
tdx := 5:
nvar := 5:
k := 3:
tdc := 5:
maxit := 10:
x := Matrix([[77.3, 13, 9.699999999999999, 1.5, 6.4], [82.5, 10, 7.5, 1.5, 6.5], [66.90000000000001, 20.6, 12.5, 2.3, 7], [47.2, 33.8, 19, 2.8, 5.8], [65.3, 20.5, 14.2, 1.9, 6.9], [83.3, 10, 6.7, 2.2, 7], [81.59999999999999, 12.7, 5.7, 2.9, 6.7], [47.8, 36.5, 15.7, 2.3, 7.2], [48.6, 37.1, 14.3, 2.1, 7.2], [61.6, 25.5, 12.9, 1.9, 7.3], [58.6, 26.5, 14.9, 2.4, 6.7], [69.3, 22.3, 8.4, 4, 7], [61.8, 30.8, 7.4, 2.7, 6.4], [67.7, 25.3, 7, 4.8, 7.3], [57.2, 31.2, 11.6, 2.4, 6.5], [67.2, 22.7, 10.1, 3.3, 6.2], [59.2, 31.2, 9.6, 2.4, 6], [80.2, 13.2, 6.6, 2, 5.8], [82.2, 11.1, 6.7, 2.2, 7.2], [69.7, 20.7, 9.6, 3.1, 5.9]], datatype=float[8], order='C_order'):
isx := Vector([1, 1, 1, 1, 1], datatype=integer[kernelopts('wordsize')/8]):
cmeans := Matrix([[82.5, 10, 7.5, 1.5, 6.5], [47.8, 36.5, 15.7, 2.3, 7.2], [67.2, 22.7, 10.1, 3.3, 6.2]], datatype=float[8], order='C_order'):
wt := Vector([], datatype=float[8]):
inc := Vector(20, datatype=integer[kernelopts('wordsize')/8]):
nic := Vector(3, datatype=integer[kernelopts('wordsize')/8]):
css := Vector(3, datatype=float[8]):
csw := Vector(3, datatype=float[8]):
NAG:-g03efc(x, isx, cmeans, inc, nic, css, csw, 'n' = n, 'm' = m, 'tdx' = tdx, 'nvar' = nvar, 'k' = k, 'tdc' = tdc, 'wt' = wt, 'maxit' = maxit):
|