Displaying 1 result from an estimated 1 matches for "memplog".
Did you mean:
mem_log
2012 Jan 20
0
fit Johnson Sb with fitdist(method="mme")
...xi)/l)*(1-((x-xi)/l))))*exp(-0.5*(g +
d*log(((x-xi)/l)/(1-((x-xi)/l))))^2)
}
pjohn <- function(x,xi,l,d,g) {
pnorm(g + d*log(((x-xi)/l)/(1-((x-xi)/l))))
}
qjohn <- function(p,xi,l,d,g) {
qnorm(xi + (l*exp((qnorm(p) - g)/d))/(1 + exp((qnorm(p) - g)/d)))
}
library(moments)
memplog <- function(x,order) {
ifelse(order==1, mean(x), moment(x,order,central=TRUE))
}
fjsb <-
fitdist(c(data1),"john",method="mme",order=c(1,2),memp="memplog",start=list
(xi=0,l=max(data1),d=1,g=1))
But it results in:
> fjsb <-
fitdist(c(data1),"...