Bonjour,
Je cherche à tracer un qq-plot avec ma distribution empirique et une loi de pareto.
qq un serait comment faire svp ?
merci
Modérateur : Groupe des modérateurs
Code : Tout sélectionner
?qqplot
Code : Tout sélectionner
require(actuar)
# ou si tu ne l'as pas : install.packages("actuar")
? rpareto
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])
Code : Tout sélectionner
Erreur dans qgpd(ppoints(data), xi = xi) :
argument(s) inutilisé(s) (xi = xi)
MARI WAFA a écrit : Je cherche à tracer un qq-plot avec ma distribution empirique et une 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", ...)
}
)
Retourner vers « Questions en cours »
Utilisateurs parcourant ce forum : Aucun utilisateur enregistré et 1 invité