Nothing
"vec2par" <-
function(vec, type, nowarn=FALSE, paracheck=TRUE, ...) {
if(type == 'aep4') {
z <- list(type = 'aep4', para = vec, source = "vec2par")
names(z$para) <- c("xi","alpha", "kappa", "h")
}
else if(type == 'cau') {
z <- list(type = 'cau', para = vec, source = "vec2par")
names(z$para) <- c("xi","alpha")
}
else if(type == 'emu') {
z <- list(type = 'emu', para = vec, source = "vec2par")
names(z$para) <- c("eta","mu")
}
else if(type == 'exp') {
z <- list(type = 'exp', para = vec, source = "vec2par")
names(z$para) <- c("xi","alpha")
}
else if(type == 'gam') {
z <- list(type = 'gam', para = vec, source = "vec2par")
if(length(vec) == 2) {
names(z$para) <- c("alpha","beta")
} else {
names(z$para) <- c("mu","sigma","nu")
#theta <- 1/(z$para[2]^2*z$para[3]^2)
#Q <- z$para[2]*z$para[3];
#d <- 1/(z$para[2]*Q)
#p <- z$para[3]*d;
#a <- exp(log(z$para[1]) - (log(d)-log(p))/p)
#mu <- a*exp(lgamma((d+1)/p)-lgamma(d/p))
#sd <- 1/(p*sqrt(d))
#z$other <- c(theta, Q, a, d, p, mu)
#names(z$other) <- c("theta", "Q", "Wiki(a)", "Wiki(d)", "Wiki(p)", "Wiki(mean)")
}
}
else if(type == 'gep') {
z <- list(type = 'gep', para = vec, source = "vec2par")
names(z$para) <- c("beta","kappa","h")
}
else if(type == 'gev') {
z <- list(type = 'gev', para = vec, source = "vec2par")
names(z$para) <- c("xi","alpha","kappa")
}
else if(type == 'gld') {
z <- list(type = 'gld', para = vec, source = "vec2par")
names(z$para) <- c("xi","alpha","kappa","h")
}
else if(type == 'glo') {
z <- list(type = 'glo', para = vec, source = "vec2par")
names(z$para) <- c("xi","alpha","kappa")
}
else if(type == 'gno') {
z <- list(type = 'gno', para = vec, source = "vec2par")
names(z$para) <- c("xi","alpha","kappa")
}
else if(type == 'gov') {
z <- list(type = 'gov', para = vec, source = "vec2par")
names(z$para) <- c("xi","alpha","beta")
}
else if(type == 'gpa') {
z <- list(type = 'gpa', para = c(vec[1],vec[2],vec[3]),
zeta=vec[4], source = "vec2par")
names(z$para) <- c("xi","alpha","kappa")
}
else if(type == 'gum') {
z <- list(type = 'gum', para = vec, source = "vec2par")
names(z$para) <- c("xi","alpha")
}
else if(type == 'kap') {
z <- list(type = 'kap', para = vec, source = "vec2par")
names(z$para) <- c("xi","alpha","kappa","h")
}
else if(type == 'kmu') {
z <- list(type = 'kmu', para = vec, diracdelta=NA, source = "vec2par")
names(z$para) <- c("kappa","mu")
}
else if(type == 'kur') {
z <- list(type = 'kur', para = vec, source = "vec2par")
names(z$para) <- c("alpha","beta")
}
else if(type == 'lap') {
z <- list(type = 'lap', para = vec, source = "vec2par")
names(z$para) <- c("xi","alpha")
}
else if(type == 'lmrq') {
z <- list(type = 'lmrq', para = vec, source = "vec2par")
names(z$para) <- c("mu", "alpha")
}
else if(type == 'ln3') {
z <- list(type = 'ln3', para = vec, source = "vec2par")
names(z$para) <- c("zeta","mu","sigma")
}
else if(type == 'nor') {
z <- list(type = 'nor', para = vec, source = "vec2par")
names(z$para) <- c("mu","sigma")
}
else if(type == 'pe3') {
z <- list(type = 'pe3', para = vec, source = "vec2par")
names(z$para) <- c("mu","sigma","gamma")
}
else if(type == 'pdq3') {
z <- list(type="pdq3", para=vec, source="vec2par")
names(z$para) <- c("xi","alpha","kappa")
}
else if(type == 'pdq4') {
z <- list(type="pdq4", para=vec, source="vec2par")
names(z$para) <- c("xi","alpha","kappa")
}
else if(type == 'ray') {
z <- list(type = 'ray', para = vec, source = "vec2par")
names(z$para) <- c("xi","alpha")
}
else if(type == 'revgum') {
z <- list(type = 'revgum', para = c(vec[1],vec[2]),
zeta=vec[3], source = "vec2par")
names(z$para) <- c("xi","alpha")
}
else if(type == 'rice') {
z <- list(type = 'rice', para = vec, source = "vec2par")
names(z$para) <- c("nu", "alpha")
}
else if(type == 'sla') {
z <- list(type = 'sla', para = vec, source = "vec2par")
names(z$para) <- c("xi","alpha")
}
else if(type == 'smd') {
z <- list(type = 'smd', para = vec, source = "vec2par")
names(z$para) <- c("xi", "a", "b", "q")
}
else if(type == 'st3') {
z <- list(type = 'st3', para = vec, source = "vec2par")
names(z$para) <- c("xi","alpha","nu")
}
else if(type == 'texp') {
z <- list(type = 'texp', para = vec, source = "vec2par")
names(z$para) <- c("xi","alpha","is.stationary")
}
else if(type == 'tri') {
z <- list(type = 'tri', para = vec, source = "vec2par")
names(z$para) <- c("min","mode","max")
}
else if(type == 'wak') {
z <- list(type = 'wak', para = vec, source = "vec2par")
names(z$para) <- c("xi","alpha","beta","gamma","delta")
}
else if(type == 'wei') {
z <- list(type = 'wei', para = vec, source = "vec2par")
names(z$para) <- c("zeta","beta","delta")
}
else {
z <- list(type = type, para = vec, source = "vec2par", byuser=TRUE)
return(z)
}
if(! paracheck) return(z)
if(! are.par.valid(z,nowarn,...)) {
ops <- options()
if(nowarn == TRUE) options(warn=-1)
warning("The parameters are invalid for the distribution")
options(ops)
return()
}
if(type == "kmu") {
if(! is.finite(vec[1])) {
names(z$para) <- c("kappa", "M")
z$diracdelta <- pdfkmu(0, z)
}
}
return(z)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.