Génération d'une distribution gamma généralisée à l'aide de scripts JMP<meta http-equiv="Content-Type" content="text/html; charset=UTF-8" />Tout,
Existe-t-il une fonction qui permettrait de générer des nombres aléatoires
d'une distribution gamma généralisée ? Le guide de script a
fonction gamma aléatoire pour générer des nombres aléatoires à partir
Distribution gamma mais je ne vois rien pour le gamma généralisé
Distribution.
Merci de votre aide.
-ASVP
Génération d'une distribution gamma généralisée à l'aide de scripts JMPMerci Matt et MS pour avoir montré un moyen de générer des nombres aléatoires à partir de
distribution gamma généralisée.
Salutations,
-ASVP
Génération d'une distribution gamma généralisée à l'aide de scripts JMPSalut:
MS nous a donné un indice.À partir de l'invite R :
> bibliothèque (VGAM)
>rgengamma
et
>?rgengamma
révèle que rgengamma renvoie une fonction rgamma de base R modifiée
y = rgamma(n, k)
échelle*y^(1/j)
Dans JMP, nous pouvons également utiliser une modification de la fonction de base "Random Gamma", par exemple
Table de données actuelle() << Nouvelle colonne( "JMP GenGamma",
Formule( 1*Gamma aléatoire( 1 )^(1/2) ) );
Table de données actuelle() << Exécuter les formules ;
(Bien sûr, il faudra bidouiller les paramètres pour obtenir une correspondance)
Faire correspondre les générateurs de nombres aléatoires, n'importe qui - peut-être en fixant des graines ?
Meilleures salutations,
-Mat
Meilleur
Génération d'une distribution gamma généralisée à l'aide de scripts JMPJe ne pense pas que vous puissiez obtenir des nombres aléatoires pour une distribution gamma généralisée à trois paramètres en utilisant les générateurs intégrés. Si c'est possible, j'aimerais aussi savoir comment faire.
Cependant, JMP 9 peut obtenir des nombres aléatoires à partir de R. Voir l'exemple JSL ci-dessous.
// Fonction qui renvoie des nombres aléatoires à partir de
// distribution gamma généralisée générée
// par package R VGAM
rgengamma = Fonction ( { n = 100 , b = 1 , ré = 1 , k = 1 } ,
R Init () ;
Soumettre ( "
bibliothèque (VGAM)
rggam <- rgengamma(100,1,2,3)
" ) ;
jeu de rôle = R Obtenir ( jeu de rôle ) ;
Terme R () ;
jeu de rôle
) ;
// Exemple d'utilisation de la fonction pour remplir un
// table de données avec n=100 nombres aléatoires
dt = Nouveau tableau ( "Gamma généralisé" , Nouvelle colonne ( "Aléatoire" ) ) ;
Colonne ( dt , "Aléatoire" ) << valeurs définies ( rgengamma ( 100 , 1 , 2 , 3 ) ) ;