「統計学関連なんでもあり」の過去ログ--- 048

No.22584　Dataを部分的に取り出す方法　　【赤岳】　2018/08/07(Tue) 22:31

いつも勉強させていただいています。
Rの基本操作だと思いますが，わかりません。
今，次のようなデータが1000個あるとします。
x<- rnorm(1000, 15.5, 3.3)

このデータの95％区間にあるデータを取り出したいのですが，どのように取り出せばよいのでしょうか。
たとえば，
> y<- quantile(x, c(0.025, 0.975))
とすると，yのlength=2となり，950個のデータが取り出せません。
どなたか，yが950個のベクトルになるようなコマンドを教えてもらますでしょうか。

No.22585　Re: Dataを部分的に取り出す方法　　【青木繁伸】　2018/08/08(Wed) 05:48

x <- rnorm(1000, 15.5, 3.3) のベクトル x から，10 以上， 15 以下のデータを取り出すときはどのようにしますか？
x[x >= 10 & x <= 15] ですね

では y <- quantile(x, c(0.025, 0.975)) から，下限と上限を取り出すのはどのようにしますか？
y[1], y[2] ですね

それをあわせれば
x[x >= y[1] & x <= y[2]]
です

No.22586　Re: Dataを部分的に取り出す方法　　【赤岳】　2018/08/08(Wed) 08:26

青木先生，
ありがとうございました。
ここ2日くらいそのような関数というかコマンドばかり探していました。
今後ともよろしくお願いします。