problème avec la "tapply" ==aide SVP

Bonjour je veux calculer la moy avec tapply,
mais voici ce que j'obtiens comme erreur:

Code:

Error in tapply(ph2$ph_2, repi, mean) : arguments must have same length

voici mon script:

Code:

ph2<- read.table (file= file.choose(), header=T )
repi<- rep((1:42), each=3)
m<- tapply(ph2$ph_2, repi,  mean)

Et je vous mets le jeu de données :
je veux calculer la moy pour toutes les valeurs de Ph_2 (n=42) dans chaque parcelle par site (1moy = 1 parcelle)

site parcelle pts ph_2
cast 9 1 4.82 ....|\
cast 9 2 4.84 ....|-==»1moy
cast 9 3 4.76 ....|/
cast 10 1 4.82
cast 10 2 4.51
cast 10 3 4.46
Cléri 9 1 4.91
cléri 9 2 4.63
cléri 9 3 4.72
cléri 10 1 4.86
cléri 10 2 4.69
cléri 10 3 4.53
dollard 9 1 4.83
dollard 9 2 5.03
dollard 9 3 4.77
dollard 10 1 5.00
dollard 10 2 5.02
dollard 10 3 5.00
dup_est 9 1 5.06
dup_est 9 2 5.18
dup_est 9 3 4.95
dup_est 10 1 5.03
dup_est 10 2 5.01
dup_est 10 3 4.86
dup_oue 9 1 4.68
dup_oue 9 2 4.66
dup_oue 9 3 4.38
dup_oue 10 1 4.80
dup_oue 10 2 4.36
dup_oue 10 3 4.15
land 9 1 5.04
land 9 2 4.75
land 9 3 4.91
land 10 1 4.31
land 10 2 4.63
land 10 3 4.65
Paradis 9 1 4.49
Paradis 9 2 4.29
Paradis 9 3 4.10
Paradis 10 1 4.22
Paradis 10 2 4.43
Paradis 10 3 4.20

bonjour,

il est normal que ça plante : tu n'as pas 42 combinaisons différentes mais 14 !

Code:

repi <- rep(1:14, each=3)

3*14= 42 alors que 3*42 = 126.

Sinon le mieux est encore de faire :

Code:

repi <- with(ph2, paste(site, parcelle, sep= "_ "))

après avoir corrigé Cléri en cléri.

c'est extraordinaire
merci
maintenant j'ai bien compris le raisonnement,
et pour ta formule c'est très pratique, parce que je savais pas comment exporter chaque moy avec parcelle-site (dans une 1 colonne)
merci bcp

droopy a écrit:bonjour,

il est normal que ça plante : tu n'as pas 42 combinaisons différentes mais 14 !
repi <- rep(1:14, each=3). 3*14= 42 alors que 3*42 = 126.

Sinon le mieux est encore de faire :
repi <- with(ph2, paste(site, parcelle, sep= "_ "))
après avoir corrigé Cléri en cléri.