3

私は R の初心者です。素朴な質問で申し訳ありません。

私のデータの一部は次のとおりです。

POS         Obs_FCT
4345     0.10049191
4484     0.08930793
4515    -0.00468725
4547    -0.00749802
4628     0.01143598
5347     0.05673895
6071     0.01143598
6449    -0.00070778
6498    -0.0109962
7320     0.00543984
7457    -0.00206247
7761     0.08018589
7875    -0.00601304
7988    -0.00070778
9459     0.01447144
9483     0.08269899
9495    -0.03353229
9552    -0.00206247
9602    -0.0269909
9701    -0.00206247
11809   -0.01952417
12593   -0.02173826
16787    0.08930793
17049   -0.06738125
17058    0.01325792

Obs_FCT(Y 軸) VS (X 軸)をプロットしたいのですPOSが、まったくプロットされていません。また、列にスライディング ウィンドウを適用したいと考えていPOSます。ウィンドウは次のようになります (つまり、ウィンドウ サイズは 1000 で、200 ずつスライドします)。

1--1000
200 ---1200
400 ---1400 
600 ----1600
--------

列の数値がPOS各ウィンドウにある場合、 の平均値を計算しますObs_FCT。(X 軸はウィンドウの中点を使用します)。

これを達成するためにRコードを使用する方法を教えてください。一般的に、私は自分が望むものを手に入れることができるかもしれないことを知っていrollapplyます。しかし、機能が必要なようです。

4

1 に答える 1

1

データをテキスト ファイルにコピーペーストし、read.table を使用して R に読み込みました。plot コマンドを使用すると、必要な結果が得られますplot

> d<-read.table("path/to/file.txt", header=TRUE, sep="")
> head(d)
   POS     Obs_FCT
1 4345  0.10049191
2 4484  0.08930793
3 4515 -0.00468725
4 4547 -0.00749802
5 4628  0.01143598
6 5347  0.05673895
>plot(d)

これにより、X 軸に POS、Y 軸に Obs_FCT を持つプロットが作成されます。

パッケージ caTools の runmean を使用すると、実行中の中央値を作成できます。ウィンドウ内のポイント数は指定できますが、ウィンドウのスライド方法は指定できません。

> md <- runmean(x=d$Obs_FCT, k=1000)

その後、次のようにプロットできます。

plot(d$POS, md)
于 2013-03-19T12:29:40.307 に答える