Nothing
### Create parameter sets for Hyperbolic
###
### require(GeneralizedHyperbolic)
mus <- c(-1,0,1,2)
deltas <- c(1,2,5,10)
xis <- rep(c(0.1,0.3,0.5,0.7,0.9), 1:5)
chis <- c(0,-0.25,0.25,-0.45,0,0.45,-0.65,-0.3,0.3,0.65,
-0.85,-0.4,0,0.4,0.85)
xiChis <- cbind(xis, chis)
lambdas <- 1
### hyperbLargeParam
maxrows <- length(mus)*length(deltas)*NROW(xiChis)*length(lambdas)
hyperbLargeParam <- matrix(nrow=maxrows,ncol=5)
rownum <- 1
for (i in 1:length(mus)){
for (j in 1:length(deltas)){
for (k in 1:NROW(xiChis)){
for (l in 1:length(lambdas)){
xi <- xiChis[k,1]
chi <- xiChis[k,2]
alpha <- (1 - xi^2)/(deltas[j]*xi*sqrt(xi^2 - chi^2))
beta <- alpha*chi/xi
param <- c(mus[i],deltas[j],round(alpha,3),round(beta,3),lambdas[l])
if (param[3] <= abs(param[4])) print("Error")
hyperbLargeParam[rownum,] <- param
rownum <- rownum + 1
}
}
}
}
hyperbLargeParam <- hyperbLargeParam[,-5]
### hyperbSmallParam
mus <- c(0,1)
deltas <- c(1,5)
xis <- rep(c(0.1,0.3,0.7), c(1,2,4))
chis <- c(0,-0.25,0.25,-0.65,-0.3,0.3,0.65)
xiChis <- cbind(xis, chis)
lambdas <- 1
maxrows <- length(mus)*length(deltas)*NROW(xiChis)*length(lambdas)
hyperbSmallParam <- matrix(nrow=maxrows,ncol=5)
rownum <- 1
for (i in 1:length(mus)){
for (j in 1:length(deltas)){
for (k in 1:NROW(xiChis)){
for (l in 1:length(lambdas)){
xi <- xiChis[k,1]
chi <- xiChis[k,2]
alpha <- (1 - xi^2)/(deltas[j]*xi*sqrt(xi^2 - chi^2))
beta <- alpha*chi/xi
param <- c(mus[i],deltas[j],round(alpha,3),round(beta,3),lambdas[l])
if (param[3] <= abs(param[4])) print("Error")
hyperbSmallParam[rownum,] <- param
rownum <- rownum + 1
}
}
}
}
hyperbSmallParam <- hyperbSmallParam[,-5]
### hyperbLargeShape
mus <- 0
deltas <- 1
xis <- rep(c(0.1,0.3,0.5,0.7,0.9), 1:5)
chis <- c(0,-0.25,0.25,-0.45,0,0.45,-0.65,-0.3,0.3,0.65,
-0.85,-0.4,0,0.4,0.85)
xiChis <- cbind(xis, chis)
lambdas <- 1
maxrows <- length(mus)*length(deltas)*NROW(xiChis)*length(lambdas)
hyperbLargeShape <- matrix(nrow=maxrows,ncol=5)
rownum <- 1
for (i in 1:length(mus)){
for (j in 1:length(deltas)){
for (k in 1:NROW(xiChis)){
for (l in 1:length(lambdas)){
xi <- xiChis[k,1]
chi <- xiChis[k,2]
alpha <- (1 - xi^2)/(deltas[j]*xi*sqrt(xi^2 - chi^2))
beta <- alpha*chi/xi
param <- c(mus[i],deltas[j],round(alpha,3),round(beta,3),lambdas[l])
if (param[3] <= abs(param[4])) print("Error")
hyperbLargeShape[rownum,] <- param
rownum <- rownum + 1
}
}
}
}
hyperbLargeShape <- hyperbLargeShape[,-5]
### hyperbSmallShape
mus <- 0
deltas <- 1
xis <- rep(c(0.1,0.3,0.7), c(1,2,4))
chis <- c(0,-0.25,0.25,-0.65,-0.3,0.3,0.65)
xiChis <- cbind(xis, chis)
lambdas <- 1
maxrows <- length(mus)*length(deltas)*NROW(xiChis)*length(lambdas)
hyperbSmallShape <- matrix(nrow=maxrows,ncol=5)
rownum <- 1
for (i in 1:length(mus)){
for (j in 1:length(deltas)){
for (k in 1:NROW(xiChis)){
for (l in 1:length(lambdas)){
xi <- xiChis[k,1]
chi <- xiChis[k,2]
alpha <- (1 - xi^2)/(deltas[j]*xi*sqrt(xi^2 - chi^2))
beta <- alpha*chi/xi
param <- c(mus[i],deltas[j],round(alpha,3),round(beta,3),lambdas[l])
if (param[3] <= abs(param[4])) print("Error")
hyperbSmallShape[rownum,] <- param
rownum <- rownum + 1
}
}
}
}
hyperbSmallShape <- hyperbSmallShape[,-5]
## save(hyperbLargeParam, hyperbSmallParam, hyperbLargeShape, hyperbSmallShape,
## file = "hyperbParam.rda")
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.