Discussion :

[Hadrien35]

Bonjour à tous.

Je possède un tableau de données de type ind*var, et je cherche à opérer des permutations aléatoires de ce tableau (nommons le tab1).
Concrètement, je souhaite fixer un seuil de permutation x de façon à ce que dans mon tableau tab2, je retrouve mon tableau tab1 avec x% de données permutées au sein de chaque variable.
Je ne souhaite pas ici conserver le fait qu'à un individus i1 correspond un profil ligne donné (c'est à dire que l'on permute sa valeur pour la variable p1 avec la valeur p1 d'un individus i2 sans permuter leurs valeurs de leurs autres variables).

J'ai cherché via la fonction sample, mais elle me parait assez obscure...
J'aurai besoin d'un coup de main, je n'attend pas, bien entendu, un joli code qui me tombe tout cuit dans la bouche, seulement un petit coup de pouce sur 1 ou 2 lignes de codes permettant la permutation de données à hauteur de x%.

Merci d'avance!

[dickoa]

Salut
je viens de lire assez rapidement ta question, et elle est semble assez obscure......
plus sérieusement, si tu veux obtenir plus rapidement de l'aide un exemple avec un jeu de donnée et le code que tu as essayé est indispensable.

[Hadrien35]

Je vais essayer d'illustrer avec un exemple fictif...
Et en fait, pour faire plus simple, je ne considère qu'un vecteur de 100 données.

Code :

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

1
2
3
4
 
modif=0,05
a_perturber<-seq(0,1,100)
NbPermutation<-round((modif*length(a_perturber))/2)

Je souhaite perturber les données par permutation à hauteur de "modif" (ici 5%, j'effectue donc NbPermutation (ici 3).

Je cherche donc à permuter 3 couples de valeurs au sein de mon vecteur a_permuter. Je sais qu'une amie utilise la fonction sample pour extraire un échantillon d'individus à permuter, mais j'ai du mal à utiliser cette fonction.

Pour la suite de mon problème, je ferai cela sur chacune de mes variables (je ne cherche pas ici à conserver la structure des données par lignes).

J'espère que c'est plus clair,
désolé si j'ai été un peu obscure quand j'ai posé mon problème.

Merci d'avance!

[A. D.]

Bonjour,

Il y a sûrement moyen de faire plus simple, mais en bidouillant à la main on peut par exemple faire :

Code :

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

1
2
3
4
5
6
7
8
9
> a_perturber<-sample(1:5,100,replace=TRUE)
> modif<-0,05
> NbPermutation<-round((modif*length(a_perturber))/2)
> NbPermutation
[1] 2
 
> valeur_a_perturber<-sample(1:length(a_perturber),NbPermutation)
> valeur_a_perturber
[1] 57 80

On va donc "perturber" les 57ème et 80ème valeurs du vecteur "a_perturber".
Après, je ne sais pas ce que vous appellez "perturber" les valeurs, vous parlez de "permutation", donc est-ce que vous voulez par exemple faire une permutation circulaire (dans mon exemple par exemple, 1 devient 2, 2 devient 3, etc...)?
Dans ce cas, on peut faire quelque chose du genre :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
> a_perturber[valeur_a_perturber[1]]
[1] 3
> a_perturber[valeur_a_perturber[2]]
[1] 1
 
> a_perturber[valeur_a_perturber[1]] <- (a_perturber[valeur_a_perturber[1]]+1)%%5
> a_perturber[valeur_a_perturber[2]] <- (a_perturber[valeur_a_perturber[2]]+1)%%5
> if(a_perturber[valeur_a_perturber[1]]==0) { a_perturber[valeur_a_perturber[1]]<-5 }
> if(a_perturber[valeur_a_perturber[2]]==0) { a_perturber[valeur_a_perturber[2]]<-5 }
 
> a_perturber[valeur_a_perturber[1]]
[1] 4
> a_perturber[valeur_a_perturber[2]]
[1] 2

Bon voilà, ce n'est pas du tout propre et je ne sais pas si ça correspond à ce que vous recherchez, mais c'est une idée qui m'est venue comme ça...

Bonne continuation


Cordialement,

A.D.

[Hadrien35]

Merci beaucoup,
j'ai pu comprendre comment fonctionnait la fonction sample du coup, et pu régler mon problème en changeant un peu la formulation de NbPermutations.

En fait mon objectif était, avec les valeurs que tu trouves par exemple, d'échanger dans un tableau et pour une variable donnée (a_perturber dans l'exemple) la valeur de la ligne 57 avec la valeur de la ligne 80 (perturbation par permutation donc). Le tout en échangeant au final 5% des valeurs par variable.

Merci encore,

Hadrien

[A. D.]

Ok, je comprends mieux le sens de "permutation" ici.
Bon eh bien, contente d'avoir pu apporter un peu d'aide alors

Bien cordialement,

A.D.