科学の箱

科学・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

factorと数値型の変換

ジニ係数にファクターと数値変換についてコメントをいただいたので調べてみた。 その結果以下の3つのことがわかった。 read.csvで桁区切り文字つきの値を読み込むとfactorになる。 区切り文字つき …

no image

rvestで読み込むcssを調べる方法

rvestで対象とするhtmlのブロックは下記のように取得できる。 http://blog.rstudio.com/2014/11/24/rvest-easy-web-scraping-with-r/ …

no image

diag

3×3の対角行列を生成する > diag(3) [,1] [,2] [,3] [1,] 1 0 0 [2,] 0 1 0 [3,] 0 0 1 3×3の対角行列を生成し、対角 …

no image

Rでのランダムデータ作成 – 上級編

平均が異なるランダムデータを作成する > x <- rnorm(12, mean=rep(1:3,each=4), sd=0.2) > y <- rnorm(12, mean= …

no image

Rによるやさしい統計学/5-統計的検定-1

5章で学ぶこと 5章では検定方法について紹介している。以下の方法について検定対象と検定する条件を理解する。 Z検定 t検定 無相関検定 独立性の検定 また検定の手順についてもまとめている。 検定方法 …

2014年4月
月 火 水 木 金 土 日
« 3月   5月 »
 123456
78910111213
14151617181920
21222324252627
282930  

side bar top



アーカイブ

カテゴリー