ellipse de confiance + doc sur treillis

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

Michel Arnaud
Messages : 2
Enregistré le : 04 Oct 2005, 06:49

ellipse de confiance + doc sur treillis

Messagepar Michel Arnaud » 16 Mar 2007, 07:19

Je cherche
1) a tracer l'ellipse de confiance a p% d'un ensemble de n points de coordonnees (xi, yi)
Avez vous des noms de fctions qui pourraient faire cela

2) un doc pedago sur le package lattice

Merci de vos reponses

Renaud Lancelot
Messages : 2484
Enregistré le : 16 Déc 2004, 08:01
Contact :

Messagepar Renaud Lancelot » 16 Mar 2007, 12:35

Bonjour Michel,

1. Voir data.ellipse dans le package car:

Code : Tout sélectionner

library(car)
set.seed(123)
x <- 1:100
y <- 1.5 * x + rnorm(100, 0, 20) + 50
plot(x, y, xlim = c(-50, 150), ylim = c(0, 250))
data.ellipse(x, y, levels = .95, add = TRUE)


2. Le meilleur doc que je connaisse est le bouquin de Paul Murrell:

Murrell, P., 2005. R graphics. Chapman & Hall / CRC, London, 301 p.

Il dépasse largement le cadre lattice mais reste très pédagogique avec bcp d'exemples. C'est aussi le support de la formation sur la construction de graphes avec R que j'ai suivie l'année dernière avec V Chevalier et E Tillard.

Si toi (et d'autres !) le souhaitent, je peux essayer de faire une fiche avec le contenu lattice de cette formation, mais je n'aurai pas le temps de le faire avant cet été.

Amicalement,

Renaud

Stéphane Laurent
Messages : 1557
Enregistré le : 05 Déc 2006, 19:07

Messagepar Stéphane Laurent » 25 Jan 2009, 15:47

Bonjour,

Il se passe quelque chose d'étrange. Le code de Renaud ci-dessus fonctionne chez moi, mais voilà ce qui se passe si j'utilise autre chose pour x et y :

Code : Tout sélectionner

> ....
> x <- SIMS$alpha[1:100]; y<-SIMS$beta[1:100]
> plot(x,y, type="n")
> data.ellipse(x,y, add=TRUE, levels=.95)
Error in data.ellipse(x, y, add = TRUE, levels = 0.95) :
  x and y must be vectors of the same length
> length(x); length(y)
[1] 100
[1] 100


Qu'est-ce que ça peut être ? x et y sont bien des vecteurs numériques de même longueur.

Lorraine Bauer
Messages : 67
Enregistré le : 18 Nov 2008, 12:06

Messagepar Lorraine Bauer » 25 Jan 2009, 16:33

Bonjour,

Si toi (et d'autres !) le souhaitent, je peux essayer de faire une fiche avec le contenu lattice de cette formation, mais je n'aurai pas le temps de le faire avant cet été.


Je pense que cela pourrait intéresser pas mal de monde ici et évidemment moi-même je suis très friande de ce genre de doc. Je dis cela d'autant plus que récemment, j'ai eu quelque démêlés avec ce paquetage... Merci à vous Renaud.

Lorraine

Pierre Casadebaig
Messages : 67
Enregistré le : 07 Déc 2006, 12:31

Messagepar Pierre Casadebaig » 27 Jan 2009, 09:18

Bonjour,

l'auteur du paquet "lattice" a aussi écrit un livre sur son logiciel.
Ce livre est plein d'exemples, est assez direct aussi. J'estime qu'il ressemble à une grosse doc illustrée et bien plus agréable que le manuel en ligne (?xyplot). Il y a même moyen d'y avoir accès en ligne selon la structure qui vous héberge !

Renaud Lancelot
Messages : 2484
Enregistré le : 16 Déc 2004, 08:01
Contact :

Messagepar Renaud Lancelot » 27 Jan 2009, 14:54

Oui, et mon mail date de 2007: pas mal d'eau a coulé sous les ponts, et pas vraiment le temps de me lancer dans la rédaction d'une fiche... Désolé !

Renaud


Retourner vers « Questions en cours »

Qui est en ligne

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