Commandes et procédures relatives aux lois Uniformes continues.
dunif(\(t\),min=,max=,log=),
punif(\(t\),min=,max=,lower.tail=,log.p=),
qunif(\(q\),min=,max=,lower.tail=, log.p=),
runif(\(n\),min=,max=) :
Soit \(X\) une v.a. de loi \({\cal U}(\rbrack\theta_1\ ;\ \theta_2\lbrack)\). Dans toutes les commandes précédentes il faut indiquer
min=\(\theta_1\) (0 par défaut) et max=\(\theta_2\) (1 par défaut). Ces commandes donnent respectivement \(f_X(t)\) ou son logarithme si
log=TRUE, \(P(X \leq t)\) si lower.tail=TRUE ou son logarithme si log.p=TRUE, le quantile d’ordre \(q\) si
lower.tail=TRUE ou son logarithme si log.p=TRUE, une simulation d’un \(n\)-échantillon de \(X\).
Les commandes suivantes donnent le Graphique de la f.r. de la loi
Uniforme \({\cal U}(\rbrack 1,5\ ;\ 3,5\lbrack)\).
plot(
function(x)
punif(x,min=1.5,max=3.5), 0, 5,
xlab="x",
ylab="y",
col="blue",
main=" Fig. 1. RĂ©partition U(]1,5 ; 3,5[)",
ylim=
c(-.25,1.25))
Les commandes suivantes donnent le Graphique des fonctions puissances des deux tests de
l’Exemple concernant l’alternative d’hypothèses unilatérales pour des lois Uniformes continues. Nous compilons («sourçons») dans R
les procédures Pupsi1 et Pupsi2,
créées pour calculer ces fonctions puissances.
plot(
function(x)
Pupsi1(x,theta0=1,alpha=.1,n=5,resul),.5,1.6,
xlab="theta",
ylab="pu",
ylim=
c(-0.02,1.02),
main="Fig. 1. Fonctions puissances.",
col="sienna4");
curve(
Pupsi2(x,1,.1,5,resul),
add=TRUE,
col="palegreen4");
legend(
x="topleft",
y=NULL,
legend=
c("Pupsi1","Pupsi2"),
text.col=
c("sienna4","palegreen4"));
x0=c(1,0);
y0=c(0,0.1);
x1=c(1,1);
y1=c(0.1,0.1);
segments(x0,y0,x1,y1,
col="red");
points(x=1,y=0.1,
col="red",
pch=16);