9694

Minitab 22 - Beispiele für Simulationsmakros: Annäherung der Konstanten d2(n) und C4(n) für die erwartungstreue Schätzung

  • Erstellt am 14.9.2022
  • Überarbeitet am 30.4.2024
  • Software: Minitab 22, 21

Ich will in Form einer Simulation nachvollziehen, warum

  • die Zusammenfassung

    Spooled=140i=140si2

    der Teilgruppenstandardabweichungen in einer Prozessfähigkeitsanalyse mit 40 Teilgruppen der Größe von 5  um die Konstante C4((5-1)*40+1)=C4(161)

  • der Mittelwert der Teilgruppenstandardabweichungen

    Squer=140i=140si

    in einer Prozessfähigkeitsanalyse mit 40 Teilgruppen der Größe von 5  um die Konstante C4(5)

  • der Mittelwert

    Rquer=140i=140ri

    der Teilgruppenspannweite in einer Prozessfähigkeitsanalyse mit 40 Teilgruppen der Größe von 5  um die Konstante d2(5)

  • der Mittelwert

    MRquer=1199i=1199ri

    der gleitenden Spannweite der Länge 2 in einer Prozessfähigkeitsanalyse mit Stichprobenumfang 200 und Teilgruppengröße 1 um die Konstante d2(2)

korrigiert werden muss, um einen erwartungstreuen Schätzer für die innere Standardabweichung zu erhalten. Bei Teilgruppengröße > 1 ist die innere Standardabweichung um die Teilgruppenmittelwerte bereinigt. Gibt es ein Makro, mit dem ich das tun kann?

apspaketDas entsprechende APS-Paket ist über unseren ADDITIVE Professional Service erhältlich. Um das Paket zu erhalten, kontaktieren Sie unseren Support per E-Mail an Diese E-Mail-Adresse ist vor Spambots geschützt! Zur Anzeige muss JavaScript eingeschaltet sein. oder per Telefon unter +49 6172 5905 20 jeweils unter Angabe der APS-Paketnummer 1052.

Erläuterung

Hintergrund

Die von Minitab verwendeten Konstanten für die erwartungstreue Schätzung basieren auf dem Buch

D. J. Wheeler and D. S. Chambers. (1992). Understanding Statistical Process Control, Second Edition, SPC Press, Inc.

Man könnte sich die Schätzung der Standardabweichung durch die Werte SpooledC4((n-1)*g+1), SquerC4(n), Rquerd2(n) beziehungsweise MRquerd2(w) grob so vorstellen, dass die Werte Spooled, Squer, Rquer und MRquer für große Stichprobenumfänge näherungsweise proportional zur (inneren) Standardabweichung sein dürften.

Die jeweilige Konstante für die erwartungstreue Schätzung müsste beispielsweise in der Rquer-Methode so gewählt sein, dass im Mittel

Rquer=d2(n)*σinnerhalb

gilt, insbesondere weil ein erwartungstreuer Schätzer für einen Parameter diesen Parameter als Erwartungswert besitzt. Wenn die Grundgesamtheit die wahre innere Standardabweichung σinnerhalb=1 besitzt, gilt also im Mittel Rquer=d2(n).

