expérience de simulationJe veux faire une expérience de simulation, générer 12 variables aléatoires, chaque variable contient 100 observations, de sorte que les variables soient soumises à un problème d'hétérogénéité et suivent la distribution normale
L'expérience est répétée 10 000 fois.
# Je ne sais pas comment faire ça avec JMPpro
Merci
Modélisation statistique avancéeles fenêtresRe : expérience de simulation"comment rendre les variables hétérogènes en variance" - cela signifie-t-il que chaque colonne a besoin d'un sigma différent ? Avez-vous besoin de 12 valeurs sigma choisies au hasard ?
"Problème d'hétéroscédasticité" - cela signifie-t-il que vous avez une fonction en tête pour changer le sigma en fonction de mu ? Et une fonction pour choisir mu, probablement au hasard ? (Je n'ai aucune expertise ici, je devine juste après 30 secondes de lecture.)
Je pense que vous voudrez le faire en utilisant des boucles for imbriquées. Ma suggestion initiale est trop simple. Vous pouvez également utiliser des formules de colonnes, mais un script JSL sera plus facile à gérer. Les boucles imbriquées ressembleraient à ceci
m = J( 100 * 10000, 12, 0 );
for(row=1, row<=nrows(m), row+=1,
for(col=1, col<=ncols(m), col+=1,
m[row,col] = HSD(row,col);
);
);<BR /><BR />dt = astable(m)<BR /><BR />...
Vous devez écrire la fonction écrite par l'utilisateur "HSD" avec ce que vous voulez réellement.
Vous souhaiterez peut-être également modifier les 10 000 répétitions en 2 ou 3 répétitions jusqu'à ce que vous décidiez ce que l'analyse doit faire. J'ai utilisé le résumé car il regroupe toutes les 10 000 analyses dans un tableau de données. Vous souhaiterez peut-être exécuter une plate-forme 10 000 fois, récupérer quelques résultats, puis fermer la plate-forme. Vous ne voulez pas que les plates-formes 10 000 soient ouvertes en même temps, car JMP manquerait de mémoire.
Re : expérience de simulationCher Monsieur
Tout d'abord, merci beaucoup pour votre aide.
En fait, votre script ci-dessus se rapproche de ce que je veux,
mais je ne sais pas comment rendre les variables hétérogènes en variance.
Je ne sais pas comment modifier votre script, tel que :
Nombre de variables aléatoires = 12
Le nombre d'observations dans chaque variable = 100
Et que les variables résultantes souffrent d'un problème d'hétéroscédasticité
Avec possibilité de répéter l'expérience 10 000 fois
Salutations
Re : expérience de simulationregarde une recherche comme celle-ci https://www.google.com/search?q=site%3Ajmp.com+Heteroscedasticity pour trouver des questions précédentes qui pourraient avoir de bonnes réponses dans la communauté. Dans l'interface graphique de recherche Google, vous pouvez le faire avec site : jmp.com pour limiter les résultats :
site:jmp.com Heteroscedasticitysite:jmp.com Heteroscedasticitysite:jmp.com Hétéroscédasticité
Je ne suis pas statisticien. Les données ci-dessus sont aléatoires. La source que vous indiquez suggère que les données devraient se répéter d'une manière ou d'une autre au fil des saisons... pas complètement aléatoire. Si 1..12 représente les mois de l'année, cela n'a peut-être pas de sens d'avoir une colonne janvier, une colonne février, ... Il pourrait être beaucoup plus logique d'avoir une colonne de date et une colonne de valeur.
| date |
valeur |
| 2020/01 |
7 |
| 2020/02 |
6 |
| 2020/03 |
8 |
J'ai levé les yeux Hétéroscédasticité et je suppose que vous aurez besoin d'une approche différente pour générer vos nombres aléatoires afin qu'ils présentent réellement ce trait. Sinon, c'est juste aléatoire sans rien d'intéressant mais le mu et le sigma sont passés à la normale aléatoire, comme le montre le graphique ci-dessus.
Juste deviner, et prendre des risques... Je pense que vous devrez générer un mu aléatoire (normal ? uniforme ?) et un sigma basé sur ce mu, puis utiliser RandomNormal(mu,sigma). Ou je pourrais me tromper complètement.
Re : expérience de simulation
la source:
Autorégression périodique à coefficient aléatoire (repec.org)
Re : expérience de simulationCher Monsieur Craige_Hales
Merci beaucoup pour votre aide,
Je voulais faire une expérience de simulation :
Nombre de variables aléatoires = 12
Le nombre d'observations dans chaque variable = 100
Et que les variables résultantes souffrent d'un problème d'hétéroscédasticité
# Avec possibilité de répéter l'expérience 10 000 fois
Re : expérience de simulationLe résumé pourrait être un bon début
</CODE></P>
<PRE><CODE class=" language-jsl">dt = As Table( J( 100 * 10000, 12, Random Normal( 0, 1 ) ) );
dt << New Column( "group", formula( Floor( (Row() - 1) / 100 ) + 1 ) );
dt << runformulas;
dt:group << deleteformula;
dt << Go To( group );
dt << Move Selected Columns( To first );
dtsummary = dt << Summary(
Group( :group ),
Mean( :Col1 ),
Min( :Col1 ),
Max( :Col1 ),
Std Dev( :Col1 ),
Mean( :Col2 ),
Min( :Col2 ),
Max( :Col2 ),
Std Dev( :Col2 ),
Freq( "None" ),
Weight( "None" )
);
dtsummary << Graph Builder(
Size( 653, 452 ),
Show Control Panel( 0 ),
Variables(
X( :group ),
Y( :"Mean(Col1)"n ),
Y( :"Min(Col1)"n, Position( 1 ) ),
Y( :"Max(Col1)"n, Position( 1 ) ),
Y( :"Std Dev(Col1)"n, Position( 1 ) ),
Y( :"Mean(Col2)"n, Position( 1 ) ),
Y( :"Min(Col2)"n, Position( 1 ) ),
Y( :"Max(Col2)"n, Position( 1 ) ),
Y( :"Std Dev(Col2)"n, Position( 1 ) )
),
Elements(
Points(
X,
Y( 1 ),
Y( 2 ),
Y( 3 ),
Y( 4 ),
Y( 5 ),
Y( 6 ),
Y( 7 ),
Y( 8 ),
Legend( 7 )
)
)
);
</CODE></P>
<P><SPAN class="lia-inline-image-display-wrapper lia-image-align-inline" image-alt="Color!" style="width: 812px;"><SPAN class="lia-inline-image-display-wrapper" image-alt="Color!" style="width: 812px;"><span class="lia-inline-image-display-wrapper" image-alt="Color!" style="width: 812px;"><img src="https://us.v-cdn.net/6038925/uploads/images/41991iFA52067AFA73F318/Untitled%204%20By%20%28group%29%20-%20Graph%20Builder.png" role="button" title="Untitled 4 By (group) - Graph Builder.png" alt="Color!" /><span class="lia-inline-image-caption" onclick="event.preventDefault();">Color!</span></span><SPAN class="lia-inline-image-caption" onclick="event.preventDefault();">Color!</SPAN></SPAN><SPAN class="lia-inline-image-caption" onclick="event.preventDefault();">Couleur!</SPAN></SPAN></P></LINGO-BODY><LINGO-SUB id="lingo-sub-481697" slang="en-US" mode="NONE">Re : expérience de simulation</LINGO-SUB><LINGO-BODY id="lingo-body-481697" slang="en-US" mode="NONE"><P>Je ne sais pas quel est le problème d'hétérogénéité, mais le reste pourrait ressembler à ceci :</P>
<P> </P>
<PRE><CODE class=" language-jsl">dt = As Table( J( 100 * 10000, 12, Random Normal( 0, 1 ) ) );
dt << New Column( "group", formula( Floor( (Row() - 1) / 100 ) + 1 ) );
dt << runformulas;
dt:group << deleteformula;
dt << Go To( group );
dt << Move Selected Columns( To first );
Faites votre analyse avec la variable de groupe.
12 million random normal numbers12 million random normal numbers12 millions de nombres normaux aléatoires
