S-ML distribution
dSml = function(x, p, log = FALSE) #log of pdf of S-ML model
{
n=count(x)
loglik <- (log(pi/2)+log(p)-log(1+p)-(2*p*x)+log((1+p)*exp(x*p)*(2*p*x)-1)+
log(sin((pi/2)*(1+exp(-p*x)*(p*x)/(1+p))*exp(-x*p))))
if ( log == FALSE)
```

```
density <- exp(loglik)
else density <- loglik
return(density)
}
theta <- maxlogL(x =st, dist = "dSml",start = 0.59)
summary(theta)