元データ分析の会社で働いていた人の四方山話スマートフォン用の表示で見る元データ分析の会社で働いていた人の四方山話この広告は、90日以上更新していないブログに表示しています。線形単回帰について、最小2乗法で確認しました。a,b,cについては、前回の拡張を利用してという形式で求められます。これを解くと、係数のa,b,cが求まるので、当てはめる2次関数が求まるというわけですね。当てはめてみることを考えます。比較の為に、Rの組み込みの非線形回帰の関数nlsを利用しています。data=test.frame)#dev.off()結果です。4.441e-16図ではのようになりました。はい、うまく動いているようです。曲線も絶妙な場所を通っています。ここまでで分かるように、この最小2乗法は一般化できます。(数式を書くのはめんどくさいので割愛)すなわち、任意のn次式の係数を推定できるということですね。どんなに複雑なモデルの当てはめもできてしまいます。それに意味があるかは別ですが・・・しかし、この最小2乗法は、ノイズに異常に弱いので気をつけてる必要があります。いうのも、上記の例で言うと、(4,-10)とかの外れ値が加わると・・・のように、2次関数の向きが逆さまになるくらい影響を受けてしまいます。このような外れ値やノイズへの耐性を持った推定方法がロバスト推定になります。次回は余裕があればロバスト推定を扱いたいところです。昨日のTokyo.Rで、判別分析とかロバスト推定とかの話を聞いて、...2010-01-05関数の時間計測さて、前回C拡張を利用すればRが爆速化するというネタを書きま...2008-11-24Numpyメモ(1):単位行列が欲しいそんなときはidentity()関数を使えばよい。impo...2008-07-04汚いコードさらし。一応コードも載せておく。Entriesはてなブログをはじめよう!rindai87さんは、はてなブログを使っています。ブログを報告するスターをつけました引用をストックしましたストック一覧を見る閉じる引用するにはまずログインしてくださいログイン閉じる引用をストックできませんでした。再度お試しください閉じる限定公開記事のため引用できません。