QQ-PLOT

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

MARI WAFA
Messages : 18
Enregistré le : 06 Sep 2011, 20:50

QQ-PLOT

Messagepar MARI WAFA » 15 Sep 2011, 20:31

Bonjour,

Je cherche à tracer un qq-plot avec ma distribution empirique et une loi de pareto.

qq un serait comment faire svp ?

merci

dicko ahmadou
Messages : 444
Enregistré le : 21 Nov 2009, 20:15

Messagepar dicko ahmadou » 15 Sep 2011, 21:36

Salut

Regarde l'aide de la fonction qqplot.
Quand on débute dans R, une des premières choses c'est de se familiariser avec le système d'aide.

Code : Tout sélectionner

?qqplot


Maintenant si tu lis bien l'aide et les exemples, tu comprendras que tu as surement besoin de générer un échantillon qui suit approximativement une distribution de Paréto (pour les paramètres j'imagine que tu sais quoi y mettre). Dans le package stats cette distribution n'est pas disponible, par contre elle existe dans d'autres packages.
Je te conseille le package actuar

Code : Tout sélectionner

require(actuar)
# ou si tu ne l'as pas : install.packages("actuar")
? rpareto
The best thing about being a statistician is that you get to play in everyone's backyard.
John Tukey

MARI WAFA
Messages : 18
Enregistré le : 06 Sep 2011, 20:50

Messagepar MARI WAFA » 17 Sep 2011, 16:29

Bonjour,

entre temps, j'ai trouvé ça :

Code : Tout sélectionner

gpd.model <- gpd(indice1, 0.0125) #0.0125 le seuil détermiré à partir d'un ME-PLOT
gpd.model$par.ests[1]
qplot(indice1, gpd.model$par.ests[1])



Mais j'ai l'erreur suivante sur le qplot :

Code : Tout sélectionner

Erreur dans qgpd(ppoints(data), xi = xi) :
  argument(s) inutilisé(s) (xi = xi)


pourquoi l'argument n'est pas utilisé ??[/quote]

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

Messagepar Renaud Lancelot » 18 Sep 2011, 06:16

qqplot, pas qplot. Regarder aussi la fonction qqmath du package lattice.
Renaud

MARI WAFA
Messages : 18
Enregistré le : 06 Sep 2011, 20:50

Messagepar MARI WAFA » 18 Sep 2011, 11:21

il y a bien une fonction qplot dans le package evir (Exploratory QQplot for Extreme Value Analysis).

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

Messagepar Renaud Lancelot » 19 Sep 2011, 05:34

Il y en a même d'autres, mais le conseil portait sur qqplot disponible dans les fonctions de base de R.
Renaud

Matthieu Lesnoff
Messages : 118
Enregistré le : 29 Nov 2004, 12:41

Re: QQ-PLOT

Messagepar Matthieu Lesnoff » 20 Sep 2011, 22:54

MARI WAFA a écrit : Je cherche à tracer un qq-plot avec ma distribution empirique et une loi de pareto

Un exemple de code avec qqmath de lattice et une fonction construite qchisq01 (mélange de lois de Chi2). Pour votre cas, il suffit de remplacer qchisq01 par la fonction définissant les quantiles de la loi de Pareto :

Code : Tout sélectionner

rchisq01 <- function(n, theta = 0.5) c(rep(0, round(theta * n)), rchisq(n - round(theta * n), df = 1))
qchisq01 <- function(p, theta = 0.5) {
  z <- rchisq01(100000, theta)
  z <- quantile(z, probs = p)
  as.numeric(z)
  }

obs <- rchisq(n = 5, df = 1) # une distribution empirique quelconque
qqmath(
  x = z,
  distribution = function(p) qchisq01(p), # la loi utilisée pour la comparaison
  panel = function(x, ...) {
    panel.qqmath(x, ...)
    panel.qqmathline(x, lty = 2, col = "red", ...)
    }
  )
ML


Retourner vers « Questions en cours »

Qui est en ligne

Utilisateurs parcourant ce forum : Bing [Bot] et 1 invité