"CcorrProf.MBparams.MBnetwork.MBdataDE.MBgibbsPars" <-
function(params,graph,data,gibbsPars,quiet=FALSE,scores.ref,contrasts=NULL,...)
{
# Compatibility checks
    isCompat1 <- compatObjects(params,graph)
    isCompat2 <- compatObjects(graph,data)
    isValid <- validObject(gibbsPars)

# Run Gibbs sampler on simulated data
    if(!quiet) cat("\nRunning gibbs sampler on simulated data...\n\n")
    graph <- calcCPDs(x=graph,data=data,contrasts)
    graph <- CgibbsSampler(x=graph,params=params,gibbsPars=gibbsPars,quiet=quiet,...)
    scores <- calcLods(graph)

# CORR analysis
    if(!quiet) cat("\nCalculating correlation...\n\n")
    #scores[is.na(scores)] <- min(scores,na.rm=TRUE)-1
    corr <- cor(scores.ref,scores,use="pairwise.complete.obs",method="spearman")

    return(corr)
}