Mercurial > repos > ecology > phylo_index
diff PhyloIndex.R @ 0:71a46f890035 draft
planemo upload for repository https://github.com/galaxyecology/tools-ecology/tree/master/tools/Phylodiversity_workflow commit 0de557d919c26eb0b5ab61504bc597d551503ac3
author | ecology |
---|---|
date | Tue, 20 May 2025 09:53:58 +0000 |
parents | |
children |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/PhyloIndex.R Tue May 20 09:53:58 2025 +0000 @@ -0,0 +1,41 @@ +#!/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() +}