view PhyloIndex.R @ 1:1d73a4af71f5 draft default tip

planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Phylodiversity_workflow commit ce74877d21915627f36eae4d2eb3851a6a858aae
author ecology
date Thu, 22 May 2025 16:43:58 +0000
parents 71a46f890035
children
line wrap: on
line source

#!/bin/Rscript
# Indice de phylodiversite


# args 
args = commandArgs(trailingOnly=TRUE)
#args = c("comm_matrix", "comm_tree.txt")

# library
library(Matrix)
library(ape)
library(phyloregion)
library(SparseArray)

# functions
write_results <- function(){
  write.table(ses_PD, output_file, sep = "\t", row.names=FALSE)
}

# verif args
if (length(args)<5){stop('Usage : input_tree_file input_matrix_file random_seed clustering_model output_file')
}else{
  # read enter file
  comm_tree <- read.tree(args[1])

  comm_matrix <- readSparseCSV(args[2] , sep="\t")
  comm_matrix <- as(comm_matrix,"dgCMatrix")

  seed <- args[3]

  model <- args[4]

  output_file <- args[5]

  # rdm + compute
  set.seed(seed)
  ses_PD <- PD_ses(comm_matrix, comm_tree, model = model, reps = 999)
  
  # write res 
  write_results()
}