2013-05-31 1 views
0

J'utilise R, et je voudrais générer un certain nombre d'observations en utilisant rweibull(n, shape, scale = 1).Besoin d'un échantillon aléatoire: Distribution de Weibull avec des temps d'inter-arrivée spécifiques

J'ai le taux d'arrivée (c'est-à-dire 1/temps interarrivée), mais je ne sais pas comment l'utiliser dans la fonction rweibull.

+0

Ce n'est pas vraiment une question sur R, ce est une question de statistiques sur la façon d'utiliser la distribution de Weibull. Vous devriez essayer de demander cela à http://stats.stackexchange.com/ –

+1

Vous devez déterminer quel taux d'arrivée correspond à; c'est probablement l'échelle (ou 1/échelle). – huon

Répondre

2

Le paramètre d'échelle est ce que vous devez travailler avec et le paramètre de forme est ce qui doit être défini sur 1 pour créer une distribution exponentielle. Le paramètre d'échelle est de 1/taux:

interT = 8 
plot(density(rexp(100, rate=1/interT))) 
with(density(rweibull(100, scale=interT, shape=1)), 
    lines(x,y, col="red")) 

(. Mais si vous utilisez le package de survie vous devez être conscient que les paramètres sont différents)

+0

Merci les gars. Ce que j'essaie de faire est de générer des nombres aléatoires comme une heure d'arrivée du trafic Web, ce qui suit Wiebull. Je connais le taux d'arrivée pour le trafic web/s. Dwin pour l'échelle = 1/taux. voulez-vous dire par le taux est l'heure d'arrivée ou le taux d'arrivée? J'ai seulement le taux d'arrivée, donc je peux dire échelle = 1/taux d'arrivée? Donc, la fonction fonctionne comme ceci: rweibull (n, shape = 1, scale = 1/taux d'arrivée)? – user2440912

+0

Ce que je voulais dire, c'est que le paramètre 'rweibull' 'scale' était équivalent à un paramètre' rexp' 1/'rate', dans le cas où le rweibull: shape = 1. Mon exemple utilisait un "temps d'arrivée" de 8 (dans des unités non spécifiées) et visait à démontrer (tirages aléatoires modulo) cette équivalence. –

Questions connexes