génération de séquence

Bonjour,

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

Code:

---

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

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:

---

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

---

Citation:

---

Envoyé par Jesk34

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.

Citation:

---

Envoyé par 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:

---

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.

Citation:

---

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:

---

tab[,1:200]

---

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

Code:

---

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:

---

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.

Jeu de données :
Pièce jointe 92917

Code:

---

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 )
Pièce jointe 92918

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.

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:

---

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.

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

Citation:

---

Envoyé par Jesk34

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:

---

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.

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.

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)?

Citation:

---

Envoyé par Jesk34

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.

Code:

---

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.

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:

---

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.

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