Graphique de régression multiple

Bonjour, j'ai un problème urgent.

Je veux faire un graphe de régression multiple (j'ai 3 facteurs). Je ne comprends pas comment je peux avoir une représentation graphique de cela.
Je cherche depuis le début de l'après-midi, je n'en peux plus, est-ce que quelqu'un pourrait m'aider s'il vous plaît?

Merci,

L.

Salut,

je comprends à moitié, tu peux poster un bout de code ?

Tu cherches une représentation graphique d'une quantité à 4D ?

++

En fait, je dispose d'une température, et je veux voir l'effet de l'obliquité, la précession et le CO2 sur cette température par une analyse de régression.

Par exemple :

Code:

---

1 2 3 4 5 6 7 8

T = [-52.4698 -52.5111 -52.5573 -52.5917 -52.5068 -52.6141 -52.5618 -52.5595 -52.5753 -52.5500];   obliquite = [0.4162 0.4169 0.4173 0.4176 0.4175 0.4173 0.4168 0.4161 0.4152 0.4141];   precession = [-0.0360 -0.0322 -0.0255 -0.0166 -0.0060 0.0053 0.0162 0.0258 0.0333 0.0378];   co2 = [246.3996 245.4632 246.2506 246.0240 242.2793 239.9962 238.3901 241.9025 235.3792 235.4883];

---

En fait j'ai fait un graphe de la température en fonction de l'obliquité et une droite de régression linéaire :

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15

figure(1) plot(obliquite,T ,'.', 'linewidth',2) xlabel('obliquité (rad)', 'Fontsize',16); ylabel('température moyenne (C°)', 'Fontsize',16); set(gca, 'LineWidth',2,'FontSize',16); hold on; num=0; den=0; for i=1:501 num=num+(obl(i)-Xobl)*(tout(i)-Xtout); den=den+(obl(i)-Xobl)^2; end b=num/den; a=XT-b*Xobl; Y=a+b.*obl; plot(obl,Y ,'.r');

---

Voyant que cela avait une certaine dispersion (en fait j'ai 500 valeurs de chaque, je n'ai donné ci-dessus qu'une partie), on m'a demandé d'effectuer une régression multiple. J'ai l'équation Y = a0+a1*obliquite+a2*precession+a3*co2, je ne sais juste pas comment faire pour représenter ça graphiquement.

Est-ce possible?

Dans le cas où je n'aurais que 2 facteur (ex : obliquite et precession), comment puis-je représenter cela sur un graphe?

Merci beaucoup, les stat, c'est pas mon fort !

L.

Tu peux faire un graph en utilisant la fonction mesh par exemple(deux facteurs et le troisième constant).

Une idée serait que tu fasses plusieurs graphes (2 variables et le troisième facteur constant que tu incrémenterait à chaque graphe), que tu les enregistres et que tu utilises la fonction movie, qui te donnerait la quatrième dimension : le temps ^^.

Salut,


Essayes avec ceci :

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15

% Données ----------------- T = [-52.4698 -52.5111 -52.5573 -52.5917 -52.5068 -52.6141 -52.5618 -52.5595 -52.5753 -52.5500]; obliquite = [0.4162 0.4169 0.4173 0.4176 0.4175 0.4173 0.4168 0.4161 0.4152 0.4141]; precession = [-0.0360 -0.0322 -0.0255 -0.0166 -0.0060 0.0053 0.0162 0.0258 0.0333 0.0378]; co2 = [246.3996 245.4632 246.2506 246.0240 242.2793 239.9962 238.3901 241.9025 235.3792 235.4883];   % regression ------------ X = [ones(size(T')) obliquite' precession' co2']; a = X\(T'); % Les coefficients [a0 a1 a2 a3]   % Fig 1 ------------------------ figure(1) plot(T,T,'.r'); hold on plot(T, X*a,'ob');

---

L'idée est de tracer la (T,T) et (T,regression(T)) pour voir si cela colle bien ;) .

++

Salut!

Merci dam2227, merci mr samurai !


J'ai utilisé ce que tu m'as conseillé samurai, cela correspond à ce que je voulais !
C'est génial !
Je ne sais pas comment j'aurais fait sans toi ! :hola:

Lucie