1

私はデータのセットを持っていますx(1000,2)

2   18.1199158925616
3   9.56472607516328
7   34.3629828698699
4   30.9865760010538
7   54.4869549747429
8   54.7696742935546
2   1.45832407805944
3   15.4231700658276
8   57.0479651194063
1   5.34909555648445
2   11.5593766410445
2   17.4392978010484
6   25.4072607558367
4   15.1248062077068
2   13.2035176450795
7   33.2507069131440
5   21.6859046133920
3   18.4458236149053
1   0.824962982951128
8   46.6322862041701
8   43.6889420657037
2   17.0930689315346
3   13.9754117410459
2   1.99260423380962
11  86.3172411757665
8   50.6244869641866
1   1.16602336459361
2   12.6604166373800
2   17.4519078987572
5   27.4799515521804

plot(x(:,1),x(:,2),'*r');私がスキャッター形式で プロットしたもの:ここに画像の説明を入力してください

normpdfしかし、x(:、1)の選択された値に 以下のようなものを使用したいと思います:ここに画像の説明を入力してください

x(:,2)に対応する列の値を見つけるにはx(:,1)==3

[x3]=x(find(x(:,1)==3),2);
mu=mean(x3);
sig=std(x3);
Y=normpdf(x3,mu,sig);

figure('color',[1,1,1]);
hhh=plot3(3*ones(length(xi)),x(find(x(:,1)==3),2),Y,'-r')
p = polyfit(x(:,1),x(:,2),1);
plin=plot(x(:,1),polyval(p,x(:,1)),'-b');

私のバージョンはパンツに見えます しかし、「-r」でジグザグの線が表示されるので、赤い点線のpdfを滑らかな曲線にするにはどうすればよいですか。何らかの方法でデータを注文する必要がありますか?また、正規曲線の中央にある赤い点を青い線まで下げるにはどうすればよいですか?そしておそらく最も重要なのは、十分な知識のある読者として、このタイプのプロットまたはおそらくもっと派手でない何かを見たいですか?

4

1 に答える 1

1

はい、x-valuesを並べ替えると、赤い曲線が表示されます。また、ここでは必要ありませんfind()。コードをテストしていませんが、これらの変更は機能するはずです。

[x3]=sort(x(x(:,1)==3,2));
mu=mean(x3);
sig=std(x3);
Y=normpdf(x3,mu,sig);

figure('color',[1,1,1]);
hhh=plot3(3*ones(length(xi)),x3,Y,'-r')
p = polyfit(x(:,1),x(:,2),1);
plin=plot(x(:,1),polyval(p,x(:,1)),'-b');

そのようなプロットの価値については、かなり雑然としていると思います。しかし、それはデータの分布の一般的な見方を与えるために使用できると思います。

于 2013-03-26T15:32:39.283 に答える