# NOT RUN {
x <- seq(1, 5, 0.1)
y <- c()
for (theta in x) y <- c(y, structural_function(theta, c(2, 3), log, sqrt, theta))
plot(x, y)
lines(x, log(x), col = "blue")
lines(x, sqrt(x), col = "red")
####
f <- function(theta) {
p <- c(1, 1)
structural_function(
theta,
c(15, 20),
function(p) CD_A(alpha = 5, Beta = c(0.6, 0.4), p),
function(p) CD_A(alpha = 15, Beta = c(0.6, 0.4), p),
p
)
}
tmp <- sapply(1:25,f)
matplot(t(tmp), type="l")
# }
Run the code above in your browser using DataLab