# NOT RUN {
#Only one group
beta <- matrix(c(1,10,5,0.1,1,0.5,0.2,2,1),nrow=3)
L_Phi <- class_group(beta)
group_types(beta) #gives c(1), NULL and NULL
#Two groups
n <- 3
beta <- diag(1,n)
beta[1,2] <- -0.32
beta[2,1] <- 1/beta[1,2]
L_Phi <- class_group(beta)
group_types(beta) #gives NULL, c(2) and c(1)
#with data
data(data_sim_RegNeg_1grpNeg1grpPos)
group_types(data_sim_RegNeg_1grpNeg1grpPos$param$beta)
# }
Run the code above in your browser using DataLab