proc(xx,yy1,yy2)=idxgraph(out,x)
xx=matrix(max(x[,1])+1)*0
yy1=xx
yy2=xx
i=1
while(i <= max(x[,1])+1)
xx[i]= i-1
yy1[i]= max(paf(out.uu'[,1], x[,1]==i-1))
yy2[i]= max(paf(out.uu'[,2], x[,1]==i-1))
i=i+1
endo
endp
library("xclust")
library("plot")
z=read("butterfly")
x=z[,2:3]
c=2
m=2
e=0.01
out=xcfcgk(x,c,m,e) //apply fuzzy gustafson-kessel clustering
out.uu'
tt=idxgraph(out,x)
tyy1=tt.xx~tt.yy1
tyy2=tt.xx~tt.yy2
tyy1 = setmask(tyy1, "line", "red", "thin")
tyy2 = setmask(tyy2, "line", "blue", "medium")
disp = createdisplay(1,1)
show(disp,1,1, tyy1, tyy2)
strm = string("%1.0f", m)
setgopt(disp,1,1, "title","Membership Function for Fuzzy GK with m="+ strm)