[r] タグが付いた質問

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

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

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

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

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

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

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

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

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

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

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

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

条件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 回答
58 閲覧数

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 回答
69 閲覧数

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

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

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
0 回答
35 閲覧数

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

以下のような縦積みのデータセットで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
0 回答
64 閲覧数

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
0 回答
50 閲覧数

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

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

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

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

Rのsegmentedパッケージの方法

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

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

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 回答
67 閲覧数

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 回答
116 閲覧数

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 ← ...
0
0 回答
64 閲覧数

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

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

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

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

回帰分析をする際に、説明変数のデータ型の指定はする必要がありますか

はじめまして。Rでアンケート調査の分析を試みております。 目的変数はダミー変数で、説明変数には5段階評価の回答結果、ダミー変数、連続変数がそれぞれ複数含まれています。 そこで、回帰をする際に、説明変数のデータ型の指定はする必要があるかどうかをお聞きしたいです。 データを読み込むと、最初はどの変数もデータ型がdbl扱いになっているのですが、 5段階評価の回答のうち、1<2<3<4<...
  • 1
1
1 回答
119 閲覧数

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

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

Rのtcplパッケージを利用して解析を行いたい

ToxCastのプログラムで、細胞実験の結果からデータ解析が可能なtcplパッケージが作成されている。 GitHubもあるが、使い方が分からない。 https://github.com/USEPA/CompTox-ToxCast-tcpl 使い方を少しでも共有して頂けると幸いです。
user avatar
0
0 回答
46 閲覧数

Rプログラムのパッケージ、dplyrのsummarise関数について

データ分析をした際に疑問が浮かび上がったので、質問させていただきます。 データ分析をする際に定番のgroupby関数からsummarise関数を適用していくと思われますが、summarise関数適用後と適用前で格納されているデータの数が違うのはなぜでしょうか。 私の予想では、SQLのdistinctなどが適用されているように感じますが、確証がないので、質問させていただきます。 また、...
1
1 回答
76 閲覧数

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 ...
  • 15
0
0 回答
48 閲覧数

RでGCNs (Graph Convolutional Networks)はできますでしょうか?

RでGCNs (Graph Convolutional Networks)はできますでしょうか?
0
1 回答
191 閲覧数

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 回答
106 閲覧数

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

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

Rのlibrary実行時のエラー

R・Rstudio・rtoolsをインストールしたあと、install.packagesでtidyverse・rgl等をインストールしたのですが、library(rgl) を実行したときのみ、以下のエラーが表示されてしまいます。 install.packages("rgl", dependencies=TRUE) ...
0
0 回答
1,125 閲覧数

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

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

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

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

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

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

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

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

Rで多変量時系列データをクラスタリングしたい

私は、150名の参加者を対象に3時点の感情状態を尋ねた時系列データを持っています。 感情状態は28種類あり、それぞれその強さについて6段階で尋ねています。 この3時点の変化パターンに基づいてクラスタリングしたいと考えているのですが、 なにかよい分析方法はありますでしょうか。 統計ソフトはRを使っています。 データの形は以下のようなものです。 PP Time1: Happy Time2: ...
0
1 回答
464 閲覧数

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

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

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

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

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

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

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

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

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-...
  • 11
1
1 回答
66 閲覧数

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

RStudioでRmarkdownファイルを作成中、日本語メッセージが文字化けする

環境 OS: Windows 10 Home 64ビット Rのバージョン: 4.0.3 RStudioのバージョン: 1.4.1106 ディレクトリ:C:\Users\xxx\rmarkdownfile\test.Rmd (xxxには自分のユーザー名が入ります) 問題 初めてこちらで質問させて頂きます。 表題の通り、RStudioでRmarkdownファイルを作成しているのですが、...
  • 11

1
2 3 4 5
7