Discussion :

[Jesk34]

Bonjour,

J'ai un autre soucis dans R.
Je veux faire une génération de séquence.

Code :

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

1
2
matrice = c()
Sequence= seq(1,200,by=1) # j'ai 200 colonnes

Le problème c'est que mes données vont de 400 à 600 et non 0 à 200.
N'y a t'il pas moyen de lui demander de me le faire de 400 à 600 parce qu'après je dois faire un graphe en fonction de "Sequence" et qu'il m'affiche forcément 0 à 200 sur l'axe des abcisses.

Tant que j'y suis, chaque nombre correspond à une autre séquence distribuée non régulièrement :
exemple :
400 = 567
401 = 571
402 = 578
403 = 687 etc....
Peut-on remplacer les chiffres d'origine par les nouveaux et respecter les échelles sur le graphe ?

Je ne sais pas si je me suis bien exprimée
Merci

[A. D.]

Bonjour,

Votre problème ne me semble pas très clair mais voici quelques éléments qui pourront peut-être vous aider...

Une séquence de 400 à 600 par pas de 1 :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
sequence<-seq(400,600,by=1)
sequence
 [1] 400 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 416 417
 [19] 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435
 [37] 436 437 438 439 440 441 442 443 444 445 446 447 448 449 450 451 452 453
 [55] 454 455 456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471
 [73] 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 489
 [91] 490 491 492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507
[109] 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525
[127] 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543
[145] 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560 561
[163] 562 563 564 565 566 567 568 569 570 571 572 573 574 575 576 577 578 579
[181] 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597
[199] 598 599 600

Tant que j'y suis, chaque nombre correspond à une autre séquence distribuée non régulièrement :
exemple :
400 = 567
401 = 571
402 = 578
403 = 687 etc....
Peut-on remplacer les chiffres d'origine par les nouveaux et respecter les échelles sur le graphe ?

Pour moi, le mieux serait de prendre comme vecteur d'abscisses le vecteur correspondant à vos valeurs (par exemple : c(567,571,578,687,...) ), et comme vecteur d'ordonnées les valeurs correspondantes.
En effet, les écarts entre les différentes valeurs seront respectés par défaut par la fonction "plot" et les valeurs affichées sur l'axe seront celles du vecteur des abscisses.

En espérant que cela puisse vous aider. Sinon le mieux serait que vous postiez un court mais représentatif extrait de vos données et du code testé.

Bonne continuation


Cordialement,

A.D.

[Jesk34]

Bonjour,

Votre problème ne me semble pas très clair mais voici quelques éléments qui pourront peut-être vous aider...

Une séquence de 400 à 600 par pas de 1 :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
sequence<-seq(400,600,by=1)
sequence
 [1] 400 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 416 417
 [19] 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435
 [37] 436 437 438 439 440 441 442 443 444 445 446 447 448 449 450 451 452 453
 [55] 454 455 456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471
 [73] 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 489
 [91] 490 491 492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507
[109] 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525
[127] 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543
[145] 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560 561
[163] 562 563 564 565 566 567 568 569 570 571 572 573 574 575 576 577 578 579
[181] 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597
[199] 598 599 600

Le problème si je fais ça c'est qu'il va me séquencer de la 400ème colonne à la 600ème colonne de mon tableau... Alors que moi j'ai bien 200 colonnes mais qui démarrent pas à 1.

[A. D.]

Le problème si je fais ça c'est qu'il va me séquencer de la 400ème colonne à la 600ème colonne de mon tableau... Alors que moi j'ai bien 200 colonnes mais qui démarrent pas à 1.

Dans ce cas, votre problématique n'est pas bien posée...

Si vous souhaitez récupérer les 200 premières colonnes d'un tableau appelé par exemple "tab", vous pouvez procéder ainsi :

Code :

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

tab[,1:200]

Si les colonnes qui vous intéressent ne démarrent pas à la première mais par exemple à la 3ème :

Code :

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

tab[,3:203]

Vous pouvez également créer un vecteur "vect" qui contient les numéros de colonnes qui vous intéressent, puis faire comme ceci pour les récupérer :

Code :

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

tab[,vect]

Je ne comprends pas bien ce que vous cherchez à faire... Le mieux serait que vous postiez un extrait de votre jeu de données (au moins pour que nous ayons une idée de la structure), et que vous précisiez votre problématique.


Cordialement,

A.D.

[Jesk34]

Jeu de données :
testforum.xlsx

Code :

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

1
2
3
4
5
6
7
8
9
matrice = c()
SeqNiv= seq(1,10,by=1) 
for (i in SeqNiv)
{
a = hist(gbl[,i],breaks=seq(0,1,by=0.01),plot=TRUE)
matrice = cbind(matrice,a$density/max(a$density))
}
dim(matrice)
image(SeqNiv,seq(0,1,0.01),t(matrice),col=gray(1-seq(0,1,length.out=100)))

Ce que je veux : [Edit : mon premier jet était sur 20 colonnes... d'où le 20, dans le fichier envoyéc'est 10 )


