I have a function XOLD() from a nearly verbatim port of legacy
FORTRAN in a package. I have remplemented this function as XNEW()
using much cleaner native R and built-in functions of R. I have
switched the package to the XNEW(), but for historical reasons would
like to retain the XOLD()
Amina,
If you can follow this code, you will see that the black dots scatter
around an L-skew of zero and L-kurtosis of 0.122602, which are the
theoretical L-moments of the standard normal distribution.
library(lmomco)
lmrdia - lmrdia()
plotlmrdia(lmrdia,ylim=c(-.2,.2),xlim=c(-.2,.2))
lms -
I am teaching grad/undergrad course in which I am introducting
students to R as a side topic; however, I am simultaneously
supporting Windows/MacOSX/Linux for the students.
Under some conditions, the Windows users can not read files
originating from MacOSX. The end-of-line differences
Paul,
Package lmomco fits generalized pareto (three parameter) using method
of L-moments. I suspect that other packages that Brian identified
use method of moments or other.
William
On Sep 3, 2006, at 10:55 AM, Paul Smith wrote:
Dear All
I am trying to fit Pareto distribution to some
CRAN et al.,
I would like to add an extented introduction or other arbitrary
sections to my package lmomco.
I have been shipping inst/doc/Introduction.Rd. I would like to have
this content inserted to the front of the PDF build for the CRAN. The
R-exts.pdf seems to be a little silent on
here.
Thanks,
William Asquith (author package lmomco)
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
I have question (curiosity) regarding returned values of R's qcauchy
() function,
for nonexceedance probability (F). It seems the ideal returned range
of cauchy distribution should be [-Inf,Inf].
For F=0
qcauchy(0)
[1] -Inf
but for F=1
qcauchy(1)
[1] 8.16562e+15
It seems to me that the