R/inference_vertical_transmission.R

In BPerezLamarque/HOME: Phylogenetic Models of HOst-Microbiota Evolution

inference_vertical_transmission <-
function(mu,name,index,sequences){
  
  host_tree <- NULL
  eig_val <- NULL
  eig_vect <- NULL
  ivp <- NULL
  propinv <- NULL
  
  load(paste("data/data_model_",name,"_",index,".RData",sep=""))
  
  n=nrow(sequences)
  N=ncol(sequences)
  
  fit <- LL_tree(mu,host_tree,sequences,n,N,eig_val, eig_vect, ivp, propinv)
  
  return(fit)
}