sie fits the breakpoint regression to the data,
richness is log(S+0.5) transformed, area is log(A)
transformed before analysis.
There is a coef, summary, print
method for fitted objects.
sieplot plots the observed (transformed)
data and the fitted line.
References
Lomolino, M. V., and M. D. Weiser. 2001.
Towards a more general species-area relationship:
diversity on all islands, great and small.
Journal of Biogeography, 28, 431--445.