library("cusp") ### ### Example 1: The Attitudes data set ### data("attitudes") data("attitudeStartingValues") fit <- cusp(y ~ Attitude, alpha ~ Orient + Involv, beta ~ Involv, data = attitudes, start = attitudeStartingValues) summary(fit, logist = TRUE) plot(fit) ### ### Example 2: The Oliva et al. based data set ### # Model I data("oliva") oliva.fit1 <- cusp(y ~ z1 + z2 - 1, alpha ~ x1 + x2 + x3 - 1, beta ~ y1 + y2 + y3 + y4 - 1, data = oliva) summary(oliva.fit1) par(mar = c(5.1,4.6,4.1,2.1)) plot(oliva.fit1, what = "bifurcation") cusp3d(oliva.fit1, theta=150, Y = 2.95, B=5.2, np=601, n.surface=70, lwd=0.1) # Model II set.seed(121) oliva.fit2 <- cusp(y ~ z1 + z2, alpha ~ x1 + x2 + x3 + y1 + y2 + y3 + y4, beta ~ x1 + x2 + x3 + y1 + y2 + y3 + y4, data = oliva, start = rnorm(19)) summary(oliva.fit2, logist = TRUE) ### ### Example 3: Zeeman catastrophe machine data ### data("zeeman1") fit.zeeman <- cusp(y ~ z, alpha ~ x + y, beta ~ x + y, data = zeeman1) summary(fit.zeeman, logist = TRUE) cusp3d(fit.zeeman, surf.hue = 40, theta=215, phi=37.5, B=5.25)