Discussion :

[magictol53]

Salut a tous,

Alors voilà, nouveau problème.Je devais prédire des points à l'aide de transformations polynomiales. Jusque là pas de problème, le script fonctionne et je calcule une erreur moyenne quadratique (rms) pour chaque point. Je souhaiterai
maintenant calculer l'erreur mais par validation croisée et plus précisément leave one out.
Je vous mets mon script pr mieux comprendre:

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
16
17
18
19
% load gcp.txt
load gcp
x = gcp(:,1:2);
y = gcp(:,3:4);
[xx Mx mx] = scale(x);
[yy My my] = scale(y);
a = yy(:,1);
b = yy(:,2);
c = a.*b;
one = ones(99,1);
T = [a b one];
% obtention de la matrice des coefficents coeff = (M*T)\(M*xx)
M = T';
H = inv(M*T)*M; 					
coeff = H*xx;
xxp=T*coeff;
xp = descale(xxp,Mx,mx);   % J'obtiens mes points prédits
scatter(xp(:,1),x(:,1));
rms = sqrt(mean(sum((xp-x).^2,2)));

Je poste mon problème ici car l'aide matlab pour ce problème est assez incomplète et il est très difficile d'appliquer les exemples pour mon cas. J'aimerais juste qu'il génère le modèle avec 98 points et qu'il le teste sur la 99e observation. Chaque point doit être testé ainsi

Merci d'avance

Pierrot

[magictol53]

Salut a tous,

je vais renouveler ma question et être plus clair car un précédent post n'a eu que peu de succès.
J'ai deux vecteurs du type:

x1 x'1
x2 x'2
x3 x'3
x4 x'4
x5 x'5
x6 x'6
J'aimerai effectuer un calcul répétitif entre ces deux vecteurs mais seulement entre cinq valeurs. Il faudrait donc que je supprime seulement la première ligne, puis seulement la seconde, puis seulement la troisième pour qu'à chaque fois il me reste bien cinq valeurs mais ce n'est jamais la même qui disparait.
Quelle est la boucle qui peut me permettre de faire cette chose?

Merci d'avance

Pierrot

[Jean Dumoncel]

Ca n'a pas l'air très compliqué, as-tu essayé avec une boucle for?

Sinon, la toolbox bioinfo possède des fonctions spécifiques à ce genre de problème comme crossvalind.

[magictol53]

Salut Magellan,

"Ca n'a pas l'air très compliqué" ne m'a que très peu aidé. Si j'en viens à poster quelque chose sur votre forum, c'est qu'auparavant j'ai cherché plusieurs heures..
Je me doute qu'il faille une boucle for et l'ensemble des fonctions traitant de la cross-validation (y compris crossvalind) ne permettent pas de faire ce que je souhaite. Je peine toujours à implémenter ma boucle donc si quelqu'un a une idée, je suis preneur.

Merci

[Jean Dumoncel]

S
"Ca n'a pas l'air très compliqué" ne m'a que très peu aidé. Si j'en viens à poster quelque chose sur votre forum, c'est qu'auparavant j'ai cherché plusieurs heures..

J'espérais que ces quelques mots te redonneraient au moins un peu d'espoir!!

Tu as donné le principe dans ton deuxième post, as-tu essayé de le mettre en oeuvre?
Il suffit de faire une boucle for qui supprimerait une ligne de tes données initiales, quelque chose comme :

Code :

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

1
2
3
4
5
for k=1:size(ma_variable,1)
     copie_de_ma_variable = ma_variable;
     copie_de_ma_variable(k,:) = [];
     traitement sur copie_de_ma_variable;
end

[magictol53]

Merci Magellan,

J'ai quelque chose de semblable mais sans résultat.
J'essaie ton code et je reviens au post pour une éventuelle solution.

A plus tard

[ol9245]

Code :

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

1
2
3
4
5
6
function [tout_moins_un, un_qui_reste] = separe(M, i)
L = size(M, 1) ;
cle = (1:L) == i ;
tout_moins_un = M(cle, :) ;
un_qui_reste = M(~cle, :) ;
end

à utiliser comme ça :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
function main
X = (1:10)' ; % un exemple
Y = 3 * X + 1 + rand(size(X)) ; % un exemple
 
for i = 1:10
    [Xa, Xb] = separe(X, i) ;
    [Ya, Yb] = separe(Y, i) ;
    % faire l'ajustement de Ya sur Xa
    % faire la vérification sur le couple (Xb, Yb)
end
end