タグ付けされた質問 「r」

これが私のRコードです。関数は次のように定義されます。 f <- function(x, T) { 10 * sin(0.3 * x) * sin(1.3 * x ^ 2) + 0.001 * x ^ 3 + 0.2 * x + 80 } g <- function(x, T, f=f) { exp(-f(x) / T) } test <- function(g=g, T=1) { g(1, T) } 実行エラーは次のとおりです。 > …

143 r 

ggplot2の凡例について質問があります。同じグラフに3本の線を描くことができ、3つの色を使用した凡例を追加したいと思います。これは使用されるコードです library(ggplot2) require(RCurl) link<-getURL("https://dl.dropbox.com/s/ds5zp9jonznpuwb/dat.txt") datos<- read.csv(textConnection(link),header=TRUE,sep=";") datos$fecha <- as.POSIXct(datos[,1], format="%d/%m/%Y") temp = ggplot(data=datos,aes(x=fecha, y=TempMax,colour="1")) + geom_line(colour="red") + opts(title="TITULO") + ylab("Temperatura (C)") + xlab(" ") + scale_y_continuous(limits = c(-10,40)) + geom_line(aes(x=fecha, y=TempMedia,colour="2"),colour="green") + geom_line(aes(x=fecha, y=TempMin,colour="2"),colour="blue") + scale_colour_manual(values=c("red","green","blue")) temp と出力 使用した3つの色と変数の名前（TempMax、TempMedia、TempMin）の凡例を追加します。私が試してみました scale_colour_manual 正確な方法を見つけることができません。 残念ながら、元のデータはリンクされたサイトから削除され、復元できませんでした。しかし、彼らはこのフォーマットの気象データファイルから来ました "date","Tmax","Tmin","Tmed","Precip.diaria","Wmax","Wmed" 2000-07-31 00:00:00,-1.7,-1.7,-1.7,-99.9,20.4,20.4 2000-08-01 00:00:00,22.9,19,21.11,-99.9,6.3,2.83 2000-08-03 00:00:00,24.8,12.3,19.23,-99.9,6.8,3.87 2000-08-04 …

143 r  ggplot2  legend  r-faq 

NAデータフレーム列の値の数を数えたい。私のデータフレームが呼び出されdf、検討している列の名前がそうだとしますcol。私が思いついた方法は次のとおりです： sapply(df$col, function(x) sum(length(which(is.na(x))))) これはこれを行うための良い/最も効率的な方法ですか？

142 r  dataframe 

