Discussion :

[orland]

Bonjour,

J'ai un code en C et je n'arrive pas à l'écrire en R

Voici le code en C :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
      cste = 1.0/(double)N ;
 
      for ( i = 2 ; i < (N-3) ; i++ )
	{
	  ang = 0.0 ;
	  m = 0 ;
	  for ( j = -1 ; j <= +3 ; j++ )
	    {
	      ang = ang + atan2 ( (arr[i+j]-arr[i+j-1]), ctse ) ;  /* normalisation */
	      m = m + 1 ;
	    }
 
	  angle[i] = ang/(double)m ;   /* moyenne */
	}

Voici le code que j'ai tapé, mais malheureusement ça ne marche pas

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
am = function(data)
{
library(dagR)
n          = length(data)
cste       = 1.0/n
angle      = NULL
for(i in 3 : n-3)
{
bb = paste("b",i,sep="") 
k = NULL 
for (j in -1:3)
    {
     m=NULL
     aa = paste("a",j,sep="")
     k = c(m,assign(aa,atan2(data[i+j] - data[i+j-1],cste), pos=1))
    m=m+1
     }
  angle = c(angle,assign(bb,(k/m),pos=1))
}
return(angle)
}

Veuillez s'il vous plaît voir et m'aider à trouver une solution.
Merci d'avance

[A. D.]

Bonjour,

Si j'ai bien compris, vous souhaitez réécrire votre algorithme C en R ?

A mon avis, vous vous compliquez la vie pour rien ! Je pense que les "assign" ne sont absolument pas indispensables, avez-vous vraiment besoin de garder ces éléments "en plus" ?

Bonne continuation


Cordialement,

A.D.

[orland]

Bonjour

Si j'ai bien compris, vous souhaitez réécrire votre algorithme C en R ?

Oui c'est ce que je veux

A mon avis, vous vous compliquez la vie pour rien ! Je pense que les "assign" ne sont absolument pas indispensables,

C'est ce que je pensais faire ou plutot l'idée que j'ai trouvée ...

avez-vous vraiment besoin de garder ces éléments "en plus" ?

Non je n'ai pas besoin des éléments "en plus"

Pouvez-vous, s'il vous plaît, m'aider à améliorer le code que j'ai tapé ou me donner le code correct, avec un grand merci d'avance

[A. D.]

Re-bonjour,

Je n'ai pas vérifié la "concordance" entre votre code R et votre code C, mais concernant votre code R, il ne manque pas grand chose pour que cela fonctionne à mon avis...

Essayez de supprimer les "assign" pour voir ce que ça donne. Il restera ensuite quelques erreurs que vous devriez pouvoir corriger par vous-même en reprenant votre code ligne par ligne.

Bonne continuation


Cordialement,

A.D.

[orland]

Ok i'll try
Et j'espère que je trouverais la solution
Je vous transmettrais les résultats que je trouverais
Cordialement B.H

[orland]

Bonjour,
j'ai essayé avec ce code :

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
am = function(data)
{
library(dagR)
n          = length(data)
cste       = 1.0/n
angle      = NULL
for(i in 3 : n-3)
{
bb = paste("b",i,sep="") 
k = NULL 
for (j in -1:3)
    {
     aa = paste("a",j,sep="")
     k = c(m,assign(aa,atan2(data[i+j] - data[i+j-1],cste), pos=1))
    }
  angle = c(angle,assign(bb,(k/5),pos=1))
}
return(angle)
}

Et çane marche toujours pas !
Je ne sais pas pourquoi vous m'avez dit de supprimer les "assign"
Je suis vraiment bloqué