Tracé de courbe + Régression linéaire

Bonsoir à tous, j'ai un problème que je n'arrive pas à résoudre.

J'ai une colonne de valeurs, appelée C, du type:

Code:

---

1 2	C'= 1 2 4 2 5 1 3 6 4 5 5 3

---

Je veux pouvoir tracer une courbe qui serait une espèce d'histogramme de cette colonne.

Mettons que j'ai y fois le nombre x. Je voudrais pouvoir placer le point de coordonnées x en abscisse et y en ordonnée (nombre d'occurrences de la valeur x), pour avoir la courbe y=f(x).

Avec un histogramme, donc, c'est faisable rapidement (moyennant quelques ajouts à hist pour rendre la chose exploitable avec les bandes), mais un plot 8O

J'ai longuement réfléchi au problème, pour stocker le nombre d'occurrences des différentes valeurs etc., mais je dois dire que là, je ne sais du tout comment faire :cry:
Et je n'ai pas trouvé de commande du genre histplot :D

Pourriez vous m'aider? Merci d'avance!

bonus: et une courbe logarithmique en log(nombre d'occurrences de x)=f(log(valeur x)), y a t-il moyen de le faire?

Salut,

Essaie:

Code:

---

1 2 3 4

x=round(10*rand(1,100)); [n,out]=hist(x,unique(x)); plot(out,n)

---

pour l'echelle logarithmique, il y a la fonction loglog:

Code:

---

loglog(out,n)

---

Bonus: 8-) pour les semilogarithmique: semilogx et semilogy.

Superbe, ça marche! Merci, merci, et merci! :P

Et allez, j'ose le demander:

Je tombe donc sur une certaine courbe qui a l'apparence d'une droite. Comment faire pour tracer la droite qui approche au mieux la courbe? Je pense qu'il doit y avoir une commande puissante toute faite, mais je ne la connais pas. Merci encore? :)

Tu peux utiliser polyfit.

[EDIT]
ou bien la division a gauche, exemple:

Code:

---

1 2 3 4 5 6 7 8

x=1:5; y=5*x+2+randn(1,5); c=[x' ones(5,1)]\y'; plot(x,y,'*') hold all yy=c(1).*x+c(2); plot(x,yy)

---

Citation:

---

Envoyé par AC88mm

Superbe, ça marche! Merci, merci, et merci! :P

Et allez, j'ose le demander:

Je tombe donc sur une certaine courbe qui a l'apparence d'une droite. Comment faire pour tracer la droite qui approche au mieux la courbe? Je pense qu'il doit y avoir une commande puissante toute faite, mais je ne la connais pas. Merci encore? :)

---

Si x et y sont des vecteurs de taille 1xn et que y=a*x+b+e avec a et b constant et e alors on trouve a et b minimisant e en faisant le code suivant :

Code:

---

1 2 3 4 5 6

A=[x(:) ones(size(x(:)))]; B=y(:); X=A\B a=X(1) b=X(2)

---