data.frame（またはdata.table）で、NAを以前の最も近い非NA値で「転送」します。（の代わりにdata.frame）ベクトルを使用する簡単な例は次のとおりです。 > y <- c(NA, 2, 2, NA, NA, 3, NA, 4, NA, NA) 私はそのようなものfill.NAs()を構築することを可能にする関数を望みますyy： > yy [1] NA NA NA 2 2 2 2 3 3 3 4 4 この操作を多くの（合計〜1 Tb）小さなサイズdata.frame（〜30-50 Mb）で繰り返す必要があります。ここで、行はNAであり、そのエントリはすべてです。問題に取り組む良い方法は何ですか？ 私が作った醜い解決策はこの関数を使っています： last <- function (x){ x[length(x)] } fill.NAs <- function(isNA){ if (isNA[1] == 1) { isNA[1:max({which(isNA==0)[1]-1},1)] <- …

141 r  data.table  zoo  r-faq 

Rイントロでのスコープについて読み終え​​たばかりで、その<<-割り当てについて非常に興味があります。 マニュアルには、の1つの（非常に興味深い）例が示さ<<-れていました。私がまだ足りないのは、これがいつ役立つかというコンテキストです。 だから私があなたから読んでもらいたいのは、の使用<<-が興味深い/役に立つ場合の例（または例へのリンク）です。それを使用することの危険性とは何か（簡単に追跡できなくなっているように見える）、および共有したいと思うかもしれないヒント。

140 r  scoping  lexical-scope  r-faq 

Rでエラーを「スロー」するにはどうすればよいですか？データフレームといくつかの列名を受け取り、それらを処理する関数があります。列が存在しない場合、関数を停止し、それに依存するすべての関数を停止します。 私が見てきたrecoverとbrowseしてtracebackも、彼らは私が探しています何近いがないように思われたが、。

140 r  error-handling 

私が持っているとしましょう： v = rep(c(1,2, 2, 2), 25) ここで、各一意の値が出現する回数をカウントします。unique(v) 一意の値は何であるかを返しますが、その数は返しません。 > unique(v) [1] 1 2 私に与えてくれるものが欲しい length(v[v==1]) [1] 25 length(v[v==2]) [1] 75 しかし、より一般的なワンライナーとして:)次のように近いものですが（かなりではありません）： #<doesn't work right> length(v[v==unique(v)])

140 r  count  unique 

（私ができる）バイオコンダクターのhexbinを使用して、（png）表示領域全体を埋めるプロットを生成します-軸なし、ラベルなし、背景なし、ヌチンなし。

139 r  ggplot2 

ggplot2でy軸とx軸の原点/遮断を設定するにはどうすればよいですか？ x軸の線は正確ににあるはずy=Zです。 とZ=0または別の指定された値。

139 r  plot  ggplot2  axes 

ユーザーがキーを押すまでRスクリプトを一時停止したいのですが。 どうすればよいですか？

139 r  keypress  wait  readline  readkey 

したがって、私たちはすべてのRの新規ユーザーに「applyベクトル化されていません。パトリックバーンズRインフェルノサークル4をチェックしてください」と言っています（引用します）。 一般的な反射は、適用ファミリーで関数を使用することです。これは ベクトル化ではなく、ループ非表示です。apply関数の定義にはforループがあります。lapply関数はループを埋めますが、実行時間は明示的なforループとほぼ同じになる傾向があります。 実際、applyソースコードをざっと見てみると、ループがわかります。 grep("for", capture.output(getAnywhere("apply")), value = TRUE) ## [1] " for (i in 1L:d2) {" " else for (i in 1L:d2) {" これまでのところ、わかりましたが、見てみるとlapply、vapply実際にはまったく異なる画像が表示されます。 lapply ## function (X, FUN, ...) ## { ## FUN <- match.fun(FUN) ## if (!is.vector(X) || is.object(X)) ## X <- as.list(X) ## .Internal(lapply(X, FUN)) ## …

138 r  performance  loops  vectorization  apply 

2つのggplotsがあり、これをと水平方向に揃えますgrid.arrange。私はたくさんのフォーラム投稿を調べてきましたが、私が試したものはすべて、更新されて別の名前が付けられたコマンドのようです。 私のデータは次のようになります。 # Data plot 1 axis1 axis2 group1 -0.212201 0.358867 group2 -0.279756 -0.126194 group3 0.186860 -0.203273 group4 0.417117 -0.002592 group1 -0.212201 0.358867 group2 -0.279756 -0.126194 group3 0.186860 -0.203273 group4 0.186860 -0.203273 # Data plot 2 axis1 axis2 group1 0.211826 -0.306214 group2 -0.072626 0.104988 group3 -0.072626 0.104988 group4 -0.072626 0.104988 …

138 r  ggplot2  legend  gridextra 

現在のところ、この質問はQ＆A形式には適していません。事実、参考文献、専門知識によって回答が裏付けられることを期待していますが、この質問は、討論、議論、投票、または拡張ディスカッションを求める可能性があります。この質問が改善され、場合によっては再開できると思われる場合は、ヘルプセンターにアクセスしてください。 8年前に閉鎖。 ロックされています。質問はトピックから外れていますが、歴史的に重要であるため、この質問とその回答はロックされています。現在、新しい回答や相互作用を受け入れていません。 多くの場合、MATLABのライセンスの高さについて不満を言う人がいます。次に、なぜOctaveやRを使用しないのか、疑問に思います。しかし、後者は正しいですか？Rを使用してMATLABを置き換えることはできますか？

137 matlab  r  scientific-computing 

Rでそれほど大きくないdata.frameを作成しましたが、ビルドにかなり時間がかかります。それをファイルとして保存したいのですが、Rで再び開くことができますか？

137 r  dataframe 

以下のサンプルのように、周辺ヒストグラムを持つ散布図を作成する方法はありggplot2ますか？Matlabではそれがscatterhist()関数であり、Rと同等のものも存在します。ただし、ggplot2の場合は見ていません。 私は単一のグラフを作成することから始めましたが、それらを適切に配置する方法がわかりません。 require(ggplot2) x<-rnorm(300) y<-rt(300,df=2) xy<-data.frame(x,y) xhist <- qplot(x, geom="histogram") + scale_x_continuous(limits=c(min(x),max(x))) + opts(axis.text.x = theme_blank(), axis.title.x=theme_blank(), axis.ticks = theme_blank(), aspect.ratio = 5/16, axis.text.y = theme_blank(), axis.title.y=theme_blank(), background.colour="white") yhist <- qplot(y, geom="histogram") + coord_flip() + opts(background.fill = "white", background.color ="black") yhist <- yhist + scale_x_continuous(limits=c(min(x),max(x))) + opts(axis.text.x = theme_blank(), axis.title.x=theme_blank(), …

137 r  ggplot2  histogram  scatter-plot