WICHTIG: Der Betrieb von goMatlab.de wird privat finanziert fortgesetzt. - Mehr Infos...

Mein MATLAB Forum - goMatlab.de

Mein MATLAB Forum

 
Gast > Registrieren       Autologin?   

Partner:




Forum
      Option
[Erweitert]
  • Diese Seite per Mail weiterempfehlen
     


Gehe zu:  
Neues Thema eröffnen Neue Antwort erstellen

Test auf Normalverteilung / sonstige Verteilung

 

Tuomas

Gast


Beiträge: ---
Anmeldedatum: ---
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 07.09.2014, 10:08     Titel: Test auf Normalverteilung / sonstige Verteilung
  Antworten mit Zitat      
Hallo,

ich habe einen Vektor mit ca. 8000 Werten zwischen -55 und +80.
Dazu möchte ich prüfen, ob für diese Werte von einer Normalverteilung oder von einer anderen Verteilung ausgenagnen werden kann.

Mein Vorgehen war bisher so:

1. Ich habe mit der Distribution Fitting Toolbox verschiedene Verteilungen angepasst. Bei der grafischen Darstellung fällt auf, dass die Werte des Vektors symmetrisch verteilt sind, jedoch im Vergleich zur Normalverteilung eine Konzentration der Werte um den Mittelwert aufweisen (d.h. die Glockenkurve ist spitzer und höher als bei der Normalverteilung).
Passt man eine t-location-scale (Student-t-Verteilung) an, so ergibt sich rein visuell eine sehr gute Abbildung.

Die Distribution Fitting Toolbox liefert folgende Outputs:

Für die Anpassung der Normalverteilung
Code:
Distribution:    Normal
Log likelihood:  -28892.6
Domain:          -Inf < y < Inf
Mean:            0.003129
Variance:        42.8931

Parameter  Estimate  Std. Err.
mu         0.003129  0.0699748
sigma       6.54928  0.0494839

Estimated covariance of parameter estimates:
       mu           sigma      
mu      0.00489647  3.94128e-20
sigma  3.94128e-20   0.00244865


Für die Anpassung der t-location-scale Verteilung
Code:
Distribution:    t Location-Scale
Log likelihood:  -27583
Domain:          -Inf < y < Inf
Mean:            -0.345743
Variance:        72.6054

Parameter  Estimate   Std. Err.
mu         -0.345743  0.0485707
sigma        3.63633  0.0574168
nu           2.44535  0.0831268

Estimated covariance of parameter estimates:
       mu           sigma        nu        
mu      0.00235911  0.000231406  0.000339635
sigma  0.000231406   0.00329669   0.00327788
nu     0.000339635   0.00327788   0.00691007


Frage 1: Lässt sich allein aufgrund von diesen Werten schon sagen, dass es sich um eine Normalverteilte oder Student-t Verteilte Werte handelt?



Weiteres Vorgehen:
Ich habe dann versucht einen Chi2-Test sowie einen KSTest zu machen um die Hypothese auf Normalverteilung bzw. Student-t Verteilung zu legegen oder zu widerlegen. Dabei habe ich jedoch einige Probleme:

1. Ich bin nicht sicher wie der Output aus Matlab interpretiert werden muss.
2. Ich bin nicht sicher, ob ich die Werte richtig in die kstest bzw. chi2gof Funktion eingebe.

Kann mir dazu jemand vielleich ein Beispiel geben wo auch zu erkennen ist, wie die eingegebenen Vektoren aussehen?

Viele Grüße
Tuomas


Harald
Forum-Meister

Forum-Meister


Beiträge: 24.499
Anmeldedatum: 26.03.09
Wohnort: Nähe München
Version: ab 2017b
     Beitrag Verfasst am: 07.09.2014, 10:21     Titel:
  Antworten mit Zitat      
Hallo,

ich habe dazu gerade was geschrieben:
http://www.gomatlab.de/kolomogorov-.....v-test-kstest-t34768.html

Zitat:
jedoch im Vergleich zur Normalverteilung eine Konzentration der Werte um den Mittelwert aufweisen (d.h. die Glockenkurve ist spitzer und höher als bei der Normalverteilung)

