🎄teratailクリスマスプレゼントキャンペーン2024🎄』開催中!

\teratail特別グッズやAmazonギフトカード最大2,000円分が当たる!/

詳細はこちら
ggplot2

ggplot2は、統計解析向けのプログラミング言語及びその開発実行環境であるR言語のグラフィック作成用のライブラリです。

R

R言語は、「S言語」をオープンソースとして実装なおした、統計解析向けのプログラミング言語です。 計算がとても速くグラフィックも充実しているため、数値計算に向いています。 文法的には、統計解析部分はS言語を参考にしており、データ処理部分はSchemeの影響を受けています。 世界中の専門家が開発に関わり、日々新しい手法やアルゴリズムが追加されています。

Q&A

解決済

1回答

5096閲覧

Rのggplot2における散布図に回帰式や決定係数、P値などを書き込みたい

tapyon23

総合スコア1

ggplot2

ggplot2は、統計解析向けのプログラミング言語及びその開発実行環境であるR言語のグラフィック作成用のライブラリです。

R

R言語は、「S言語」をオープンソースとして実装なおした、統計解析向けのプログラミング言語です。 計算がとても速くグラフィックも充実しているため、数値計算に向いています。 文法的には、統計解析部分はS言語を参考にしており、データ処理部分はSchemeの影響を受けています。 世界中の専門家が開発に関わり、日々新しい手法やアルゴリズムが追加されています。

0グッド

0クリップ

投稿2020/12/19 23:07

編集2020/12/19 23:43

前提・実現したいこと

Rのggplot2で作成した散布図に、回帰式や決定係数、P値などを表記したい
対処方法などを、ご教授いただけると幸いです。

発生している問題・エラーメッセージ

ggpmisc にて表記したところ、数式などの文字が以下のような状態で表記される。
正常に表記するためにはどうすれば良いか?

italic(y)~'='~6+0.8*~italic(x) ※画像が添付できないので状態を書き込みました などと表記される。 表記したいのは y= 6+0.8x

該当のソースコード

library(ggplot2) library(ggpmisc) library(tidyverse) windowsFonts(Cen = windowsFont("Century")) fit <- lm(M ~ F, data = d) new_mtcars <- cbind(d, predict(fit, interval = "prediction")) head(new_mtcars) g <- ggplot(data = new_mtcars, aes(x = F, y = M)) g + geom_ribbon(aes(ymin = lwr, ymax = upr), alpha = 0.2, fill = 'green') + geom_point() +geom_smooth(method = "lm", formula = y ~ x) + stat_poly_eq(formula = y ~ x,aes(label = paste(..eq.label..,rr.label,p.value.label,"N=", nrow(d))))

試したこと

ggpmiscを読み込んだ時点で、「annotate」 がマスクされていると表示されます。
annotate:: で試してみたが、上手くいかなかった。

補足情報(FW/ツールのバージョンなど)

R4.0.3
パッケージのアップデートは最新にしています

気になる質問をクリップする

クリップした質問は、後からいつでもMYページで確認できます。

またクリップした質問に回答があった際、通知やメールを受け取ることができます。

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

guest

回答1

0

ベストアンサー

stat_poly_eq()とpaste()を使う場合,区切り文字の指定方法が
ややトリッキーになるようです.

cranによくわかる例があります!
Cranによれば,
空白スペースの場合は,sep=""のように(チルダ)を指定するようです.

今回のように,sepを使わずに空白や
文字列で"N="等を挿入したい場合は,以下のように実装可能なようです.

R

1d <- cars 2fit <- lm(speed ~ dist, data = d) 3new_mtcars <- cbind(d, predict(fit, interval = "prediction")) 4g <- ggplot(data = new_mtcars, aes(x = dist, y = speed)) 5g<-g + geom_ribbon(aes(ymin = lwr, ymax = upr), alpha = 0.2, fill = 'green') + 6 geom_point() +geom_smooth(method = "lm", formula = y ~ x) + 7 stat_poly_eq(formula = y ~ x,aes(label = 8 paste(stat(..eq.label..),"*\" \"*", 9 stat(rr.label),"*\" \"*", 10 stat(p.value.label),"*\" N= \"*", 11 nrow(d),"*\" \"", 12 sep="")),parse=T) 13print(g)

投稿2020/12/24 11:53

編集2020/12/24 11:54
yyki

総合スコア38

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

tapyon23

2020/12/29 07:06

ご回答ありがとうございました。 Cranを見ればよかったのですね!しかもご丁寧に記述までして頂き、ありがとうございました。 助かりました。
guest

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

15分調べてもわからないことは
teratailで質問しよう!

ただいまの回答率
85.36%

質問をまとめることで
思考を整理して素早く解決

テンプレート機能で
簡単に質問をまとめる

質問する

関連した質問