[r] タグが付いた質問
Rはフリーかつオープンソースなプログラミング言語とソフトウェア環境で、統計解析、生物情報学、グラフ描画などに用いられます。質問の際は再現性のある最小限の例を添えていただけると助かります。
316
質問
2
票
2
回答
125
閲覧数
modelsummary関数で数値を右に並べたい
標準誤差を推定値の下ではなく、右に表示させたいです。
画像では推定値の下に標準誤差が出力されています。
試したコード: rstudioで実行し、バージョンはR version 4.1.3 (2022-03-10)になります。
library(modelsummary)
model_summary <- modelsummary(model, output = "default&...
0
票
0
回答
31
閲覧数
fastLinkを用いたProbabilistic data Linkageについて
fastLinkパッケージを用いて、二つのデータセットのリンクを試みています。リンクに用いる項目日付、時間、年齢、性別で全て数値に変換しています(時間:2/11→0211、11:15→1115、男性矢印、以下のようなスクリプトを作りました。
result <- fastLink(dfA = data_2, dfB = data_1, varnames = c("date",...
0
票
0
回答
60
閲覧数
行列の逆行列が計算できない(Rの系統樹の処理)
系統樹は自分で以下のように作成しています。
どうやってもエラーが解消されなくて困っています。作成方法がいけないのでしょうか??
> newick_tree<-"(A:4,(B:2,(C:2.953469601,D:1.443532938):1.136357167):1.104667235);"
> phylo_tree <- read.tree(text ...
1
票
2
回答
131
閲覧数
R言語でデータを年代別にし、箱ひげ図を作成したい
データは年齢と幸福度になります。
年代別の箱ひげ図をRで作成したいです。
私の想定としては、以下の二つが必要だと思っています。
・年齢を20~29などに設定し、まず年代別にする。
・箱ひげ図をそれぞれの年代でプロットするようにする。
#hensuu1 <- age, hensuu2 <- kouhukudo
ggplot(data, aes(x = hensuu1, y = ...
0
票
0
回答
100
閲覧数
アーカイブから"mvpart"パッケージをダウンロードできません
いくつかのドキュメントを参考にしましたが、アーカイブから"mvpart"パッケージをダウンロードしようとすると下記の通りエラーになってしまいます。
両方ともunknown type name 'Sint'が原因のようですが、何のことだか見当がつきません。
ご教示いただけますと幸いです。
どうぞよろしくお願いいたします。
試行したコード①
install.packages(&...
0
票
0
回答
49
閲覧数
Rselenium でブラウザが起動しない
Rseleniumを使用したいのですがwebブラウザを開く段階でエラーが出てしまいました。
エラーメッセージ
Selenium message:Unable to create new service: ChromeDriverService
Build info: version: '3.141.59', revision: 'e82be7d358', time: '2018-11-14T08:...
0
票
0
回答
96
閲覧数
Rでのwebスクレイピング
Rでwebスクレイピングしたいと考えています。
抜き取りたい情報はM-1グランプリの情報で
(https://www.m-1gp.com/schedule/detail.html?id=539)におけるグループ名、コンビ名、結果の情報とコンビ名をクリックした先にある結成年、所属、前年度における結果等を抜き出そうと考えています。
初めのグループ名、コンビ名、...
0
票
2
回答
205
閲覧数
Rでのwebスクレイピングで、データが抽出できない
Rでwebスクレイピングを行っているのですが, 抽出したいデータが変数内に入りません.
抽出したいのは, テーブル型のデータです.
サイトは, 東京都教育委員会の都立高校検索のサイトで、アクセスのところの駅・路線・所要時間のデータを抽出したいと考えております.
使用したライブラリ:rvest
Rのバージョン:R 4.2.1
使用したOS:macOS
以下にコードを添付します.
library(...
0
票
1
回答
93
閲覧数
ロジスティック回帰の最尤推定におけるヘッセ行列の計算ができない
ロジスティック回帰の最尤推定におけるヘッセ行列の計算についてご教示ください。
以下、R言語での記述ですが、
hessian <- function(beta, X) {
p <- exp(X %*% beta) / (1 + exp(X %*% beta))
W <- diag(p * (1 - p))
return(t(X) %*% W %*% X) # ...
0
票
0
回答
164
閲覧数
Rでnnetがダウンロードできません
install.packages("nnet")
Warning in install.packages :
リポジトリー https://cran.ism.ac.jp/src/contrib に対する索引にアクセスできません :
cannot open URL 'https://cran.ism.ac.jp/src/contrib/PACKAGES'
WARNING:...
1
票
1
回答
91
閲覧数
Rで0.17 * 10000 == 1700がFALSEになる理由と解消法
タイトルの通りなのですが,Rで0.17 * 10000 == 1700の真偽を判定させると,FALSEが返ってきます。以下はTRUEと判定されます(当然ですが…)。
1700 * 1 == 1700
170 * 10 == 1700
17 * 100 == 1700
1.7 * 1000 == 1700
Rが正しく計算できていないように見えるにははなぜでしょうか?また解消法はあるでしょうか?
追記:...
0
票
0
回答
43
閲覧数
Rで傾向スコアマッチングがうまくいきません
matched_data<-matchit(ope_timing~age+height+sex+weight+flag_anticoag_antiplatelet+flag_pulse+flag_icu+flag_hcu+BarthelIndex_adm_cat+cci_quanindex+jcsadm_1, data = H_pre, method = "nearest"...
0
票
1
回答
89
閲覧数
RでNAを含む連続変数がカテゴリー変数と判断されます
Rで-100から∔100までの連続変数(ただし計算上は±5毎に区切られます)の平均、SD値を、条件A/B群でそれぞれ比較し、TableOneで解析をしています。
この連続変数(BI_dif)の中にNAが多数あるからでしょうが、TableOneを作成すると、5毎のカテゴリー変数別に割合が表示され、平均・SD値が表示できません。
pt_01_excel %>% replace_na((list(...
0
票
1
回答
179
閲覧数
Rをバッチファイルのみで実行する方法
Rをインストール済みのWindowsPCでバッチファイルからRを実行したいと考えております。
2行の短いプログラムなので、プログラムファイルは別に置かずに、バッチファイル内にコードを記述する方法を探していますが見つかりません。
(バッチファイルからRプログラムを実行する方法は見つかるのですが。)
対応できる方法はありますでしょうか。ご教示いただけますと幸いです。
どうぞよろしくお願いいたします。
1
票
2
回答
154
閲覧数
PythonまたはRでクロス集計をしたい
PythonまたはRで以下のCSVファイルを読み込んで、クロス集計を行いたいのですが、うまくいきません。
男性┗肉のようにクロス軸を二次元にする事は出来るのですが、一次元にする事が出来ません。
CSVデータイメージ
クロス集計イメージ
ChatGPTで質問したところ以下の回答でしたが、主力イメージが意図した結果になりませんでした。
ご教示いただけますと幸いです。宜しくお願いいたします。
# ...
0
票
1
回答
65
閲覧数
forで回して作成した複数のヒストグラムにサブタイトルを追加したい
下記のスクリプトだと、サブタイトルに subtitle_file[i, 2] のすべてが結合された状態で追加されてしまいます。私は、ヒストグラムが作成されるたびに、 subtitle_file[i, 2] の1つが順番に追加されるスクリプトを作成したいです。
ヒストグラム作成に使用するファイル数とID_Name.txtの行数は一致しています。
# サブタイトルを格納するための変数を初期化
...
0
票
0
回答
53
閲覧数
Missing required option: '--preset <name>'というエラーで困っています。
”Platypus”というRのプラットホームでレパトア解析を行いたいと思っています。まず開発元の Vignettes を用いて使い方を確認しているのですが、7.3のMIXCRを使うところで"Missing required option: '--preset '"のメッセージがでてきます。
どなたかわかる方がいましたらご教示いただけるとありがたいです。よろしくお願い致します。
...
1
票
2
回答
166
閲覧数
Boxplotを条件ごとに作成する。
ggplot2を用いて、経時的なサンプルの変化を箱ひげ図で示したいと思っています。
以下のような手順で進めると、添付画像のように横軸のDayが反映されない箱ひげ図が出力されてしまいます。
Dayごとに箱ひげ図をまとめるには、どのようにすれば良いでしょうか?
Boxplotに詳しい方、コメントを頂けないでしょうか?
library(tidyverse)
library(easystats)
#...
0
票
1
回答
99
閲覧数
Rstudio起動時に前回エラーを起こした作業スペースが読み込まれて作業ができない
RStudioでコードを書いていたところ、以下のエラーが表示され、その後スクリプトに書き込みはできるもの計算ができなくなりました。
fatal error; unexpected exception: string too long
ファイルを消去すれば解決するだろうと思い、エラーが出たファイルを削除してRStudioを再び起動すると、...
0
票
1
回答
65
閲覧数
条件A、条件Bの数値を元に条件Cを除したい
Rでの解析において、条件A、条件Bの数値を元に、条件Cを除するプログラムを書きたいのですが、うまくいきません。
(例)
Sample = c("A", "A", "A", "B", "B", "B")
Day = c(0, 1, 2, 0, 1, 2)
Value = c(100,...
0
票
1
回答
169
閲覧数
Rにてtidyverseが読み込めません
version R version 3.5.3 (2019-03-11)
os Windows 10 x64 (build 19045)
system x86_64, mingw32
ui RStudio
rstudio 2022.07.2+576 Spotted Wakerobin (desktop)
以上の条件下で実行したいと思ってますが、...
1
票
1
回答
139
閲覧数
Rでlmeを用いた解析をする際のエラー;invalid type (list) for variable 'data'
Rでlmeモデルを使った解析をしようとしていますが、以下のように解析しようとするとエラーが出てしまいます。原因を調べてみましたがわかりません。
データの形はdataframeであり、列の型はcharacterとintが混在しています。
コード:
data <- read.table("data.txt",header=T,sep="\t")
...
0
票
1
回答
170
閲覧数
Rを使用して列の合計値を計算、各列の合計値が100以上の列のみを抽出する方法
Rを使用して、列の合計値を計算、列の合計値が100以上の列のみを抽出してデータフレームにする方法を教えていただけますでしょうか。
df <- data.frame(id=1:3, score1 = c(200, 30, 100), score2=c(0, 30, 10), score3=c(90, 40, 10), score4=c(200, 220, 300))
...
0
票
1
回答
71
閲覧数
グループ列があるデータセットの相関行列の作成
以下のような縦積みのデータセットでgroup別に相関行列をggpairs関数で作成したいのですが、ご教示いただけますと幸いです。
ID value group
1 2 A
1 4 A
1 1 B
1 7 B
2 3 A
2 9 A
2 7 B
2 2 B
0
票
1
回答
120
閲覧数
1つの列の2つのグループ(男女)から四分位の変数を作成→新しい1つの変数を作成→正確に分類できているかの確認まで教えてください。
下記のようなデータフレームで、男女別にBPの四分位変数(BPq)を作成し、それぞれの下位25%を新変数(lowBP)= 1(その他は0)として作成したいです。
df <- data.frame(
ID = 1:4, age = c(43, 62, 54, 55), sex = c(0, 1, 1, 0), group = c(1, 2, 3, 3)
BP = c(120, ...
0
票
1
回答
164
閲覧数
RStudioで Markdown の出力時、警告が出て止まってしまう
RStidioでパッケージをロードすると、マークダウンの出力画面では警告が出て(しかもその警告の文が途中までで)、以後の処理は行われず止まってしまいます。
左上画面でチャンクごとに実行したところ、以下の警告が表示されますが、以後の処理もチャンクごとに実行することはできています。
Warning: package ‘ggplot2’ was built under R version 4.1.3
...
1
票
1
回答
54
閲覧数
2条件のグループ別の平均値
以下のデータフレームに対してgroup1とgroup2それぞれの要約統計量(平均、最大値、最小値など)を出力したいです。ご教示いただけますと幸いです。
例)
MenのAグループの平均値=
MenのBグループの平均値=
WomenのAグループの平均値=
WomenのBグループの平均値=
df <- data.frame(age=1:100, group1=c(rep("Men&...
1
票
1
回答
267
閲覧数
複数の列の条件から新しい列を作成する方法
sexが1、かつ、groupが3の場合に、新しい変数として「category」=1、その他(sex=0, group=1,2)は0と作成したいです。ご教示お願い致します。
df <- data.frame(
ID = 1:4, age = c(43, 62, 54, 55), sex = c(0, 1, 1, 0), group = c(1, 2, 3, 3)
BP = c(...
1
票
1
回答
105
閲覧数
Rstudioでグループ別のglmを教えてください
下記のデータフレーム(df)で
glm(BP ~ age + sex + BMI, data = df) を実施しています。
データフレームをそのままでグループ別(sex)に分析する方法はありますか?
対象のデータフレーム:
ID age sex BP BMI
1 43 0 120 21
2 62 1 130 26
3 54 1 132 23
4 55 0 110 19
0
票
1
回答
260
閲覧数
Rを使用して列の値を他の列にあるNAに置き換えたい
Rを使用してデータをクリーニングしています。
下記に3つの列(ID, s1, s2)があり、s1には2つの欠測があります。
s1に欠測がある場合、s2の同じ行をNAに置き換えた列を作成したいのですが、どなたかおしえていただけますか?
ID s1 s2
1 2 1
2 3 1
3 2 3
4 NA 2 ← ここにs1のNAを入れたい
5 NA 2 ← ...
2
票
0
回答
88
閲覧数
R ggplot2で100%積み上げグラフ
ある疾患の届出数の推移(2010~2021年)を患者の年齢階層別に100%積み上げグラフで作成したいと思っています。
以下のコードでグラフを作成したのですが、X軸の年が2010から4年間隔で2022年まで表示されてしまいます(2010、2013、・・・2022)。この2022を消すか、最後のみ2021を表示させるためにはどうしたらよいでしょうか。
よろしくお願いします。
Fig3dat<-(...
1
票
1
回答
266
閲覧数
Rで文字列の日付データをDate型に変換したい
文字列の日付データをDate型に変換したいのですが、以下のようなエラーが出ました。
エラーコードの指している事がどういう事なのかがわからずに困っています。
実行コード
asp <- asp %>%
select(date,genre_name,amount,device_type)
asp$date <-
as.Date(asp$date)
エラーコード
...
1
票
1
回答
90
閲覧数
Rのleft_joinでキーが~未満で一番大きい値になったものをjoinさせる方法
Rにおいて、以下のようなdfがあったとします。
AAA 200903
AAA 201003
BBB 200903
BBB 201003
これにleft_joinで以下のdfを右側に結合させ、
AAA 200901 1
AAA 200902 2
AAA 201003 3
AAA 201004 4
AAA 201005 5
BBB 200902 6
BBB 200912 7
BBB 201011 8
...
0
票
1
回答
447
閲覧数
R でhist(x) 関数でヒストグラムを描けない( 'x' は数値でなければなりません とエラーを返す)
Rで,以下のようなコードで 関数を作り,その結果をヒストグラムにしたいのですが, 'x' は数値でなければなりません とエラーを返してきてうまくいきません. 助けてください.
func01 = function(n){
sample1 <- 100
df1 = data.frame()
for(i in 1:n){
sample1 <- sample1+1000
df1 <- ...
0
票
1
回答
171
閲覧数
Rのクロス集計から条件をつけて削除
Rでクロス集計(table(X,Y))をした結果、例えばX*Yが10未満だった行を削除したいのですが、一気にやる方法はありますでしょうか?
0
票
0
回答
2,112
閲覧数
Rのパッケージインストール時のエラー
Rstudioでパッケージを何かインストールしようとすると毎回
Error in install.packages : error reading from connection
というエラーが表示されます。
何度かR, R studio, rtoolsをダウンロードしなおしました。
原因不明です・・既にほかのアプリも入れているPCなので、あまりクリーニングはしたくありません。
Rで、...
0
票
1
回答
81
閲覧数
R言語で土日祝日を除いた日数の計算結果をデータフレームへ
Rで土日祝日を除いた開始日と終了日の間の日数を計算しようとしましたが、上手くいきません。
ベクトルの型が上手くいかないものと思いますが、何か上手くやる方法はありませんでしょうか。
library(dplyr)
p1 <- c(1:2)
p2 <- as.Date(c("2013-01-01","2013-01-22"))
p3 <- ...
0
票
1
回答
584
閲覧数
Rの色々なエラーで悩まされています
rglのライブラリを読み込もうとすると下記のエラーが出ます。
既にxfunがwin-libraryにもRのlibraryフォルダの中にあることも確認済みだし、何回かxfunとrglの再インストールもしていて、原因がわかりません。
rdsのインストールも一応してあります。
Error: package or namespace load failed for ‘rgl’ in library....
1
票
1
回答
64
閲覧数
R 文字型データから日時データへの変換
"2022.04.30 午前 5:55"
のような文字列データをRの日時データに変換したいです。
文字列を地道に置換していけば出来るような気もしなくはないですが、スマートなやり方をご存知でしたら教えて頂きたく。
0
票
1
回答
809
閲覧数
Rでパッケージのインストールでエラーが表示されてしまう
Rx64 4.1.2でSteel-Dwass検定をする際、
「パッケージ」→「パッケージのインストール」→「Japan(Tokyo)」→「NSM3」
の指示通りにするが、Japanを選択以降、以下のメッセージが表示されてしまい、NSM3を選択出来ません。
utils:::menuInstallPkgs()
--- このセッションで使うために、CRAN のミラーサイトを選んでください ---
...
0
票
1
回答
165
閲覧数
rpy2に関して、Rのディレクトリの場所をどのように特定しているのでしょうか
<環境>
os:windows10
実行環境:anaconda3 - jupyter notebook
主要ライブラリ:rpy2 3.7
〈本文〉
pythonでRを使用するためのrpy2というライブラリについて質問です。python環境下でrpy2を利用してRを使用することに成功しました。
一つ疑問なのですが、タイトルにもあるとおり、...
0
票
2
回答
147
閲覧数
R ヒストグラム x軸が名義変数の場合の作図方法を知りたいです。
Rを使って、ヒストグラムを描きたいです。
対象とする変数は、数値の大小に意味がある順序変数ではなく、「名義変数」に分類される業種コード(ex.製造業:01,金融:19,食品:38など)です。
最終的に業種コードごとに度数が表示される図にしたいのですが、現状としてx軸が順序変数としてみなされているため、業種コードのなかでとんでいる(もともと存在しない)番号が数値として残ってしまう状態です。
...
0
票
1
回答
259
閲覧数
fixest::etable()を使ってクロスリファレンスをする方法はありますか?
英語版でも同様の質問 を投稿しています。
質問
fixest::etable() を使ってクロスリファレンスをする方法はありますか?etable は,簡易的なメソッドですが,とても便利です.もしできる方法があるととても助かります.
環境は,Windows 10で,R 4.1.2を使っています.RStudioでbookdown::render_book()をしています.ファイル構成と内容は以下です.
...
0
票
2
回答
1,595
閲覧数
R エラーメッセージ 「`.rowNamesDF<-`(x, value = value) でエラー: 重複した 'row.names' は許されません」
使用するデータは以下のようなパネルデータ(unbalanced_pdf)です。pdata.frame()で、plm用の形式に変換しています。(一番左列: 企業id-年)
////////name price PBR
1-2009 会社a 1890 0.8
1-2010 会社a 1890 0.9
1-...
1
票
1
回答
99
閲覧数
データフレームからあるベクトルの要素を含む行だけを抜き出し、新しいデータフレームを作成したい
使用するデータフレームは以下のような形です。
name year price
A社 1999 200
A社 2000 202
A社 2001 199
A社 2002 400
A社 2003 207
B社 1999 300
B社 2000 500
...
1
票
1
回答
661
閲覧数
RMeCabが日本語ではうまくいきません。
質問:日本語に対してRMecCabをうまく動かせるようにしたい。
RMeCabパッケージの関数を日本語のテキストデータに対して実行すると以下のようになり失敗します。
極楽not foundとなっていますが、実際には極楽は存在します。ちなみに英語のみのテキストデータに対してはうまく動きます。
> library("RMeCab")
> r <- ...
0
票
2
回答
84
閲覧数
完全パネルデータの作成にあたり、観測期間(年次データ)のすべてが揃った個体のデータのみを、新しいデータフレームに保存したいです
使用するデータは多少簡略化していますが、以下のようなものです。
name year ...
A社 1999 ...
A社 2000 ...
A社 2001 ...
A社 2002 ...
B社 1999 ...
B社 2000 ...
C社 2001 ...
C社 2002 ...
D社 ...
1
票
2
回答
297
閲覧数
同じ構成のfor文でも、処理数を増やすとエラーが出るようになる
3つ目のfor文がなぜエラーになるのか、わかる方いらっしゃいますか?
下記のコードは各オブジェクト作成からfor文の処理まで一部始終のコードと、途中からの出力結果となります。
とてもシンプルに書いているつもりなのですが、同じ構成の3つ目のfor文でエラーがでます。
ご存知の方いらっしゃいましたらご教授のほどお願いいたします。
環境は Windows 10 です。
書いたコード
av <- 1:...
2
票
1
回答
668
閲覧数
Stata の dta ファイルを開くには?
現在イギリス大学院課程をオンラインにて履修中の者です。
統計学の授業内で R と RStudio、および指定のdtaファイルをダウンロードするように指示されました。
RStudioのダウンロードまでは問題なく進んだのですが、dtaファイルをダウンロードすると、添付の画像のように文字化けして表記されます。よく見るとところどころ英語が確認できるのですが…
こちらはなぜ文字化けしているのか、...
0
票
1
回答
278
閲覧数
Rを使ってエクセルの列のデータをループさせオブジェクトに格納したい
R言語を用いて、エクセルのデータを列ごとにオブジェクトに格納したいです。
この時、オブジェクト数をエクセルのデータの列数に応じて用意することができるようなプログラムを考えたいです。
以下の画像のようなデータがエクセルのブックに格納されています。(sam_data245.xlsx)
ファイルの形式は.xlsx形式で、データはブック内sheet1に入っています。
1行目にはデータのラベル(データ名)、...