Je me demande comment puis-je simuler un échantillon de n durées de vie de distribution Weibull qui incluent des observations de type I censurées à droite. Par exemple, ayons le n = 3, la forme = 3, l'échelle = 1 et le taux de censure = .15, et le temps de censure = .88. Je sais comment générer un...
Je cherche à générer 450 points de données en R. Il y a trois ensembles distincts 150 chacun répartis dans une bande circulaire avec des rayons différents (à 1, 2,8 et 5). En particulier, je cherche à reproduire le premier graphique à la p546 de The Elements of Statistical Learning. Je serais très...
Je sais comment générer une séquence avec une moyenne de . Par exemple, dans Matlab, si je veux générer une séquence de longueur , c'est:0 ± 1 10000±1±1\pm 1000±1±1\pm 1100001000010000 2*(rand(1, 10000, 1)<=.5)-1 Cependant, comment générer une séquence avec une moyenne de , c'est-à-dire avec...
Dans ce document , qui concerne la commande "set seed", les gens de Stata discutent des problèmes liés à la définition des graines lors de la génération de nombres pseudo-aléatoires. Un "non" notable est "ne pas utiliser en série la séquence de nombres naturels comme graines, car cela a un modèle...
Pourquoi les générateurs de nombres aléatoires comme runif()dans R ne génèrent-ils pas le même résultat à chaque fois? Par exemple: X <- runif(100) X génère des sorties différentes à chaque fois. Quelle est la raison de générer des sorties différentes à chaque fois? De quelles fonctionnalités...
Supposons un fichier de données avec plus de 80 millions de uns et de zéros, générés aléatoirement. A partir de ce fichier, nous voulons créer une liste d'entiers décimaux aléatoires. C'est le plan pour faire cette conversion. Divisez les 80 millions de chiffres en groupes de 4 chiffres binaires....
Fermé. Cette question est hors sujet . Il n'accepte pas actuellement les réponses. Voulez-vous améliorer cette question? Mettez à jour la question afin qu'elle soit sur le sujet pour la validation croisée. Fermé il y a 2 ans . Mon but: J'aimerais avoir une fonction qui prend une adresse e-mail et...
Comment dois-je échantillonner efficacement la distribution suivante? x ∼ B ( α , β) , x > k x∼B(α,β), x>k x \sim B(\alpha, \beta),\space x > k Si n'est pas trop grand, l'échantillonnage de rejet peut être la meilleure approche, mais je ne sais pas comment procéder lorsque k est grand....
J'exécute une simulation sur R et un cluster d'ordinateurs et j'ai le problème suivant. Sur chacun des ordinateurs X, je lance: fxT2 <- function(i) runif(10) nessay <- 100 c(mclapply(1:nessay, fxT2), recursive=TRUE) Il y a 32 ordinateurs, chacun avec 16 cœurs. Cependant, environ 2% des...
Dans mon programme, je dois exécuter N threads séparés chacun avec leur propre RNG qui est utilisé pour échantillonner un grand ensemble de données. J'ai besoin de pouvoir semer tout ce processus avec une seule valeur afin de pouvoir reproduire les résultats. Suffit-il d'augmenter simplement...
J'ai besoin de créer des vecteurs aléatoires de nombres réels a_i satisfaisant aux contraintes suivantes: abs(a_i) < c_i; sum(a_i)< A; # sum of elements smaller than A sum(b_i * a_i) < B; # weighted sum is smaller than B aT*A*a < D # quadratic multiplication with A smaller than D where...
J'ai besoin de tirer des nombres aléatoires à partir d'une distribution log-cauchy qui a la densité: Quelqu'un peut-il m'aider ou me diriger vers un livre / papier qui pourrait me montrer comment?F( x ; μ , σ) = 1x πσ[ 1 + ( l n ( x ) -
J'ai un GLMM du formulaire: lmer(present? ~ factor1 + factor2 + continuous + factor1*continuous + (1 | factor3), family=binomial) Lorsque j'utilise drop1(model, test="Chi"), j'obtiens des résultats différents de ceux que j'utilise à Anova(model, type="III")partir du package de voiture ou...
Étant donné une pièce avec un biais inconnu , comment puis-je générer des variations - aussi efficacement que possible - qui sont distribuées par Bernoulli avec une probabilité de 0,5? Autrement dit, en utilisant le nombre minimum de flips par variable
J'essaie d'adapter un modèle à temps discret dans R, mais je ne sais pas comment le faire. J'ai lu que vous pouvez organiser la variable dépendante dans différentes lignes, une pour chaque observation de temps, et utiliser la glmfonction avec un lien logit ou cloglog. En ce sens, j'ai trois...
J'utilise des séquences à faible écart depuis un certain temps pour les distributions uniformes, car j'ai trouvé leurs propriétés utiles (principalement en infographie pour leur apparence aléatoire et leur capacité à couvrir densément [0,1] de manière incrémentielle). Par exemple, des valeurs...
J'ai besoin de générer une liste de variables aléatoires soumises à des contraintes qui peuvent être exprimées sous la forme où est une matrice m \ fois n si \ bf {x } a n entrées. Dans tous les cas que je traite, n >> m , par exemple n sera environ 14 000 et m sera 50. Je ne sais pas quelle...
J'écris actuellement un algorithme de confidentialité différentielle en utilisant le mécanisme de Laplace. Malheureusement, je n'ai aucune expérience en statistique, donc beaucoup de termes me sont inconnus. Alors maintenant, je trébuche sur le terme: bruit de Laplace . Pour rendre un différentiel...
Je n'ai jamais suivi de cours de statistique officiel, mais en raison de mes recherches, je tombe constamment sur des articles qui appliquent plusieurs concepts statistiques. Souvent, je vais voir une description d'un processus de Monte Carlo appliqué à une situation donnée, et pour ce que je peux...
Supposons que vous ayez une pièce de monnaie équitable que vous pouvez lancer autant de fois que vous le souhaitez (éventuellement infiniment). Est-il possible de générer la distribution uniforme discrète sur( 1 , 2 , . . . , K )(1,2,...,k)(1,2,...,k), où kkkn'est pas une puissance de 2? Comment...