123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242 |
- #!/usr/bin/Rscript
- args <- commandArgs(TRUE)
- dist <- read.table(args[1])
- taille = args[2]
- taille = as.integer(taille)
- tailleTotale = length(dist$V1)
- #cat("taille demandée = ",taille,"\n");
- #cat("taille totale = ",tailleTotale,"\n");
- if (taille>tailleTotale){
- taille = tailleTotale;
- }
- png("fig.png", width = 400, height = 300)
- #pdf("fig.pdf")
- #x11()
- #options(devices="X11")
- #vdist = dist$V1[1:100]
- vdist = dist$V1[1:taille]
- #vdist = dist$V1
- #v = ks.test(vdist,pnorm,mean(vdist),sd(vdist))
- n = length(vdist)
- #skew = mean((vdist-mean(vdist))^3)/sd.pop(vdist)^3
- m3 <- sum((vdist - mean(vdist))^3)/length(vdist)
- s3 <- sqrt(var(vdist))^3
- skew = m3/s3
- cat("taille de l'échantillon = ",n,"\n\n")
- cat("moyenne = ",mean(vdist)," ")
- cat("mediane = ",median(vdist),"\n\n")
- cat("min = ",min(vdist)," ")
- cat("max = ",max(vdist),"\n\n")
- cv = 100*sd(vdist)/mean(vdist)
- if (mean(vdist) == 0) cv=0
- cat("écart-type",sd(vdist),"\n")
- cat("coef. de variation = ",cv)
- if ( cv < 5 ) {
- cat (" pas de bruit\n\n")
- } else if (cv<33) cat (" bruit faible\n\n") else cat (" bruit important\n\n")
-
- quant <- quantile(vdist)
- cat("quartiles ",quant,"\n");
- cat("écart semi-interquartiles ",(quant[4]-quant[2])/2,"\n");
- #skew = mean((vdist-mean(vdist))^3)/sd.pop(vdist)^3
- m3 <- sum((vdist - mean(vdist))^3)/length(vdist)
- s3 <- sqrt(var(vdist))^3
- skew = m3/s3
- cat("skewness = ",skew,"\n");
- hist(vdist,col = "blue",breaks = 50,freq = FALSE)
- den <- density(vdist)
- lines(den, col = "red")
- dev.off()
|