Rで時系列スペクトル分析 - nokunoの日記
id:yokkunsさんが主催する第4回R勉強会のid:hamadakoichiさんの発表で議論になったのですが、スペクトル分析（フーリエ変換）において0に近い低周波の成分が出るときは、全体のトレンドを表しています、という話。 データはRに付属のUKgasを使って、spectrum関数にかけました。まずは、普通にデータ点とそのスペクトル、そしてARモデルで平滑化したスペクトルを見てみます。以下のように、低周波付近と1付近にピークがあることがわかります。 > plot(UKgas) > plot(spectrum(UKgas)) > plot(spectrum(UKgas),method="ar") そして、diff関数を使って差分を取ると、以下のようにトレンドが取り除かれ、周波数領域で見ると低周波のピークがなくなることが分かります。残ったピークは、1年周期なので季節変動を表します。 > p
R-Source
group1 <- c(0.7,-1.6,-0.2,-1.2,-0.1,3.4,3.7,0.8,0.0,2.0) # グループ 1 の睡眠時間の増加 group2 <- c(1.9, 0.8, 1.1, 0.1,-0.1,4.4,5.5,1.6,4.6,3.4) # グループ 2 の睡眠時間の増加 t 検定を行う関数 t.test() を用いてデータ sleep に関する検定を行う手順を紹介する． 検定例(1) ： 一標本 t 検定 関数 t.test() で 1 標本 t 検定を行う場合の書式と引数を挙げる．一番目の引数にデータを与えればとりあえず検定が出来る．全ての引数を指定する必要は無い．を指定すれば検定が行われる．
Rコマンダーのインストール
Rコマンダーという便利なパッケージがあります。これを使えば、コマンド型のソフトを使うのが苦手な方も気楽にRを使えるようになると思います。ぜひぜひインストールしてみましょう。 Last modified: Thu, 23 Apr 2009 10:25:41 +0900 Rの追加パッケージ Rをインストールしては見たものの、コマンドを自分で打ち込んで使うことに慣れていない人はすぐに難しいと考えてしまいがちのようです（実際にはそんなことはないのですが…）。しかし、心配は無用です。Rにはたくさん追加パッケージが用意されていて、これらを追加することでより便利に使うことができるようになっています。今回は、そのパッケージの中でも特に最初に使う人に便利だと考えられるRコマンダーのインストールについて説明します。 パッケージのインストール法 インストールの際にはインターネットへの接続が必要です。ご注意くださ
1
リリース、障害情報などのサービスのお知らせ
最新の人気エントリーの配信
処理を実行中です
j次のブックマーク
k前のブックマーク
lあとで読む
eコメント一覧を開く
oページを開く
R-Source
Rのお仕事:信頼区間 - livedoor Blog（ブログ）
