# NOT RUN {
library(magrittr)
# For the categorical logit model
depoDeltaB <- c(0.05, 0.08)
plot3logit:::DeltaB2pc_cat3logit(depoDeltaB) %>%
plot3logit:::pc2p0(depoDeltaB, flink = list(
linkfun = plot3logit:::linkfun_cat3logit,
linkinv = plot3logit:::linkinv_cat3logit))
# For the ordinal logit model
depoDeltaB <- 0.08
depoalpha <- c(-0.4, 0.4)
plot3logit:::DeltaB2pc_ord3logit(depoDeltaB, depoalpha) %>%
plot3logit:::pc2p0(depoDeltaB, flink = list(
linkfun = function(x) plot3logit:::linkfun_ord3logit(x, depoalpha),
linkinv = function(x) plot3logit:::linkinv_ord3logit(x, depoalpha)))
# }
Run the code above in your browser using DataLab