゙ソフトウェア開発ツールの便利な使い方を紹介。この広告は、90日以上更新していないブログに表示しています。2013-07-09Rのグラフに回帰分析の結果である数式を書き込むR言語データをプロットしたグラフに、推定された予測値や、推定されたパラメーターを数式ごと書き込みたいときはありませんか?そう難しい話では無いのですが、幾つかコツがあるので例を作ってみました。外気温(冷却材は氷でゼロ度だが、室温はもっと高い模様。推定を行う推定を行い、予測値を計算してみましょう。モデルの制約によりpredict関数は使えません。プロットする座標の範囲を決めるグラフを重ねる前に縦横の座標の範囲を決めておかないと、ズレます*1。データをプロットするまずは観測値をプロットしてみます。推定値を重ね書きする推定値をグラフに書き込みます。ここはRの紹介で良く触れられますね。lty=2)するべきだろと突っ込まれました。直線なのが気に入らない人は、plotの代わりに以下のようにcurveを使ってください。推定された式を凡例として表示するまず、数式をsprintf関数で作って、parse関数で表現式に変えます。expression関数は、代入が手間なので使いません。plotmathをすると一覧が見られます。coef(r)[1]))次に、legend関数で凡例を表示します。推定値はlty、観測値はpchになりますが、使わない方には-1を入れていることに注意してください。データの出所を明記するデータを拝借しているので、出所を明記しておきましょう。相関係数などを書いても良かったかも知れません。文字が化ける場合は、以下のようにpar()関数で指定するグラフィックス・パラメーターでフォントを指定してください。フォントはMeiryoを指定していますが、お使いのOSにあわせましょう。family="GraphFont")画面のグラフを保存したい場合は、以下のような感じで。PNGファイルを閉じる*1:推定値のプロットにpar(new=T)とplot()を使うからズレるのであって、低水準関数のlines()を使えば座標の範囲を決めなくてもいいと言うツッコミがありましたが、ここを決めておかないと軸目盛などの位置もズレると言うか、あわせるのが困難になるので、やはり決めておいた方がいいです。ニュースの社会科学的...2015-01-29時系列データの季節調整をしてみよう時系列データには、季節バイアスが入りがちです。»プロフィールid:uncorrelatedソフトウェア開発ツールの便利な使い方を紹介。読者になる検索リンクはてなブログブログをはじめる週刊はてなブログはてなブログPro最新記事Rでガブリエル比較区間を計算する方法クロス集計表のセルごとの有意確率を出してみよう行列を使った計算で、SURの感じを掴んでみよう相関行列から一定以上の大きさの変数の組み合わせを抽出よくある操作変数法(IV)のための練習データとOLSによる推定バイアス月別アーカイブはてなブログをはじめよう!uncorrelatedさんは、はてなブログを使っています。ブログを報告するスターをつけました引用をストックしましたストック一覧を見る閉じる引用するにはまずログインしてくださいログイン閉じる引用をストックできませんでした。再度お試しください閉じる限定公開記事のため引用できません。