R/ggm.estimate.pcor.R

In GeneNet: Modeling and Inferring Gene Networks

### ggm.estimate.pcor  (2006-06-02)
###
###     Small-samples estimate of GGM partial correlation coefficients
###
### Copyright 2003-06 Juliane Schaefer, Rainer Opgen-Rhein, and Korbinian Strimmer
###
###
### This file is part of the `GeneNet' library for R and related languages.
### It is made available under the terms of the GNU General Public
### License, version 2, or at your option, any later version,
### incorporated herein by reference.
### 
### This program is distributed in the hope that it will be
### useful, but WITHOUT ANY WARRANTY; without even the implied
### warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
### PURPOSE.  See the GNU General Public License for more
### details.
### 
### You should have received a copy of the GNU General Public
### License along with this program; if not, write to the Free
### Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
### MA 02111-1307, USA


# estimate partial correlation coefficients
ggm.estimate.pcor <- function(x, method=c("static", "dynamic"), ...)
{
  method <- match.arg(method)

  ########
 
 
  if (method == "dynamic")
  {  
     return( dyn.pcor(x, ...) )  
  } 
   
  if (method == "static")
  {  
     return( pcor.shrink(x, ...) )  
  } 
 
 
  ########
  
  #if (method == "observed.pcor")
  #{    
  #  return( pcor.shrink(x, lambda=0) )
  #} 

  ########
  
  #if (method == "partial.bagged.cor")
  #{
  #   cor.bag <- cor.bagged(x, R=R, ...)   
  #   return( cor2pcor(cor.bag) )
  #} 
  
  ########

  #if (method == "bagged.pcor")
  #{  
  #   return( pcor.bagged(x, R=R, ...) )
  #} 
}