科学の箱

科学・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. […] この図についてこちらで詳細な手順について記載している。 […]

Wilcoxonの符号付順位決定 | 科学の箱 へ返信する コメントをキャンセル

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

関連記事

no image

data.frameについて注意すること

以前predict()について検証を進めていた際に、不可解なエラーが出てしまい、途中で調査をやめてしまったことがあった。 具体的なエラーは以下である。 > newGalton$child &lt …

no image

Rオブジェクト指向のメモ

whichではwhich.minとwhich.maxが用意されているという指摘を受けたので調べ見てた。 which.minとwhich.maxはRでいうところのオブジェクト指向を利用した構造にみえる。 …

no image

conjoint分析の資料

マニュアル conjointパッケージ caFactorialDesign Rでconjointパッケージを利用した例 Rでコンジョイント分析 Rでコンジョイント分析 |極めて個人的なメモ コンジョイ …

no image

ジニ係数

計量経済学の本を読んでいたら、ジニ係数の利用について記述されていたのでRで実データを使って計算してみる。 経済学の勉強でジニ係数というのは知っていたが、その利用方法としては所得の分配がどの程度平等であ …

no image

データフレームで層別に操作をする

Rでデータベースに似たデータを取り扱う際には、データフレームを使うことがおおい。JavaはPHPなどのプログラミング言語では、層別に集計作業をする際にはデータベース上のSQLか、データを読み込んでから …

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

side bar top



アーカイブ

カテゴリー