Discussion :

[leovenus]

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.

[mr_samurai]

Salut,

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

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

++

[leovenus]

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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.

[Dam2227]

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 ^^.

[mr_samurai]

Salut,


Essayes avec ceci :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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 .

++

[leovenus]

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 !

Lucie