1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37
|
# creation d'un jeu de donnees pour tester le code
x<-100*rnorm(20)
y<-sample(c("canne","mais","blé"),20,replace=TRUE)
z<-paste("Etude",1:20,sep="")
M<-data.frame(CO2=x,Culture=y,Etude=z) # si les donnees ne sont pas dans un data.frame, soit les passer en data.frame soit adapter un peu le code
# code utilise
tab<-table(M$Culture)
Cult<-names(tab)
lvs<-levels(M$Etude)
MinMax<-rep(0,6*2)
dim(MinMax)<-c(2,6)
colnames(MinMax)<-paste(c("Min","Max"),sort(rep(Cult,2)),sep="-")
k<-0
for(i in 1:length(tab))
{
S<-M[which(M$Culture==Cult[i]),]
MinMax[2,i+k]<-lvs[S$Etude[which(S$CO2==min(S[,1]))]]
MinMax[2,i+k+1]<-lvs[S$Etude[which(S$CO2==max(S[,1]))]]
MinMax[1,i+k]<-min(S[,1])
MinMax[1,i+k+1]<-max(S[,1])
k<-k+1
} |
Partager