Das APS-Paket Nr. 1052 besteht aus dem lokalen Makro ADD_sup_Prozessfaehigkeit_Konstanten_Erwartungstreu_Simulation. Dieses erstellt,

  • wenn es mit dem Unterbefehl MRQuer aufgerufen wird, für

    • einen Stichprobenumfang n (const.1)
    • eine Länge w (const.2) der gleitenden Spannweite zwischen 2 und 50 und
    • eine Anzahl anzit (const.3) der Iterationen

    in jeder Iteration n standardnormalverteilte Zufallszahlen (Mittelwert 0, Standardabweichung 1), berechnet den Mittelwert der gleitenden Spannweite der Länge w (kurz MR-quer)

  • wenn es mit dem Unterbefehl SPooled, SQuer oder RQuer aufgerufen wird, für

    • eine Teilgruppengröße ngrp (const.1)
    • eine Anzahl anzgrp (const.2) an Teilgruppen und
    • eine Anzahl anzit (const.3) der Iterationen

    in jeder Iteration n = ngrp*anzgrp in anzgrp Teilgruppen mit je ngrp Werten unterteilte standardnormalverteilte Zufallszahlen (Mittelwert 0, Standardabweichung 1) und berechnet

    • im Fall des Unterbefehls SPooled die noch nicht erwartungstreu korrigierte Zusammengefasste Standardabweichung (kurz S-pooled)
    • im Fall des Unterbefehls SQuer den Mittelwert der anzgrp Teilgruppenstandardabweichungen (kurz S-quer)
    • im Fall des Unterbefehls RQuer den Mittelwert der anzgrp Teilgruppenspannweiten (kurz R-quer)

Anschließend stellt es die MR-quer- bzw. S-pooled- bzw. S-quer- bzw. R-quer-Werte aller Iterationen in einem Histogramm dar, zusammen einer roten Referenzlinie an deren Mittelwert und zum Vergleich einer blauen Referenzlinie am Wert

  • d2(w) im Fall des Unterbefehls MRQuer.
  • C4((ngrp-1)*anzgrp+1) im Fall des Unterbefehls SPooled.
  • C4(ngrp) im Fall des Unterbefehls SQuer.
  • d2(ngrp) im Fall des Unterbefehls RQuer.

Die Idee ist, dass die beiden Referenzlinien bei vielen Iterationen nahe beieinander liegen müssten. Bitte legen Sie das Makro in dem Verzeichnis ab, dass Sie in Minitab unter Datei: Optionen als Speicherort für Makros festgelegt haben.

Beispiel 1

Die Simulation soll 50000 Iterationen mit je 100 Zufallszahlen enthalten, und die Konstante d2(2) für die erwartungstreue Schätzung soll mit dem Mittelwert der Mittelwerte der gleitenden Spannweiten der Iterationen approximiert werden.

Der Befehl für den Makroaufruf ist

%ADD_sup_Prozessfaehigkeit_Konstanten_Erwartungstreu_Simulation 100 2 50000;
  MRQuer.

Simulation für d2(2) mit 50000 Iterationen

 

Für jede Iteration werden 200 standardnormalverteilte Zufallszahlen erstellt.  

Prozessfaehigkeit_Konstanten_fuer_die_Erwartungstreue_Schaetzung_Simulation_01

Beispiel 2

Die Simulation soll 50000 Iterationen mit je 5*40 Zufallszahlen (40 Teilgruppen mit je 5 Werten) enthalten, und die Konstante C4(5) für die erwartungstreue Schätzung soll mit dem Mittelwert der S-quer-Werte approximiert werden.

Der Befehl für den Makroaufruf ist

%ADD_sup_Prozessfaehigkeit_Konstanten_Erwartungstreu_Simulation 5 40 50000;
  SQuer.

Simulation für C4(5) mit 50000 Iterationen

 

Für jede Iteration werden standardnormalverteilte Zufallszahlen erstellt, die in 40 Gruppen von jeweils 5 Werten unterteilt sind.                            

Prozessfaehigkeit_Konstanten_fuer_die_Erwartungstreue_Schaetzung_Simulation_02

Dieses Makro ist ein Beispiel für die Automatisierungsmöglichkeiten für Minitab. Trotz aller Sorgfalt übernehmen wir keine Gewährleistung für die Richtigkeit der Berechnungen und Ergebnisse.

Question?

Diese E-Mail-Adresse ist vor Spambots geschützt! Zur Anzeige muss JavaScript eingeschaltet sein., wenn Sie eine Frage zu diesem Artikel haben.

  • Download-ID: 1052