Discussion :

[Fibulavie]

Bonjour,

Débutante en R, je souhaite réaliser une boucle for dans laquelle le nom de mes variables dépend de mon incrément, et dans la boucle je réalise des opérations avec ces variables. (par ex, pour ce cas l'utilisation de paste ne fonctionnera pas)
données:
une db.class appelée grid, dans laquelle j'ai par ex 6 variables (dans la réalité j'en ai 300, d'où l'utilisation de ma boucle for):
simu1, simu2, simu3, result1, result2, result3
je souhaite créer des nouvelles variables dans cette db.class telles que:
test1=simu1+result1
en utilisant la fonction db.add

d'avance merci bcp pour votre aide!!

[tototode]

Bonjour,

Pourquoi dans le cas que tu indiques la fonction paste ou sprintf ne pourrait pas marcher ?

cdlt

[Fibulavie]

j'ai essayé la fonction paste mais dans le db.add elle ne fonctionne pas car elle ne reconnaît pas les noms de variables comme étant des variables justement, sur lesquelles on peut faire des opérations.

[tototode]

est-ce que tu peux nous donner la commande que tu utilises et la librairie dont tu te sers ?

après tu peux te servir de fonction comme eval, parse, subsitute pour recréer une ligne de commande et l'évaluer :

Code :

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

eval(parse(text = sprintf("add.db(test%d = simud%d + result%d)", 1, 1, 1)))

cdlt

[Fibulavie]

Bonjour,
Oui alors je travaille sous RGeostat, qui utilise le code de R.
ma fonction est la suivante:

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
20
21
22
23
24
25
26
27
28
test=db.create(nx=c(200,100)) # je créé ma db.class test
panel_test=db.create(nx=c(20,10),dx=c(10,10)) # je crée une 2ème db.class
test=db.add(test,Calorific) # je créé mes deux variables (dans mon il s'agissait de test; sauf qu'ici j'en ai deux)
test=db.add(test,Accu)
 
# ici je vais récupérer depuis 3 autres db.class des variables; toutes incrémentées de 1 à 4
for (i in 1:4){
test=db.grid.copy(simu100_thickness,test,names=paste("Simu.Thickness.S",i,sep="")) 
}
for (i in 1:4){
test=db.grid.copy(simu100_A0,test,names=paste("A0.Accumulation.S",i,sep=""))
}
for (i in 1:4){
test=db.grid.copy(simu100_A1,test,names=paste("A1.Accumulation.S",i,sep=""))
}
 
# là, je réalise mes calculs en écrasant à chaque fois mes variables Accu et Calorific
for (i in 1:4){
test[,5]=test[,i+9]+test[,i+13]
 
test[,4]=test[,5]%/%test[,i+5]
 
# sauf qu'avant de les écraser j'aimerais les copier dans une variable incrémentée et c'est là que mon paste ne fonctionne pas...
test_panel=db.stat.grid(test,panel_test,fun="mean",names="Calorific") # ça c'est juste une fonction qui me permet de récupérer pour chaque panel composé de 100 points la moyenne de Calorific de ces 100 points. Ma variable de sortie s'appelle Stats.Calorific
 
test_panel=db.add(test_panel,paste("MeanCalorific",1,sep="")=Stats.Calorific) # et là c'est le bout de ma fonction qui ne fonctionne pas... qu'à la base j'aurais bien utilisée dès le début pour avoir des Accu1, Accu2, Calorific1, Calorific2, etc...
 
}

D'avance merci beaucoup.

[tototode]

bonjour,

tu peux remplacer la ligne de commande qui ne va pas par celle là :

Code :

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

eval(parse(text = sprintf("test_panel=db.add(test_panel,%s=Stats.Calorific)", paste("MeanCalorific",1,sep=""))))

Après il faudra peut-être jouer avec les notions d'environnement :

Code :

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

eval(parse(text = sprintf("test_panel=db.add(test_panel,%s=Stats.Calorific)", paste("MeanCalorific",1,sep=""))), envir = .GlobalEnv)

La logique du truc c'est que tu utilises les fonctions paste, sprintf etc. pour écrire la ligne de commande que tu souhaites mais sous forme de chaines de caractères, ensuite tu transformes cette chaîne de caractère en expression avec la fonction parse et enfin tu executes cette expression avec la fonction eval.

Code :

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

1
2
3
com <- sprintf("test_panel=db.add(test_panel,%s=Stats.Calorific)", paste("MeanCalorific",1,sep=""))
expr <- parse(text = com)
eval(expr)

cdlt

[Theta]

Je pense qu'il y a moyen de faire quelque chose de plus propre qu'avec eval avec la fonction do.call().

[tototode]

J'aime bien l'expression plus 'propre', je ne savais pas l'eval sale :-)

Je suis preneur d'un exemple, parce que j'utilise fréquemment do.call mais dans ce cas précis je n'y arrive pas.
Merci

cdlt

[tototode]

Est-ce que ceci peut fonctionner ?

Code :

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

1
2
3
l1 <- list(test_panel, quote(Stats.Calorific))
names(l1)[2] <- paste("MeanCalorific",1,sep="")
test_panel <- do.call(add.db, l1)

C'est plus "propre" en terme de script mais pas d'objet parce qu'on est obligé de passer par une liste intermédiaire.
Bref j'ai pas mieux.

cdlt

[Theta]

@tototode: C'est plus où moins le code que j'avais en tête, je n'avais pas le temps de l'écrire au moment où j'ai posté. En espérant que ça marche.

Au niveau de la "propreté" c'est surtout que dans d'autres langages (js, php, C#...) l'utilisation d'eval est généralement déconseillée quand on peut faire autrement, pour des raisons de performances, sécurité, comportement inattendu, difficulté de débogage etc. En R ces problèmes se posent un peu moins.