[r] タグが付いた質問
Rはフリーかつオープンソースなプログラミング言語とソフトウェア環境で、統計解析、生物情報学、グラフ描画などに用いられます。質問の際は再現性のある最小限の例を添えていただけると助かります。
323
件の質問
0
票
0
回答
37
閲覧数
mutate()で追加中の値を参照する方法
csvデータを読み込んで,mutate()とif_else()を使って新しい列を追加しようとしています.
一つ前の行の値を参照してmutete()を行う方法を,こちらでお教えいただきました.
こちらの続きなのですが,mutate()で追加中の値を参照する方法はあるのでしょうか.
csvデータは20行あります.
ex.csv
name value
hoge 5
huga 7
piyo 4
puka ...
1
票
1
回答
76
閲覧数
dplyr::mutate()におけるif_else()で, "size" に関するエラー
csvデータを読み込んで,mutate()とif_else()を使って新しい列を追加しようとしています.
csvデータは20行あります.
ex.csv
name value
hoge 5
huga 7
piyo 4
puka 2
...
以下コード
library('tidyverse')
d <- read.csv('ex.csv', sep = '\t', header = T)
d2 ...
0
票
1
回答
48
閲覧数
Rのkernlabで、kernel matrixからSVMを行ったときの不具合について
ガウシアンカーネルを用いていますが、説明変数からSVMを行った場合と、カーネル行列からSVMを行った場合で結果が異なります。どうしてこのような違いが生じるかわかりません。
data(spam)
index <- sample(1:dim(spam)[1])
spamtrain <- spam[index[1:2000],]
spxx <- as.matrix(spamtrain[,...
0
票
0
回答
27
閲覧数
ローカルのNexusインスタンスをAzure Container Registryに接続できない
私はプロジェクトのために、Azure Container Registry (ACR)に接続するように設定されたDockerプロキシで、3つの別々のローカルNexusインスタンス(各データセンターに1つ)をセットアップしようとしています。目標は、各データセンターのNexusリポジトリをACRに保存されたイメージのローカルキャッシュとして機能させ、...
0
票
0
回答
26
閲覧数
R packageのtidychatmodelsを用いる際にtime_outを管理する方法について
R言語のtidychatmodelsパッケージを用いてR環境内にchatbotの作成に取り組んでいます。
modelでgpt-4-turboを選択し、コードを実行すると以下のエラーが発生します。
add_modelの引数をgpt-3.5-turboとすると実行はされるのですが、hallucinationが出現し、意味の無い回答となってしまったため、gpt-4-...
0
票
0
回答
62
閲覧数
pythonまたはRで共分散構造分析をする際、2値データを扱う方法を教えてください。
アンケートデータを使って共分散構造分析を行いたいのですが、回答がはい/いいえの2値の場合、テトラコリック相関を使う、という情報にたどり着きました。
pythonかRの共分散構造分析用ライブラリ(手元ではsemopyを試しています)で、2値データを使って分析を行う場合、どのようにコードを書けばよいでしょうか。
0
票
0
回答
40
閲覧数
shpファイルとcsvファイルを読み込みモラン分析を行った。その分析結果をコロプレス図に表したいが、どのような手順が必要か?
shpファイルとcsvファイルを読み込みモラン分析を行った。その分析結果をコロプレス図に表したいが、どのような手順が必要か?
-3
票
1
回答
241
閲覧数
csvを読み込みモラン散布図を作成したいが手順が分からない。
csvデータをもとに、モラン散布図を作成したいのですが、データを読み込みそこから何をすればモラン散布図が作れるのか理解ができていない。しなければいけない手順やコードとそのコードの意味を教えてほしい。
0
票
1
回答
68
閲覧数
e_mark_point で min、avg、max の異なる値を返す方法は?
このデータがあります:
library(echarts4r)
x <- c(10, 20, 30, 45, 46, 50)
y <- c(12, 25, 36, 46, 49, 59)
matrixinc1 <- data.frame(x, y)
near_line <- function(matriz) {
diferenca <- abs(matriz[,...
2
票
1
回答
174
閲覧数
modelsummary関数で数値を右に並べたい
標準誤差を推定値の下ではなく、右に表示させたいです。
画像では推定値の下に標準誤差が出力されています。
試したコード: rstudioで実行し、バージョンはR version 4.1.3 (2022-03-10)になります。
library(modelsummary)
model_summary <- modelsummary(model, output = "default&...
0
票
0
回答
41
閲覧数
fastLinkを用いたProbabilistic data Linkageについて
fastLinkパッケージを用いて、二つのデータセットのリンクを試みています。リンクに用いる項目日付、時間、年齢、性別で全て数値に変換しています(時間:2/11→0211、11:15→1115、男性矢印、以下のようなスクリプトを作りました。
result <- fastLink(dfA = data_2, dfB = data_1, varnames = c("date",...
0
票
0
回答
85
閲覧数
行列の逆行列が計算できない(Rの系統樹の処理)
系統樹は自分で以下のように作成しています。
どうやってもエラーが解消されなくて困っています。作成方法がいけないのでしょうか??
> newick_tree<-"(A:4,(B:2,(C:2.953469601,D:1.443532938):1.136357167):1.104667235);"
> phylo_tree <- read.tree(text ...
1
票
2
回答
278
閲覧数
R言語でデータを年代別にし、箱ひげ図を作成したい
データは年齢と幸福度になります。
年代別の箱ひげ図をRで作成したいです。
私の想定としては、以下の二つが必要だと思っています。
・年齢を20~29などに設定し、まず年代別にする。
・箱ひげ図をそれぞれの年代でプロットするようにする。
#hensuu1 <- age, hensuu2 <- kouhukudo
ggplot(data, aes(x = hensuu1, y = ...
0
票
0
回答
215
閲覧数
アーカイブから"mvpart"パッケージをダウンロードできません
いくつかのドキュメントを参考にしましたが、アーカイブから"mvpart"パッケージをダウンロードしようとすると下記の通りエラーになってしまいます。
両方ともunknown type name 'Sint'が原因のようですが、何のことだか見当がつきません。
ご教示いただけますと幸いです。
どうぞよろしくお願いいたします。
試行したコード①
install.packages(&...
0
票
0
回答
136
閲覧数
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
回答
109
閲覧数
Rでのwebスクレイピング
Rでwebスクレイピングしたいと考えています。
抜き取りたい情報はM-1グランプリの情報で
(https://www.m-1gp.com/schedule/detail.html?id=539)におけるグループ名、コンビ名、結果の情報とコンビ名をクリックした先にある結成年、所属、前年度における結果等を抜き出そうと考えています。
初めのグループ名、コンビ名、...
0
票
2
回答
328
閲覧数
Rでのwebスクレイピングで、データが抽出できない
Rでwebスクレイピングを行っているのですが, 抽出したいデータが変数内に入りません.
抽出したいのは, テーブル型のデータです.
サイトは, 東京都教育委員会の都立高校検索のサイトで、アクセスのところの駅・路線・所要時間のデータを抽出したいと考えております.
使用したライブラリ:rvest
Rのバージョン:R 4.2.1
使用したOS:macOS
以下にコードを添付します.
library(...
0
票
1
回答
123
閲覧数
ロジスティック回帰の最尤推定におけるヘッセ行列の計算ができない
ロジスティック回帰の最尤推定におけるヘッセ行列の計算についてご教示ください。
以下、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
回答
364
閲覧数
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
回答
96
閲覧数
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
票
1
回答
132
閲覧数
RでNAを含む連続変数がカテゴリー変数と判断されます
Rで-100から∔100までの連続変数(ただし計算上は±5毎に区切られます)の平均、SD値を、条件A/B群でそれぞれ比較し、TableOneで解析をしています。
この連続変数(BI_dif)の中にNAが多数あるからでしょうが、TableOneを作成すると、5毎のカテゴリー変数別に割合が表示され、平均・SD値が表示できません。
pt_01_excel %>% replace_na((list(...
0
票
1
回答
396
閲覧数
Rをバッチファイルのみで実行する方法
Rをインストール済みのWindowsPCでバッチファイルからRを実行したいと考えております。
2行の短いプログラムなので、プログラムファイルは別に置かずに、バッチファイル内にコードを記述する方法を探していますが見つかりません。
(バッチファイルからRプログラムを実行する方法は見つかるのですが。)
対応できる方法はありますでしょうか。ご教示いただけますと幸いです。
どうぞよろしくお願いいたします。
1
票
2
回答
176
閲覧数
PythonまたはRでクロス集計をしたい
PythonまたはRで以下のCSVファイルを読み込んで、クロス集計を行いたいのですが、うまくいきません。
男性┗肉のようにクロス軸を二次元にする事は出来るのですが、一次元にする事が出来ません。
CSVデータイメージ
クロス集計イメージ
ChatGPTで質問したところ以下の回答でしたが、主力イメージが意図した結果になりませんでした。
ご教示いただけますと幸いです。宜しくお願いいたします。
# ...
0
票
1
回答
81
閲覧数
forで回して作成した複数のヒストグラムにサブタイトルを追加したい
下記のスクリプトだと、サブタイトルに subtitle_file[i, 2] のすべてが結合された状態で追加されてしまいます。私は、ヒストグラムが作成されるたびに、 subtitle_file[i, 2] の1つが順番に追加されるスクリプトを作成したいです。
ヒストグラム作成に使用するファイル数とID_Name.txtの行数は一致しています。
# サブタイトルを格納するための変数を初期化
...
1
票
2
回答
262
閲覧数
Boxplotを条件ごとに作成する。
ggplot2を用いて、経時的なサンプルの変化を箱ひげ図で示したいと思っています。
以下のような手順で進めると、添付画像のように横軸のDayが反映されない箱ひげ図が出力されてしまいます。
Dayごとに箱ひげ図をまとめるには、どのようにすれば良いでしょうか?
Boxplotに詳しい方、コメントを頂けないでしょうか?
library(tidyverse)
library(easystats)
#...
0
票
1
回答
233
閲覧数
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
回答
256
閲覧数
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
回答
253
閲覧数
Rでlmeを用いた解析をする際のエラー;invalid type (list) for variable 'data'
Rでlmeモデルを使った解析をしようとしていますが、以下のように解析しようとするとエラーが出てしまいます。原因を調べてみましたがわかりません。
データの形はdataframeであり、列の型はcharacterとintが混在しています。
コード:
data <- read.table("data.txt",header=T,sep="\t")
...
0
票
1
回答
222
閲覧数
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
回答
90
閲覧数
グループ列があるデータセットの相関行列の作成
以下のような縦積みのデータセットで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
回答
132
閲覧数
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
回答
318
閲覧数
RStudioで Markdown の出力時、警告が出て止まってしまう
RStidioでパッケージをロードすると、マークダウンの出力画面では警告が出て(しかもその警告の文が途中までで)、以後の処理は行われず止まってしまいます。
左上画面でチャンクごとに実行したところ、以下の警告が表示されますが、以後の処理もチャンクごとに実行することはできています。
Warning: package ‘ggplot2’ was built under R version 4.1.3
...
1
票
1
回答
65
閲覧数
2条件のグループ別の平均値
以下のデータフレームに対してgroup1とgroup2それぞれの要約統計量(平均、最大値、最小値など)を出力したいです。ご教示いただけますと幸いです。
例)
MenのAグループの平均値=
MenのBグループの平均値=
WomenのAグループの平均値=
WomenのBグループの平均値=
df <- data.frame(age=1:100, group1=c(rep("Men&...
1
票
1
回答
369
閲覧数
複数の列の条件から新しい列を作成する方法
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
回答
142
閲覧数
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
回答
420
閲覧数
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
回答
101
閲覧数
R ggplot2で100%積み上げグラフ
ある疾患の届出数の推移(2010~2021年)を患者の年齢階層別に100%積み上げグラフで作成したいと思っています。
以下のコードでグラフを作成したのですが、X軸の年が2010から4年間隔で2022年まで表示されてしまいます(2010、2013、・・・2022)。この2022を消すか、最後のみ2021を表示させるためにはどうしたらよいでしょうか。
よろしくお願いします。
Fig3dat<-(...
1
票
1
回答
449
閲覧数
Rで文字列の日付データをDate型に変換したい
文字列の日付データをDate型に変換したいのですが、以下のようなエラーが出ました。
エラーコードの指している事がどういう事なのかがわからずに困っています。
実行コード
asp <- asp %>%
select(date,genre_name,amount,device_type)
asp$date <-
as.Date(asp$date)
エラーコード
...
1
票
1
回答
103
閲覧数
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
回答
716
閲覧数
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
回答
214
閲覧数
Rのクロス集計から条件をつけて削除
Rでクロス集計(table(X,Y))をした結果、例えばX*Yが10未満だった行を削除したいのですが、一気にやる方法はありますでしょうか?
0
票
0
回答
3,018
閲覧数
Rのパッケージインストール時のエラー
Rstudioでパッケージを何かインストールしようとすると毎回
Error in install.packages : error reading from connection
というエラーが表示されます。
何度かR, R studio, rtoolsをダウンロードしなおしました。
原因不明です・・既にほかのアプリも入れているPCなので、あまりクリーニングはしたくありません。
Rで、...
0
票
1
回答
105
閲覧数
R言語で土日祝日を除いた日数の計算結果をデータフレームへ
Rで土日祝日を除いた開始日と終了日の間の日数を計算しようとしましたが、上手くいきません。
ベクトルの型が上手くいかないものと思いますが、何か上手くやる方法はありませんでしょうか。
library(dplyr)
p1 <- c(1:2)
p2 <- as.Date(c("2013-01-01","2013-01-22"))
p3 <- ...
0
票
1
回答
746
閲覧数
Rの色々なエラーで悩まされています
rglのライブラリを読み込もうとすると下記のエラーが出ます。
既にxfunがwin-libraryにもRのlibraryフォルダの中にあることも確認済みだし、何回かxfunとrglの再インストールもしていて、原因がわかりません。
rdsのインストールも一応してあります。
Error: package or namespace load failed for ‘rgl’ in library....
1
票
1
回答
67
閲覧数
R 文字型データから日時データへの変換
"2022.04.30 午前 5:55"
のような文字列データをRの日時データに変換したいです。
文字列を地道に置換していけば出来るような気もしなくはないですが、スマートなやり方をご存知でしたら教えて頂きたく。
0
票
1
回答
1,233
閲覧数
Rでパッケージのインストールでエラーが表示されてしまう
Rx64 4.1.2でSteel-Dwass検定をする際、
「パッケージ」→「パッケージのインストール」→「Japan(Tokyo)」→「NSM3」
の指示通りにするが、Japanを選択以降、以下のメッセージが表示されてしまい、NSM3を選択出来ません。
utils:::menuInstallPkgs()
--- このセッションで使うために、CRAN のミラーサイトを選んでください ---
...
0
票
1
回答
219
閲覧数
rpy2に関して、Rのディレクトリの場所をどのように特定しているのでしょうか
<環境>
os:windows10
実行環境:anaconda3 - jupyter notebook
主要ライブラリ:rpy2 3.7
〈本文〉
pythonでRを使用するためのrpy2というライブラリについて質問です。python環境下でrpy2を利用してRを使用することに成功しました。
一つ疑問なのですが、タイトルにもあるとおり、...
0
票
2
回答
191
閲覧数
R ヒストグラム x軸が名義変数の場合の作図方法を知りたいです。
Rを使って、ヒストグラムを描きたいです。
対象とする変数は、数値の大小に意味がある順序変数ではなく、「名義変数」に分類される業種コード(ex.製造業:01,金融:19,食品:38など)です。
最終的に業種コードごとに度数が表示される図にしたいのですが、現状としてx軸が順序変数としてみなされているため、業種コードのなかでとんでいる(もともと存在しない)番号が数値として残ってしまう状態です。
...
0
票
1
回答
293
閲覧数
fixest::etable()を使ってクロスリファレンスをする方法はありますか?
英語版でも同様の質問 を投稿しています。
質問
fixest::etable() を使ってクロスリファレンスをする方法はありますか?etable は,簡易的なメソッドですが,とても便利です.もしできる方法があるととても助かります.
環境は,Windows 10で,R 4.1.2を使っています.RStudioでbookdown::render_book()をしています.ファイル構成と内容は以下です.
...