##---- Should be DIRECTLY executable !! ----##-- ==> Define data, use random,##-- or do help(data=index) for the standard data sets.## The function is currently defined asfunction (x,bend=1.28,op=2){n<-length(x)mestse<-sqrt(sum((ifmest(x,bend,op=2)^2))/(n*(n-1)))mestse}