[r] タグが付いた質問

Rはフリーかつオープンソースなプログラミング言語とソフトウェア環境で、統計解析、生物情報学、グラフ描画などに用いられます。質問の際は再現性のある最小限の例を添えていただけると助かります。

絞り込み
並べ替え
タグ付けされた
0
2 回答
103 閲覧数

Rでのwebスクレイピングで、データが抽出できない

Rでwebスクレイピングを行っているのですが, 抽出したいデータが変数内に入りません. 抽出したいのは, テーブル型のデータです. サイトは, 東京都教育委員会の都立高校検索のサイトで、アクセスのところの駅・路線・所要時間のデータを抽出したいと考えております. 使用したライブラリ:rvest Rのバージョン:R 4.2.1 使用したOS:macOS 以下にコードを添付します. library(...
ぴっぴ's user avatar
0
1 回答
82 閲覧数

ロジスティック回帰の最尤推定におけるヘッセ行列の計算ができない

ロジスティック回帰の最尤推定におけるヘッセ行列の計算についてご教示ください。 以下、R言語での記述ですが、 hessian <- function(beta, X) { p <- exp(X %*% beta) / (1 + exp(X %*% beta)) W <- diag(p * (1 - p)) return(t(X) %*% W %*% X) # ...
Yuki Ino's user avatar
0
0 回答
105 閲覧数

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:...
K. Doi's user avatar
  • 1
1
1 回答
86 閲覧数

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が正しく計算できていないように見えるにははなぜでしょうか?また解消法はあるでしょうか? 追記:...
YAGT's user avatar
  • 11
0
0 回答
38 閲覧数

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"...
K. Doi's user avatar
  • 1
0
1 回答
79 閲覧数

RでNAを含む連続変数がカテゴリー変数と判断されます

Rで-100から∔100までの連続変数(ただし計算上は±5毎に区切られます)の平均、SD値を、条件A/B群でそれぞれ比較し、TableOneで解析をしています。 この連続変数(BI_dif)の中にNAが多数あるからでしょうが、TableOneを作成すると、5毎のカテゴリー変数別に割合が表示され、平均・SD値が表示できません。 pt_01_excel %>% replace_na((list(...
K. Doi's user avatar
  • 1
0
1 回答
119 閲覧数

Rをバッチファイルのみで実行する方法

Rをインストール済みのWindowsPCでバッチファイルからRを実行したいと考えております。 2行の短いプログラムなので、プログラムファイルは別に置かずに、バッチファイル内にコードを記述する方法を探していますが見つかりません。 (バッチファイルからRプログラムを実行する方法は見つかるのですが。) 対応できる方法はありますでしょうか。ご教示いただけますと幸いです。 どうぞよろしくお願いいたします。
user87562's user avatar
  • 923
1
2 回答
138 閲覧数

PythonまたはRでクロス集計をしたい

PythonまたはRで以下のCSVファイルを読み込んで、クロス集計を行いたいのですが、うまくいきません。 男性┗肉のようにクロス軸を二次元にする事は出来るのですが、一次元にする事が出来ません。 CSVデータイメージ クロス集計イメージ ChatGPTで質問したところ以下の回答でしたが、主力イメージが意図した結果になりませんでした。 ご教示いただけますと幸いです。宜しくお願いいたします。 # ...
PG初心者's user avatar
0
1 回答
63 閲覧数

forで回して作成した複数のヒストグラムにサブタイトルを追加したい

下記のスクリプトだと、サブタイトルに subtitle_file[i, 2] のすべてが結合された状態で追加されてしまいます。私は、ヒストグラムが作成されるたびに、 subtitle_file[i, 2] の1つが順番に追加されるスクリプトを作成したいです。 ヒストグラム作成に使用するファイル数とID_Name.txtの行数は一致しています。 # サブタイトルを格納するための変数を初期化 ...
さかな's user avatar
0
0 回答
51 閲覧数

Missing required option: '--preset <name>'というエラーで困っています。

”Platypus”というRのプラットホームでレパトア解析を行いたいと思っています。まず開発元の Vignettes を用いて使い方を確認しているのですが、7.3のMIXCRを使うところで"Missing required option: '--preset '"のメッセージがでてきます。 どなたかわかる方がいましたらご教示いただけるとありがたいです。よろしくお願い致します。 ...
Sam's user avatar
  • 1
1
2 回答
129 閲覧数

Boxplotを条件ごとに作成する。

ggplot2を用いて、経時的なサンプルの変化を箱ひげ図で示したいと思っています。 以下のような手順で進めると、添付画像のように横軸のDayが反映されない箱ひげ図が出力されてしまいます。 Dayごとに箱ひげ図をまとめるには、どのようにすれば良いでしょうか? Boxplotに詳しい方、コメントを頂けないでしょうか? library(tidyverse) library(easystats) #...
user56816's user avatar
0
1 回答
66 閲覧数

Rstudio起動時に前回エラーを起こした作業スペースが読み込まれて作業ができない

RStudioでコードを書いていたところ、以下のエラーが表示され、その後スクリプトに書き込みはできるもの計算ができなくなりました。 fatal error; unexpected exception: string too long ファイルを消去すれば解決するだろうと思い、エラーが出たファイルを削除してRStudioを再び起動すると、...
c o Q's user avatar
  • 1
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,...
user56816's user avatar
0
1 回答
127 閲覧数

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) 以上の条件下で実行したいと思ってますが、...
mtaya33's user avatar
1
1 回答
121 閲覧数

Rでlmeを用いた解析をする際のエラー;invalid type (list) for variable 'data'

Rでlmeモデルを使った解析をしようとしていますが、以下のように解析しようとするとエラーが出てしまいます。原因を調べてみましたがわかりません。 データの形はdataframeであり、列の型はcharacterとintが混在しています。 コード: data <- read.table("data.txt",header=T,sep="\t") ...
金光由夏's user avatar
0
1 回答
158 閲覧数

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)) ...
user54886's user avatar
0
1 回答
68 閲覧数

グループ列があるデータセットの相関行列の作成

以下のような縦積みのデータセットで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
user54886's user avatar
0
1 回答
118 閲覧数

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, ...
user54886's user avatar
0
1 回答
131 閲覧数

RStudioで Markdown の出力時、警告が出て止まってしまう

RStidioでパッケージをロードすると、マークダウンの出力画面では警告が出て(しかもその警告の文が途中までで)、以後の処理は行われず止まってしまいます。 左上画面でチャンクごとに実行したところ、以下の警告が表示されますが、以後の処理もチャンクごとに実行することはできています。 Warning: package ‘ggplot2’ was built under R version 4.1.3 ...
Otaku Rounin's user avatar
1
1 回答
49 閲覧数

2条件のグループ別の平均値

以下のデータフレームに対してgroup1とgroup2それぞれの要約統計量(平均、最大値、最小値など)を出力したいです。ご教示いただけますと幸いです。 例) MenのAグループの平均値= MenのBグループの平均値= WomenのAグループの平均値= WomenのBグループの平均値= df <- data.frame(age=1:100, group1=c(rep("Men&...
user54886's user avatar
0
0 回答
107 閲覧数

Rのsegmentedパッケージの方法

Rのsegmentedパッケージを使用してセグメント回帰の図を出力をしたいと考えています。ブレークポイントは散布図を視覚的に確認して任意で設定するのでしょうか?もしくはRが自動で検出してくれるのでしょうか?コード含めご教示いただけますと幸いです。
user54886's user avatar
1
1 回答
223 閲覧数

複数の列の条件から新しい列を作成する方法

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(...
user54886's user avatar
1
1 回答
88 閲覧数

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
user54886's user avatar
0
1 回答
206 閲覧数

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 ← ...
user54886's user avatar
0
0 回答
86 閲覧数

Rstudio 遺伝子発現解析のt検定に際してのエラーが解決できません。

Rについて質問です。 初心者で申し訳ございません。 現在、腫瘍組織と正常組織を比較して、各遺伝子の発現量に差があるかを、t検定しております。 それに際してのエラーが解決できず、質問させていただきます。 func<-function(x){ result<-t.test(tumor[,x],normal[,x]) pvalue<-result$p.value ...
Tak Ta's user avatar
  • 1
2
0 回答
77 閲覧数

R ggplot2で100%積み上げグラフ

ある疾患の届出数の推移(2010~2021年)を患者の年齢階層別に100%積み上げグラフで作成したいと思っています。 以下のコードでグラフを作成したのですが、X軸の年が2010から4年間隔で2022年まで表示されてしまいます(2010、2013、・・・2022)。この2022を消すか、最後のみ2021を表示させるためにはどうしたらよいでしょうか。 よろしくお願いします。 Fig3dat<-(...
Lisa's user avatar
  • 21
1
1 回答
229 閲覧数

Rで文字列の日付データをDate型に変換したい

文字列の日付データをDate型に変換したいのですが、以下のようなエラーが出ました。 エラーコードの指している事がどういう事なのかがわからずに困っています。 実行コード asp <- asp %>% select(date,genre_name,amount,device_type) asp$date <- as.Date(asp$date) エラーコード ...
timon's user avatar
  • 11
1
1 回答
86 閲覧数

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 ...
r_panda's user avatar
  • 15
0
1 回答
330 閲覧数

R でhist(x) 関数でヒストグラムを描けない( 'x' は数値でなければなりません とエラーを返す)

Rで,以下のようなコードで 関数を作り,その結果をヒストグラムにしたいのですが, 'x' は数値でなければなりません とエラーを返してきてうまくいきません. 助けてください. func01 = function(n){ sample1 <- 100 df1 = data.frame() for(i in 1:n){ sample1 <- sample1+1000 df1 <- ...
takashikun's user avatar
0
1 回答
148 閲覧数

Rのクロス集計から条件をつけて削除

Rでクロス集計(table(X,Y))をした結果、例えばX*Yが10未満だった行を削除したいのですが、一気にやる方法はありますでしょうか?
r_panda's user avatar
  • 15
0
0 回答
1,870 閲覧数

Rのパッケージインストール時のエラー

Rstudioでパッケージを何かインストールしようとすると毎回 Error in install.packages : error reading from connection というエラーが表示されます。 何度かR, R studio, rtoolsをダウンロードしなおしました。 原因不明です・・既にほかのアプリも入れているPCなので、あまりクリーニングはしたくありません。 Rで、...
sheep_apprentice's user avatar
0
1 回答
72 閲覧数

R言語で土日祝日を除いた日数の計算結果をデータフレームへ

Rで土日祝日を除いた開始日と終了日の間の日数を計算しようとしましたが、上手くいきません。 ベクトルの型が上手くいかないものと思いますが、何か上手くやる方法はありませんでしょうか。 library(dplyr) p1 <- c(1:2) p2 <- as.Date(c("2013-01-01","2013-01-22")) p3 <- ...
r_panda's user avatar
  • 15
0
1 回答
532 閲覧数

Rの色々なエラーで悩まされています

rglのライブラリを読み込もうとすると下記のエラーが出ます。 既にxfunがwin-libraryにもRのlibraryフォルダの中にあることも確認済みだし、何回かxfunとrglの再インストールもしていて、原因がわかりません。 rdsのインストールも一応してあります。 Error: package or namespace load failed for ‘rgl’ in library....
sheep_apprentice's user avatar
1
1 回答
64 閲覧数

R 文字型データから日時データへの変換

"2022.04.30 午前 5:55" のような文字列データをRの日時データに変換したいです。 文字列を地道に置換していけば出来るような気もしなくはないですが、スマートなやり方をご存知でしたら教えて頂きたく。
Yuki Ino's user avatar
0
1 回答
733 閲覧数

Rでパッケージのインストールでエラーが表示されてしまう

Rx64 4.1.2でSteel-Dwass検定をする際、 「パッケージ」→「パッケージのインストール」→「Japan(Tokyo)」→「NSM3」 の指示通りにするが、Japanを選択以降、以下のメッセージが表示されてしまい、NSM3を選択出来ません。 utils:::menuInstallPkgs() --- このセッションで使うために、CRAN のミラーサイトを選んでください --- ...
user51545's user avatar
0
1 回答
140 閲覧数

rpy2に関して、Rのディレクトリの場所をどのように特定しているのでしょうか

<環境> os:windows10 実行環境:anaconda3 - jupyter notebook 主要ライブラリ:rpy2 3.7 〈本文〉 pythonでRを使用するためのrpy2というライブラリについて質問です。python環境下でrpy2を利用してRを使用することに成功しました。 一つ疑問なのですが、タイトルにもあるとおり、...
shihatu's user avatar
  • 11
0
2 回答
137 閲覧数

R ヒストグラム x軸が名義変数の場合の作図方法を知りたいです。

Rを使って、ヒストグラムを描きたいです。 対象とする変数は、数値の大小に意味がある順序変数ではなく、「名義変数」に分類される業種コード(ex.製造業:01,金融:19,食品:38など)です。 最終的に業種コードごとに度数が表示される図にしたいのですが、現状としてx軸が順序変数としてみなされているため、業種コードのなかでとんでいる(もともと存在しない)番号が数値として残ってしまう状態です。 ...
mochi's user avatar
  • 11
0
1 回答
226 閲覧数

fixest::etable()を使ってクロスリファレンスをする方法はありますか?

英語版でも同様の質問 を投稿しています。 質問 fixest::etable() を使ってクロスリファレンスをする方法はありますか?etable は,簡易的なメソッドですが,とても便利です.もしできる方法があるととても助かります. 環境は,Windows 10で,R 4.1.2を使っています.RStudioでbookdown::render_book()をしています.ファイル構成と内容は以下です. ...
ryoto's user avatar
  • 115
0
2 回答
1,412 閲覧数

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-...
mochi's user avatar
  • 11
1
1 回答
87 閲覧数

データフレームからあるベクトルの要素を含む行だけを抜き出し、新しいデータフレームを作成したい

使用するデータフレームは以下のような形です。 name year price A社 1999 200 A社 2000 202 A社 2001 199 A社 2002 400 A社 2003 207 B社 1999 300 B社 2000 500 ...
mochi's user avatar
  • 11
1
1 回答
582 閲覧数

RMeCabが日本語ではうまくいきません。

質問:日本語に対してRMecCabをうまく動かせるようにしたい。 RMeCabパッケージの関数を日本語のテキストデータに対して実行すると以下のようになり失敗します。 極楽not foundとなっていますが、実際には極楽は存在します。ちなみに英語のみのテキストデータに対してはうまく動きます。 > library("RMeCab") > r <- ...
Camford Oxbridge's user avatar
0
2 回答
78 閲覧数

完全パネルデータの作成にあたり、観測期間(年次データ)のすべてが揃った個体のデータのみを、新しいデータフレームに保存したいです

使用するデータは多少簡略化していますが、以下のようなものです。 name year ... A社 1999 ... A社 2000 ... A社 2001 ... A社 2002 ... B社 1999 ... B社 2000 ... C社 2001 ... C社 2002 ... D社 ...
mochi's user avatar
  • 11
1
2 回答
278 閲覧数

同じ構成のfor文でも、処理数を増やすとエラーが出るようになる

3つ目のfor文がなぜエラーになるのか、わかる方いらっしゃいますか? 下記のコードは各オブジェクト作成からfor文の処理まで一部始終のコードと、途中からの出力結果となります。 とてもシンプルに書いているつもりなのですが、同じ構成の3つ目のfor文でエラーがでます。 ご存知の方いらっしゃいましたらご教授のほどお願いいたします。 環境は Windows 10 です。 書いたコード av <- 1:...
haru's user avatar
  • 81
2
1 回答
603 閲覧数

Stata の dta ファイルを開くには?

現在イギリス大学院課程をオンラインにて履修中の者です。 統計学の授業内で R と RStudio、および指定のdtaファイルをダウンロードするように指示されました。 RStudioのダウンロードまでは問題なく進んだのですが、dtaファイルをダウンロードすると、添付の画像のように文字化けして表記されます。よく見るとところどころ英語が確認できるのですが… こちらはなぜ文字化けしているのか、...
プログラミング超初心者's user avatar
0
1 回答
268 閲覧数

Rを使ってエクセルの列のデータをループさせオブジェクトに格納したい

R言語を用いて、エクセルのデータを列ごとにオブジェクトに格納したいです。 この時、オブジェクト数をエクセルのデータの列数に応じて用意することができるようなプログラムを考えたいです。 以下の画像のようなデータがエクセルのブックに格納されています。(sam_data245.xlsx) ファイルの形式は.xlsx形式で、データはブック内sheet1に入っています。 1行目にはデータのラベル(データ名)、...
milk's user avatar
  • 73
0
1 回答
112 閲覧数

ステップワイズ法によるBCsFt集団のQTL解析

ステップワイズ法によるQTL解析を行うときに以下のエラーが出てしまいます。 No default penalties available for cross type bcsft cross type bcsftでモデルを構築するにはどうすればよいのでしょうか。 調べてみましたがわかりませんでした。よろしくお願いいたします。 mapFunc <- "kosambi" ...
user48327's user avatar
2
1 回答
888 閲覧数

Rでロジスティック回帰を行ったときのdeviance残差の計算方法

Rでロジスティック回帰を行った際のresidual devianceを取得したいと考えています。 glm.fit関数を用いたところ、以下のような完全分離の警告が出ました。 glm.fit: fitted probabilities numerically 0 or 1 occurred そこで、こちらを参考にlogistf関数を利用してみました。 cont <- logistf....
katht's user avatar
  • 21
2
1 回答
394 閲覧数

データフレームを表示するとき、列名と各列の数値の水平位置を揃える方法

画像のように、出力時に列名称が各列の数値と位置が揃わない現象をなんとかしたいと考えております。 画像はオブジェクトrr データフレームの一部です。 揃える方法をご存知の方いらっしゃいましたら ご教授のほどお願いいたします。 環境は Windows 10 です。 ※視認性向上のため質問情報を整理しました。
haru's user avatar
  • 81
1
1 回答
196 閲覧数

modelsummaryのmodelplot()を使って文字化けしない図をpdfで出力するためにはどうすればよいですか?

R Markdownを使ってpdf文書を作成しています。 modelsummaryパッケージに含まれるmodelplot()関数を使ってモデルの推定値と標準誤差を図で示したいですが、図の凡例に日本語が含まれる場合、文字化けを起こします。公式のドキュメントでは、日本語を扱う方法は触れられていませんでした。この問題を解決するためにはどうすればよいでしょうか? ...
ryoto's user avatar
  • 115
2
1 回答
1,444 閲覧数

library(Rcmdr) でエラー: ‘Rcmdr’ という名前のパッケージはありません

プログラミングに関して素人もいいところなのですが 腸内細菌の研究をやる必要がありqiime2を使用しています。 解析途中で躓いてしまい教えていただけたらと思います。 denoiseの作業のところで、以前は問題なくできていたのですが 別のことでezRをダウンロードしてから?(関係あるかわかりませんが) ストップするようになってしまいました。 --verbose でエラーをみると (qiime2-...
yokmok's user avatar
  • 21

1
2 3 4 5
7