[r] タグが付いた質問
Rはフリーかつオープンソースなプログラミング言語とソフトウェア環境で、統計解析、生物情報学、グラフ描画などに用いられます。質問の際は再現性のある最小限の例を添えていただけると助かります。
81
未回答やベストアンサーなし質問
3
票
0回答
75 件の閲覧数
Rで中央値medianを計算する際のエラー
CSVで読み込んだデータのうち一つのカラムの中央値を求めようとするとエラーが出てしまいます。
データは2カラムのみで、12行、該当カラムはdbl型式です。
原因がもしわかればご教示いただければ幸いです。
mean(df$y)
#error: the condition has length > 1 and only the first element will be used[1] &...
3
票
1回答
189 件の閲覧数
Rで同じidのカラム内の変数の移行のチェック
R3.3.2とパッケージtidyverseを用いて前処理中です。
以下のようなデータなのですが、同じidについて複数回調査(chousaが調査回数)しています。
同じ調査項目(今回はQ1)が0から1に変わってから1回目、2回目、…4回目までのデータがほしいです。
どのようにmutateすればよいでしょうか?
よろしくお願いします。
ID Q1 Q6 chousa …
1 0 3 ...
3
票
0回答
346 件の閲覧数
R {hexbin} libraryを使用してplotする時のx-y rangeの設定について
Rを使って、図を作成する時に、
x-y座標のサイズを設定する方法について質問させていただいております。
先に使用するサンプルデータ及び問題があるコードを書かせていただきます。
以下のデータをsample.dataとして使用しています。
1247211284201945,A-vector-1,0.028382835,0.327685692,0.288425494
...
2
票
1回答
119 件の閲覧数
長期記憶性のある時系列データの解析に関する質問
現在、表題のようなデータに対して分析を実行しています。
データの内容に関しては規約があるため詳しく説明することができませんが、以下のような方法で採取されたデータです。ここで、被説明変数をy,説明変数をxとします。また、時刻をtとします。
時刻tにy(t)が観測される。直近の説明変数x(T|t-100<T<t),(直近100イベントのレコードx)が収集される。また、...
2
票
0回答
64 件の閲覧数
Rでデータフレームから特定のデータフレームをフィルタリングする方法
形態素解析の結果であるデータフレームdata1(1523行4列)から、別のデータフレームdata2(150行4列)の値と一致する行を削除したいのですが、
どのような関数を用いたらいいのかがわかりません。
以下が試してみた結果になります。
コード①
data1[!(data1$Term == data2$Term),]
結果①
Term Info1 Info2 Freq
...
2
票
1回答
119 件の閲覧数
Jupyter notebook上でRパッケージが読み込めない
プログラミング初心者です。
minicondaを経由してJupyter notebook上でRを使った解析を行おうとしています。
Rのカーネルの追加までは上手くいったのですが、Rパッケージを読み込もうとすると下記のようなエラーが出ます。
何が原因でしょうか?
OSはmacOS10.15.2を使っています。
pyenv1.2.15を使ってminiconda3-4.3.30を導入しています。
...
1
票
1回答
87 件の閲覧数
reticulateを使ってRstudio上でPythonを書こうとすると,R Session Abortedになってしまいます。
お世話になります。初めて質問させていただきます。
RstudioでPythonを使いたいのですが,毎回強制終了になってしまいます。
いくつか対処を試みたのですが,自分では解決できなかったので教えていただきたいです。
これまでに行ったことを記載していますが,必要な情報があれば,ご指摘ください。
恐れ入りますが,よろしくお願いいたします。
● 実施したいこと
...
1
票
0回答
31 件の閲覧数
{ggspatial}のannotation_map_tile()ができない
OpenStreetMapを背景にした地図を作りたいと思い、{ggspatital}のannotation_map_tile()を使って描写しようとしているのですが、できません。
library(tidyverse)
library(sf)
library(ggspatial)
dat <- data.frame(
lat= runif(10,120,140),
lng= ...
1
票
0回答
42 件の閲覧数
nlmeの表記について教えてほしいです
Rでの解析を始めたばかりの初心者なのですが、nlmeパッケージのnlme関数について質問があります。
以下はnlmeのhelpからとってきた例ですが、
fixed = Asym + R0 + lrc ~ 1,
random = Asym ~ 1,
のところの意味が分かりませんので教えていただけないでしょうか。
...
1
票
0回答
41 件の閲覧数
AtomでHydrogenを用いてRをコンパイルする場合に日本語がでない
AtomエディタでHydrogenパッケージ、jupyter、IRkernelを用いてRを動かそうとしており、
以下のサイトなどを参考に構築をしました。
https://nussilovesdata.eu/post/setting-up-atom-for-r/
https://gomiba.co.in/blog/archives/831
しかし、...
1
票
1回答
172 件の閲覧数
ggplot2で図が表示されたりされなかったりする。
Rでtidyverseとqiime2Rを用いて作図を試みています。
データを入力してスクリプトを実行すると、
1回目は作図されるのですが、図を保存して図ウィンドウ(Quartz)を閉じ、再度同じスクリプトを入力すると、
UseMethod("depth") でエラー:
'depth' をクラス "NULL" ...
1
票
1回答
90 件の閲覧数
リストをfor文で代入したい
6つの既存のリストを list[1],list[2],list[3],list[4],list[5],list[6] というような形で for 文で代入したいのですが、どうすればいいかわかりません。
ベクトルで for 文を用いて、f[1],f[2],f[3],f[4],f[5],f[6] に対して代入するときは、
以下のようになるのはわかるのですが、リストを同じようにするとうまく行きません。
...
1
票
1回答
152 件の閲覧数
gtfファイルの読み込みで ”R Session Aborted" 、RNA-seqデータ解析
R初心者でゲノム発現解析を始めました。
環境:MacOS、メモリ32GB使用、R version 4.0.0、Rtudio利用
HiSeqで解析したRNA-seq readカウント後のデータで、geneIDをgeneSymbolにするため、RパッケージrefGenomeを使用。
DAVIDでうまく変換できなかったりと他にも色々試しました。
アノテーションのgtfファイル(1.46G)...
1
票
0回答
57 件の閲覧数
正規分布のグラフに色を塗る
以下のようにJRIを使って正規分布のグラフを作り画像として保存しています。
import org.rosuda.JRI.REXP;
import org.rosuda.JRI.Rengine;
public class JriTest
{
public static void main(String[] args)
{
Rengine engine = new ...
1
票
1回答
624 件の閲覧数
anacondaターミナルを開くとerror reportが表示される
anacondaターミナルを開くと以下のエラーレポートが表示されます。
errorの解消方法が分かる方教えてください。
/Users/username/.anaconda/navigator/a.tool ; exit;
# >>>>>>>>>>>>>>>>>>>>>&...
1
票
1回答
87 件の閲覧数
R言語 オブジェクトの中身のアクセス
オブジェクトの中身のアクセス方法を教えていただきたいです。
具体的には,,,
sex<-c("F","M","M","M","F")
number<-c(1,2,3,4,5)
x<-data.frame(SEX=sex,NUMBER=number)
y<-c(1,2,4)
このような時,,,
c<-subset(x,NUMBER!=y)
...
1
票
1回答
216 件の閲覧数
java r interfaceの有効桁数
【質問】
java r interface(JRI)を使用して、Rの結果を取得しています。
JRIで返す桁数とRのGUIで返す桁数に差異があり、どこでJRIの桁数を指定したら良いのか分かりません。
どなたかご教示お願いします。
【事象】
RのGUI
1.getOption(digits)を実行して、デフォルトの7であることを確認。
2.「t.test」コマンドを実行して、tが-1....
1
票
0回答
152 件の閲覧数
Rによるデータテーブルの作成について
行が1900、列が15のデータフレームを読み込み、パッケージによる分析のためにTableもしくはmatrixへ変換しようと試みました。ですが、要素が上限を超えることから、エラーが生じます。
Data_table <- as.table(data)
Error in as.table.default(data) : table に強制変換できません
2^31 ...
1
票
0回答
103 件の閲覧数
オブジェクトの読み込みについて
私はSeuratを使用していてGEOからインストールしたデータセットを分析しようとしています。 しかし、CreateSeuratObjectに関するエラーメッセージが表示されました。
Warning messageの解決方法を教えていただけたらと思います。
よろしくお願いします。
> library(dplyr)
> library(Seurat)
> library(...
1
票
0回答
186 件の閲覧数
jupyter notebook で保存したグラフが正常に表示されない
jupyter notebookのR言語で、plotを使ってグラフを出力しました。
File>Download as>HTML(.html)で結果を保存したところ、以下のようにグラフが表示されません。何か解決方法はありますでしょうか。
環境はMacbook,ブラウザはChromeです。よろしくお願いいたします!!
ダウンロードする前
ダウンロードした後
1
票
0回答
209 件の閲覧数
vscodeの自動整形による文字化け
vscodeを使用してRのプログラミングを行っているのですが,vscodeの機能である自動整形(Format On Save)をONにしたところ文字化けしてしまいました.
print(paste("標本平均 : ", mean(y), seq = ""))
これを保存すると
print(paste("<e6><a8><99><e6><9c&...
1
票
0回答
680 件の閲覧数
Rでの特定データの順位付け
今,データフレームをとりこんで,例をあげるとこのようなかんじです。
このデータの「v2」に関して,aとbでグルーピングをし,それぞれのグループの中で順位をつけたいと思っています。
v1 v2 v3 v4
1 H29年 a 120 m
2 H29年 a 130 m
3 H30年 ...
1
票
1回答
91 件の閲覧数
Rで組み分けをその中から特定の項目をカウントしたいのですが
組 得意科目
1 A 国語
2 A 数学
3 B 国語
4 C 数学
5 C 数学
6 C 国語
7 C 国語
このような表から国語が得意な人間の数を組ごとに出したいのですが、
A 1
B 1
C 2
のような形で回答が出るようにしたいです。
どのように書けばよいでしょうか。
1
票
0回答
61 件の閲覧数
reticulateパッケージのiterate関数の使い方を教えてください。
library(reticulate)
x <- r_to_py(1:10)
x
[1, 2, 3, 4, 5, 6, 7, 8, 9, 10]
class(x)
[1] "python.builtin.list" "python.builtin.object"
iterate(x)
Error in iterate(x) : iterate function ...
1
票
1回答
359 件の閲覧数
Rで日本語ツイートの取得ができない
Rで以下のコードでツイート取得しようとしたらエラーが出ます。
OSはWindows10で、Rstudioを使ってます。
ソースコード
searchTwitter(searchString = "#100歳までの計画",
n = 500,
lang = "ja",
...
1
票
0回答
246 件の閲覧数
Rを使ったローリング回帰がうまくワークしません。(PerformanceAnalyticsno)
以下のコードでローリング回帰を行いたいのですが、
↓のようなエラーがでしまいます。。。。対応方法を教えて頂きたいです。
coefs= rollapplyr(dataset1, size.window, function(x){
+ coef(lm(y~x1+x2+x3+x4,data=x))}, by.column = FALSE)
Error in model.frame.default(...
1
票
0回答
77 件の閲覧数
Rの行指定について
Rのプログラミングについて質問です。
データフレームから行を抽出するときに、
ランダムな数が入った変数で指定したいのですが、
下記のコードで行なったところ1行しか出力できません。
data[[a : a + b], "列名"]
a,b はともにfloor(runif)で取得した
整数が入っているのですが、
a行目しか出力されません。。
解決策をご存知の方教えていただけませんか?
1
票
0回答
996 件の閲覧数
RStudio 起動時にフリーズし、操作不能のまま動作が停止する
何度アンインストールをし、再インストールしても件名の症状が出て動作停止してしまいます。
あるRmdファイルを開いたまま終了q()してしまったため、起動時に自動的にそれを読み込む際にエラーが生じているのかと思い、削除したり戻したりや、.RDataを同様にしたりしてもうまくいきません。
つきましては、
1)動作停止を回避する手段(エラー内容は文末に記載)
または、
2)起動時に.rや....
1
票
0回答
113 件の閲覧数
Rでデータフレームの日本語を0, 1に置き換える関数を自作した時のトラブル
よろしくお願いします。
Rで
のようなデータフレームの「あり」、「なし」を1, 0 に変更しようと思って以下のような自作の関数を書きましたが、もとのデータフレームに変更が反映されておらず、要素は「あり」、「なし」のままでした。
nihongosayonara <- function(df) {
for(i in 1:length(names(df))){
df[i][df[...
1
票
1回答
304 件の閲覧数
直交回帰直線から残差平方和を求める方法
直交回帰直線から残差平方和を求めたいです。
deviance(lm( d[,1] ~ d[,2] ))で線形回帰の残差平方和は求められるようなのですが、直交回帰の場合はどのように求めればよいでしょうか?
直交回帰直線は以下のように求めました。
r <- prcomp( ~ d[,1] + d[,2] )
slope <- r$rotation[2,1]/r$rotation[1,1]...
1
票
1回答
2,390 件の閲覧数
R言語におけるループでの複数のファイルの入力および処理
たびたび初歩的な質問で申し訳ありません。
下記質問では大変お世話になりました。
R言語におけるデータフレームから欠損値を除去する方法、欠損値しかない列の除去(もしくは0への変換)
現在、csv作成→Rへ入力→Mann-whitney計算→csv出力を、様々なご指導を元に下記にように実行しています。
Data <- read.csv("sample.csv", header=TRUE)
...
1
票
0回答
137 件の閲覧数
RcppでDataFrameのベクターの型を判定する方法
R,C++ 共に初心者です。
RcppでDataFrame内のベクターの型を判定する方法がわからずに躓いています。
的外れなコードで申し訳ございませんがご指摘をお願いします。
【環境】
Windows 10
RStudio 1.044
Rcpp 0.12.5
実験コード
// [[Rcpp::export]]
StringVector test(DataFrame df)
{
...
1
票
0回答
2,060 件の閲覧数
二つのデータフレームの値の違いを一覧で表示してくれる R関数はありますか?
Linuxで言うところの diff関数のようなものを探しています。データフレーム1 とデータフレーム2 は同じデータ構造ですが、一部の値が異なります。この異なる場所を R で特定したいと思います。良い関数/パッケージがあれば紹介して下さい。
1
票
0回答
249 件の閲覧数
R/Rscript/Rio コマンドライン上でのパッケージインストールについて
『コマンドラインではじめるデータサイエンス』
https://www.oreilly.co.jp/books/9784873117416/
の本のp127において
< data/iris.csv Rio -e 'skewness(df$sepal_length)'
とあります。skewnessはmoments パッケージに含まれてているのですが、
...
1
票
1回答
5,630 件の閲覧数
RStudio Presentationで文字化け
RStudioのPresentation機能を使ってスライドを作成しています。
OS Windows7
RStudio ver. 0.99.491(最新)
Global/Project OptionのEncodingはUTF8です。
文字コードの選択につき、以下のような問題に悩んでいます。
1.ソースのrpresファイルでSJISを選ぶと、Preview画面でフォルダ名しか表示されない
2....
1
票
0回答
908 件の閲覧数
Rで{e1071}パッケージを用いるとpredictで(subscript) 論理値添え字が長すぎますというエラーが出る件
表題の件で、ご質問させてください。
コードは下記の要領です。
いつもpredictの箇所でエラーになります。
どうしたらよいか、ご教授願います。
install.packages("e1071")
library(e1071)
data.train <-read.csv("predict_shop4_7_1.csv",header=T)
head(data.train)
data....
1
票
0回答
1,444 件の閲覧数
Rのエラーについて:以下にエラーeval(expr, envir, enclos) : オブジェクト “cv” は存在しません
> d.svm<-svm(cv~.,data)
以下にエラーeval(expr, envir, enclos) : オブジェクト "cv" は存在しません
上記のエラーはどのようにすればよいでしょうか?
ご教授願います。
1
票
0回答
383 件の閲覧数
R言語にアンドゥ・リドゥはございますでしょうか?
お世話になっております。
現在Windows7(64Bit)、VisualBasic(.NetFramework4.0)、R.NETを介してRとやり取りをするアプリケーションを作ろうとしています。
表題の件につきまして、所謂Rに対して行った作業について取消しややり直しをウィンドウズフォームアプリケーションで実現したく思っております。
R言語のコマンドやR....
0
票
0回答
39 件の閲覧数
分散が等しい時、1標本t検定と2標本t検定の違いについて
例えば、以下データのように、分散は等しく、平均値が違う場合です。
1標本t検定(0.002111)と2標本t検定(0.001053)の検定結果から、単純に1標本t検定の結果の方が、2倍高いと考えてよいのでしょうか。
どなたかアドバイスをお願いしてもよろしいでしょうか。
x1<-c(1,2,3,4,5)
x2<-c(6,7,8,9,10)
(m1<-mean(x1)) # 3
(...
0
票
0回答
30 件の閲覧数
0
票
0回答
51 件の閲覧数
入れ子構造になったデータのクラスター分析がしたい
ただいま、国勢調査等のデータを使って地域差の分析をしています。
そのなかで、以下のように入れ子になったデータがセットを作成しました。ほとんどのデータは町丁目まであるのですが、統計によっては一部市区町村までしかないデータ(財政力指数など)があるため入れ子構造になっており、このデータを2層そのまま使ってクラスター分析を行いたいです。
...
0
票
1回答
77 件の閲覧数
R でパッケージのインストールでエラーが表示されてしまう
パッケージをインストールしようとすると、以下の文が出てきます。
どうすればインストールできるようになるのでしょうか?
行った手順:
Rを開く→パッケージ→パッケージのインストール→このメッセージが出てくる
といった流れです。
パッケージの読み込みなどでも同じメッセージが出てきます。
エラーメッセージ:
utils:::menuInstallPkgs()
警告: リポジトリー https://...
0
票
0回答
39 件の閲覧数
Rで2つのデータフレームに含まれる緯度経度を同じグループID同士の中から近接結合させるには?
2つのデータフレームがあります(ddd, env0)。これらを緯度経度の近さに基づいて空間結合させたいと思っています。ただし、近い緯度経度は同じstormから選ばれるようにしたいです。
head(ddd)
storm lat long Mesh3rd
1 270 34.88750 137.2562 52372260
2 270 34.88750 137.2688 ...
0
票
0回答
31 件の閲覧数
Rvestを用いた複数のページからのtableデータの取り出し
Rでwebページからサッカーのtableデータを取り出し、csvファイルに出力したいと考えています。
一日だけのデータならば、以下のコードで成功しました。
html <- read_html()
data <-html %>%
html_table(fill=TRUE)
しかし、一ヶ月分のデータを得るために以下のコード
for (i in 1:31)
を用いた場合、...
0
票
1回答
35 件の閲覧数
Rstudioによるパネルデータ作成について
(従属変数)
A:1人あたりのゴミ排出量(日/g)
(説明変数)
a:人口 b:人口密度 c:平均世帯人員 d:核家族率 e:公営借家率 f:民営借家率
g:1人あたりの住宅延べ面積(㎡) h:対可住地農用地割合 i:対可住地宅地割合
j:財政力指数 k:100人当たりの事業所数 l:農林漁業事業所数 m:鉱業事業所数
n:卸売・小売業・飲食店事業所割合 o:...
0
票
0回答
35 件の閲覧数
Rstudioの共分散構造分析で、「parameterTableでない」というエラーが解消されません
やろうとしていること
パネルデータを共分散構造分析(SEM)で分析しようとしています(交差遅延効果モデル)
起こっているエラー
parameterTableとはなんでしょうか?
lavaan ERROR: model is a list, but not a parameterTable?
lavaan NOTE: missing column(s) in parameter table: [...
0
票
1回答
268 件の閲覧数
Rstudioにてlmer関数を使用したいが,lme4,lmerTest,mlmRevパッケージをインストールしても使用できません
Rstudioにて,線形階層モデルを行おうとしています。
昔,同様の処理をLme4パッケージのインストール,lmerライブラリの呼び出し,lmer関数の使用でできたと記憶しているのですが,lmerライラブリの呼び出しの時点でうまくいきません。
library("lmer")
library("lmer") でエラー: ‘lmer’ ...
0
票
1回答
191 件の閲覧数
Rstudioの共分散構造分析のエラーが解消されません(lavaan関数)
やろうとしていること
Rstudioのlavaanパッケージを用いて共分散構造分析をしたいです(交差遅延効果モデル)
用いているデータは2時点のパネルデータです
現在のコード
model2009 <-'
#internal, institutionという2つの潜在変数でINFLU, COMPLEX, PARTY, ELECTION, CONGRESSという5つの観測変数を予測しています(...
0
票
1回答
116 件の閲覧数
パスの\(エスケープ)を/に置換する方法を教えてください
パスが指定されたときに指定されたパスの(エスケープ)を/に置換したいと思ってます。
例えば
C:\R\data\covid
を
C:/R/data/covid/
というように置換したいです。
以下のプログラムのように、gsub関数を使ってパスの(エスケープ)を/に置換しようとしましたが、エラーになってしまいました。
gsub("\\","/","C:...
0
票
1回答
110 件の閲覧数
Xgboost 回帰でエラーが発生する
以下のコードを実行したところ、エラーが発生してしまいます。
エラーメッセージ
xgb.iter.update(bst$handle, dtrain, iteration - 1, obj) でエラー:
[23:11:58] amalgamation/../src/objective/regression_obj.cu:61: Check failed: preds.Size() == ...