|
|
|
Durchsichtiger 3D-Gitterplot |
|
| hukd321 |
Forum-Newbie
|
 |
Beiträge: 4
|
|
|
|
Anmeldedatum: 22.02.16
|
|
|
|
Wohnort: CGN
|
|
|
|
Version: R2013a
|
|
|
|
 |
|
Verfasst am: 22.02.2016, 00:50
Titel: Durchsichtiger 3D-Gitterplot
|
 |
| |
 |
 |
Hallo Community 
ich beschäftige mich gerade mit multivariaten Konfidenzregionen und habe nach googlen ein schönes Beispielbild gefunden, welches ich jetzt mal fix nachprogrammieren möchte:
Mir gefällt dabei ganz gut, dass man in der X-Y-Ebene die Konfidenzellipse sieht und dann mit etwas Abstand in der Z-Ebene die zugehörige Wahrscheinlichkeitsdichte.
Dass das in Matlab vermutlich nicht so leicht umsetzbar ist, damit habe ich mich bereits abgefunden, allerdings würde es mir schon reichen, hätte die Ellipse schön sichtbar in der X-Y-Ebene und ein durchsichtiges Gitter (ähnlich wie auf dem Bild) in auf der Z-Achse.
Ich hab dazu ein kurzes Toybeispiel geschrieben:
- Daten mittels y = 2x + 2 generieren und mit N(0,1) verrauschen
- Modell mittels nlinfit dranfitten (modelf = beta(1)*x + beta(2)), Startwerte beta0 = (1,1) --> liefert beta_opt und Cov_beta_opt
- approximierte Konfidenzellipse (90%) berechnen über Stichprobe von 20.000 aus N(beta_opt, Cov_beta_opt)
die mit Wkt < 90% werden rausgeschmissen und neu gezogen
- Ellipse wird geplottet (Kasten könnt ihr ignorieren)
- Wkt Verteilung über Ellipse wird geplottet (ACHTUNG, HIER ANDERE SKALIERUNG BEI X- UND Y-Achse)
bzw. "Draufsicht"
Wie krieg ich jetzt einen Plot hin, der in etwa so aussieht wie der da oben? Also mit dem "durchsichtigen" Gitter?
Anbei findet ihr das Toybeispiel. Den Kasten um die Ellipse könnt ihr dabei ignorieren (im Code alles, was mit '-k' und plot zu tun hat  ), da will ich noch was anderes versuchen.
Ich hoffe, ihr könnt mir helfen
LG hukd321
| Beschreibung: |
|
Download |
| Dateiname: |
ConfEllipsoid_toy.m |
| Dateigröße: |
1.7 KB |
| Heruntergeladen: |
388 mal |
|
|
|
|
|
|
| Harald |
Forum-Meister
|
|
Beiträge: 24.495
|
|
|
|
Anmeldedatum: 26.03.09
|
|
|
|
Wohnort: Nähe München
|
|
|
|
Version: ab 2017b
|
|
|
|
|
|
|
Verfasst am: 22.02.2016, 08:51
Titel:
|
|
Hallo,
versuchs mal mit mesh, oder mit Setzen der Eigenschaft 'FaceAlpha'.
Grüße,
Harald
|
|
|
|
| hukd321 |
Themenstarter
Forum-Newbie
|
|
Beiträge: 4
|
|
|
|
Anmeldedatum: 22.02.16
|
|
|
|
Wohnort: CGN
|
|
|
|
Version: R2013a
|
|
|
|
|
|
|
Verfasst am: 22.02.2016, 10:21
Titel:
|
|
Hallo Harald,
vielen Dank für deine Antwort.
Der Befehl:
erzeugt bei mir den Plot:
Der ist an sich ziemlich gut, nur leider viel zu dunkel. Das liegt ja wahrscheinlich an der 'Engmaschigkeit' des Grids.
Gibt es da irgendeinen Klugen weg, das zu ändern? Wenn ich das Grid nur halb so groß wähle, wird der Plot nicht mehr so schön:
Viele Grüße und vielen Dank
hukd
|
|
|
|
| Harald |
Forum-Meister
|
|
Beiträge: 24.495
|
|
|
|
Anmeldedatum: 26.03.09
|
|
|
|
Wohnort: Nähe München
|
|
|
|
Version: ab 2017b
|
|
|
|
|
|
|
Verfasst am: 22.02.2016, 10:29
Titel:
|
|
Hallo,
das sieht bei mir gut aus:
Mit welchem Release arbeitest du denn? Das sieht nach einem älteren Release aus?
Grüße,
Harald
|
|
|
|
| hukd321 |
Themenstarter
Forum-Newbie
|
|
Beiträge: 4
|
|
|
|
Anmeldedatum: 22.02.16
|
|
|
|
Wohnort: CGN
|
|
|
|
Version: R2013a
|
|
|
|
|
|
|
Verfasst am: 22.02.2016, 10:39
Titel:
|
|
Hi Harald,
ich arbeite mit R2013a.
so sieht das bei mir aus
:/
|
|
|
|
| Winkow |
Moderator
|
|
Beiträge: 3.842
|
|
|
|
Anmeldedatum: 04.11.11
|
|
|
|
Wohnort: Dresden
|
|
|
|
Version: R2014a 2015a
|
|
|
|
|
|
|
Verfasst am: 22.02.2016, 14:14
Titel:
|
|
|
|
|
| hukd321 |
Themenstarter
Forum-Newbie
|
|
Beiträge: 4
|
|
|
|
Anmeldedatum: 22.02.16
|
|
|
|
Wohnort: CGN
|
|
|
|
Version: R2013a
|
|
|
|
|
|
|
Verfasst am: 22.02.2016, 20:41
Titel:
|
|
das ist nicht schlecht
Jetzt muss ich das gitter nur noch besser hinkriegen. Auf der einen Seite nicht so engmaschig, auf der anderen seite aber die Mode glatter.
Und dann die Ellipse (bzw. die Punktwolke) von oben in die X-Y Ebene plotten...
|
|
|
|
|
|
|
Einstellungen und Berechtigungen
|
|
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
|  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.
|
|
