Conclusions
L'été a été sec et chaud. Les
tensiomètres ont très vite atteint le maximum de mesure.
Néanmoins les seules différences observées ont
été :
- le couvert Ecosem dessèche plus vite que le couvert
Fétuque
- le rang désherbé chimiquement dessèche
plus vite et se réhumecte moins vite que le rang sous bâche
tissée.
Sinon la fauche n'a pas d'influence sur la
disponibilité azotée ou hydrique. Le type de couvert (Ecosem
contre Fétuque) n'a pas d'influence sur la disponibilité
azotée.
Qu'en est-il de la concurrence racinaire ? La
compétition se ferait-elle plutôt au niveau de l'espace racinaire,
plutôt que de la disponibilité directe en nutriment ? Des profils
racinaires sont prévus pour l'hiver 2015, afin d'observer l'interface
racinaire pommier/interrang.
Si concurrence accrue il y a, résulte-t-elle d'une
différence d'espèces (biomasse produite, profondeur racinaire,
période de croissance préférentielle) ou du mode de
végétation (vivace contre annuelle) ?
13
MORGANE FOURNIER
ETUDE DE LA COMPETITION INTERRANG/POMMIER ET
RANG/POMMIER
Scripts
tensio_test <- function(data) j
#Lorsque toutes les valeurs sont à 199, ne fait pas de
test
if (length(unique(data$Mesure))!=1)
j
# On fait un test de normalité et des tests
d'homogénéité des variances
norm=shapiro.test(data(,'Mesure'])
homo=bartlett.test((data(,'Mesure'])~(data(,'Couvert']))
a=''
testrealise=''
if (norm(2]>0.05 & homo(3]>0.05)j
#MODIFIER ICI LE TEST SOUHAITE : LIGNE ou RANG
lm.tensio=lm(Mesure~Couvert+Bloc, data=data)
anov=Anova(lm.tensio)
#on récupère la plue value du test de l'effet
modalité pvalue=anov((4]](1] pvalue2=anov((4]](2] testrealise="Test
realise : ANOVA"
if (pvalue<=0.05) j
a="montre une difference significative entre couvert"
b=1
} else j
a="ne montre pas de difference significative"
b=0
}
} else j # Si le test de normalité n'est
pas ok, on fait un test de
Kruskal Wallis
kw=kruskal.test(list(data(data$Couvert=='ecosem','Mesure'],data(data$Couver
t=='fetuque et rg','Mesure']))
pvalue2='pas de test'
# Si le facteur a un effet
if(kw(3]<=0.05)j
testrealise="Test réalisé : Kruskal-Wallis"
a="montre une difference significative entre couvert"
b=1
} else j
testrealise="Test réalisé : Kruskal-Wallis"
a="ne montre pas de difference significative"
b=0
}
}
# Quand valeurs à 199 } else j
testrealise='aucun' a=''
b=''
pvalue2=''
}
meandata=tapply(data(,'Mesure'],data(,'Couvert'],mean,na.rm=T)
print(list(paste(testrealise, a, " avec des moyennes de (en cbar) :", sep=' '),
meandata),quote=F) sortie=data.frame('Date'=unique(data$Date), 'Test'=b,
'ECOSEM'=meandata('ecosem'], 'FETUQUE et RG'=meandata('fetuque et rg'], 'effet
bloc X'=pvalue2) return(sortie)
}
1
MORGANE FOURNIER
| 
