科学の箱

科学・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 – ElemStatLearn – Coutries

データの概要 国ごとの相違を距離で表しているクロス表 > dim(countries) [1] 12 12 > names(countries) NULL > str(countri …

no image

Rでのラインの表示

Rでラインを引くために関数がいくつか用意されている。 abline() 画像全体に線を引く segments() 始点と終点を指定してその間に線を引く points() “点” …

no image

Rによるやさしい統計学/3-2つの変数の記述統計

散布図を作成 Rで散布図を作成する。データはrunifで10個ずつ作成する。データのレンジは50~100. x_test <- as.integer(runif(10, 50, 100)) y_ …

no image

{psych}

psychパッケージは性格分析、精神分析、心理統計に役立つ関数をまとめたパッケージになる。 モジュールとしては以下のように分かれている。 データ入力および記述統計 データクレンジングおよびクラスタリン …

no image

クラスタリング

Rによるクラスタ分析の実行 県別の持家率と自家用車普及率からクラスタ分析を実行し似た県を探してみる。。 データソース 持家と自家用車普及率のデータは下記から取得する 持家普及率のデータはこちらから取得 …

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

side bar top



アーカイブ

カテゴリー