R/network_gen.R

In networkABC: Network Reverse Engineering with Approximate Bayesian Computation

#' @title Random scale-free network generation. This function is used intensively in the abc function.
#' @description  Generate random network topology
#' @param number_genes A number
#' @param clust_coef A number
#' @return A list with the number of of genes, the targeted clustering coefficient and the resulting network
#' @examples
#' network_gen(10,1)
#' @export
#' @useDynLib networkABC C_run



network_gen<-function(number_genes,clust_coef){
  result=.C(C_run,number_genes= as.integer(number_genes), clust_coef=as.numeric(clust_coef), network=as.integer(1:(number_genes*number_genes)))
      #No longer needed since the symbol is registered in the NAMESPACE
      #,PACKAGE="networkABC")
  
  result$network<-matrix(result$network,result$number_genes,result$number_genes,byrow=TRUE)
return(result)
}