wissi stress inversion implemented

Author: tobiste

NAMESPACE

@@ -385,6 +385,7 @@ export(slip_inversion_hansen_boot)
 export(slip_inversion_michael)
 export(slip_inversion_misfit)
 export(slip_inversion_simple)
+export(slip_inversion_wissi)
 export(slip_inversion_yamaji_sato)
 export(slip_inversion_yamaji_sato_boot)
 export(slip_tendency)
@@ -439,6 +440,7 @@ export(watson_MLE)
 export(watson_inference)
 export(wellner)
 export(wellner_inference)
+export(wissi_polyphase)
 export(woodcock_plot)
 import(data.table)
 import(ggplot2)

R/stress_inversion.R

@@ -67,24 +67,27 @@
 #'   res_yamaji <- slip_inversion(x, method = "yamaji")
 #'
 #'   res_hansen <- slip_inversion(x, method = "hansen", type = "6d")
+#'   
+#'   res_wissi <- slip_inversion(x, method = 'wissi')
 #'
 #'   stereoplot(guides = FALSE)
 #'   fault_plot(x, col = "gray80")
 #'   points(res_michael$principal_axes, pch = 1:3, col = 2)
 #'   points(res_angelier$principal_axes, pch = 1:3, col = 3)
 #'   points(res_yamaji$principal_axes, pch = 1:3, col = 4)
 #'   points(res_hansen$principal_axes, pch = 1:3, col = 5)
+#'   points(res_wissi$principal_axes, pch = 1:3, col = 6)
 #'   legend("topleft",
 #'     pch = 1,
-#'     legend = c("Michael (1984)", "Angelier (1990)", "Yamaji & Sato (2006)", "Hansen (2013)"),
-#'     col = 2:5
+#'     legend = c("Michael (1984)", "Angelier (1990)", "Yamaji & Sato (2006)", "Hansen (2013)", "WISSI"),
+#'     col = 2:6
 #'   )
 #'   legend("bottomright",
 #'     pch = 1:3,
 #'     legend = c("S1", "S2", "S3")
 #'   )
 #' }))
-slip_inversion <- function(x, method = c("michael", "angelier", "hansen", "yamaji"), ...) {
+slip_inversion <- function(x, method = c("michael", "angelier", "hansen", "yamaji", "wissi"), ...) {
   method <- match.arg(method)
 
   if (method == "angelier") {
@@ -93,6 +96,8 @@ slip_inversion <- function(x, method = c("michael", "angelier", "hansen", "yamaj
     slip_inversion_michael(x, ...)
   } else if(method == "yamaji") {
     slip_inversion_yamaji_sato(x, ...)
+  } else if(method == "wissi"){
+    slip_inversion_wissi(x, ...)
   } else {
     slip_inversion_hansen(x, ...)
   }