Das kann auch lediglich heißen, dass die Standardabweichung klein ist.

Zitat:
Frage 1: Lässt sich allein aufgrund von diesen Werten schon sagen, dass es sich um eine Normalverteilte oder Student-t Verteilte Werte handelt?

Meines Erachtens kann man lediglich durch Vergleich der Log Likelihoods sagen, welche Verteilung besser ist.

Zitat:
1. Ich bin nicht sicher wie der Output aus Matlab interpretiert werden muss.
2. Ich bin nicht sicher, ob ich die Werte richtig in die kstest bzw. chi2gof Funktion eingebe.

Dann stelle doch bitte deinen Code zur Verfügung, und man kann dir sagen, ob das so sinnvoll ist.
Beispiele finden sich in der Dokumentation.

Grüße,
Harald
Private Nachricht senden Benutzer-Profile anzeigen
 
Tuomas

Gast


Beiträge: ---
Anmeldedatum: ---
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 07.09.2014, 11:58     Titel:
  Antworten mit Zitat      
Ich habe mal 2 Grafiken aus der Distibution Fitting Toolbox angehängt. Angepasst habe ich die Werte an eine Normalverteilung und eine T-Verteilung.

Danach habe ich versucht eine kstest zu machen. Dabei bin ich äquivalent zu dem Beispiel in der Matlab Doku vorgegangen. Der Vektor mit den Werten ist x. Die Variable tloc ist ein Export aus der Distribution Fitting Toolbox und enthält u.a. die Verteilungsparameter:

Eingabe:

Code:
testcdf_tloc = [x,cdf('tlocationscale',x,pd_tloc.mu,pd_tloc.sigma,pd_tloc.nu)];
[h,p,k,c]=kstest(x,'CDF',testcdf_tloc)


Ausgabe:
Code:
h = 1

p = 9.9428e-05

k = 0.0238

c =  0.0145


Wenn ich das Ergebnis richtig interpretiere bedeutet dass, dass man nicht davon ausgehen kann, dass x t-Verteilt ist. Dieses Ergebnis ist für mich aufgrund der Grafiken aus der Distribution Fitting Toolbox schwer nachvollziehbar. Daher stelle ich mir die Frage ob die Vorgehensweise korrekt ist.

Für die Normalverteilung sieht das ganze wie folgt aus:
Input:
Code:

testcdf_norm = [x,cdf('norm',x,pd_norm.mu,pd_norm.sigma)];
[h,p,k,c]=kstest(x,'CDF',testcdf_norm);


Output:
Code:
h =  1

p = 2.2976e-76

k = 0.0998

c = 0.0145


dfit2.jpg
 Beschreibung:

Download
 Dateiname:  dfit2.jpg
 Dateigröße:  24.66 KB
 Heruntergeladen:  465 mal
dfit1.jpg
 Beschreibung:

Download
 Dateiname:  dfit1.jpg
 Dateigröße:  22.08 KB
 Heruntergeladen:  370 mal
 
Neues Thema eröffnen Neue Antwort erstellen



Einstellungen und Berechtigungen
Beiträge der letzten Zeit anzeigen:

Du kannst Beiträge in dieses Forum schreiben.
Du kannst auf Beiträge in diesem Forum antworten.
Du kannst deine Beiträge in diesem Forum nicht bearbeiten.
Du kannst deine Beiträge in diesem Forum nicht löschen.
Du kannst an Umfragen in diesem Forum nicht mitmachen.
Du kannst Dateien in diesem Forum posten
Du kannst Dateien in diesem Forum herunterladen
.





 Impressum  | Nutzungsbedingungen  | Datenschutz | FAQ | goMatlab RSS Button RSS

Hosted by:


Copyright © 2007 - 2025 goMatlab.de | Dies ist keine offizielle Website der Firma The Mathworks

MATLAB, Simulink, Stateflow, Handle Graphics, Real-Time Workshop, SimBiology, SimHydraulics, SimEvents, and xPC TargetBox are registered trademarks and The MathWorks, the L-shaped membrane logo, and Embedded MATLAB are trademarks of The MathWorks, Inc.