subintensity_matrix <- matrix(c(-1.5, 1.5, 0,
0, -1, 1,
0, 0, -0.5),
ncol = 3,
byrow = TRUE)
ph1 <- PH(subintensity_matrix)
print(ph1)
Run the code above in your browser using DataLab