R (プログラミング言語)

はじめに R でグラフを描画していたのですが、legend（凡例）だけ出力したい、という場面がありました。 検索ワードが悪かったようで、答えに辿り着くのに少し時間がかかったため、メモとして残しておきます。 方法 以下の Stack Overflow の方法をそのまま使いました。 ggplot2 - How can I create a legend without a plot in R? - Stack Overflow 空のグラフを描画して、その凡例として出力したい凡例を設定すればいいようです。 コード例： 出力： 出力された画像 このままだと余白や画像サイズが不自然なので、png() や…

本稿ではRを用いてグループ内平均を計算する方法をご紹介いたします。 パッケージの読み込み 使用するデータセット グループ内平均の計算 最後に パッケージの読み込み library(dplyr) インストールされていない方はinstall.packages("dplyr")をご利用下さい。 使用するデータセット 今回は、男女2人ずつのテストのスコア("Tesy_score")が含まれるデータセットを作成しました。 男性("male")と女性("female")ごとの平均を計算していきましょう。 data <- data.frame(name = c("Taro", "Jiro", "Hanako…

書籍の紹介 2023年9月13日に、技術評論社より『数値シミュレーションで読み解く統計のしくみ 〜Rでためしてわかる心理統計』という書籍を上梓します。 gihyo.jp 本書は、確率分布や、確率分布が関係する様々な定理、帰無仮説検定や信頼区間、そしてサンプルサイズ設計（研究実施前に、取得すべきデータのサイズを見積もる研究実践）など様々な内容を、プログラミング言語Rを用いたシミュレーションにより理解することを目指した本です。 僕は3章（乱数生成シミュレーションの基礎）・4章（母数の推定のシミュレーション）の執筆を担当しました。 「ベルヌーイ分布・二項分布・正規分布・多変量正規分布・χ2乗分布・t…

はじめに ggpubrで箱ヒゲ図とバイオリンプロットを書いてみます。 私としては箱ひげ図の方が好きです。バイオリンプロットはなんとなく、KDEでの推定があっているのかわからない時があって、誤解を与えそうな気がするからです。 データの紹介 今回使うのはToothGrowth というデータセットです。 モルモットにビタミンCを投与した際の歯の長さについてのデータセットだそうです。 投与方法が2つあって、オレンジジュースとして投与したデータと、精製したアスコルビン酸で投与したデータがありますが、今回はまとめて表示しています。 処理と結果 箱ヒゲ図 まずは何も考えずに箱ヒゲ図を作ってみます。 data…

はじめに 最近少しRを使うようになりました。 これまでIgor ProかPythonを使ってグラフを描いていたのですが、Igorは有料、PythonのMatplotlibは使い勝手がイマイチということで、晴れてRデビューです。 ちょっと調べるとRにはggplot2というのがあり、こちらで綺麗にグラフが書けるということでしたが、面倒くさがりな私にはもっとお似合いのggpubrというのがありました。 rpkgs.datanovia.com ggpubrはggplot2ベースの、簡単にPublication qualityのグラフが作れるモジュールのようです。試しに対応のあるデータで遊んでみました。…

何となくで使ってきたロジスティック回帰モデルを勉強しなおした。 ロジスティック回帰モデルとは ざっくり言うと、ある出来事が「起こる」、「起こらない」ことを予測する方法で、また、起こる/起こらないは、二値分類でも多クラス分類にも使える。 この確率の計算は、データをロジット（logit）関数に渡し、特徴量と重みパラメータを計算して、ロジスティック関数に当てはめている。また、ロジットを対数オッズ（log-odds）、ロジスティック関数をシグモイド関数と呼ぶこともある。 このとき、ロジスティック関数は、起こる/起こらないを0から1の確率で示し、0.5以上なら「起こる」と0.5未満なら「起こらない」と予…

公開せずに下書きに残っていたので、 TokyoRのLT ライフログとしてLTを記録。 ライフログ tokyor.connpass.com 発表者 タイトル 資料 LT1: @bob3bob3 Tokyo.Rも100回目なので、そろそろPythonと仲良くしたい。 LT2: @kazutan Rに管理されてみる LT3: @eitsupi dtplyr鬼はええ！このままdplyrパイプライン全部data.tableで計算しようぜ！ LT4: @zaiya グラフで可視化する話 LT5: @dropout009 回帰分析ではlm()ではなくestimatr::lm_robust()を使おう LT6…

はじめに 無料公開データセット10選 おわりに 【最終更新日：2023年7月10日】 はじめに これらのデータセットは、PowerBIを使用したビジュアライゼーション、レポート作成、分析の練習に最適です。また、各データセットの利用規約とプライバシーポリシーを確認し、適切に使用してください。 無料公開データセット10選 Available datasets • Rdatasets このリンク先では、Rプログラミング言語用の大量のデータセットが提供されています。データ分析の練習やテストに非常に適しています。 オリジナル統計データ体験版｜ダウンロードサービス｜日本統計センター 日本の各種統計情報を網…