科学の箱

科学・IT・登山の話題

R

splineを利用する

投稿日:2014年4月17日 更新日:

データが連続していない時には、となりあうデータを使ってその間の値を予測できる。簡単な例では1と2という値があり、この間を一つの値で補完すると1.5となる。補間する値を増やしていけば滑らかな曲線になる。

複数のデータについて補間する方法のひとつがスプライン補間である。他に利用できる方法としてはラグランジェ補間がある。

Rでスプライン補間を実行する際にはspline関数が使える。

spline関数についてはこちらのリンクがヘルプになっている。http://stat.ethz.ch/R-manual/R-patched/library/stats/html/smooth.spline.html

boneデータを使って男女別に補間をしてみる。

plot(spnbmd ~ age, data=bone, col=as.numeric(gender) + 1)
bone.spline.male <- with(subset(bone,gender=="male"), smooth.spline(age, spnbmd,df=12))
bone.spline.female <- with(subset(bone, gender=="female"), smooth.spline(age, spnbmd, df=12))
lines(bone.spline.male, col="blue")
lines(bone.spline.female, col="red2")

spline

メタ情報

inarticle



メタ情報

inarticle



-R
-

執筆者:


  1. […] この図についてこちらで詳細な手順について記載している。 […]

comment

メールアドレスが公開されることはありません。 * が付いている欄は必須項目です

関連記事

no image

Rによるやさしい統計学/6-2つの平均値を比較する

(1) 統計が好きか嫌いかで統計テスト1の得点について有意な差はあるか。 指導法データは score.csvに作成する。 ID,名前,性別,数学,統計,心理学テスト,統計テスト1,統計テスト2,指導法 …

no image

Wilcoxonの符号付順位決定

boneのspnbmdを使って年齢別に男女差があるかを確認してみる。検定ではWilcoxonの符号付順位決定を利用する。 boneは特定個人について骨密度を計測している。データにはIDと骨密度以外に性 …

no image

Rによるやさしい統計学/7-分散分析

(1)学部によりテストの結果に有意な差があるかを調べる。 データはcsv形式で作成し、7_1_test_score.csvとする。以下はデータの内容。 法学部,文学部,理学部,工学部 75,62,65 …

no image

主成分分析

Rで提供されているirisを使って主成分分析を実施する。 主成分分析の手順 データが分離できるか検証する 主成分負荷量と寄与率を求める 新しい変数に名前を付ける 変量プロット 主成分プロット 今回はR …

no image

dataanalysis-002-week3

exploratory analysis グラフの目的 データのおおよそを理解する パターンを探す モデリングを探す デバッグ コミュニケーション boxplot 値のレンジを確認する 他の値とレンジ …

2014年4月
« 3月   5月 »
 123456
78910111213
14151617181920
21222324252627
282930  

side bar top



アーカイブ

カテゴリー