Construction d'une table de référence

Postez ici vos questions, réponses, commentaires ou suggestions - Les sujets seront ultérieurement répartis dans les archives par les modérateurs

Modérateur : Groupe des modérateurs

Remy Duchemin
Messages : 63
Enregistré le : 21 Mai 2011, 18:55

Construction d'une table de référence

Messagepar Remy Duchemin » 28 Jan 2014, 21:05

Bonsoir,

Mon objectif est de construire la table de la loi normale centrée réduite et d'intégrer le graphique au dessus de la table dans une sortie pdf, par exemple.
Pour ce faire, première étape construire le graph :

Code : Tout sélectionner

x=seq(-4,4,length=200)
y=dnorm(x)
plot(x,y,type="l", lwd=2, col="blue")
x=seq(-1,1,length=100)
y=dnorm(x)
polygon(c(-1,x,1),c(0,y,0),col="gray")

On peut bien entendu l'arranger mais ce n'est pas mon objet ici.
Deuxième étape, la construction de la table :
Pour avoir une aire particulière, voilà comment je procède : ici entre 0 et 1

Code : Tout sélectionner

 pnorm(1,mean=0,sd=1)-pnorm(0,mean=0,sd=1)

mais c'est pour construire le tableau final où là je bloque.
L'idée est d'obtenir un résultat de cette forme :
http://probastat.over-blog.com/article-acceder-a-la-table-de-la-loi-normale-55117187.html
Est-ce qu'il possible d'obtenir un peu d'aide s'il vous plaît ?
Merci par avance

Logez Maxime
Messages : 3138
Enregistré le : 26 Sep 2006, 11:35

Messagepar Logez Maxime » 29 Jan 2014, 08:11

Bonjour,

une possibilité :

Code : Tout sélectionner

fun <- function(x, y) pnorm(x+y)
tab <- outer(x<-seq(0,2.9, by=0.1), y<-seq(0,0.09, by=0.01), fun)
colnames(tab) <- y
rownames(tab) <- x
tab <- round(tab,4)
Cordialement,
Maxime

Remy Duchemin
Messages : 63
Enregistré le : 21 Mai 2011, 18:55

Messagepar Remy Duchemin » 29 Jan 2014, 09:02

Merci beaucoup Maxime!
Je ne connaissais pas "outer" :?

Remy Duchemin
Messages : 63
Enregistré le : 21 Mai 2011, 18:55

Messagepar Remy Duchemin » 29 Jan 2014, 09:17

Un petit complément, si je veux sortir dans un pdf le graph et le tableau que dois-je utiliser ?

Code : Tout sélectionner

pdf(file="sortie.pdf")
...
dev.off()


Retourner vers « Questions en cours »

Qui est en ligne

Utilisateurs parcourant ce forum : Aucun utilisateur enregistré et 1 invité