Displaying 2 results from an estimated 2 matches for "qjohn".
Did you mean:
john
2010 Oct 03
1
Johnson Distribution Fit
...me the
follwing errors. Any help would be appreciated
#xi = xi
#lambda =l
#delta =d
#gamma = g
djohn = function(x,xi,l,d,g)
(d/(l*sqrt(2*pi)*((x-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) xi + (l*exp((qnorm(p) - g)/d))/(1 +
exp((qnorm(p) - g)/d))
f1c <- fitdist(data2,"john",start=list(xi = 0.5 ,l = 50, d = 1, g = 1))
Error in fitdist(data2, "john", start = list(xi = 0.5, l = 50, d = 1, :
the function mle failed to estimate the par...
2012 Jan 20
0
fit Johnson Sb with fitdist(method="mme")
...bility density function of
Johnson SB distribution (source:
http://www.ntrand.com/johnson-sb-distribution/)
(d/(l*sqrt(2*pi)*((x-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=...