【R言語】統計解析フリーソフトＲ 第４章【GNU R】
■ このスレッドは過去ログ倉庫に格納されています
0001１３２人目の素数さん
2011/01/09(日) 17:31:49統計計算で重宝するデータ型や、複数要素を処理する演算や関数、
解析結果を表示するグラフィックなど、多彩な機能を提供します。
●関連サイト
The R Project
http://www.r-project.org/
RjpWiki
http://www.okada.jp.org/RWiki/
リンク集
http://www.okada.jp.org/RWiki/?%A5%EA%A5%F3%A5%AF%BD%B8
0304１３２人目の素数さん
2012/01/15(日) 12:02:12.630305１３２人目の素数さん
2012/01/15(日) 19:09:49.95まぁ
r-project
r-script
r 統計
このあたりですね
0306１３２人目の素数さん
2012/01/16(月) 12:50:23.28まぁまずはRSiteSearch()だろ。
そんで、ttp://www.rseek.org/ とか。
最後に>>305の検索語でgoogleを使えばいいんじゃない？
○○を実装したパッケージを探すなら、「cran ○○」で検索する。
0307１３２人目の素数さん
2012/01/17(火) 00:15:09.670308１３２人目の素数さん
2012/01/17(火) 12:46:55.11言いたいことは、こういうこと？
正方行列Aが逆行列を持ち、A^{-1} の固有値が m,n なら、
A^{-1/2} = (A^{-1}+ mnI)/(m+n)
違ってたら、ごめん。
0309308
2012/01/17(火) 12:48:44.46A^{-1/2} の固有値が m,n なら、
A^{-1/2} = (A^{-1}+ mnI)/(m+n)
かな。
0310308
2012/01/17(火) 20:39:11.48>>308,309 が、しっちゃかめっちゃかだったので捲土重来を目指して…
A^t 転置、A^* 複素共役 として、
A が、A^t* A = A A^t* を満たし、A のすべての固有値が正だったら、
(U^t*)U = I を満たす行列 U と対角行列 D を使って、
A = UDU^{-1} と書けるから、
A^{-1/2} = UD^{-1/2}U^{-1}
（ただし、D^{-1/2} は、D の各対角成分d_{i,i}に対し、
D^{-1/2} の各対角成分 x_ {i,i} = 1/√d_{i,i} としたもの。）
という理解であってるのかな。
0311308
2012/01/17(火) 22:19:59.38R-Sources ttp://cse.naro.affrc.go.jp/takezawa/r-tips/r/20.html
見たけど、これで行けるんじゃね？
> A <- matrix(c(1,2,0,3), 2, 2)
> B <- eigen(A)$vectors %*% diag(1/sqrt(eigen(A)$values)) %*% solve(eigen(A)$vectors)
> solve(A)
[,1] [,2]
[1,] 1.0000000 0.0000000
[2,] -0.6666667 0.3333333
> B %*% B
[,1] [,2]
[1,] 1.0000000 0.0000000
[2,] -0.6666667 0.3333333
> A %*% B %*% B
[,1] [,2]
[1,] 1.000000e+00 0
[2,] 2.217736e-16 1
最後のがキッチリ単位行列にならないのは、
丸め誤差の所為だから仕方ないと思われ。
0312308
2012/01/17(火) 22:26:17.13> B <- eigen(A)$vectors %*% diag(1/sqrt(eigen(A)$values)) %*% solve(eigen(A)$vectors)
> eigen(A)$values
[1] 6.3722813 3.0000000 0.6277187
> solve(A)
[,1] [,2] [,3]
[1,] 0.8333333 -0.75 0.5
[2,] -0.5000000 0.75 -0.5
[3,] 0.1666667 -0.25 0.5
> B %*% B
[,1] [,2] [,3]
[1,] 0.8333333 -0.75 0.5
[2,] -0.5000000 0.75 -0.5
[3,] 0.1666667 -0.25 0.5
0313１３２人目の素数さん
2012/01/18(水) 10:57:06.79> A %*% B %*% B
[,1] [,2]
[1,] 1.000000e+00 0
[2,] 3.330669e-16 1
> zapsmall(A %*% B %*% B)
[,1] [,2]
[1,] 1 0
[2,] 0 1
0314308
2012/01/18(水) 12:32:51.26おお、そんな関数があったとは。
0315308
2012/01/18(水) 12:33:01.61おお、そんな関数があったとは。
0316308
2012/01/18(水) 22:59:15.60おお、そんな関数があったとは。
0317308
2012/01/18(水) 22:59:44.42おお、そんな関数があったとは。
0318308
2012/01/18(水) 23:01:26.21おお、そんな関数があったとは。
0319308
2012/01/18(水) 23:08:57.17おお、そんな関数があったとは。
0320308
2012/01/18(水) 23:35:40.76おお、そんな関数があったとは。
0321308
2012/01/19(木) 08:20:42.49おお、そんな関数があったとは。
0322308
2012/01/19(木) 08:22:06.35おお、そんな関数があったとは。
0323308
2012/01/19(木) 08:23:29.36おお、そんな関数があったとは。
0324１３２人目の素数さん
2012/01/19(木) 08:24:58.31おお、そんな関数があったとは。
0325１３２人目の素数さん
2012/01/19(木) 08:25:39.88おお、そんな関数があったとは。
0326308
2012/01/19(木) 08:26:49.60おお、そんな関数があったとは。
0327１３２人目の素数さん
2012/01/19(木) 16:46:26.30おお、そんな関数があったとは。
0328１３２人目の素数さん
2012/01/19(木) 18:16:15.47おお、そんな関数があったとは。
0329308
2012/01/19(木) 19:21:49.25おお、そんな関数があったとは。
0330１３２人目の素数さん
2012/01/19(木) 21:24:53.73あと、R sourceとかね。
0331１３２人目の素数さん
2012/01/19(木) 21:40:30.51おお、そんな関数があったとは。
0332１３２人目の素数さん
2012/01/19(木) 22:15:43.34おお、そんなsourceがあったとは。
0333１３２人目の素数さん
2012/01/19(木) 22:16:42.63おお、そんなsourceがあったとは。
0334１３２人目の素数さん
2012/01/19(木) 22:17:29.69おお、そんなsourceがあったとは。
0335１３２人目の素数さん
2012/01/19(木) 22:18:01.01おお、そんなsourceがあったとは。
0336１３２人目の素数さん
2012/01/19(木) 22:18:19.65おお、そんなsourceがあったとは。
0340１３２人目の素数さん
2012/01/19(木) 23:14:29.360341308
2012/01/20(金) 04:32:17.580342１３２人目の素数さん
2012/01/20(金) 05:22:51.540343313
2012/01/20(金) 20:36:44.68でも、実際の話、zapsmall()を知らないで、
条件分岐における数値型の同値の比較を
どうやって今まで書いていたのだろうと不思議に思う。
0344１３２人目の素数さん
2012/01/20(金) 20:42:23.02こいつ自分で質問の意味わかってのかいな？
0345１３２人目の素数さん
2012/01/20(金) 23:50:31.94他の言語と同じように書いてたよ。
0346１３２人目の素数さん
2012/01/20(金) 23:57:15.16abs(x-y) < 10^-14 とか
0347１３２人目の素数さん
2012/01/23(月) 09:35:12.7910*3の配列があり，1列目と2列目，1列目と3列目の計算し，
元の配列に戻したいばあい，どのように配列指定をしてやれば良いのでしょうか？
たとえば
x <- matrix(1:30, ncol=3)
で，
4列目に
x[,1]+x[,2]
5列目に
x[,1]*x[,2]
おそらく
x <- matrix(x[,1]+x[,2], x[,1]*x[,3], ???????)
で，「???????」を指定すればいいと思うのですが，，，，
よろしくお願いします
0348１３２人目の素数さん
2012/01/23(月) 09:52:39.90> cbind(x,x[,1]+x[,2],x[,1]*x[,2])
[,1] [,2] [,3] [,4] [,5]
[1,] 1 11 21 12 11
[2,] 2 12 22 14 24
[3,] 3 13 23 16 39
[4,] 4 14 24 18 56
[5,] 5 15 25 20 75
[6,] 6 16 26 22 96
[7,] 7 17 27 24 119
[8,] 8 18 28 26 144
[9,] 9 19 29 28 171
[10,] 10 20 30 30 200
こういうこと？
0349347
2012/01/23(月) 10:02:32.95>348
matrix()を使わなければいけないと悩んでいた私が馬鹿でした。
matrizだと横展開にしかならず苦悩しておりました。
0350348
2012/01/23(月) 12:02:36.25いや、cbind()かmatrix()かは好みの問題だと思うぞ。
> matrix(c(x,x[,1]+x[,2],x[,1]*x[,2]),ncol=5)
[,1] [,2] [,3] [,4] [,5]
[1,] 1 11 21 12 11
[2,] 2 12 22 14 24
[3,] 3 13 23 16 39
[4,] 4 14 24 18 56
[5,] 5 15 25 20 75
[6,] 6 16 26 22 96
[7,] 7 17 27 24 119
[8,] 8 18 28 26 144
[9,] 9 19 29 28 171
[10,] 10 20 30 30 200
0351１３２人目の素数さん
2012/01/23(月) 13:13:39.14RStudioご利用の方いましたら、ご一報お願いします
0352１３２人目の素数さん
2012/01/23(月) 13:28:06.82$ dpkg -l rstudio r-base |tail -2
ii r-base 2.14.1-1oneiric0 GNU R statistical computation and graphics system
ii rstudio 0.94.110 RStudio
問題なし。
0353１３２人目の素数さん
2012/01/23(月) 18:42:28.25河童さんが、par(mfrow=c(2,1))を使わずに、layout(matrix(1:2, 2))を使う理由って何だろう。
layout()だとなんかいいことがあるのかな？
0354308
2012/01/23(月) 21:15:31.430355１３２人目の素数さん
2012/01/25(水) 18:23:35.62例えば
1 2 3 4 5 6 7 8 9 10 11 12....
を
3 2 4 ....
という空白ごとに区切って
1 2 3
4 5
6 7 8 9
10 11 12....
としたいのです。
もし可能ならこのような場合，最初の文字列を配列に組み込むべきなのか，文字列として正規表現として処理するべきなのかどちらが好ましいのでしょうか？
よろしくお願いします
0356１３２人目の素数さん
2012/01/25(水) 19:06:09.47> f355 <- function(x,y){
+ s <- unlist(strsplit(x," "))
+ res <- vector()
+ j <- cumsum(y)
+ res[1] <- paste(s[1:j[1]],collapse='')
+ res[2] <- paste(s[(j[1]+1):j[2]],collapse='')
+ res[3] <- paste(s[(j[2]+1):length(s)],collapse='')
+ return(res)
+ }
>
> f355("1 2 3 4 5 6 7 8 9 10 11 12",c(3,2,4))
[1] "123" "45" "6789101112"
こんな感じか。
あとは、res[1]からres[3]のところをfor文と条件分岐で処理して、
汎用化すればよいと思う。
0357356
2012/01/25(水) 19:08:12.54> f355 <- function(x,y){
+ s <- unlist(strsplit(x," "))
+ res <- vector()
+ j <- cumsum(y)
+ res[1] <- paste(s[1:j[1]],collapse='')
+ res[2] <- paste(s[(j[1]+1):j[2]],collapse='')
+ res[3] <- paste(s[(j[2]+1):j[3]],collapse='')
+ res[4] <- paste(s[(j[3]+1):length(s)],collapse='')
+ return(res)
+ }
>
> f355("1 2 3 4 5 6 7 8 9 10 11 12",c(3,2,4))
[1] "123" "45" "6789" "101112"
0358356
2012/01/25(水) 19:36:19.47> f355.2 <- function(x,y){
+ s <- unlist(strsplit(x," "))
+ n <- as.integer(cut(as.integer(s),c(0,cumsum(y),Inf)))
+ sapply(1:length(unique(n)),function(x){s[n==x]})
+ }
>
> f355.2("1 2 3 4 5 6 7 8 9 10 11 12",c(3,2,4))
[[1]]
[1] "1" "2" "3"
[[2]]
[1] "4" "5"
[[3]]
[1] "6" "7" "8" "9"
[[4]]
[1] "10" "11" "12"
多分、識者ならsapply()部分をもっと効率よくできる気がする。
0359１３２人目の素数さん
2012/01/26(木) 01:27:31.67{
STR.SPLIT <- strsplit(x," ")[[1]]
## y = c(3, 5, 1) で c(1,1,1,2,2,2,2,2,3)を返す
GROUP <- inverse.rle(list(lengths = y,
values = seq_along(y)))
## 上記のGROUPにあわせてsplit()
split(x = STR.SPLIT[seq_along(GROUP)], f = GROUP)
}
## これなら数字の順番がめちゃくちゃでも、文字列でも行けるかも。spaceが続くとおかしくなるけど
f355.3(x = "2 ch R 2012 01 26 2 0 1 2 0 1 2 6 R", y = c(3, 5, 1, 2, 2))
0360355
2012/01/26(木) 07:03:57.07ご教授ありがとうございます。
配列に組み込んだものをループで処理することしか考えてませんでしたので，
以外な回答です。
プログラムに関しては勉強中なのですがなかなか上手くいかず悩んでおります。
ありがとうございます
0361355
2012/01/26(木) 08:28:33.03これまで統計処理に関してRを利用してきましたが，この度プログラミングもしっかり勉強したいと思うようになりました。
先にも記したようにプログラミングは初学者です。
CやPerlなど多様な言語がありますが，Rの場合どの系統のプログラミング教本で勉強を始めれば良いのでしょうか？
オススメ教科書があればよろしくお願い致します。
0362１３２人目の素数さん
2012/01/26(木) 10:23:38.00>>356だけど、正統なプログラミングの勉強をしたことがない。
プログラミングの勉強って、語学の勉強と似ていると思う。
一通り、文法を勉強したら、あとはひたすら経験を積むだけ。
人の表現を聞いて、あぁこういう場面ではこういう風に言えばいいのか、
と人まねを積み重ねる。
一通りの文法レベルの話なら、私はまだ未読だけど、
間瀬先生のRプログラミングマニュアルはどうかな。
間瀬先生の本なら、間違いはないと思う。
>>359
おぉ素晴らしい。rle()の存在を失念していました。
0363359
2012/01/26(木) 11:44:04.18Rデータ自由自在(Data Manipulation with R)がおすすめ。
Rプログラミングマニュアル(間瀬)は辞書的な感じ。
Rの基礎とプログラミング技法(石田訳)、プログラミングR(高階)は読み物的な感じ。
レイアウト的に後者の方が読みやすいか。
海外では最近出たThe Art of R Programmingがよいとされている様子。
0364１３２人目の素数さん
2012/01/27(金) 14:53:34.000365364
2012/01/27(金) 15:34:02.11すみません
0366１３２人目の素数さん
2012/01/28(土) 12:45:02.86sum(vector.name[p:q])
が一番普通か
0367１３２人目の素数さん
2012/01/30(月) 00:42:00.17(x,y,z)=
(0,0,0)
(1,1,2)
(3,2,4)
(0,9,2)
(2,0,4)
みたいに不規則な点列で構成された場合のコンタープロットを作成することは
出来ませんでしょうか？
Ｚの値を0〜１に正規化して、(x,y)の位置にgray(z)で塗りつぶした円を二次元プロットする
っていう手法でそれらしいものを作ってはみたものの、やっぱりイマイチなのでｗ
0368１３２人目の素数さん
2012/01/30(月) 11:00:34.81ぱっと思いつくのは、グリッドを作成して、(x,y,z)の値から、スプラインやkrigingなど
空間平滑法で各グリッド点における値を推定し、それで等値線を描く。
0369名無しさん
2012/01/30(月) 21:51:46.12ID:???krigingでググると
http://web.sfc.keio.ac.jp/~maunz/wiki/index.php?%C2%E807%B2%F3%B1%E9%BD%AC%A5%DE%A5%CB%A5%E5%A5%A2%A5%EB
っていうのが見つかりました。
ちょっとやってみます。どもでした。
0370名無しさん
2012/01/30(月) 23:07:03.68ID:???http://blog.livedoor.jp/dqnplus/archives/1692812.html
自分にはこんなのが精一杯でした．．．
library(ggplot2)
DF <- data.frame(age=c("10〜20代", "30代","40代","50代"), value=c(97,78,78,94))
ggplot(DF, aes(x=age, y=value, fill = age)) +
geom_bar(width = 1) +
scale_y_continuous(breaks = 0:100) +
coord_polar() +
labs(x = "", y = "") +
opts(legend.position = "none", axis.text.y = theme_blank(),axis.ticks = theme_blank())
0372名無しさん
2012/01/31(火) 00:39:21.68ID:???できるとしたら地図関係のパッケージだろうね。
0373名無しさん
2012/01/31(火) 10:47:16.46ID:???こんな感じじゃないの
i <- (1:360)/360*2*pi
cntr <- c(0.1,0.3)
v <- c(94,78,78,97)
arc <- cut(1:360,c(0,round(cumsum(v)*360/sum(v)),Inf))
plot(cos(i),sin(i),type="n",asp=1,axes=FALSE,xlab="",ylab="")
l <- levels(arc)
cols <- c("deepskyblue","dodgerblue1","dodgerblue3","brown1")
lbs <- c("50代\n94人","40代\n78人","30代\n78人","10〜20代\n97人")
cexs <- c(1.5,1.5,1.5,3)
for (j in 1:4){
k <- c(i[arc==l[j]],max(i[arc==l[j]])+1/360*2*pi)
polygon(c(cntr[1],cos(k)),c(cntr[2],sin(k)),
col=cols[j])
text(mean(cos(k))/1.5,mean(sin(k))/1.5,lbs[j],cex=cexs[j])
}
0374名無しさん
2012/01/31(火) 12:12:08.71ID:???ウジのグラフは、センターがずれてるだけじゃなくて「弧の長さの割合」も意図的に変えてるから、完全再現は無理みたいよｗ
参考 http://blog.kzfmix.com/entry/1327735061
そもそも、あんなおぞましいグラフを好んで描く必要ないと思うけど、こんなんでどう?
library(grid)
quartz(height=4, width=4)
#png(filename = "fuji-pie.png", width = 360, height = 360, units = "px")
num <- c(94, 97, 78, 78)
distortion <- c(1, 1.2, 1, 1) # <- !!
dist_num <- num * distortion
cum <- c(0, cumsum(dist_num))
rad <- cum / sum(dist_num) * 360 / 180 * pi
colors <- c("#0000FFAA", "#FF0000AA", "#000055AA", "#0000AAAA")
labels <- paste(c("50ies¥n", "10~20ies¥n", "30ies¥n", "40ies¥n"), num, " ps", sep="")
lab_cex <- c(1, 1.5, 1, 1)
lab_col <- c("white", "yellow", "white", "white")
div <- 200 # 円弧部分の分割数
dx <- 0 # 半径=0.5 として、x 方向への中心のズレ (-0.1 とかにすると、より強烈ｗ)
dy <- 0.1 # 半径=0.5 として、y 方向への中心のズレ
for (i in 1:length(num)) {
x <- c(0.5 + dx, 0.5 + 0.5 * sin(rad[i]), 0.5 + 0.5 * sin(seq(rad[i], rad[i+1], length.out=div)), 0.5 + 0.5 * sin(rad[i+1]))
y <- c(0.5 + dy, 0.5 + 0.5 * cos(rad[i]), 0.5 + 0.5 * cos(seq(rad[i], rad[i+1], length.out=div)), 0.5 + 0.5 * cos(rad[i+1]))
grid.polygon(x=x, y=y, gp=gpar(fill=colors[i]))
grid.text(x=0.5 + 0.4 * sin((rad[i]+rad[i+1])/2), y=0.5 + 0.4 * cos((rad[i]+rad[i+1])/2), label=labels[i], gp=gpar(col=lab_col[i], cex=lab_cex[i]))
}
#dev.off()
RjpWikiにも置いときました。
http://www.okada.jp.org/RWiki/?R%20%A4%C7%A5%A8%A5%F3%A5%BF%A5%E1#nf026bb6
0376373
2012/01/31(火) 14:45:06.35ID:???>>374さんは河童さんかと一瞬思ったけど、
そういえば河童さんはgridが嫌いだったから、別な人ですね。
0377375
2012/01/31(火) 23:04:16.60ID:???裏Wikiの方でも、がっつりとダメ出しされましたｗ
でも「裏」氏のファンなので、ちょっと嬉しかったりもしますｗ
0378374
2012/01/31(火) 23:16:34.60ID:???0379名無しさん
2012/01/31(火) 23:52:13.20ID:???0380１３２人目の素数さん
2012/02/02(木) 15:19:45.320〜6ぐらいまでの整数を要素にもつサイズが50のベクトルを作りたい。
できればポアソン分布に従って。
ここまでなら、
> rpois(50,2.5)
[1] 1 1 1 4 4 3 0 3 3 3 0 1 1 0 4 5 2 2 2 2 4 0 4 4 2 1 3 1 1 3 3 2 1 1 3 4 3 3
[39] 3 3 4 1 1 1 5 5 1 0 3 6
こんな感じでできるんだけど、このベクトルの合計をきっちり150にしたい。
固定したいのは、ベクトルサイズ50とベクトル合計値150だけで、lambdaとかは適当でいい。
当たり前だけど、合計値は変動する。
> sum(rpois(50,2.5))
[1] 132
> sum(rpois(50,2.5))
[1] 131
> sum(rpois(50,2.5))
[1] 118
こういうときはどうすればいいんだろうか。
よろしくお願いします。
0381１３２人目の素数さん
2012/02/02(木) 16:52:09.32無理やりの処理だが、こんなんじゃ駄目？
mypois <- function(n=50,sum=150,adjtimes=10){
myvec <- rep(0,times=n)
mysum <- 0
for(i in 1:adjtimes){
mysum <- sum(myvec)
myvec <- myvec + sign(sum-mysum)*rpois(n,abs(sum-mysum)/n)
myvec[myvec < 0] <- 0
if( abs(sum-mysum) == 0 ){ break }
}
return(myvec)
}
> x <- mypois(50,150)
> x
[1] 3 2 3 2 5 3 4 5 4 0 1 1 5 4 6 3 5 5 0 4 4 3 4 5 4 7 1 4 1 2 1 2 3 2 1 3 7 2
[39] 0 2 2 3 0 1 3 6 3 4 2 3
> sum(x)
[1] 150
0382381
2012/02/02(木) 17:00:53.33mypois <- function(n=50,sum=150,max=6,adjtimes=10){
myvec <- rep(0,times=n)
mysum <- 0
for(i in 1:adjtimes){
myvec <- myvec + sign(sum-mysum)*rpois(n,abs(sum-mysum)/n)
myvec[myvec < 0] <- 0
myvec[myvec > max] <- max
mysum <- sum(myvec)
if( abs(sum-mysum) == 0 ){ break }
}
return(myvec)
}
> x <- mypois(50,150,6)
> x
[1] 1 6 1 0 5 3 5 1 2 6 3 3 6 3 3 3 3 3 1 4 2 3 3 0 5 4 3 6 3 5 2 4 4 0 2 1 3 3
[39] 2 4 0 2 1 3 6 4 4 2 1 6
> sum(x)
[1] 150
0383381
2012/02/02(木) 17:33:39.030〜6からはみ出たデータを丸めたからかな。上手い方法がないものかねぇ。
> x <- mypois(50,150,6)
> x
[1] 5 2 1 3 4 4 3 4 3 0 5 5 3 3 1 2 4 4 4 2 5 0 3 4 2 2 3 4 4 1 6 1 1 5 4 4 2 4
[39] 2 5 1 3 2 0 5 2 5 4 1 3
> poissondist(hist(x,plot=F)$counts)
ポアソン分布への適合度の検定
data: hist(x, plot = F)$counts
X-squared = 9.7821, df = 4, p-value = 0.04426
sample estimates:
n lambda
50.00 2.06
0384１３２人目の素数さん
2012/02/02(木) 18:40:11.25おぉ、十分です。最大は「6ぐらい」なので７とかでもOKですから、
丸めはマイナスになるときだけ必要です。
for文のところはとても勉強になりました。
ありがとうございました。
0385１３２人目の素数さん
2012/02/06(月) 23:36:20.10一通り調べたのですが、分からない事があるので質問させてください。
○行ごとのplotを行いたい場合、どのようなコマンドを打てば良いのでしょうか？
例えば、
y1 y2...y10
x1
x2
.
.
.
x100
このような100行×10列のデータフレーム名Xの2列目と4列目をplotしたい場合、
plot(X$y2,X$y4)で1行ごとのplotは出来るのですが、
これをx1,x5,...x100と5行おきにplotしたいのです。
どうかご教授ください。
0386１３２人目の素数さん
2012/02/07(火) 00:25:49.86X1 <- X[(0:20)*5+2, ]
0387１３２人目の素数さん
2012/02/07(火) 08:46:05.39＞このような100行×10列のデータフレーム名Xの2列目と4列目をplotしたい場合、
＞plot(X$y2,X$y4)で1行ごとのplotは出来るのですが、
＞これをx1,x5,...x100と5行おきにplotしたいのです。
plot(X[seq(1, 100, 5), c(2, 4)])
でどう？
[行番号, 列番号]
で要素を指定
seq(from, to, by)で1,6,11.16...,100をつくる
0388１３２人目の素数さん
2012/02/07(火) 22:43:27.831,5,...,100 は、最初だけ4行の間隔ですね。
5,10,15,...,100 、あるいは、1,6,11,...,96 ですか？
でしたら、私のお薦めも >>387 と同じ方法。
0389１３２人目の素数さん
2012/02/08(水) 13:38:50.77なんであれがデフォルトなんだろう
0390１３２人目の素数さん
2012/02/08(水) 17:34:58.860391１３２人目の素数さん
2012/02/08(水) 20:31:18.89アンダースコア2度打ちで、アンダースコアを入力できるというのを
知った上で迷惑と行っているのだろうね。
そもそもRではアンダースコアは使用禁止だったので、
空いているアンダースコアにその機能が割り当てられたという歴史的経緯がある。
デフォルトの理由は、>>390の言うとおり便利で楽だかから。
いやな人は設定を変更すればいいし。
逆にターミナルでRを使っていて、<-のつもりでアンダースコアを入力して
イラッと来ることがある。
0392１３２人目の素数さん
2012/02/08(水) 22:46:01.26設定変えるためにググってたらrjpに飛んで、そこで初めて二度打ちやC-q _ でアンダースコアを入力できることを知った
アンダースコアを打つためにわざわざテキストモードに変更してた
よく調べもしないで愚痴言ってすいません
0393１３２人目の素数さん
2012/02/13(月) 12:25:46.70うまくいってるかたいますか？
0394393
2012/02/13(月) 14:44:06.36emacsのR側で一行でも改行がコマンドを実行すればコマンダー側のフリーズが直りました。
フリーズ自体は何台かのPCで試したところ起こる場合と起こらない場合がありました。
環境に依存するようです。
0395１３２人目の素数さん
2012/02/13(月) 22:00:37.21emacs ess で Rコマンダーを使うメリットって、何？
利点がまったく分からないのだけど。
0396393
2012/02/14(火) 00:04:31.33ご指摘のとおりで。
研究仲間がコマンダー使ってて、両方同時に使えると便利だから試しました。
まぁそうじゃなくても使えないよりは使えたほうが少しくらい便利かなと、、、
0397１３２人目の素数さん
2012/02/15(水) 23:54:50.83せっかくベクトル演算言語なので数え上げも高階関数で出来ればと思ったんですが。
0398１３２人目の素数さん
2012/02/16(木) 09:49:58.30http://stat.ethz.ch/R-manual/R-patched/library/base/html/funprog.html
プログラミング言語はよくわからないけど、これは参考になるだろうか？
0399１３２人目の素数さん
2012/02/16(木) 09:57:29.11たぶん Reduce()がfoldに相当するんじゃね
0402１３２人目の素数さん
2012/02/17(金) 05:50:16.64切片を取り除くには「-1」を指定するということはわかりますが、ここを好きな値に固定したいのですが・・・
0403１３２人目の素数さん
2012/02/17(金) 13:00:44.62ttp://r.789695.n4.nabble.com/R-Forcing-the-intercept-td812422.html
0404397
2012/02/17(金) 19:52:56.72おもしろいかどうかわかりませんが、
data <- as.integer(runif(10,1,21))
のようなデータがある時に、
predicate <- function(x) { x < 5 }
が真になるような項目を数え上げたいとします。
for文を使うと
ret <- 0
for(i in 1:length(data)){ if (predicate(data[i])) { ret <- ret + 1 } }
のような形になりますが、高階関数を使うと、
length(Filter(predicate, data))
と書けるようになります。
これは関数プログラミングの技法で、もう少し詳しく学びたいようでしたら
programming in haskellなんかが薄くていい本なんじゃないかと思います。
■ このスレッドは過去ログ倉庫に格納されています