get_partial_correlation_from_precision = function(Theta,selfEdges=FALSE)
{
# by default, does not return self edges (diagonal is set to zero)
ggm = -cov2cor(Theta)
if(!selfEdges)
ggm[diag(ggm)] = 0
return(ggm)
}
ggm[diag(ggm)] = 0 needs to be changed to diag(ggm)=0
@katehoffshutta
ggm[diag(ggm)] = 0
needs to be changed todiag(ggm)=0
Thanks, Chen