## Fitting a 4-dimensional HR MST tree
my_graph <- igraph::graph_from_adjacency_matrix(
rbind(c(0, 1, 0, 0),
c(1, 0, 1, 1),
c(0, 1, 0, 0),
c(0, 1, 0, 0)),
mode = "undirected")
n <- 100
Gamma_vec <- c(.5,1.4,.8)
complete_Gamma(Gamma = Gamma_vec, graph = my_graph) ## full Gamma matrix
set.seed(123)
my_data <- rmpareto_tree(n, "HR", tree = my_graph, par = Gamma_vec)
my_fit <- mst_HR(my_data, p = NULL, cens = FALSE)
Run the code above in your browser using DataCamp Workspace