Sauf que j'aimerais dans un premier temps qu'il y ait 400 à 410 en abcisse.
En fait ce sont des carottages faits dans le sol, ce sont des profondeurs et chaque profondeur correspond à une date 14C. Si je pouvais avoir ces dates en abcisse à la place de la profondeur ce serait top.

[A. D.]

Bon, je ne suis toujours pas sûre de bien comprendre : ce que vous voudriez c'est simplement qu'au lieu des valeurs 0 à 10 vous ayez 400 à 410 en abscisse sur votre graphe?
Si oui, vous pouvez par exemple faire comme ceci :

Code :

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

1
2
image(SeqNiv,seq(0,1,0.01),t(matrice),col=gray(1-seq(0,1,length.out=100)),xaxt="n")
axis(1,at=0:10,labels=seq(400,410,1)) # vous pouvez remplacer l'argument "labels" par un vecteur contenant vos dates par exemple

En espérant que ceci vous aide à avancer...


Cordialement,

A.D.

[Jesk34]

il semble que la fonction xaxt ne fonctionne qu'avec la fonction plot et pas image.

[A. D.]

il semble que la fonction xaxt ne fonctionne qu'avec la fonction plot et pas image.

Etrange, je l'ai testé sous Windows et sous linux (R 2.7 et R 2.14) et cela semble fonctionner très bien.

J'ai testé avec le code suivant (exemple obtenu dans l'aide de la fonction "image") :

Code :

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

1
2
3
4
x <- y <- seq(-4*pi, 4*pi, len=27)
r <- sqrt(outer(x^2, y^2, "+"))
image(z = z <- cos(r^2)*exp(-r/6), col=gray((0:32)/32),xaxt="n")
axis(1,at=0:1,labels=c("A","B"))

Après c'est peut-être au niveau de l'utilisation de la fonction "axis" qu'il faut ajuster les paramètres en fonction de votre cas précis.


Cordialement,

A.D.

[Jesk34]

Ca fonctionne pour mes profondeurs avec votre aide et un peu de bidouille.
Merci
Je galère toujours avec mes dates, j'ai essayé de créer un vecteur ( la flemme d'écrire 200 dates) mais impossible de l'intégrer.

[A. D.]

Bonjour,

Afin que nous puissions vous aider au mieux, il serait bien que vous essayiez de développer au maximum vos posts sur le forum : quelles sont vos données? de quoi sont-elles composées? que cherchez-vous à faire précisément? qu'avez-vous déjà essayé (exemple de code)?

Je galère toujours avec mes dates, j'ai essayé de créer un vecteur ( la flemme d'écrire 200 dates) mais impossible de l'intégrer.

Dans quoi sont contenues / où sont les dates dont vous parlez (je ne les vois pas dans le fichier que vous avez mis en pièce jointe)?
A quel endroit souhaitez vous intégrer le vecteur contenant ces dates? Au niveau de l'axe des abscisses du graphe? Quel code avez-vous utilisé pour le moment?


Cordialement,

A.D.

[Jesk34]

Code :

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

1
2
3
vect<-c(1580.561642,1565.721486,[j'ai copié les 200 dates],-10834.50975,-10988.19226)
image(SeqNiv,seq(0,1,0.01),t(matrice),col=gray(1-seq(0,1,length.out=100)),xaxt="n")
axis(1,at=0:221,labels=vect)

Le problème c'est qu'elles ne sont pas distribuées selon un pas de temps régulier et le graphe que ça me sort ne respecte pas l'échelle.

[A. D.]

Ok, donc si j'ai bien compris, "vect" est un vecteur des "dates" que vous voudriez voir s'afficher en abscisse. De plus, ces dates sont des valeurs réelles et vous voudriez que l'échelle de l'axe des abscisses respecte les écarts entre deux "dates" consécutives. C'est bien cela?

Pour cela, il va falloir remplacer l'élément correspondant aux abscisses dans la fonction "image".

Si votre vecteur de "dates" est ordonné (mais il ne semble malheureusement pas l'être dans votre exemple), vous pouvez utiliser directement ce vecteur comme paramètre correspondant aux abscisses dans la fonction image.
Sinon, je crains fort qu'il faille créer vous même le vecteur à utiliser. Une idée pour vous servir de base (mais il y a sûrement moyen de faire mieux) :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
vect<-c(1580.561642,1565.721486,-10834.50975,-10988.19226)
 
ecart<-abs(diff(vect))
 
ecart
[1]    14.84016 12400.23124   153.68251
 
 
SeQ<-c(0,cumsum(ecart))
 
SeQ
[1]     0.00000    14.84016 12415.07139 12568.75390

Puis utiliser "SeQ" comme paramètre pour les abscisses dans la fonction "image", et également pour at dans la fonction "axis".


Bonne continuation


Cordialement,

A.D.

[Jesk34]

Merci beaucoup pour votre aide.
Ça marche bien !
Désolé, j'ai mis le temps mais j'y suis arrivée !