Eaux1<-read.table("Eaux1.txt",h=T, row.names=7) #loading of the data in the R software in an object named Eaux1. #This object is a dataframe and not a matrix boxplot(Eaux1) Dist.forme<-function(x) { par(mfrow=c(2,2)) x.naomit<-na.omit(x) hist(x.naomit, col="gray", prob=TRUE,xlim=c(min(x.naomit), max(x.naomit)), main="") curve(dnorm(x,mean=mean(x.naomit),sd=sd(x.naomit)),add=TRUE,lwd=2,col="red") boxplot(x.naomit) iqd<-summary(x.naomit)[5] - summary(x.naomit)[2] points(mean(x.naomit), col = "orange", pch = 18) plot(density(x.naomit,width=2*iqd),xlab="x",ylab="",type="l",main="") qqnorm(x.naomit) qqline(x.naomit) } Dist.forme(Eaux1[,1]) Dist.forme(Eaux1[,2]) Eaux1[,1] rank(Eaux1[,1]) sort(Eaux1[,1]) median(Eaux1[,1]) mean(Eaux1[,1]) mean(Eaux1[,1],trim=0.5) mean(Eaux1[,1],trim=0.05) mean(Eaux1[,2]) mean(Eaux1[,2],trim=0.5) mean(Eaux1[,2],trim=0.05) var(Eaux1[,1],na.rm=T) sd(Eaux1[,1]) var(Eaux1[,2]) sd(Eaux1[,2]) summary(Eaux1) library(e1071) for (i in 1:6 ) { print(sd(Eaux1[, i],na.rm=T)) print(skewness(Eaux1[, i],na.rm=T)) print(kurtosis(Eaux1[,i],na.rm=T)) } ceram18<-read.table("ceram18.txt",h=T, row.names=1,sep="\t") dim(ceram18) ; monnaie=t(ceram18[1,]) summary(na.omit(monnaie)) monnaie<-na.omit(monnaie)/100 par(mfrow=c(2,2)) hist(monnaie,nclass=4,xlab= "4 classes",prob=TRUE, col="lightblue", border="pink") hist(monnaie,nclass=7,xlab= "7 classes",prob=TRUE, col="lightblue", border="pink") hist(monnaie,nclass=10,xlab= "10 classes",prob=TRUE, col="lightblue", border="pink") plot(x=c(3,16),y=c(0,0.7),type="n",bty="l",xlab="Annee/100", ylab="Estimation de la densite ") rug(monnaie) lines(density(monnaie,bw=0.15),col = "blue") lines(density(monnaie,bw=0.25),col = "green") lines(density(monnaie,bw=0.5),col = "red") lines(density(monnaie,bw=1),col = "pink") panel.hist <- function(x, ...) { usr <- par("usr"); on.exit(par(usr)) par(usr = c(usr[1:2], 0, 1.5) ) h <- hist(x, plot = FALSE,col="lightblue") breaks <- h$breaks; nB <- length(breaks) y <- h$counts; y <- y/max(y) rect(breaks[-nB], 0, breaks[-1], y, col="cyan", ...) } panel.cor <- function(x, y, digits=2, prefix="", cex.cor, ...) { usr <- par("usr"); on.exit(par(usr)) par(usr = c(0, 1, 0, 1)) r <- abs(cor(x, y)) txt <- format(c(r, 0.123456789), digits=digits)[1] txt <- paste(prefix, txt, sep="") if(missing(cex.cor)) cex.cor <- 0.8/strwidth(txt) text(0.5, 0.5, txt, cex = cex.cor * r) } panel.lm<-function(x,y){ points(x,y) abline(lm(y~x)) #lines(lowess(x,y),col="red") } pairs(Eaux1, diag.panel=panel.hist, cex.labels = 2, font.labels=2, upper.panel=panel.cor,lower.panel=panel.lm )