workflow4metabolomics / anova Goto Github PK

Jean-François relevently suggested to add graphical outputs for interactions (currently there are only plots for individual effects). He also highlighted that means of interactions are not provided in the variableMetadata file.

I am not fully satisfied in the way manova() works in case of NA.
Maybe we could consider using directly aov() in a loop or an "apply" instead of using manova().

Just to say that on my Firefox there is
System Message: ERROR/3 (, line 69)
in the help section of the wraper

cheers

@melpetera in an email - 15/09/15

Point 1 : paramètre Mode

Ce paramètre est prévu a priori pour pouvoir faire des anova au choix sur les lignes et colonnes (good point d'un point de vue purement stat). Mais du coup, comme maintenant il ne prend en entrée comme metadata que celles des samples (vu qu'orienté métabo), ben on ne peut plus les faire que pour les variables (une anova par variable). En metabo c'est ce qu'on fait, donc no problem, mais dans l'absolu du coup il faudrait :

• soit supprimer le paramètre "Mode" (car du coup useless si on ne fait que des anova par variable) => easy to be done (c'est ce que je recommande perso)
• soit garder la possibilité d'aussi faire des anova par sample (pour des gens qui voudraient utiliser l'outil pour autre chose que des données métabo) => plus de dev en perspective, et perso j'aurais plutôt tendance à dire que si quelqu'un veut faire ça sur des données non métabo, puisque ce n'est pas de la métabo il peut très bien mettre ses variables dans la sampleMetada et faire concorder sa datamatrix (après tout c'est son choix de venir utiliser des outils designés pour la métabo !)

Point 2 : Separators

A priori, vues nos conventions actuelles de format, il n'y a plus besoin des paramètres "Separator of columns" et "Decimal separator".

Point 3 : Outputs

Actuellement, il y a deux fichiers en sortie :

• un fichier correspondant à la dataMatrix avec en plus les colonnes de p-value
• un fichier correspondant à la dataMatrix avec seulement les variables significatives

Sauf que du coup, ces deux fichiers ne permettent pas de chaîner ensuite avec d'autres outils. Et du coup, voici deux propositions :

• Option minimaliste : au lieu de mettre les colonnes de pvalue dans la dataMatrix, les mettre dans la variableMetadata (nécessite de l'ajoute en input de l'outil), et supprimer la sortie filtrée. Il faut également ajouter une indicatrice supplémentaire en cas d'intéraction pour l'option Selection method (histoire que l'utilisateur puisse filtrer selon se critère avec Generic_Filter ensuite).
  ¤ Point fort : un seul dataset de sortie
  ¤ Point faible : il faut faire une étape de Generic Filter après, et ça fait un peu plus de dev à faire derrière (à noter qu'il faudrait alors aussi regrouper les paramètres Threshold et Selection Method sous la condition du N-way, donc un peu de changement dans l'agencement des paramètres)
• Option tout compris : il faut également mettre les colonnes de pvalue dans le variableMetadata (à ajouter en input donc), mais on ajoute également en plus de la datamatrix filtrée déjà présente un variableMetadata filtré (pour pouvoir chainer ensuite avec d'autres outils).
  ¤ Point fort : pas besoin d'utiliser Generic_filter derrière ni de rajouter la création d'une indicatrice supplémentaire dans variableMetadata, et moins de modif d'interface
  ¤ Point faible : plus de datasets en sortie

Point 4 : Working Example

Les sorties du Output files ne sont plus à jour.
Perso je recommande l'abandon de cette section, mais ce n'est que mon point de vue.