Outils pour utilisateurs

Outils du site


public:specs_distribution_specificites

Différences

Ci-dessous, les différences entre deux révisions de la page.

Lien vers cette vue comparative

Les deux révisions précédentes Révision précédente
public:specs_distribution_specificites [2013/10/17 17:20]
matthieu.decorde@ens-lyon.fr
public:specs_distribution_specificites [2014/03/27 16:05] (Version actuelle)
matthieu.decorde@ens-lyon.fr
Ligne 18: Ligne 18:
  
 <​code>​ <​code>​
-x=14 +library(textometry) 
-F=200 +`specif_distrib` <​-function(xF, t, T) { 
-T=10000 +
-t=500+
 xlimit=(F/​4) xlimit=(F/​4)
 f=0:xlimit f=0:xlimit
Ligne 39: Ligne 38:
 pfsumx=pfsum[x+1] pfsumx=pfsum[x+1]
  
-plot(f,​pfsum,​ type="​l",​ col="​blue"​)+title <- sprintf("​phyper(%s,​ F=%d, T-F=%d, t=%d) avec f=%d", paste("​[",​f[1],​ ",",​ f[xlimit+1],"​]"​),​ F, T-F, t, x) 
 +plot(f,​pfsum,​ type="​l",​ col="​blue"​, main=title)
 lines(f, pf, type="​l",​ col="​green"​) lines(f, pf, type="​l",​ col="​green"​)
  
 al=max(pfsum)/​10 al=max(pfsum)/​10
 ad=1 ad=1
-if (m > x) +if (m > x) {ad = -ad
- ad = -ad+tpos=4 
 +if (m > x) {tpos = 2}
 arrows(m, al*0.7,m,0, length=0.1) arrows(m, al*0.7,m,0, length=0.1)
-arrows(x+ad,​ al*0.7+px,​x,​px,​ length=0.1) +arrows(x+ad,​ al*0.7+px, x, px, length=0.1) 
-arrows(x+ad,​ al*0.7+pfsumx,​x,​pfsumx,​ length=0.1)+arrows(x+ad,​ al*1.4+pfsumx, x, pfsumx, length=0.1)
  
-text(m,​al,​paste("​mode=",​round(m))) +smode <- sprintf("mode = %d", round(m)) 
-text(x,​al+px,​paste("p(",​x,"​) = ", format.default(px,​ sci = TRUE, digits = 4))) +spf <- sprintf("p(f'=%d) = %s", x, format.default(px,​ sci = TRUE, digits = 4)) 
-text(x,​al+pfsumx,​paste("p(n>​=",​x,"​) = ", format.default(pfsumx,​ sci = TRUE, digits = 4)))+if (m > x) { 
 + sintpf <- sprintf("p(f'>=%d) = %s", x, format.default(pfsumx,​ sci = TRUE, digits = 4)) 
 +} else { 
 + sintpf <- sprintf("p(f'<​=%d) = %s", x, format.default(pfsumx,​ sci = TRUE, digits = 4)) 
 +}
  
-s1 <- sprintf("​phyper(f=%sF=%dT-F=%dt=%d)"paste("​["​,f[1]","f[xlimit+1],"​]"​),F,T-F,t+text(m,al,smode,xpd=TRUE) 
-s2 <- sprintf("​integral ​%s", s1+text(x, (al*0.7)+px,spfpos=tpos,xpd=TRUE) 
-legend("​topright",​ c(s1, s2),lty = 1, title = "​Légende",​ col=c("​green",​ "​blue"​),​+text(x(al*1.4)+pfsumx,sintpfpos=tpos,xpd=TRUE) 
 + 
 +s1 <- sprintf("​P(f'​ = f)"
 +s2 <- sprintf("​integral ​de P(f' <= f) et P(f' >= f)") 
 + 
 +legend("​topright",​ c(s1, s2, spf, sintpf),lty = 1, title = "​Légende",​ col=c("​green",​ "blue", "​white",​ "white"),
        inset = .02, cex=0.9)        inset = .02, cex=0.9)
 +}
 +specif_distrib(x=15,​F=200,​t=500,​T=10000)
 </​code>​ </​code>​
  
public/specs_distribution_specificites.txt · Dernière modification: 2014/03/27 16:05 par matthieu.decorde@ens-lyon.fr