Durante il tentativo di capire come viene calcolata la soglia T^2 per l'estrapolazione dal aiuto sezione, la mia conclusione è che la soglia T^2 per un dato livello sigma dovrebbe essere costante per qualsiasi set di dati con lo stesso numero di colonne (p) e righe (n) senza dati mancanti. Questo perché la formula: T^2 soglia = {(n+1)(n-1)p/[n(np)}*F(p, np); è solo una funzione di "n" e "p".
Le due tabelle allegate di seguito hanno sia 10 colonne che 100 righe con diverse matrici di covarianza. In entrambe le tabelle, la formula per il rispettivo T^2 (non regolato) viene salvata come 11a colonna tramite Multivariata>Analisi anomala>T^2 (script salvato).In entrambe le tabelle, ho usato il profiler di questa formula per chiamare il controllo di estrapolazione (script salvato). Nella tabella Simulazione 1, la soglia è 18,20937. Nella tabella Simulazione 2, è 17.95714.
Dato che entrambe le tabelle hanno gli stessi "n" e "p", come può la loro soglia T^2 essere diversa? Oppure "n" e "p" si riferiscono a qualcos'altro diverso da righe e colonne?
Altro: nella stessa sezione della guida, si dice che il limite K-sigma di una distribuzione Normale viene utilizzato per definire la soglia. Se K=3, il valore predefinito, significa che il quantile F a 0,9973 viene utilizzato in F(p, np) poiché 3-sigma ~ 99,73%?
Mi sono appena reso conto che la formula: Soglia T^2 = {(n+1)(n-1)p/[n(np)}*F(p, np); deve essere utilizzato per meno di 10 valori non mancanti, altrimenti è T^2 threshold = T^2-hat + K*sigma(T^2), che spiega le diverse soglie nelle due tabelle.
