Je souhaiterais calculer et regrouper dans un tableau les pvalues d’un test de Wilcoxon pour toutes les comparaisons 2 à 2 possibles.
Je pense qu’il faut utiliser pairwise.wilcox.test.
J’ai un dataframe dans lequel chaque ligne correspond à une cellule.
Je souhaite comparer les moyennes de fluo (Mean) pour chaque Cond1.
Mais je butte sur « Error in factor(g) : objet 'Cond1' introuvable ».
J’ai essayé de « convertir » Cond1 en ClassNb (Numeric, factor ou character) Même message d’erreur qu’avec Cond1 …
En fait, je ne comprends pas à quoi correspond g (grouping vector or factor)
Voici mon script :
Code : Tout sélectionner
df3 <- read.csv("df3.csv",header=T, stringsAsFactors = F)
pairwise.wilcox.test(Mean, Cond1 , p.adj = "bonf")
et mon dataframe :
Code : Tout sélectionner
> dput(df3)
structure(list(Cond1 = c("OP_L2_J04", "OP_L2_J04", "OP_L2_J04",
"OP_L2_J04", "OP_L2_J04", "OP_L2_J04", "OP_L2_J04", "OP_L2_J04",
"OP_L2_J04", "OP_L2_J04", "OP_L2_J04", "OP_L2_J04", "OP_L2_J04",
"OP_L2_J04", "OP_L2_J04", "OP_L2_J04", "OP_L2_J04", "OP_L2_J04",
"OP_L2_J04", "OP_L2_J04", "OP_L2_J04", "OP_L2_J04", "OP_L2_J04",
"OP_L2_J04", "OP_L2_J04", "OP_L2_J04", "OP_L2_J04", "OP_L2_J04",
"OP_L2_J04", "OP_L2_J04", "OP_L2_J01", "OP_L2_J01", "OP_L2_J01",
"OP_L2_J01", "OP_L2_J01", "OP_L2_J01", "OP_L2_J01", "OP_L2_J01",
"OP_L2_J01", "OP_L2_J01", "OP_L2_J01", "OP_L2_J01", "OP_L2_J01",
"OP_L2_J01", "OP_L2_J01", "OP_L2_J01", "OP_L2_J01", "OP_L2_J01",
"OP_L2_J01", "OP_L2_J01", "OP_L2_J01", "OP_L2_J01", "OP_L2_J01",
"OP_L2_J01", "OP_L2_J01", "OP_L2_J01", "OP_L2_J01", "OP_L2_J01",
"OP_L2_J01", "OP_L2_J01", "OP_L2_J14", "OP_L2_J14", "OP_L2_J14",
"OP_L2_J14", "OP_L2_J14", "OP_L2_J14", "OP_L2_J14", "OP_L2_J14",
"OP_L2_J14", "OP_L2_J14", "OP_L2_J14", "OP_L2_J14", "OP_L2_J14",
"OP_L2_J14", "OP_L2_J14", "OP_L2_J14", "OP_L2_J14", "OP_L2_J14",
"OP_L2_J14", "OP_L2_J14", "OP_L2_J14", "OP_L2_J14", "OP_L2_J14",
"OP_L2_J14", "OP_L2_J14", "OP_L2_J14", "OP_L2_J14", "OP_L2_J14",
"OP_L2_J14", "OP_L2_J14"), ClassNb = c(19L, 19L, 19L, 19L, 19L,
19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L,
19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L,
19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L,
19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L, 19L,
19L, 19L, 19L, 21L, 21L, 21L, 21L, 21L, 21L, 21L, 21L, 21L, 21L,
21L, 21L, 21L, 21L, 21L, 21L, 21L, 21L, 21L, 21L, 21L, 21L, 21L,
21L, 21L, 21L, 21L, 21L, 21L, 21L), Mean = c(3826.67, 3344.29,
3407.1, 3534.19, 3299.52, 3426.14, 3270, 3587, 3246, 4820.25,
4821.56, 5175.2, 4814, 4799.75, 4749.67, 4990.92, 4743, 4744.67,
5004.33, 6785.59, 5803.29, 6358.91, 34211.12, 27182, 26775.33,
32392.54, 1803, 2129.9, 2032.4, 1808, 2005.91, 1918.75, 1941.06,
1835.5, 3129, 3220.5, 3175, 3162, 4833, 4830, 4911.47, 5628.67,
4364.5, 4523.36, 4739.65, 4346.76, 4145, 6064.6, 4730, 4673,
5369.44, 5216.62, 6299.1, 6329.84, 5829.26, 5882, 5854, 6216.22,
7468.46, 4863.86, 13039.77, 11975.23, 10076.5, 12738.8, 16307.05,
14469.84, 11440.64, 16341.18, 8936, 8522, 9361.63, 9300.52, 8457,
8493.83, 8421.91, 8694, 8920.93, 9293.31, 8662.5, 9379.71, 12223.58,
8600, 12682.33, 9269.21, 14209.79, 8794.5, 12094.5, 13122.72,
13350.33, 9342)), class = "data.frame", row.names = c(NA, -90L
))
>
Question subsidiaire, je ne sais pas comment écrire que je souhaite enregistrer en csv le tableau de pvalues généré.
Je vous remercie par avance de votre aide et conseils.
Patricia