je cherche à afficher sur un graphique représentant sur une grille l'expression de différent gene en x et 2 calcul de cette expression en y.
Pour chaque sous graphique je cherche à exprimer un calcul anova un facteur sans test posthoc dans un premier temps. Voici un premier ecm qui permet d'avoir le graphique de base je cherche à comparer les valeurs de chaque carré indépendamment.
Code : Tout sélectionner
structure(list(Traitement = c("controle", "controle", "controle", "controle",
"controle", "controle", "controle", "controle",
"A-D1 ", "A-D1 ", "A-D1 ", "A-D1 ", "A-D1 ",
"A-D1 ", "A-D1 ","A-D1 ", "A-D2 ", "A-D2 ",
"A-D2 ", "A-D2 ", "A-D2 ", "A-D2 ", "A-D2 ",
"A-D2 ", "tt", "tt", " tt", " tt", " tt", " tt",
" tt", " tt", "A-D1 tt", "A-D1 tt", "A-D1 tt",
"A-D1 tt", "A-D1 tt", "A-D1 tt", "A-D1 tt",
"A-D1 tt", "A-D2 tt", "A-D2 tt", "A-D2 tt",
"A-D2 tt", "A-D2 tt", "A-D2 tt", "A-D2 tt", "A-D2 tt"),
individus = c(1, 1, 1, 1, 2, 2, 2, 2, 1, 1, 1, 1, 2, 2, 2, 2, 1, 1, 1, 1, 2, 2, 2, 2, 1,
1, 1, 1, 2, 2, 2, 2, 1, 1, 1, 1, 2, 2, 2, 2, 1, 1, 1, 1, 2, 2, 2, 2),
calcul = c("p2", "p2", "rq", "rq", "p2", "p2", "rq", "rq", "p2", "p2", "rq", "rq", "p2",
"p2", "rq", "rq", "p2", "p2", "rq", "rq", "p2", "p2", "rq", "rq", "p2", "p2",
"rq", "rq", "p2", "p2", "rq", "rq", "p2", "p2", "rq", "rq", "p2", "p2", "rq",
"rq", "p2", "p2", "rq", "rq", "p2", "p2", "rq", "rq"),
gene = c("Ab1", "AE", "Ab1", "AE", "Ab1", "AE", "Ab1", "AE", "Ab1", "AE", "Ab1", "AE",
"Ab1", "AE", "Ab1", "AE", "Ab1", "AE", "Ab1", "AE",
"Ab1", "AE", "Ab1", "AE", "Ab1", "AE", "Ab1", "AE",
"Ab1", "AE", "Ab1", "AE", "Ab1", "AE", "Ab1", "AE",
"Ab1", "AE", "Ab1", "AE", "Ab1", "AE", "Ab1", "AE",
"Ab1", "AE", "Ab1", "AE"),
exemple = c(1.19155298160405, 0.0724818202710361, 2.02156962843567, 0.384370027854773, 0.54991364994453,
0.00103756547537115, 0.579273470531843, 1.94371801242812, 1.70745867059473,1.6153066453706,
3.30113402776652, 3.25558460945979, 0.597367907304677,
0.76807566193499, 0.977831588376713, 0.811469434599755, 3.33840383376832,
0.0223708709402352, 0.425528356825245, 0.62844516251807, 1.95106168399692,
2.0153249055041, 1.18796291908844, 0.0530971721093183, 0.912160997006902,
1.23409393372274, 0.497099437885367, 0.552044980045594, 1.13474370257865,
1.05750011086696, 2.22045583334724, 0.215237661268395, 1.66412954997914,
2.19434943809382, 2.53727405265649, 2.35800460808463, 1.89055181006916,
0.0243139030902235, 0.595713335097338, 0.536213734244759, 1.91057989778762,
1.31894320754891, 1.83627997732475, 0.604549881851668, 1.25781122754147,
1.50537479363401, 3.02598020999327, 0.065350196112971)),
row.names = c(NA,-48L),
class = c("tbl_df", "tbl", "data.frame")) %>%
ggplot() +
aes(x = Traitement, y = exemple, col = calcul) +
facet_grid(calcul~gene) +
stat_summary(fun.data = mean_se, fun.args = list(mult = 1), geom = "pointrange", position = position_dodge(0.5)) +
geom_point(position=position_dodge(0.5), alpha =0.3) +
stat_compare_means(mapping = aes(x = Traitement, y = exemple),
method = "anova", label.y = 3)+
theme(axis.text.x = element_text(angle = 45))
Par ailleurs pensez vous qu'il soit possible de mettre la valeur Controle en début de graphique (cela tourne dans un script et cette valeur est stockée dans une variable et j'ai un peu de mal à jouer avec facteur pour se faire...