初級Ｑ＆Ａ アーカイブ(7)
をテンプレートにして作成
[
トップ
] [
新規
|
一覧
|
検索
|
最終更新
|
ヘルプ
]
開始行:
COLOR(green){SIZE(20){初心者のための R および RjpWiki に...
新規投稿はできません
----
-[[初級Ｑ＆Ａ アーカイブ(6)]] (元記事が 2006-08-24 より 2...
-[[初級Ｑ＆Ａ アーカイブ(5)]] (元記事が 2005-03-13 より 2...
-[[初級Ｑ＆Ａ アーカイブ(4)]] (元記事が 2005-11-09 より 2...
-[[初級Ｑ＆Ａ アーカイブ(3)]] (元記事が 2005-05-02 より 2...
-[[初級Ｑ＆Ａ アーカイブ(2)]] (元記事が 2004-12-13 より 2...
-[[初級Ｑ＆Ａ アーカイブ(1)]] (元記事が 2004-08-03 より 2...
----
#contents
----
**凡例中に数式と数値変数を繰り返し書きたいのですが [#db8d...
>[[TAG]] (2007-06-27 (水) 21:50:40)~
~
使用環境:R version 2.5.0 (2007-04-23), i386-apple-darwin8...
~
#ref(legendmiss.png)~
凡例として~
「1 ミュー g／2 ミュー g／3 ミュー g」~
のように数字と数式（ギリシャ文字）を混ぜて書きたいと思い~
cols <- rainbow(3)
plot(sin, xlim = c(0, pi), ylim = c(-1, 1), ann=F, col=c...
a <- NULL
b <- NULL
a <- expression(paste(1, mu, "g"))
b <- c(b, a)
a <- expression(paste(2, mu, "g"))
b <- c(b, a)
a <- expression(paste(3, mu , "g"))
b <- c(b, a)
legend(0.3, -0.2, b, col = cols, lwd=1)
上のようにコードを書いてみたのですが（プロットの部分はダ...
面倒なので凡例に関する部分を~
a <- NULL
b <- NULL
for(i in 1:3)
{
a <- expression(paste(i, mu, "g"))
b <- c(b, a)
}
legend(0.3, -0.2, b, col = cols, lwd=1)
とすると、「iミューg／iミューg／iミューg」のように変数iが...
paste関数の中なので当然なのかもしれませんが何か良い方法は...
~
RjpWiki内のグラフィックス参考実例集：数式のプロットの例2に~
「数式と数値変数の結合」とあるので参考にしてexpressionの...
a <- NULL
b <- NULL
for(i in 1:3)
{
a <- substitute(paste(x, mu, "g"), list(x=i))
b <- c(b, a)
}
legend(0.3, -0.2, b, col = cols, lwd=1)
このようにすると「paste(1, mu, "g")／paste(2, mu, "g")／p...
数値変数は反映されたのですがその他が駄目になるという結果...
解決策はあるのでしょうか？~
どなたか対策をご存知でしたらご教授願えますでしょうか。~
//
- 今の場合なら，以下のようにしたらいかが？ -- &new{2007-...
cols <- rainbow(3)
plot(sin, xlim = c(0, pi), ylim = c(-1, 1), ann=F, col=c...
b <- paste(1:3, "μg")
legend(0.3, -0.2, b, col = cols, lwd=1)
- なるほど！全角文字を使ってしまえば良いのですね。盲点で...
- legend(0.3, -0.2, sapply(1:3,function(i)(eval(substitut...
- legend(0.3, -0.2, parse(text=sapply(1:3,function(i)past...
- さっそく試してみました。見事に上手くいきました！ありが...
- 一応付け加えておくと、こんなにややこしいのは、legendにv...
- 確かにtext, mtextだと何とかできたのですが、今回は自分の...
**Condor & R [#c0745713]
>[[HPC]] (2007-06-27 (水) 00:31:04)~
~
Condor というPCクラスター上での Rの動作例ってありますか...
//
**NAを含む行列を関数distで計算 [#c99db955]
>[[su]] (2007-06-25 (月) 13:34:43)~
~
使用環境：R version 2.5.0 (2007-04-23) i386-apple-darwin8...
NAを含む行列xyで行xと行yとの間の市街距離を計算しました。~
x = c(1,1,NA)
y = c(1,2,NA)
xy = rbind(x,y)
dist(xy[,-3], method="manhattan") # 3列目を除いた距離は...
dist(xy[,3], method="manhattan") # 3列目だけの距離は、NA
dist(xy, method="manhattan") # 3列目を含めた距離は、1.5
dist(xy, method="manhattan")の値が、1でもなく、NAでもなく...
//
- If some columns are excluded in calculating a Euclidean...
- 意訳すれば、NAを含む列は計算に使われないけど、その分距...
- ありがとうございました -- [[su]] &new{2007-06-25 (月) 2...
**データフレームからデータが取り出せない [#z9b0bbf4]
>[[K]] (2007-06-22 (金) 10:03:40)~
~
度数分布表の値をもとに処理をしようとしているのですが、デ...
> dbh_table <- table(data$DBH.Cls0)
> dbh_df <- data.frame(dbh_table)
> dbh_df
DBH.Cls0 Freq
1 18 5
2 22 10
3 24 35
4 26 45
5 28 30
(中略)
20 58 5
21 60 5
22 72 5
としてデータフレームdbh_dfを作成したのですが、2列目のデー...
> dbh_df[1,2]
[1]5
として欲しい"5"だけが得られるのですが、1列目のデータにつ...
> dbh_df[1,1]
[1] 18
22 Levels: 18 22 24 26 28 30 32 34 36 38 40 42 44 46 48 ...
となってしまい、必要な"18"の他に余分な1行が表示されてしま...
x <- dbh_df[1,1]; y <- dbh_df[1,2]
としても、xに余分な1行が入ってしまい、処理できずに困って...
どなたか、原因と対策をご教授いただけますでしょうか。
よろしくお願いします。
~
//
- 処理が出来ないとは具体的にはどういうことでしょう? 1列目...
- コメント、ありがとうございます。「処理ができない」と書...
- 自己レスです。as.vector(x)とすることで、期待する値を得...
- 18だけを取り出せているんですよ。ただし，18は数値ではな...
data.frame を作った後に， dbh_df[,1] <- as.numeric(as.vec...
- なるほど、Rはfactorのリストを知らせてくれていたわけです...
- ちなみに ?as.factor によれば、因子化された数値ベクトル...
**平行座標プロットの各ラインの選択 [#i0b82e1d]
>[[EDA]] (2007-06-21 (木) 19:01:15)~
~
以前、マウスによるポイントの選択方法が掲載されていまし...
//
- どのようなグラフで，どのように選択して，どのような情報...
parcoord の中を見て，マウスでクリックされた点から一番近い...
結構なプログラム量になりそう。 -- &new{2007-06-21 (木) 1...
**sub()で文字列を置き換えると，特定の文字で文字化けする [...
>[[奥村泰之]] (2007-06-21 (木) 13:38:26)~
~
御世話になっております。
下記のコードに示したように，sub()で文字列を置き換えようと...
文字化けする現象が出て困っております。
対応策を御存じの方，お教え頂けると幸いです。
> sessionInfo()
R version 2.5.0 Patched (2007-04-27 r41346)
i386-pc-mingw32
locale:
LC_COLLATE=Japanese_Japan.932;LC_CTYPE=Japanese_Japan.93...
LC_NUMERIC=C;LC_TIME=Japanese_Japan.932
> sub("AA", "山本", "AA氏")
[1] "山本氏"
> sub("AA", "五十嵐", "AA氏")
[1] "五女虫?x81"
//
- sub("AA", "五十??嵐", "AA氏"); でいいんですが、内部エ...
- ただのバグです. CP932環境で, sub("A","?u8868","A") がう...
- takahashi様，なかま様，ご回答を頂きありがとうございまし...
- 先程パッチ付きレポートを送りました. -- [[なかま]] &new{...
- 2.5.1以降は大丈夫. -- [[なかま]] &new{2007-06-27 (水) 1...
- なかま様。R-develにレポートまでして頂き，誠に感謝してお...
- version 2.5.1 で対応されていることを確認致しました。ま...
**散布図：同じような値を重ならることなく、横並びで表示し...
>[[ちょろ]] (2007-06-20 (水) 14:29:14)~
~
散布図を加工としております。~
しかし、plot関数ではデータが一直線に並んでしまい、同じよ...
一つ一つのドットが重ならないようにしたいのですが、どのよ...
イメージ的には以下のURLの一番上にある画像の左上部のような...
http://www.ne.jp/asahi/room/kuro/ScatterMakerReadMe.html~
~
ご教授いただければ幸いです。~
//
- library(gplots); ?space -- [[takahashi]] &new{2007-06-2...
- 宣伝しておきます。http://aoki2.si.gunma-u.ac.jp/R/dot_p...
- takahashi様、青木先生、アドバイスありがとうございます。...
- 似たようなものとしてはstripchart(Y~X,vert=T,method="sta...
**3次元のデータをX-Yと色分けでplotする方法 [#ga03839a]
>[[ぼう]] (2007-06-19 (火) 17:59:04)~
~
お世話になります。XYZの3次元のデータがあり、Zを色分けして...
~
mat <- matrix(rnorm(30),nrow=10)
plot(mat[,1],mat[,1],col=
ifelse(mat[,3]>0.6,2,
ifelse(mat[,3]>0.3,3,
ifelse(mat[,3]>0,4,
ifelse(mat[,3]>-0.3,5,
ifelse(mat[,3]>-0.6,6,
7
)
)
)
)
)
)
//
- データ数＝色数でよいのなら -- [[g]] &new{2007-06-19 (火...
~
mat <- matrix(rnorm(30),nrow=10)
plot(mat[order(mat[,3]),][,1],mat[order(mat[,3]),][,2],c...
- ありがとうございます。orderでzの順番にするとは妙技です...
- 値を階級に分けて，その階級に色を適用するのでしょうから...
mat <- matrix(rnorm(30), nrow=10)
m = 5 # 階級の個数
# suf は データの値により 1〜m になる
# 例ではいい加減に分類したけど findInterval の第二引数で...
suf <- findInterval(mat[,3], seq(min(mat[,3]), max(mat[...
plot(mat[,1], mat[,2], col=rainbow(m)[suf]) # rainbow で...
- じゃあ，無駄に大げさにして -- [[g]] &new{2007-06-19 (火...
~
plot(mat[order(mat[,3]),][,1],mat[order(mat[,3]),][,2],
col=rep(c("red","orange","yellow","green","blue","violet...
hist(mat[,3],breaks=c(min(mat[,3]),-.6,-.3,0,.3,.6,max(m...
$counts))
- お二方ありがとうございます。分野によっては頻繁に使いそ...
**plotでx軸とy軸を任意の座標で交わらせることはできますか...
>[[くに]] (2007-06-19 (火) 10:24:13)~
~
plotでx軸、y軸を(0,0)などの任意の座標で交わらせることはで...
#ref(kuni.png)
のように。よろしくお願いします。~
//
- axis 関数のヘルプを読むべし -- &new{2007-06-19 (火) 10...
> x <- runif(20, min=40, max=100)
> y <- runif(20, min=10, max=40)
> plot(x, y, axes=FALSE, xlim=c(0,100), ylim=c(0, 80))
> axis(1, pos=0)
> axis(2, pos=0)
- ほんとだ、helpにちゃんと書いてありますね。質問する前に...
**等高線グラフィックのスケールについて [#cc03ca46]
>[[V9]] (2007-06-13 (水) 14:31:14)~
~
たまたま等高線の話題なので、続けさせていただきます。~
contourで描いた等高線に、二次元の点をプロットしたく思って...
単なるcontourではスケールをぴったり合わせられるのですが、...
スケール（軸のメモリ）を合わせるのに、何かうまい方法はあ...
//
- 大変すみませんでした。具体的には、グラフィックス参考実...
- 説明不足で大変申し訳ございませんでした。実は、例３の中...
- 大変すみません。削除の方法を教えて頂けないでしょうか。 ...
- 実は、以下のようなものを書いて実行したら、エラーが出て...
- 上の方の「編集」をクリックして，編集画面の下の方に「テ...
ところで，何回も聞くけど ? filled.contour は読んだの？ -...
-以下のようなものを書いて実行したら、エラーが出てしまいま...
引数の長さが０とは、どういう意味なのでしょうか？
以下にエラーif (axes) { : 引数の長さが0です
filled.contour(
x,y,z,
levels=c(0,10,30,50,70,90,100),
xlim=c(-1,1),ylim=c(-100,100),key.axes=axis(4,seq(0,100,...
plot.title=title(main="main-title",xlab="x",ylab="y"),
key.title=title("ex"),
axes={axis(1);axis(2);points(exdata,pch=16,col="black")},
)
-exdataは、以下のような型式です。
x y
1 -0.2 20
2 0.8 40
3 0.1 30
・・・
- 描画できました。どうもありがとうございました。大変お手...
- ちなみに、プロットした点は、値によって色を変えることは...
- ちなみに、プロットした点は、値によって色を変えることは...
- plot.axis の { } の中を，目的に従って書けば良いだけ。 ...
- ありがとうございました。点の色分けはできたのですが、判...
- 判例の件、以下のようなものでうまく表示されませんでした...
filled.contour(
x,y,z,
levels=c(0,10,30,50,70,90,100),
xlim=c(-1,1),ylim=c(-100,100),key.axes=axis(4,seq(0,100...
plot.axes={axis(1); axis(2); points(exdata, pch=16, col=...
"deeppink", "black", "deepskyblue", "b...
legend(0, 0, paste("value", c("0-20%", "20-40%",
"40-60%", "60-80%", "80-100%")), bg="g...
col=c("blue", "deepskyblue", "black", ...
"red"), pch=16),
plot.title=title(main="main-title",xla...
key.title=title("ex")
)
- そこで、plot.axexの中で修文したところ、うまく表示できま...
**等高線用の外積計算 [#q55775cf]
> (2007-06-13 (水) 13:43:54)~
~
等高線を描くために外積を計算しようとしていますが、length...
何かうまい設定法がございましたら、ご教示頂けると幸いです。
> x<-seq(-4,4,length=20)
> y<-seq(-10000,4000,length=20)
> z<-outer(x,y,function(xi,yj) 100/(1+exp(-(0.279-1.0344...
+ )
以下にエラーouter(x, y, function(xi, yj) 100/(1 + exp(-(...
dim<- : dims [product 400] は object [20] の長さ...
//
- 基本的な間違いで、恥ずかしい限りです、、申し訳ありませ...
**無題 [#x0c3115a]
>[[あ]] (2007-06-10 (日) 20:00:15)~
~
sourceで時系列オブジェクトを読み込もうとしたら、下記のエ...
以下にエラーattributes(.Data) <- c(attributes(.Data), at...
不正な時系列パラメータが指定されました
このエラーはどのような意味なのでしょうか。~
読み込もうとしているデータは以前に読み込みを成功した事が...
//
- 無題とはいやはや・・・
読み込もうとしているデータも，変更の内容もわからないのに...
使用環境，問題の再現可能なデータセット，ロードしているラ...
ないと，誰も答えようがないと思いますが． -- &new{2007-06...
- 理由？不正な時系列パラメータを指定したからです！~
あ さん とは。如何に匿名とはいえ，もう少し考えたら？ -- ...
**ヒストグラムのｙ軸目盛（スケール？）の変更は可能でしょ...
>[[くに]] (2007-06-07 (木) 12:14:30)~
~
ヒストグラムのｙ軸目盛（スケール？）の変更は可能でしょ...
面積が0.3ヘクタールの調査区で樹木の胸高直径を測定しまし...
そのようなことは可能でしょうか。よろしくお願いします。~
//
- 元のデータを 1/0.3 倍して、そのヒストグラムを描けばいい...
- 早速のコメントありがとうございます。元のデータをxとする...
- 余計な御世話かもしれませんが、0.3ha毎の観測値が100個あ...
- ahaha,上のコメントで目が覚めた。私の提案していたのは，...
set.seed(77789)
x <- rnorm(100, mean=20, sd=5)
par(mar=c(4,4,1,4))
hist(x)
axis(4, at=0:4*10, labels=round(0:4*10*(1/0.3),1))
- 色々試行錯誤の末、少し求めているものに近づいてきました。
2cm刻みの直径階別本数を先に求めておいて、そのベクトルを...
#ref(upload.png)
でも、これではもう一息と言うところ。table関数で得られた...
>area=0.3
>layout(matrix(1:2),2)
>DBHclass <- round((data$x+0.05)/2,0)*2
>#15cm以上17cm未満は18cmという、2cm刻みの直径階にするた...
>brks <- seq(14,74,by=2)
>hist(DBHclass, breaks=brks,xlab="胸高直径階(cm)", ylab=...
>
>DBH_table <- table(DBHclass)
>plot(DBH_table/area,type="h",xlab="胸高直径階(cm)", yl...
>
>dbh<-data.frame(DBH_table/area)
>dbh
DBHclass Freq
1 18 5
2 22 10
3 24 35
4 26 45
5 28 30
6 30 55
7 32 45
8 34 60
9 36 55
10 38 90
11 40 75
12 42 65
13 44 85
14 46 30
15 48 10
16 50 10
17 52 10
18 54 15
19 56 5
20 58 5
21 60 5
22 72 5
このデータフレームから、forループか何かで c(18,18,18,18,1...
- x<-rnorm(1000);hist(x,axes=F);axis(2,0:4*50,0:4*500) --...
- x軸も一旦退避しておいてaxis()で書き込む -- &new{2007-0...
set.seed(77789)
x <- rnorm(100, mean=20, sd=5)
temp <- hist(x, breaks=seq(5,40,2), axes=FALSE, xlab="...
axis(2, 0:5*3, round(0:5*3*10/3,0))
axis(1, temp$breaks)
#ref(hist.nag.png)
- parを活用すれば, -- [[なかま]] &new{2007-06-08 (金) 10:...
set.seed(77789)
hist(x,yaxt="n")
axis(2,seq(par()$yaxp[1],par()$yaxp[2],5))
固定値は限りなく少なくなる筈.
- なるほど、
axis(2, 0:5*3, round(0:5*3*10/3,0))
としたところ、欲しいグラフが得られました。皆さん、ありが...
**和や部分リストを得る方法は [#k9c7b21b]
>[[kd]] (2007-06-06 (水) 10:16:48)~
~
リストを与えて，和や部分リストを得たいのですが方法はある...
> dwt<-dwt(x,n.levels=8, boundary="reflection")
> str(dwt@W)
List of 8
$ W1: num [1:300, 1] -0.188 0.166 -0.150 0.170 0.115 ...
$ W2: num [1:300, 1] 0.6416 0.5408 0.0322 -0.3626 -0....
:
:
次に，この dwt@W に対して和や部分リストを得たい思っていま...
dwt@W[1:4]
に対して
dwt@W[[1]]+dwt@W[[2]]+dwt@W[[3]]+dwt@W[[4]]
を得たり，
dwt@W[1:4]
に対して
list(dwt@W[[1]][1:150],dwt@W[[2]][1:150],dwt@W[[3]][1:15...
を得たいと思っています．~
//
- dwt ってのは，どのパッケージに入っているのか？x, n には...
そんなこと教えなくても答えられる人だけ答えてくれというの...
上の方に長々と書いてある注意事項は読んでないんだろうね。...
- 一個目はx<-list(1:3,4:6,9:12,13:15);do.call("+",x[1:2])...
- ただし一個目は"+"が二項関数なので３項以上あると問題だっ...
- そもそも，
dwt@W[[1]]+dwt@W[[2]]+dwt@W[[3]]+dwt@W[[4]]
ってやって正しい答えが出るならそれでいいじゃんか -- &new...
- ご教示いただきありがとうございます．forループで回す関数...
**指定データを除いたプロット [#o1a81cf5]
>[[kd]] (2007-06-03 (日) 15:05:17)~
~
下記のように，指定データを除いてプロットすることは可能で...
~
最初に下記を行いました．
> newx<-modwt(x,n.level=8,boundary="reflection")
> str(newx2)
Formal class 'modwt' [package "wavelets"] with 11 slots
..@ W :List of 8
.. ..$ W1: num [1:3061, 1] -1.82e-06 3.58e-05 -1.12e-...
:
:
..@ V :List of 8
.. ..$ V1: num [1:3061, 1] -2.14e-06 3.91e-05 1.48e-...
:
:
> plot(newx)
このplot(newx) を実行すると，x,newx@W[1],newx@W[2],...,ne...
次にnewx@V[8] だけを除いたプロットを行いたいので，下記の ...
> newx1<-????(newx)
> plot(newx1)
またnewx@W[8] だけをプロットすることも可能でしょうか？
> newx2<-????(newx)
> plot(newx2)
//
- まず str(newx2) は str(newx) の、x,newx@W[1],... は n...
- どうもありがとうございます．計算結果 y1, y2, ... を単に...
- par(mar=c(1, 1, 1, 1)) -- &new{2007-06-05 (火) 15:49:...
**locfit パッケージのロードに失敗 [#b56fb3c5]
>[[kd]] (2007-06-02 (土) 21:23:34)~
~
locfit を v. 2.5 for MacOSX 上でロードしようとすると，下...
要求されたパッケージ akima をロード中です
エラー： パッケージ 'akima' をロードできませんでした
追加情報: Warning message:
'akima' という名前のパッケージはありません in:
library(pkg, character.only = TRUE, logical = TRUE, lib...
locfit は下記からダウンロードしたものです．~
ftp://ftp.u-aizu.ac.jp/pub/lang/R/CRAN/bin/macosx/univers...
ちなみにwavelets等の他のパッケージは正常にロードできて動...
ftp://ftp.u-aizu.ac.jp/pub/lang/R/CRAN/bin/macosx/univers...
//
- 見たところ locfit パッケージは akima パッケージに依...
- この下の方にある「Macでパッケージvcdがロードできません ...
- どうもありがとうございます．「パッケージの依存関係」な...
- ずぼらしないで、まずは akima をインストールしたらど...
- 重ね重ねご教示をいただき本当にありがとうございます．aki...
**カレントディレクトリのファイル名一覧を取得 [#z6e16a93]
>[[mori]] (2007-06-01 (金) 10:41:10)~
~
カレントディレクトリに存在するファイル名の一覧をベクトル...
//
- dir <- dir(path=getwd(), pattern="file") のこと？ -- &...
- ありがとうございます。dir()でできたんですね・・・ -- [[...
- list.files() http://takenaka-akio.cool.ne.jp/doc/r_auto...
- dir と list.files は同じものです。 -- &new{2007-06-04 ...
> identical(dir, list.files)
[1] TRUE
**TIFFファイルをreadLinesで読み込みたい [#yb1adda9]
>[[akira]] (2007-05-31 (木) 20:22:00)~
~
測定器が出力するTiffファイルの1行目に測定条件が記入されて...
> x <- readLines("test.tiff", n=1)
Warning message:
'test.tif' に関する readLines で不完全な最終行が見つか...
画像ファイルをテキストで読み込む方法を教えてください。~
~
OS:WinXPSP2~
R version 2.5.0 Patched (2007-05-04 r41429)
i386-pc-mingw32
です。~
//
- そうですか。MacOS X 10.4 では読めました（何の追加情報も...
こんなファイルを
foo [9] > hexdump -c yellow.tiff
0000000 M M ?0 * ?0 ?0 002 220 301 ?0 301 ?0 ...
0000010 030 001 9 023 ?0 ] ] 037 ?0 253 z ) ...
読んでみました。
R での結果は以下のようになりました。
> x<-readLines("yellow.tiff", n=1)
> x
[1] "MM"
これがあなたの望む結果なんでしょうか。あなたが読んだ tif...
> sessionInfo()
R version 2.5.0 Patched (2007-05-18 r41630)
powerpc-apple-darwin8.9.1
- exifならそんなパッケージがあったような気がしますし, バ...
- ありがとうございます。LINUXでもうまくいきませんでした。~
ubuntu Linux~
R version 2.4.1 (2006-12-18) i486-pc-linux-gnu
本当なら、実際に使用している画像をお見せしたいのですが、...
[[デジカメ画像:http://cweb.canon.jp/camera/ixyd/810is/dat...
rgdalやrtiffのパッケージを用いるとtiffの数値情報をとれま...
> x <- readLines(dir()[1])
Warning message:
'810IS_l0.jpg' に関する readLines で不完全な最終行が見つ...
> x
[1]エラー：マルチバイト文字列が不正です
なお、測定器のTiffファイルにも同じように測定条件の文字列...
なかま様、これはバイナリなのでしょうか？コーディングの問...
- んと, JPEGやTIFFには[[EXIF(ぐぐれば沢山あるでしょう):ht...
- なかまさま、ありがとうございます。バッチ処理をしたいの...
> x <- readBin("test.tiff", what="character", endian="sw...
> x
[1] "II*"
で、scan()を試してみたら、"エラー：マルチバイト文字列が不...
> x <- scan("test.tiff", what="character")
Read 26373 items
> x
[1] "II*"
[2] ""
[3]エラー：マルチバイト文字列が不正です
これならいけそうです。~
これって、「[[ファイルTips大全の特定の文字コードのファイ...
- readBin("test.tiff", what="", n=1) でどうなりますか？わ...
- ありがとうございます。今はOSをubuntuで計算させているの...
- バイナリをバイナリで読むには, -- [[なかま]] &new{2007-0...
con<-file("binary.file",open="rb")
raw<-readBin(con,"raw",8)
raw # Hexで出る
# もし, rawの中がEUCで自分の環境がUTF-8なら"
iconv(rawToChar(raw),"EUC-JP","utf-8")
# めでたし,めでたし
close(con)
などとしないとダメでしょう.
- なかまさまありがとうございます。教えていただいた通り実...
test.tiffの1行目の"II*"から目的とする文字列の間にある制御...
- IIはTiffのフォーマットがLittleEndianと申しておりますか...
- ありがとうございます。生まれがwinなものでfileとかpipeと...
- [[libtiffライブラリ:http://www.libtiff.org/]]の[[tiffin...
libtiffはrtiffパッケージが使用しているのですが、tiffinfo...
- とりあえず、fileとscanだけでもほぼ望む結果が得られるこ...
con <- file(description=i, open="rb", blocking = TRUE)
x <- scan(con, what="logical", n=20, sep=NULL, comment.c...
どの辺が「ほぼ」かというと、^@のコードの直後の文字が取得...
**オンボードRAIDは使ってますか? [#te4364d3]
>[[通りすがりの佐々木]] (2007-05-31 (木) 16:22:29)~
~
R用PC(Athlon64,3200,2G,WXPSP2)をMBオンボードのRAIDを使っ...
//
- 昔家のパソコンで、シミュレーションを夏休み中ぶっつづけ...
- みなさんのRAID環境をしりたかったりする。Rを使うためには...
- 私は，RAID なんていうのは，今の今まで知らなかった。私は...
- RAIDなんてものを信用してはいけない。外付けHDに、こまめ...
- 性能的な物だと, 廉価版RAIDはOS+CPUに処理を丸無げするの...
- やはりオンボードRAIDは負荷がかかるのですね。実はRAIDに...
- その構成ならeSATAが使えればいいですね。NASだとEther自体...
- GentooのHPにSOFTraidについての翻訳がありました。Linux R...
- 実は、先日のクラッシュでHUBもお亡くなりになったので、ギ...
- GigaEtherの欠点はjumbo packetにNIC&HUBが対応していない...
**２変数関数の等高線プロットについて [#n6dc8ef2]
>[[GG]] (2007-05-31 (木) 11:55:06)~
~
２変数（ｘ、ｙ）関数の値（ｚ）を求めることなく、~
式の形から等高線を直接描ける方法はあるのでしょうか？~
例えば、z=1/(1+exp(-(10x-3y)))のようなグラフを直接３次元...
適当に等高線で色分けした後に、~
さらに重ねがきで幾つか点もプロットしたいと思っています。~
//
- 例示のようなきれいな関数なら，z=constant として，y=f(x,...
- 実際には例示よりももう少し複雑な式を考えておりますが、...
- y=もなくて，z も計算しなくて，コンピュータの中にこびと...
**日本語入力してもグラフが表示されない [#u4f683a9]
>[[のり]] (2007-05-27 (日) 20:39:34)~
~
Rインストール後、コマンド（R console）で日本語入力しても...
（例）
x=0:50
y=dbinom（x,50,0.25)
plot（x,y,type=’ｈ’,xlab=’x’,ylab=’ｙ’,main=’二項分布’）
＋
（注、type=ｈはヒストグラム表示、xlab、ylabはｘ、ｙ軸の名...
（使用環境 windowsXP）~
この解決方法をご存じの方がおりましたら、ご教示いただけな...
（尚、この質問は最初、初級者コースが閉鎖されていると認識...
//
- 張り付けられた文字と「日本語を入力しない他のコマンドで...
- type='h' の h までもが，全角じゃ〜〜ないですか。なんで...
x=0:50
y=dbinom(x,50,0.25)
plot(x,y,type='h',xlab='x',ylab='ｙ',main='二項分布')
#ref(baz.png)
- 有り難うございました。日本語に注意を取られすぎ、いろい...
- 本当の正解はこれなのでは？ -- [[okinawa]] &new{2007-05-...
x <- c(0:50)
y <- dbinom(x,50,0.25)
plot(x,y,type='h',xlab='x',ylab='ｙ',main='二項分布')
- あと、sessionInfo( )とコンソールから実行してみてくださ...
- c(0:50) は冗長でしょう。0:50 でよい。= を使って悪いわけ...
- c(0:50)は冗長だとは思いますが、R言語の作法に沿ってコー...
x <- 0:50
y <- dbinom(x,50,0.25)
plot(x,y,type="h",xlab="x",ylab="y",main="二項分布")
- 前にも投稿記事がありましたが、付値演算子として = を使う...
> order(x <- 1:10)
[1] 1 2 3 4 5 6 7 8 9 10
> x
[1] 1 2 3 4 5 6 7 8 9 10
> order(X = 1:10)
[1] 1 2 3 4 5 6 7 8 9 10
> X
エラー： オブジェクト "X" は存在しません
- order(x <- 1:10) がコードの簡略化に役立つとは思いません...
? c で，This is a generic function which combines its arg...
> identical(c(1:10), 1:10)
[1] TRUE
**Macでパッケージvcdがロードできません [#w5c8d728]
>[[ピーター]] (2007-05-25 (金) 18:46:30)~
~
Mac OS X 10.4.9, R 2.5.0 GUI 1.19（4308）を使っています~
「パッケージとデータ」のCRANからvcdをインストールし、「R...
要求されたパッケージ colorspace をロード中です
エラー： パッケージ 'colorspace' をロードできませんでした
追加情報: Warning message:
'colorspace' という名前のパッケージはありません
in: library(pkg, character.only = TRUE, logical = TR...
とあります。また、library(vcd)でも同様でした。~
一体なぜののでしょうか。教えてください。お願いします。~
//
- colorspaceというパッケージをインストールする必要がある...
- パッケージインストーラでインストールするときに，いつも...
- なるほど！colorspaceをインストールしたら、vcdもロードで...
**三次元プロットについて [#g5c55634]
>[[dmc]] (2007-05-23 (水) 13:42:22)~
~
scattaerplot3d関数で三次元プロットした点に、persp関数で二...
点と曲面をうまくひとつのグラフ内に描く方法がございました...
//
- perspだけでやっちゃう。 -- [[takahashi]] &new{2007-05-2...
x <- seq(-10, 10, length= 30)
y <- x
f <- function(x,y) { r <- sqrt(x^2+y^2); 10 * sin(r)/r }
z <- outer(x, y, f)
z[is.na(z)] <- 1
pmat<- persp(x, y, z, theta = 30, phi = 30, expand = 0.5...
ltheta = 120, shade = 0.75, ticktype = "det...
xlab = "X", ylab = "Y", zlab = "Sinc( r )")
p<-list(x=runif(100,-10,10),
y=runif(100,-10,10),
z=runif(100,0,5))
points(trans3d(p$x,p$y,p$z, pmat), col = 2, pch =16)
- ご回答頂き、どうもありがとうございました。無事に描画で...
- 追伸です。曲面を塗らず、透過させるにはどのようなコマン...
- perspの引数colを0とかNULLとか"transparent"とかってこと...
- 引数を０やtransparentにしてみたら、うまく描画できました...
**データがマイナスとなる個数 [#g2824fdf]
>[[port]] (2007-05-23 (水) 12:02:52)~
~
Ｒの超初心者です。~
10000個くらいのデータがあって，そのデータの中でマイナスと...
色々調べたのですが，よくわかりませんでした。~
どのたかご教授ください。~
//
- dat<-rnorm(10^4); sum(dat<0) -- [[takahashi]] &new{2007...
- どうもありがとうございます。助かりました。 -- [[port]] ...
- dat<-rnorm(10^4); length(which(dat<0))、でしょう。 -- [...
**Excelからのコピペが上手くいきません [#cc7998f8]
>[[みかん]] (2007-05-14 (月) 04:58:23)~
~
「Tips紹介」→「R にファイルを読み込む tips 集(暫定版)」→~
「Excel からコピーしたセル範囲を R に読み込む (何でも掲示...
に記載されております方法が上手くいきません。~
記載事項の通り、~
> from.excel <- function(nc)
+ {
+ matrix(scan(""), byrow=TRUE, nc=nc)
+ }
> from.excel(3)
としたのち、
エクセルで
1 F 4
2 F 3
3 M 2
という3×3のセル範囲をコピーし、貼り付けると、~
> from.excel(3)
1: 1F4
1: 2M3
以下にエラーscan("") : scan 関数は '$s' を期待したのに、...
> 3F2
エラー："3F2" に構文エラーがありました
>
というメッセージが出てきます。~
一方同じExcelのコピー範囲を一度テキストエディタに貼り付け...
> x <- read.table(stdin())
の後に貼り付けると、上手くいきます。~
しかしこの方法でも、Excelでコピーしたそのままを貼り付ける...
> x <- read.table(stdin())
0: 1F4
1: 2M3
2: 3F2
3:
>
となってしまいます。~
どうもExcelでコピーした列セル間の区切り子が、R上で認識さ...
どうしてこうなってしまうのか、ご教授をお願いいたします。~
当方WindowsXP上で、ver.2.2.1を使用しております。Excelのバ...
よろしくお願いいたします~
//
- エクセルで範囲選択してコピー後、read.table("clipboard")...
- Windowsね。x をプリントしてみたらどうなっているの？mac ...
> x <- read.table(stdin(),header=FALSE)
0: 1 F 4
1: 2 F 3
2: 3 M 2
3:
> x
V1 V2 V3
1 1 F 4
2 2 F 3
3 3 M 2
- 青木先生の方法ではmatrix関数を使用しているため~
すべての値がすべて数値あるいは文字列になっていないとうま...
~
そこで直接の回答ではありませんが、僕自身は
from.excel2 <- function(h = 0)
{
dat <- (read.delim("clipboard", header=h))
}
のような関数を別に定義してExcelのデータを読み込んでいます...
タイトル行がある場合は header=T とします。~
この場合にはデータフレームとして読み込まれるのでうまくい...
~
また、僕自身はstdin()を知らないです。-- [[あら]] &new{200...
- dat に代入するのではなく（dat って，関数外のオブジェク...
- 教えるつもりが逆に教えていただく型になってしまいました...
**直交実験を計画してくれる関数 [#v29e7c5d]
>[[TH]] (2007-05-13 (日) 09:21:44)~
~
直交実験を計画してくれる関数はどこかのRパッケージに含まれ...
//
- CRAN にあるパッケージ crossdes あたりは使えませんか。 -...
- [[Rseek:http://www.rseek.org/?cx=010923144343702598753%...
- ありがとうございます．でもcrossdesはラテン方格と釣り合...
- oa.design とそれに関する関数を用意して，その段階ではど...
- 移植しようと昨日からいろいろやっているうちに，全体の移...
**確率微分方程式 [#fa2c1ef1]
>[[SDE]] (2007-05-10 (木) 19:02:49)~
~
確率微分方程式で検索してもヒットしなかったのですが、~
解を求めるなど、確率微分方程式を取り扱ったパッケージは無...
Rの守備範囲だと思うのですが…~
//
- どうやって検索したのか？？site: r-project.org "Stochast...
- ありがとうございます。検索が雑でした。CRAN.packagesで一...
**Fizz-Buzz問題 [#x4e95cdc]
>[[X Jr.]] (2007-05-09 (水) 22:20:18)~
~
どのカテゴリに書こうか迷ったのですが、とりあえず初級Q&Aに...
「[[どうしてプログラマに・・・プログラムが書けないのか?:h...
//
- [どうしてプログラマに・・・プログラムが書けないのか?]そ...
> x <- 1:100
> print(ifelse(x%%15==0,"FizzBuzz", ifelse(x%%3==0,"Fizz...
[1] "1" "2" "Fizz" "4"
[5] "Buzz" "Fizz" "7" "8"
[9] "Fizz" "Buzz" "11" "Fizz"
[13] "13" "14" "FizzBuzz" "16"
[17] "17" "Fizz" "19" "Buzz"
:
[93] "Fizz" "94" "Buzz" "Fizz"
[97] "97" "98" "Fizz" "Buzz"
- やり方は幾つもあるけど，どれが一番「プログラマが書いた...
x <- 1:100
x[1:33*3] <- "Fizz"
x[1:20*5] <- "Buzz"
x[1:6*15] <- "FizzBuzz"
print(x)
- せめて，これを拡張して，1000までの素数をエラトステネス...
- cat の方がいいかな。それと ==0 の代わりに ! を使う -- ...
> x<-1:100;cat(ifelse(!x%%15,"FizzBuzz",ifelse(!x%%3,"Fi...
1 2 Fizz 4 Buzz Fizz 7 8 Fizz Buzz 11 Fizz 13 14 FizzBuz...
Fizz 19 Buzz Fizz 22 23 Fizz Buzz 26 Fizz 28 29 FizzBuzz...
Fizz 34 Buzz Fizz 37 38 Fizz Buzz 41 Fizz 43 44 FizzBuzz...
Fizz 49 Buzz Fizz 52 53 Fizz Buzz 56 Fizz 58 59 FizzBuzz...
Fizz 64 Buzz Fizz 67 68 Fizz Buzz 71 Fizz 73 74 FizzBuzz...
Fizz 79 Buzz Fizz 82 83 Fizz Buzz 86 Fizz 88 89 FizzBuzz...
Fizz 94 Buzz Fizz 97 98 Fizz Buzz
2文字だけ短いもの（！？）
x<-1:100;cat(ifelse(x%%15,ifelse(x%%3,ifelse(x%%5,x,"Buz...
- データを考える言語の立場から, 条件分岐はどうでしょうね....
M<-expression(1:100)
x<-eval(M)
x[!eval(M)%%3]<-"Fizz"
x[!eval(M)%%5]<-"Buzz"
x[!eval(M)%%3&!eval(M)%%5]<-"FizzBuzz"
- Rでのベクトルと行列の扱いなどの特徴を活かして(^_^;) -- ...
x<-1:105
dim(x)<-c(3,105/3)
x[3,]<-"Fizz"
dim(x)<-c(5,105/5)
x[5,]<-"Buzz"
dim(x)<-c(15,105/15)
x[15,]<-"FizzBuzz"
cat(x[1:100])
これも，ほんの少しだけ短くなる
a<-function(i,s){dim(x)<<-c(i,105/i);x[i,]<<-s}
x<-1:105
a(3,"Fizz")
a(5,"Buzz")
a(15,"FizzBuzz")
cat(x[1:100])
- いろいろと面白い例をありがとうございました。Rはループを...
- 人に聞く前に，あなたのプログラム例をまな板に載せてから...
- 初心者コーナーとあれば一番Rらしくない（笑）のもあげてお...
> for (i in 1:100) {
+ if(i%%15==0) cat("FizzBuz,")
+ else if (i%%3==0) cat("Fizz,")
+ else if(i%%5==0) cat("Buz,")
+ else cat(i,",")}
- そうですね。順序が逆になりましたが、初心者が書くとこん...
- いや，実際，挙げられた一番単純なものが一番確実に動く，...
- 私が採用係なら直前のコードを書いた人物（=> 実は私本人）...
[1] 2.744 0.000 2.757 0.000 0.000
[1] 2.737 0.000 2.748 0.000 0.000
[1] 2.712 0.000 2.723 0.000 0.000
[1] 2.676 0.004 2.698 0.000 0.000
[1] 2.628 0.000 2.633 0.000 0.000
[1] 2.632 0.000 2.649 0.000 0.000
- 私が採用担当なら, プログラマーになりたいと言う人は, 思...
- > 15,3,5 の順で条件判断すべきではなく、3,5,15の順で~
注意しないと，15 の倍数を3の倍数だと判断して（間違いじゃ...
- 実は，一番起こりやすいのは 3,5いずれの倍数でもないとい...
3,5,15の順で正しく判定するためには，かえって手間が掛かる...
ちょっと計算時間が短くなっても，バグ持ちプログラムでは意...
> a <- function(){
+ for (i in 1:1e7) {
+ if(i%%15==0) 15
+ else if (i%%3==0) 3
+ else if (i%%5==0) 5
+ else 0 }
+ }
> system.time(a())
ユーザ システム 経過
34.522 0.116 34.522
> b <- function() {
+ for (i in 1:1e7) {
+ if (i%%3==0 && i%%5!=0) 3
+ else if(i%%3!=0 && i%%5==0) 5
+ else if(i%%15==0) 15
+ else 0 }
+ }
> system.time(b())
ユーザ システム 経過
46.421 0.133 46.393
> bb <- function() { # こうすれば，無駄な比較をしないで...
+ for (i in 1:1e7) {
+ if(i%%3==0) {
+ if (i%%5==0) 15 else 3
+ }
+ else {
+ if (i%%5==0) 5 else 0}
+ }
+ }
> system.time(bb())
ユーザ システム 経過
29.060 0.138 29.072
- one-liner -- &new{2007-05-12 (土) 19:23:10};
replace(replace(replace(1:100,1:33*3,"Fizz"),1:20*5,"Buz...
**.RData の怪 [#b974396e]
>[[昔の名前で出ています]] (2007-05-09 (水) 16:13:25)~
~
R 2.5.0 をインストールして、使おうとしている途中です。~~
ディレクトリの設定をして、q() で「作業スペースを保存」し...
それでは、と、「作業スペースの保存」をして、できた .RData...
> sessionInfo()
R version 2.5.0 (2007-04-23)
i386-pc-mingw32
locale:
LC_COLLATE=Japanese_Japan.932;LC_CTYPE=Japanese_Japan.93...
;LC_NUMERIC=C;LC_TIME=Japanese_Japan.932
attached base packages:
[1] "stats" "graphics" "grDevices" "utils" "dat...
[7] "base"
> # この後、ワーキングディレクトリを設定
> getwd() # 設定の確認
[1] "C:/Documents and Settings/Owner/デスクトップ/r"
> # 作業スペースの保存
> save.image("C:??Documents and Settings??Owner??デスク...
> #ここで、いったん終了
> q()
デスクトップ??r??.RData をダブルクリックすると、なんと、R...
> sessionInfo()
R version 2.4.0 (2006-10-03)
i386-pc-mingw32
locale:
LC_COLLATE=Japanese_Japan.932;LC_CTYPE=Japanese_Japan.93...
;LC_NUMERIC=C;LC_TIME=Japanese_Japan.932
attached base packages:
[1] "methods" "stats" "graphics" "grDevices" "uti...
[7] "base"
//
- んと. explorerを起動して「ツール：フォルダオプション」...
- お返事ありがとうございました。インストール時の「Ｒを.RD...
- q()で終了したとき、イメージはそのときのワーキングディレ...
- q()で終了したとき イメージは ワーキングディレクトリの ....
単に作業スペースの変更の保存だけでなく，関数の定義なども...
私の Windows が，おかしいのかな。 -- &new{2007-05-10 (木...
- .RDataにはワーキングディレクトリは保存されません。.RDat...
- コメントありがとうございました。以下のようにまとめて宜...
(1) 現在の作業ディレクトリに .RData がない場合には，q() ...
(2) .RData がないときには，明示的に「作業スペースの保存」...
(3) 作業ディレクトリは .RData には保存されない。~
(4) .RData をダブルクリックして R を立ち上げると，立ち上...
(5) 既に .RData がある場合には，q() で「作業スペースを保...
- (1)は違います。(2)は(1)が違うので、別に明示的に作らなく...
- ああ，「一つの例外」に陥っていたのですね。 -- &new{200...
**hist で「'breaks' の数が無効です」とは？ [#waf79fcd]
>[[norion]] (2007-05-06 (日) 21:50:20)~
~
hist(5,x)というカキコミをして、（ｘにはあらかじめ数値列が...
'breaks' の数が無効です というエラーメッセージがでます。
5のところを500とかに変えても同じです。データの数は全部で2...
よろしくお願いします。~
//
- ?hist してみましょう。正しい引数を正しい順で与えていま...
- ?histは見ましたが、breaksというものの意味は（英語が得意...
-- [[norion]] &new{2007-05-06 (日) 22:17:27};
- hist(5,x) というようにしても，仰るようなエラーにはなり...
- 1行目に文字があったことに気づかなかっただけで、これを取...
- breaks とはデータ範囲をどのように分割（簡単にいえば棒と...
**複数ウインドウを開いて、描き出すウインドウを選択したい ...
>[[ななし]] (2007-05-05 (土) 13:35:36)~
~
x11() #ウインドウ一つめ
plot(1:10)
x11() #ウインドウ二つめ
plot(2:11)
このとき一つ目のウインドウにグラフを追記したいのですが可...
ひとつのウインドウを分割する場合、split.screenで分割して...
御存じの方教えてください。~
R version 2.4.0~
os linux-gnu vine/linux4.0~
//
- こんな感じの回答がお望みでしょうか。 -- [[aa]] &new{200...
x11() #ウインドウ一つめ
plot(1:10, main="1つめ")
x11() #ウインドウ二つめ
plot(2:11, main="2つめ")
dev.set(dev.prev()) ##1つめに移動
lines(1:10) #1つめに追記
dev.set(dev.next()) ##2つめに移動
lines(2:11) #2つめに追記
graphics.off()
- まさしくそれです。ありがとうございました。 -- [[ななし]...
**garchFit で「引数 "formula" がありません」というエラー ...
>[[たく]] (2007-05-01 (火) 12:46:50)~
~
gaarchFitで解析したいものがあり、試しにUKgasでやってみた...
どうすれば使えるようになるのでしょうか。
バージョンは2.4.1、OSはXPです。
> data(UKgas)
> library(fSeries)
> UKg.d<-diff(UKgas)
> UKg.m <-garchFit(formula.mean =~arma(2,0),formula.var=...
以下にエラー.modelSeries(fake = FALSE, lhs = TRUE) :
引数 "formula" がありませんし、省略時既定値もあり...
//
- その使い方はどこから探し出してきたのか？UKg.m <-garchFi...
- http://www1.doshisha.ac.jp/~mjin/R/0606_35.pdf
このサイトを参考にしました。ちなみにおっしゃるようなやり...
- どこかのpdfは見るけど，オンラインヘルプは見ていないので...
「おっしゃるようなやり方でやっても同じエラー」おかしいな...
> UKg.d<-diff(UKgas)
> UKg.m <-garchFit( formula = ~arma(2,0)+aparch(1,1),dat...
[1] "arma(2, 0)" "aparch(1, 1)"
Series Initialization:
ARMA model: arma
Formula mean: ~ arma(2, 0)
GARCH model: aparch
Formula var: ~ aparch(1, 1)
ARMA Order: 2 0
Max ARMA Order: 2
途中省略
Hessian Matrix:
mu ar1 ar2 omega
mu 0.4038139 -14.661207 -4.759985 0.005485700
ar1 -14.6612069 3721.675424 502.178082 -0.430102036
ar2 -4.7599849 502.178082 4486.891725 -1.201257840
途中省略
gamma1 -9.5579802 57.61092815 -4.6013946 -8.7055680
beta1 85.2086937 -4.60139464 240.7460409 -101.1148148
delta -22.7605912 -8.70556795 -101.1148148 58.0273266
--- END OF TRACE ---
Warning message:
計算結果が NaN になりました in: sqrt(diag(fit$cvar))
> summary(UKg.m)
Title:
GARCH Modelling
Call:
garchFit(formula = ~arma(2, 0) + aparch(1, 1), data = U...
途中省略
Information Criterion Statistics:
AIC BIC SIC HQIC
-9.893399 -9.693561 -9.903576 -9.812388
Description:
Tue May 1 15:54:32 2007 by user:
- 上のをコピペしてみたらできました・・・。さっきはどこか...
**INDSCAL(個人差MDS，３元データのMDS)の関数 [#v19fdb12]
>[[TH]] (2007-04-29 (日) 08:13:46)~
~
私が探した範囲ではRでは，多次元尺度法を行う関数として，以...
計量的多次元尺度法のための関数
パッケージ'stats'の関数 cmdscale
非計量多次元尺度法のための関数
パッケージ'MASS'の関数 isoMSD ，関数 sammon
パッケージ'vegan'の関数 metaMSD
しかし個人差を扱う，3元データ用のMDSモデルのINDSCALを実行...
//
- SensoMineR に indscal : Construct the Indscal model fo...
- 私の知識不足で，Napping data typeというのがいまひとつ分...
- 上切れました（失礼）私の知識不足で，Napping data typeと...
- 投稿は編集できるんですけどねえ。下のR/qtl も，自分でち...
- [[Rと多次元尺度法:http://www1.doshisha.ac.jp/~mjin/R/27...
- 便乗ですみません。MDSのALSCALって、Rでサポートされてい...
- こういう時は RsiteSearch でキーワード検索します。結果は...
- ALSCAL MDS R でググって，出てくるものを調べれば，http:/...
**R/qtl が load できません [#u95b6b2a]
>[[saoki]] (2007-04-27 (金) 23:29:18)~
~
R/qtlがloadできません。初めてRを使う初心者です。どなたか...
Mac OSX バージョン10.4.7のintel MacにR-2.50をインストール...
> library(qtl)
Error in dyn.load(x, as.logical(local), as.logical(now))...
共有ライブラリ
'/Library/Frameworks/R.framework/Versions/2.5/Resources/...
/qtl/libs/i386/qtl.so' を読み込めません
dlopen(/Library/Frameworks/R.framework/Versions/2.5/Re...
Library not loaded: /Library/Frameworks/R.framework/Vers...
Referenced from: /Library/Frameworks/R.framework/Versi...
Reason: image not found
以下にエラーlibrary(qtl) : .First.lib は 'qtl' に対し...
よろしく御願いいたします。~
//
- intel Mac が今ここにないので確かめられないのですが，iBo...
- ありがととうございます。iBook G4があったので（OSXバージ...
- 以下で解決しました。現在のR/qtl (qtl_1.05-2.tar)はR-2.4...
- おかしいですね。少なくとも iBookG4 の R2.5.0 で qtl_1.0...
> sessionInfo()
R version 2.5.0 Patched (2007-04-26 r41343)
powerpc-apple-darwin8.9.1
locale:
ja_JP.UTF-8/ja_JP.UTF-8/ja_JP.UTF-8/C/ja_JP.UTF-8/ja_JP....
attached base packages:
[1] "stats" "graphics" "grDevices" "utils" "dat...
[6] "methods" "base"
other attached packages:
qtl
"1.05-2"
> library(qtl)
> example(bayesint)
baysnt> data(hyper)
baysnt> ## Don't show:
baysnt> hyper <- subset(hyper, chr=c(1,4))
baysnt> ## End Don't show
baysnt> hyper <- calc.genoprob(hyper, step=0.5)
baysnt> out <- scanone(hyper, method="hk")
baysnt> bayesint(out, chr=1)
chr pos lod
c1.loc33 1 36.3 2.454921
c1.loc44.5 1 47.8 3.562391
c1.loc80.5 1 83.8 2.511073
以下略
- 再度行ったところ前述のエラーがやはり出て参りました。ses...
- おかしいですね. R2.6.0でも試しましたが，動きましたよ。...
R version 2.6.0 Under development (unstable) (2007-04-26...
Copyright (C) 2007 The R Foundation for Statistical Comp...
ISBN 3-900051-07-0
中略
'q()'と入力すればRを終了します。
[Workspace restored from /Users/***/***/***/.RData]
> library(qtl)
> library(help=qtl)
> 要するに，ここまでで何のエラーもないし，この後もちゃん...
- intel Mac でも，ok でしたよ。 -- &new{2007-04-29 (日) ...
- 手動で入れた、[/Library/Frameworks/R.framework/Versions...
**サンプルについて [#d5464f51]
>[[goto]] (2007-04-25 (水) 17:57:25)~
~
どなたか教えてください。~
以下のサンプルでは要素（球体）やノードを大きくしたり加工...
どのようにしたらこのようになるのでしょうか。~
gplot関数のパラメータを変更するのでしょうか。
http://erzuli.ss.uci.edu/R.stuff/sna/screenshots/sna_sam...
#ref(sna_sample3.png)
//
- タイトルが不適切。引用で画像のURLを書いたら表示されると...
- 説明不足の件は大変失礼致しました。snaパッケージを利用し...
- ?gplotとして、ヘルプファイルを熟読する。vertex.cex: exp...
- 引用される図があまりにも大きすぎて不都合なので，ダウン...
- 図の縮小ありがとうございます。gplotのパラメータを弄って...
- 自己レスです。gplot3という関数が用意されていました。解...
**Ｒでお絵かきツール [#e66bd6ae]
>[[ロトリング]] (2007-04-24 (火) 21:24:31)~
~
Ｒで出力したグラフに、ポリゴンやテキストをマウスで描画...
//
- > ポリゴンやテキストをマウスで描画~
その程度のことなら，テキストを指定して，描画位置をマウス...
- ご返答ありがとうございます。返事が遅れて申し訳ございま...
- ポリゴンにもよるけど，ある点がポリゴンの内か外かを判定...
- splancsパッケージにinpip()を使えば楽だと思います。マル...
- みなさん、情報ありがとうございます。ところで、最初の質...
- 単にこういうことを知りたかったとか？ -- &new{2007-04-2...
poly.xy <- locator(type = "l")
**95％信頼楕円 [#y1d9fea6]
>[[suz]] (2007-04-21 (土) 11:23:24)~
~
95％信頼区間楕円を描く方法として以下の操作は正しいでしょ...
x = rnorm(100, sd=1)
y = rnorm(100, sd=3)
plot(x, y)
library(ellipse)
r = cor(x, y)
stdev = c(sd(x), sd(y))
centre = c(mean(x), mean(y))
polygon(ellipse(r, scale=stdev, centre=centre, level=0.9...
//
- （元記事質問者と二番目のコメント投稿者は同じだと仮定し...
> V <- rbind(c(var(x), var(x,y)),c(var(x,y),var(y)))
> a <- ellipse(V, centre=centre, level=0.95)
> b <- ellipse(r, scale=stdev, centre=centre, level=0.95)
> identical(a,b)
[1] TRUE
- 注意：コメントがすでにあった場合、元記事を書き直すのは...
- 再質問します。-- &new{2007-04-22 (日) 12:47:38};
95％信頼区間楕円を描く方法として以下の操作は正しいでしょ...
set.seed(101)
x = rnorm(100, sd=1)
y = rnorm(100, sd=3)
xy = cbind(x,y)
library(ellipse)
stdev = c(sd(x), sd(y))
centre = c(mean(x), mean(y))
plot(x,y)
polygon(ellipse(var(xy), scale=stdev, centre=centre, l...
polygon(ellipse(var(xy), scale=c(1,1), centre=centre, l...
polygon(ellipse(var(xy), centre=centre, l...
polygon(ellipse(cor(x,y), scale=stdev, centre=centre, l...
polygon(ellipse(cor(x,y), scale=c(1,1), centre=centre, l...
polygon(ellipse(cor(x,y), centre=centre, l...
95％信頼区間楕円を描く方法として#1から#6のどれが正しいの...
#3と#4は同じ楕円を描きます。
scaleのデフォルトはc(1,1)ですが、#2と#3は同じ楕円を描きま...
- まあ，基本的には，その関数のオンラインヘルプをちゃんと...
- > ご指摘ありがとうございます。元記事です~
厳しく言えば，Wiki においては編集は無限の自由度があるのだ...
というか，直せば直すほど修正箇所が増えて，どうしようもな...
- ellipseのhelpを読み、exampleを実行してみました。答えの...
- 確率楕円で検索してみたのでしょうか？http://aoki2.si.gun...
- 「統計学関連なんでもあり」の過去ログ //
教えて下さいDensity Ellipse http://aoki2.si.gunma-u.ac.jp...
等確率楕円の長軸 http://aoki2.si.gunma-u.ac.jp/lecture/mb...
これは，確率楕円ですか？ http://aoki2.si.gunma-u.ac.jp/le...
二次元正規分布の確率楕円 http://aoki2.si.gunma-u.ac.jp/le...
散布図（各種の描画機能付き）散布図を描き，棄却楕円（確率...
**mfrowの値とplotmeansの動作 [#b8a84b1f]
>[[suzu]] (2007-04-20 (金) 16:32:43)~
~
グラフィックパラメータpar(mfrow=c(行数, 列数))の行数によ...
library(vegan)
data(BCI)
sp.n = specnumber(BCI)
group = gl(5, 10)
library(gplots)
par(mfrow=c(7,1)) # 行数によってplotmeansの動作が異なる
plotmeans(sp.n ~ group)
において
par(mfrow=c(6,1))
の場合、正常に動作しますが、
par(mfrow=c(7,1))
の場合、group 4のエラーバーが表示されず、以下のWarning me...
Warning messages:
1: 長さゼロの arrow は角度が不定ですので，スキップされま...
2: 長さゼロの arrow は角度が不定ですので，スキップされま...
エラーバーを正常に表示させるためには、どうすればよろしい...
使用環境 R version 2.4.0 (2006-10-03) i386-apple-darwin8....
参考文献 http://cse.naro.affrc.go.jp/takezawa/r-tips/r/68...
//
- 描画領域を大きくする。 -- [[takahashi]] &new{2007-04-20...
- アドバイスありがとうございます。「描画領域を大きくする...
**表形式でないファイルの読み込み [#qbd8961d]
>[[mori]] (2007-04-19 (木) 15:03:13)~
~
ファイルからデータを読み込む方法に関して質問があります。
[test.txt]
milk bread
bread butter
beer
milk bread butter
bread butter
上のようなtest.txtからデータを読み込んで下のように作った...
> x <- list(c("milk","bread"),
c("bread", "butter"),
c("beer"),
c("milk", "bread", "butter"),
c("bread", "butter"))
//
- read.table の fill パラメータを調べましょう。~
スマートではないが，以下のようなものも -- &new{2007-04-1...
x <- as.matrix(read.table("test.txt", header=FALSE, fill...
dimnames(x) <- NULL
apply(x, 1, function(a) a[a != ""])
こんな別解も
a <- sapply(readLines("test.txt"), function(s) unlist(st...
names(a) <- NULL
- もしくは類似解には -- &new{2007-04-19 (木) 16:28:59};
> sapply(as.list(readLines("test.txt")), function(i) str...
[[1]]
[1] "milk" "bread"
[[2]]
[1] "bread" "butter"
[[3]]
[1] "beer"
[[4]]
[1] "milk" "bread" "butter"
[[5]]
[1] "bread" "butter"
Warning message:
'test.txt' に関する readLines で不完全な最終行が見つかり...
注：上の warning は，ファイルの最後に EOF がないために出...
ちゃんとしたファイルなら warning は出ない。
- 類似解ですが，文字数が少ない(^_^;) -- &new{2007-04-19 ...
mapply(strsplit, as.list(readLines("test.txt")), split="...
- readLines()というのを知らなかったので大変勉強になりまし...
**garchFit(パッケージfSeries)のエラーについて [#d83a26c3]
>[[ケンゴ]] (2007-04-14 (土) 16:05:43)~
~
こんにちは、garchFit(パッケージfSeries)関数を使用し時系列...
おこないたいと思い以下のような処理を致しました。~
結果「エラー：!missing(data) is not TRUE」と表示されうま...
どこか間違っているのでしょうか？教えて下さい。~
使っているのはR2.4.1 OSはＸＰです。
> library(fSeries)
> UKgas
Qtr1 Qtr2 Qtr3 Qtr4
1960 160.1 129.7 84.8 120.1
1961 160.1 124.9 84.8 116.9
1962 169.7 140.9 89.7 123.3
1963 187.3 144.1 92.9 120.1
1964 176.1 147.3 89.7 123.3
1965 185.7 155.3 99.3 131.3
1966 200.1 161.7 102.5 136.1
1967 204.9 176.1 112.1 140.9
1968 227.3 195.3 115.3 142.5
1969 244.9 214.5 118.5 153.7
1970 244.9 216.1 188.9 142.5
1971 301.0 196.9 136.1 267.3
1972 317.0 230.5 152.1 336.2
1973 371.4 240.1 158.5 355.4
1974 449.9 286.6 179.3 403.4
1975 491.5 321.8 177.7 409.8
1976 593.9 329.8 176.1 483.5
1977 584.3 395.4 187.3 485.1
1978 669.2 421.0 216.1 509.1
1979 827.7 467.5 209.7 542.7
1980 840.5 414.6 217.7 670.8
1981 848.5 437.0 209.7 701.2
1982 925.3 443.4 214.5 683.6
1983 917.3 515.5 224.1 694.8
1984 989.4 477.1 233.7 730.0
1985 1087.0 534.7 281.8 787.6
1986 1163.9 613.1 347.4 782.8
> UKg.d<-diff(UKgas)
> UKg.m<-garchFit(formula=~arma(2,0)+garch(1,1),series=U...
エラー：!missing(data) is not TRUE
※条件付き平均（mean)=arma(2,0)~
※条件付き分散(variance)=garch(1,1)~
//
- 使ったことがないから適当なコメントですが、引数 formula ...
?garchFit によれば：
garchFit(formula.mean = ~arma(0, 0), formula.var = ~...
series = x, init.rec = c("mci", "uev"), delta = ...
cond.dist = c("dnorm", "dsnorm", "dged", "dsged"...
include.mean = TRUE, include.delta = NULL, inclu...
include.shape = NULL, leverage = NULL, trace = T...
algorithm = c("sqp", "nlminb", "lbfgsb", "nlminb...
control = list(), title = NULL, description = NU...
- アドバイスありがとうございます。 -- [[ケンゴ]] &new{200...
- で、アドバイスは役に立ったんですか、的はずれだったんで...
- どうも，質問の出所というか参考にしているところは同じら...
UKg.m<-garchFit(formula=~arma(2,0)+garch(1,1),data=UKg.d)
のようですよ（？） -- &new{2007-05-01 (火) 18:10:52};
**Ｒにインストールは必要ない？！ [#l0575daa]
>[[豊田秀樹]] (2007-04-14 (土) 07:03:17)~
~
信じられないようなことですが，ちゃんとインストールした~
Ｒのフォルダを単純にＵＳＢフラッシュメモリにコピーして~
おくと，そのＵＳＢフラッシュメモリを持ち歩くだけで，Ｒを...
していない計算機でＲを使えます．（そのフォルダは，全ての~
パッケージがインストールされており1ギガくらいの容量です）~
更にそのＵＳＢフラッシュメモリのＲののフォルダを単純に~
別の計算機にコピーするとその計算機でもＲが使えるようにな...
ここで「コピーで済むなら何のためのインストールだったん~
だろう？」という疑問が沸いてきました．きっとこの方法で~
コピーを繰り返すと「何か困ったことがそのうち起きる」ある...
「何か使えない高度な機能がある」としか思えません．~
文科系の学部なのでインストールはハードルが高い学生もいま...
コピーだけならだれでもでき，各学生に任せられます．~
とても便利なので，つい学生にも薦めて見たくなっていますが，~
あとで変なことが起きると困るので，どなたか「インストール~
しないと，こういう点で後で困るよ」というアドバイス~
がありましたら教えてください．それともインストールなんて~
必要ないのでしょうか？そんなことないですよね．~
//
- すみません環境・バージョンを書き忘れました．OSは Windo...
- 面白いですね。いわれてみればなるほどな、ですね。最近もV...
2.6 RをCDやUSBから起動できますか？
注意すればできます。基本的なRインストールは再配置可能（...
ですから、インストールしたRのイメージをハードディスクに...
フラッシュメモリUSBドライブ等のリムーバブル記憶装置にイ...
できます。（もしパッケージを私的なライブラリ（フォルダー...
いれば、それらの絶対パスがHTMLパッケージリストに記録され...
Rの実行には書き込み可能な一時ディレクトリとホームディレ...
で、もしそうしたものがなければ現在のディレクトリが使われ...
編成されたNTベースのWIndowsのバージョンでは何の問題も起...
もし問題があれば書き込み可能なフォルダーへのショートカッ...
ことはできないかも知れません。
- R は, 全てのOSで再配置可能です. 自分のトップ(R_HOME)が...
- 早速のアドバイスを，どうも有難うございました．了解いた...
- Rguiの引数に, R_LIBS=c:/temp 等を加えれば, パッケージの...
- 貴重な情報をありがとうございます．勉強になりました．計...
- なかまさんが言っているような「R解析キット」が分野別にあ...
- うーん, "「R解析キット」が分野別=>パッケージの追加が面...
- そうそう、私が言いたかったのは「R解析キット分野別作成キ...
- よいことを聞いたと、うれしくなり、とりあえずデスクトッ...
Warning messages:
1: 'lib.loc' 中に如何なるライブラリー木も見つかりません
in: library(package, lib.loc = lib.loc, character.only...
2: package "methods" in options("defaultPackages") was n...
起動準備中です ー 12 件の警告がありました (警告を見るに...
何かが足りないか（レジストリとかをいじらないといけないよ...
- .Rdataをダブルクリックしたまたは, コピーが不完全じゃな...
- お返事ありがとうございます。binディレクトリのRgui.exeを...
- 私の思い付く限りは, せいぜいみょーにフォルダ名の長いと...
- デスクトップの r というディレクトリに移しただけで，prof...
- 日本語ディレクトリはなるべくなら避けてほしいなぁ(希望)....
- 避けて欲しいというか，それが原因のようですね。Cドライブ...
**Windows版コンソールで使えるフォント [#sf37c914]
>[[mori]] (2007-04-12 (木) 14:32:35)~
~
Windows版のR2.4.1なのですが、GUIプリファレンスで指定でき...
//
- たとえば，Terminal とか FixedFont とか，これで，日本語...
- あれはリストボックスじゃなくって, コンボボックスだった...
- コンボボックスに直接指定すればダイジョブ。これは「HGS行...
#ref(ss2.png)
- アドバイスありがとうございます。リストにないフォントで...
**2変量散布図とヒストグラム図 [#uc13cad0]
>[[teru]] (2007-04-11 (水) 19:56:49)~
~
xとyの2変量散布図と，xおよびyのヒストグラムを書きたいので...
Ｒではどのようにすればいいのでしょうか？
**| oo
***| ooo
****| oooo
***| oo
**| oo
*| o
-------------
*********
*******
****
**
//
- 出来合いのものはないと思われるので，自分で関数を書くし...
pairs 関数の example を（実行して）見れば，望むものにちょ...
#ref(pairs.png)
- ヒストグラムではないですが、chplotパッケージがレイアウ...
- [[グラフィックス参考実例集：箱型図]]の「箱型図を軸の装...
- [[ヒストグラムと密度の推定]]の「二変量データのヒストグ...
- par()でmfrowでやるより，layout(matrix(c(1,1,4,2,2,3,2,2...
- 大変助かりました．みなさん，どうもありがとうございます ...
**ハッシュ [#x662ce31]
>[[JR]] (2007-04-11 (水) 14:24:09)~
~
Rにはperlなどのスクリプト言語でいうハッシュに相当するデー...
a => 1
b => 2
…
z => 26
というハッシュと~
v=c("a","b",…"z")
というベクトルがあったとき、vの各要素をハッシュのキーとし...
//
- こういうことでしょうか -- &new{2007-04-11 (水) 14:48:3...
> a <- 1:26
> names(a) <- letters
> a
a b c d e f g h i j k l m n o p q r s...
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19...
> a['b']
b
2
- a[names(a) == "b"] でないと最初に見付かったものしか選ば...
> x <- sample(1:100, 1e6, rep=TRUE)
> y <- sample(letters, 1e6, rep=TRUE)
> xx <- x
> names(xx) <- y
> xxx <- data.frame(x=x, y=y)
> xx["a"]
a
81
> length(xx[names(xx) == 'a'])
[1] 38482
> system.time(xx[names(xx)=="a"])
[1] 0.164 0.000 0.166 0.000 0.000
> system.time(xxx$x[xxx$y =="a"])
[1] 0.248 0.000 0.254 0.000 0.000
- ありがとうございます。ドンピシャなものは無いけど、ベク...
**区切り文字が1文字以上の空白で、#を含むNAに当たる文字が...
>[[TT]] (2007-04-07 (土) 17:30:52)~
~
こんにちは、以下のような、テキストファイルに入ったテーブ...
1 2 3 4 5
11 12 13 -1.#IND 15
これを、
tt<-read.table("hoge.txt")
で読み込もうとすると
'2' 行目には，5 個の要素がありません
というエラーが出ます。どうも2行目の4個目の要素の-1.#INDの...
tt<-read.table("hoge.txt",sep=" ")
とすると、今度は、1行目の1個目と2個目の要素の間の2個の空...
'2' 行目には，6 個の要素がありません
とエラーが出ます。何とか、1文字以上の連続した空白を1個の...
//
- hoge というのはどうも。以下のようにしたらいかが？さらな...
> x <- read.csv("foo.dat", header=FALSE, sep=" ", na.str...
> x
V1 V2 V3 V4 V5
1 1 2 3 4 5
2 11 12 13 NA 15
- read.csv は，read.table にシュガー・コートを着せたもの...
> x <- read.table("foo.dat", sep=" ", na.string="-1.#IND...
- 回答ありがとうございます。read.tableのオンラインマニュ...
> x <- read.csv("foo.dat", header=FALSE, sep=" ", na.str...
> x
V1 V2 V3 V4 V5 V6
1 1 NA 2 3 4 5
2 11 12 13 NA 15 NA
となってしまい、意図した結果になりません。もう少しご教授...
- 見れども見えず。read help. -- &new{2007-04-07 (土) 21:...
sep the field separator character. Values on each line ...
separated by this character. If sep = "" (the d...
the separator is “white space”, that is one or ...
newlines or carriage returns.
ーーーーーーーーーーーーーーーーーーーーーーーーーーーー...
> read.csv("foo.dat", header=FALSE, sep="", na.strings="...
V1 V2 V3 V4 V5
1 1 2 3 4 5
2 11 12 13 NA 15
- もしNA値に置き換えるべき項目いつも"-1.#IND" ではなく、...
> x0 <- scan("temp.txt", what=character(1)) Read 10 items
> x0
[1] "1" "2" "3" "4" "5" "...
[8] "13" "-1.#IND" "15"
> x1 <- gsub(".+#.+", "NA", x0)
> x1
[1] "1" "2" "3" "4" "5" "11" "12" "13" "NA" "15"
> x2 <- as.numeric(x1) # 直後のコメントとからもわかるよ...
> x2
[1] 1 2 3 4 5 11 12 13 NA 15
> x3 <- matrix(x2, ncol=5, by=TRUE)
> x3
[,1] [,2] [,3] [,4] [,5]
[1,] 1 2 3 4 5
[2,] 11 12 13 NA 15
> x4 <- as.data.frame(x3)
> x4
V1 V2 V3 V4 V5
1 1 2 3 4 5
2 11 12 13 NA 15
- そこまでやる必要はないのでは？~
要するに，数値を表す文字列以外を NA にするなら，以下のよ...
> x <- read.csv("foo.dat", header=FALSE, sep="", as.is=T...
> y <- as.data.frame(matrix(as.numeric(unlist(x)), nrow(...
Warning message:
強制変換により NA が生成されました
> y
V1 V2 V3 V4 V5
1 1 2 3 4 5
2 11 12 13 NA 15
無意味な空行はインデントを狂わすのでご注意
- ありがとうございました。オンラインマニュアル、まさしく...
**パッケージのインストール [#u6790519]
>[[tyamada]] (2007-04-06 (金) 18:13:53)~
~
教えてください。Rのパッケージ（たとえば、vegan)をインスト...
使っているのはR2.4.1 OSはビスタです。~
~
chooseCRANmirror()~
と打って、適当なミラーを選択し（日本の3つとも試しました）~
install.packages(c("vegan"))~
と打ちますと、はじめは調子よくいっているように見えるので...
~
URL 'ftp://ftp.u-aizu.ac.jp/pub/lang/R/CRAN/bin/windows/c...
ftp data connection made, file length 1176434 bytes~
開かれた URL~
downloaded 1148Kb~
~
以下にエラーzip.unpack(pkg, tmpDir) : ファイル 'C:/Progra...
~
などというメッセージが返されて、止まってしまいます。~
~
また、ローカルにあるZIPファイルからのインストールでも、同...
どこが間違っているのでしょうか？教えてください。~
//
- Vistaではユーザの権限管理がunix系並みに厳しくなっている...
- ありがとうございます。ご教示していただいた方法で、解決...
**ファイルタイプの識別 [#o10d8165]
>[[cam]] (2007-04-02 (月) 11:33:16)~
~
Rではsave関数でXDR形式のバイナリデータを作ることができま...
データファイルを読み込む時に、まずはファイル名からこのフ...
//
- ?data しかしこれはファイル拡張子でファイルタイプを識別...
- system 関数を使えばよろしいでしょう -- &new{2007-04-02...
> system("file test.XDR")
test.XDR: gzip compressed data, from Unix
> system("file test.data")
test.data: UTF-8 Unicode text
- ありがとうございました。参考になります。 -- [[cam]] &ne...
**作業スペースの日本語が文字化けしてしまいます [#kbd6c17f]
>[[akiba]] (2007-03-31 (土) 21:40:55)~
~
http://cran.md.tsukuba.ac.jp/から、R-2.4.1-win32.exeをイ...
文字化け部分をテキストに貼り付けると日本語としてちゃんと...
どうすれば、文字化けを直すことができるでしょうか。~
ご教示をお願いします。~
//
- まず [[日本語化掲示板]] の過去記事を読みましょう。そ...
- とりあえず直すだけなら「編集」メニューの「GUIプリファレ...
- GUIプリファレンスのFontを変えたら、直りました。ありがと...
**パッケージのインストールができません [#e02b12af]
>[[奥井正俊]] (2007-03-29 (木) 13:02:38)~
~
Rをインストールしたのちに、パッケージのインストールを試み...
//
- 現在、Aizu＆Tukubaにうまくつながっていないようです。CER...
- たぶん質問された方の環境だとプロキシ設定をしないとイン...
**Ｒの2軸プロット [#y1916d05]
>[[nick]] (2007-03-28 (水) 19:52:19)~
~
Ｒのプロットで，y軸の主軸，2軸の設定は可能でしょうか？~
色々調べたのですが，わかりません。~
どなたかご存知の方お願いします。~
//
- 例えば、plot(0,yaxt="n");axis(2,-2:2*0.1);axis(4,-4:4*0...
- お礼が遅くなりました。どうも有難うございます。 -- [[nic...
**プロットしたい行を読みとばしたい [#a23168a6]
>[[内藤]] (2007-03-28 (水) 15:41:16)~
~
Ｒの初心者です。~
プロットしたいデータが10000行あって，行を100行おきとかに...
gnuplotでいうとevery 100みたいな。~
以上，どうぞよろしくお願いします。~
//
- 例えば、 dat<-rnorm(10005); plot(dat[1:100*100]) -- [[t...
- どうもありがとうございます。助かりました。 -- [[内藤]] ...
**行列の1列目の中からある項目だけを取り出す [#z3fdd13e]
>[[norion]] (2007-03-27 (火) 23:56:46)~
~
他で作った表をRに読み込んで作業しようとしています。表全体...
行列の1列目には色々な名前が書いてあるんですが、この中から...
指定した名前に1列目が該当している行の全ての列を取り出した...
どのような関数を書けばそれが可能でしょうか。~
~
不適切な質問だったらすみません。~
//
- なんとなく行列ではなくて、データフレームではないかとい...
> M # 行、列ラベルつき行列
C1 C2 C3
R1 1 4 7
R2 2 5 8
R3 3 6 9
> M[,"C2"] # 列ラベルで第二列を取り出す
R1 R2 R3
4 5 6
> M["R2",] # 行ラベルで第二行を取り出す
C1 C2 C3
2 5 8
> str(M) # 行列なら次のように表示されるはず
int [1:3, 1:3] 1 2 3 4 5 6 7 8 9
- attr(*, "dimnames")=List of 2
..$ : chr [1:3] "R1" "R2" "R3"
..$ : chr [1:3] "C1" "C2" "C3"
> str(M) # データフレームなら次のように表示されるはず
`data.frame': 3 obs. of 3 variables:
$ C1: int 1 2 3
$ C2: int 4 5 6
$ C3: int 7 8 9
- 例えば，
> A
[,1] [,2] [,3]
[1,] "a" "1" "7"
[2,] "b" "2" "8"
[3,] "c" "3" "9"
[4,] "b" "4" "10"
[5,] "c" "5" "11"
[6,] "a" "6" "12"
> A[A[,1]=="b",]
[,1] [,2] [,3]
[1,] "b" "2" "8"
[2,] "b" "4" "10"
こんな感じでいいのでは？ -- [[g]] &new{2007-03-28 (水) 10...
- ありがとうございました。 -- [[norion]] &new{2007-03-28 ...
**関数が効かない [#peb71b5f]
>[[sunny]] (2007-03-26 (月) 23:01:23)~
~
初めて書き込みします。~
Rを使い始めてわずかの、全くの初心者です。~
barplotと入れても、表示をしてくれなくなってしまいました。~
どうやら、何かの弾みにbarplot関数が壊れてしまったようです...
壊れた関数を再設定するにはどうしたらよいのでしょう。~
//
- ls() としてみてください。もし barplot と表示されたら、...
- ls()をやってみましたが、その中にはbarplotはありませんで...
- > 表示をしてくれなくなってしまいました~
とありますが，エラーメッセージも何も表示されないのですか？~
引数なしで，barplot とだけ入力したら何か表示されますか？ ...
- Rをもう一度インストールしてみたらいかがですか。 -- &ne...
- エラーメッセージも何も表示されなかったんですが、インス...
**グラフのX,Yの刻み指定と目盛線の引き方 [#j05cfc33]
>[[のっち]] (2007-03-23 (金) 00:45:59)~
~
目的としては、CPUの負荷傾向を日々のCPU負荷状態(%)およびJO...
Rを入れて色々と調べているのですが、10, 20,30... といったY...
教えて君で申し訳ありませんが、ヒントなどいただけませんで...
//
- 目盛線については、plot(...,tck=1)で引くことが出来ました...
- abline 関数，および線種ならばその引数として渡せる lty ...
- 任意の刻みはaxisで指定しましょう。 -- &new{2007-03-23 ...
- 例として、~
plot(1:10, runif(10)*100, ylim = c(0, 100));abline(h=c(1...
で意図する刻みでの線が引けました。ありがとうございます。...
- axis のヘルプ見ました？第一引数は何だったか分かりません...
c(10,20,30,40,50,60,70,80,90)は1:9*10でよいね。 -- &new{...
- 線を引くのに悩みましたが、axis( 2, 1:9*10, lty="dotted"...
**小数以下が0を明示したいとき [#h7cb08dc]
>[[X Jr.]] (2007-03-21 (水) 21:22:39)~
~
例えば、~
y <- 1:10/2; plot(1:10, y); text(1:10, y, y, pos = 1)
とすると、小数以下が0のものは整数部分のみが表示されますが...
//
- ?format として, nsmall をみませう. -- [[なかま]] &new{2...
-
y <- 1:10/2; plot(1:10, y); text(1:10, y, format(y, nsma...
でできました。ありがとうございました。-- [[X Jr.]] &new{2...
- y <- 1:10/2; plot(1:10, y); text(1:10, y, sprintf("%.1f...
- なるほど。上の例で y <- 1:10/3 とした場合を試してみまし...
- ?formatとして, nsmallを見ると, digitsというのも(略 -- [...
**nlsで出力される回帰係数の読み込みについて [#q7aaab43]
>[[myuhe]] (2007-03-14 (水) 18:55:45)~
~
はじめまして。Rのド素人です。次のことについて、ご教示い...
~
nls()を実行すると、例えば次のような結果が得られます
Nonlinear regression model
model: inverseHeight ~ 1/(A * DBH^B) + 1/C
data: testtrees
A B C
0.2689986 3.1630292 24.6279807
residual sum-of-squares: 0.0003014474
この中の回帰係数（上の例では、0.2689986 3.1630292 24.6...
続けて行いたい処理は、出力された回帰係数で任意の独立変...
~
最初は、この結果をデータフレーム等に変換して、そこから...
~
predict()を使った方法も検討しましたが、任意の独立変数の...
~
初歩的なものかもしれませんが、よろしくお願いします。~
//
- 例えば、nls()が出力したnlsオブジェクトがhogeだとすると...
as.vector(hoge$m$getAllPars())
- coef(hoge) -- [[takahashi]] &new{2007-03-14 (水) 19:25:...
- 何事によらず R オブジェクトの構造を見たければ str() 関...
- coef()を使えば良かったのですね。str()もこれまでデータの...
- coef()を使った方法でできました。多くのアプローチまで教...
**作成したモデル（予測式）の保存・呼び出し方法 [#q17a4b19]
>[[ADMET modeler]] (2007-03-09 (金) 15:37:11)~
~
こんにちは。初めて投稿します。超初心者で非常に基本的な質...
~
仕事の関係で判別予測モデル（手法としてはBayesian, Random ...
~
たとえば、Random Forestの場合、
hlm <- read.csv('C:/work/data.csv',head=TRUE)
trainset <- hlm[1:1952,]
testset <- hlm[1953:2439,]
model <- randomForest(Class ~ ., data = trainset)
というふうにモデルを作成して、~
そのままpredictすると別のデータを予測できますよね。
pred <- predict(model, testset)
しかし、できればここで作成したモデルを保存しておいて、後...
~
この例では、Random Forestですが、一般的な回帰モデルなども...
//
- Rのオブジェクトをそのままディスク上に保存しておくには、...
- あるいは、終了時にそのまま「ワークスペースを保存」して...
- [[ファイルを読み込む tips 集(暫定版)]] で save.image() ...
**連続する整数をまとめる方法 [#s2982f8c]
>[[きゃらっと]] (2007-03-08 (木) 08:31:34)~
~
初めて投稿させていただきます．初心者です．~
(1:3)が(1,2,3)やとなる変換の逆を行わせる方法を教えてくだ...
例~
(1,2,...100, 250,251,...500, 10001, 10002,... 10005)~
変換後~
(1:100, 250:500,10001:10005)~
よろしくお願いします．~
//
- 関数を書くしかないのではないでしょうか？その関数は，以...
vecpack<-function(x) {
S<-x[-1]
E<-x[-length(x)]
TF<-E-S!=-1
s<-c(x[1],S[TF])
e<-c(E[TF],x[length(x)])
myfun<-function(x) {
if(s[x]==e[x])
as.character(s[x])
else
paste(as.character(s[x]),as.char...
}
vlist<-lapply(1:length(s),myfun)
paste("c(",paste(vlist,collapse=","),")")
}
早速付け加えてくれたようですね~
以下のように考えたのですが，ちょっと美しくないけど，せっ...
> set.seed(836478)
> ( x <- unique(sort(sample(100, 90, replace=TRUE))) ) #...
[1] 2 3 5 6 7 8 12 14 15 16 18 19 20 24 28 29 30 ...
[25] 40 42 45 49 50 52 53 54 55 59 61 62 63 65 66 68 69 ...
[49] 82 83 84 86 88 90 91 92 93 97 98 99
> n <- length(x) # ベクトルの長さ
> s <- c(TRUE, x[-n]+1 != x[-1]) # 新しいシークエンスが...
> res <- "" # 結果を得る文字列
> for (i in 1:n) { # いやらしいループで
+ if (s[i]) { # 新しいシークエンスが始まると文字列に「,...
+ res <- paste(res, x[i], sep=", ")
+ last <- x[i] # 最後に加えた整数値を覚えておく
+ }
+ if (i < n) { # シークエンスが終わるときは文字列に「：...
+ if (s[i+1] && x[i] != last) { # ただし，k:k のような...
+ res <- paste(res, x[i], sep=":")
+ }
+ }
+ else if (x[i] != last) { # 最後の要素の処理
+ res <- paste(res, x[i], sep=":")
+ }
+ }
+ res <- paste(res, ")") # 最後に「)」を付けて
> ( res <- sub(",", "c(", res) ) # 先頭が「,」で始まるの...
[1] "c( 2:3, 5:8, 12, 14:16, 18:20, 24, 28:31, 33, 35:40...
59, 61:63, 65:66, 68:74, 77, 80, 82:84, 86, 88, 90:93, 9...
vecpack の vlist lapply を使わなくてもできるので，より分...
vecpack2 <- function(x) {
S <- x[-1]
E <- x[-length(x)]
TF <- E-S != -1
s <- c(x[1], S[TF])
e <- c(E[TF], x[length(x)])
vlist <- ifelse(s==e, s, paste(s, e, sep=":")) # これだ...
paste("c(",paste(vlist,collapse=","),")")
}
巧い!
vecpack3 <- function(x)
{
s <- x[c(TRUE, x[-n]+1 != x[-1])]
e <- x[c(x[-n]+1 != x[-1], TRUE)]
return(paste("c(",paste(ifelse(s == e, s, paste(s, e, s...
collapse=","),")"))
}
マイナーな修正ですが，未定義の n とスペースを除いて
vecpack4 <- function(x)
{
s <- x[c(TRUE, x[-length(x)]+1 != x[-1])]
e <- x[c(x[-length(x)]+1 != x[-1], TRUE)]
return(paste("c(",paste(ifelse(s == e, s, paste(s, e, s...
collapse=","),")", sep=""))
}
- 早速のご返答ありがとうございました．早く皆さんのように...
**実験計画法 [#ia70d912]
>[[ちゃい]] (2007-03-05 (月) 21:56:05)~
~
Rで実験計画法ができるのでしょうか~
//
- 「実験計画法」という言葉で，あなたは何を意味しています...
実験計画を立ててくれというなら，それはたぶん無理なんじゃ...
- CRAN の貢献パッケージリストを最初から読む、または RSit...
- 有難うございました -- [[ちゃい]] &new{2007-03-06 (火) 2...
**integrate()で一様関数を積分できますか？ [#ie859bd7]
>[[藤巻十三]] (2007-02-28 (水) 17:52:22)~
~
こんにちは。~
表題の一様関数の積分は例えで、具体的にやりたいことは他に...
f<-function(x) 1
integrate(f=f,lower=0,upper=1)
とすると答えとして「１」を期待したいところでしたが、
以下にエラーintegrate(f = f, lower = 0, upper = 1) :
evaluation of function gave a result of wrong len...
となります。この原因と回避する方法を教えて頂けませんか？~
尚、「f<-function(x) ceiling(x)」なら答えは１が返ってきま...
「f<-function(x) max(x,1)」では駄目です。~
//
- すみません。わかりました。 -- [[藤巻十三]] &new{2007-02...
- なぜエラーになるかはおわかり（ヘルプ文章に書いてあるよ...
- そうじゃないです。fが返すベクトルの長さの問題です。f<-f...
- 失礼しました．なるほど、離散点での関数値を同時にサンプ...
**リスト内の各ベクトルから特定の要素のみを抽出する方法 [#...
>[[生物系大学院生]] (2007-02-26 (月) 06:11:57)~
~
こんにちは、初めて質問させていただきます。~
下記のようにして作成したブートストラップデータのリスト(z)...
> x<-c(1:100)
> y<-function(){sample(x,100,replace=TRUE)}
> z<-lapply(1:100,function(i)try(y()))
塩基配列データに基づくパラメータの推定値の分散をブートス...
~
各ベクトルについては、
> z[[1]][z[[1]]<11]
として11未満の数字のみを抜き出す方法や、
~
> sort(z[[1]])
として先頭にinformative siteを集める方法にはたどり着いた...
どなたかいい方法をご存知ないでしょうか？
//
- もう一歩。初心者ということならまずは for ループを使う...
> z <- lapply(1:100, function(i) sample(1:100, 100, rep=...
> A <- list(rep(NULL,100)) # 空リスト
> for (i in 1:100) A[[i]] <- z[[i]][z[[i]] < 11]
> B <- lapply(1:100, function(i) z[[i]][z[[i]] < 11])
> identical(A, B)
[1] TRUE
- 早速解法を教えていただきどうもありがとうございました。
新しいリストを構築するということまで考えがおよびませんで...
教えていただいた方法で思い描いていた結果が得られました。
どうもありがとうございました。 -- [[生物系大学院生]] &new...
**陰関数の表示について [#a6580f67]
>[[もも]] (2007-02-24 (土) 17:27:16)~
~
単発の質問で申し訳ないのですが~
陰関数を表示する関数、あるいはパッケージをご存知の方教え...
//
- 「陰関数を表示する」とはどういう意味ですか。 -- &new{2...
- すみません。自己解決しました。等高線＝０に気づいていま...
- そのノウハウを[[知っているといつか役にたつ(?)関数達(2)]...
**他言語で作成したプログラムのコンパイル方法について。 [#...
>[[Y.Takenaka]] (2007-02-18 (日) 16:52:55)~
~
mt19937の乱数の改良版が発表されていましたので、ラッパーを...
~
altivecやsse2と行ったものに対応していないものでは使えるよ...
http://sun.s167.xrea.com/Archives/HomePage.html (r_sfmt19...
//
- もう少し状況が分かるように説明した方がよろしいのでは？ ...
- コンパイルオプションを付けるなら,
R CMD COMPILE CFLAGS=-msse2 a.c b.c c.c
R CMD SHLIB -o abc.so a.o b.o c.o
こんな感じになります.
R CMD COMPILE -h
R CMD SHLIB -h
を見てください.-- [[なかま]] &new{2007-02-18 (日) 19:05:2...
- 情報ありがとうございます。2つめの方法でうまくいったです...
- Rで使えるようにと言う意味ならば, help(Random.user)する...
- なるほど、おっしゃるように関数のオーバヘッドは大きい事...
**自己回帰係数の算出方法 [#fb0da746]
>[[千成]] (2007-02-16 (金) 13:17:45)~
~
ARモデルの自己回帰係数の算出方法として、~
Yule-Walker法、最小二乗法、最尤法、Burg法などがありますが...
~
Yule-Walker法が他の方法に比べてどの点で秀でているのでしょ...
助言をお願い致します。~
//
- ここで尋ねるような話では無いですね。時系列の本を（何冊...
- やはりそうですか…。ありがとうございます。 -- [[千成]] &...
**整数の割り算について [#qa6ada38]
>[[廣瀬_敏之]] (2007-02-15 (木) 02:16:56)~
次のような結果をどう考えればいいのでしょう？
98-BASICの整数の割り算の結果とあまりに違うのでとまどって...
どちらが正しいか、という問題でない事はわかっていますし、...
各々のround(x,0)をとってから割り算をしていると考えても少...
この計算の定義までマニュアルに書いてあったでしょうか？
いわゆる『Ｒ言語のマニュアル』と、過去の『初心者向けＱ＆...
使用環境はＤＯＳ／Ｖ系のパソコンで、『ＸＰ』と『Ｖｉｓｔ...
R version 2.2.0, 2005-10-06, i386-pc-mingw32
attached base packages:
[1] "methods" "stats" "graphics" "grDevices" "uti...
[7] "base"
です。本当は複数のバージョンでチェックすべきなのでしょう...
なお『Ｒ』自体のソースコードをどうぞ、と言われましても、...
> -2.3%/%1
[1] -3
> test<-c(-5.99999,-5.00001,5.00001,5.99999)
> test%/%1
[1] -6 -6 5 5
> test%/%2
[1] -3 -3 2 2
> test%/%2.5
[1] -3 -3 2 2
> test%/%2.2
[1] -3 -3 2 2
> test%/%2.7
[1] -3 -2 1 2
> test%/%-1
[1] 5 5 -6 -6
> test%/%-2
[1] 2 2 -3 -3
> test%/%-2.5
[1] 2 2 -3 -3
> test%/%-2.2
[1] 2 2 -3 -3
> test%/%-2.7
[1] 2 1 -2 -2
- もうしわけありませんでした。helpで分かるそうです。 -- [...
- 演算記号に対してまでhelpが使えるとは思いませんでした。...
- %/% はいわゆる整数商ですが、98-basic （懐かしい響き！）...
**plot()のy軸ラベルと目盛りの表示について [#ff2f6052]
>[[kokko]] (2007-02-14 (水) 18:22:11)~
~
x <- 1:5
y <- c(100,1000,10000,100000,1000000)
plot(x,y, log="y", xlab="X_num", ylab="Y_num", las=1)
こんにちは。このようにy軸を対数目盛にしたグラフを作製した...
また目盛りの値が "1 e+05" のように1とeの間にスペースが入...
どなたかよろしくお願いします。~
//
- ＞ y軸のラベル "Y_num"を時計回りに45°回転させてx軸と平行~
時計回りに 90 度回転ではないですか？それはともかく，text ...
＞ "1 e+05" のように1とeの間にスペースが~
お使いの R のバージョンと，OS の種類およびバージョンは？~
私は R 2.5.0 で Mac ですが，スペースが入っているようには...
もし，どうしても気になるということなら，axis 関数で，文字...
図や表の体裁も，自分ならこのように描く（書く）のにどうし...
- ご回答ありがとうございます。余白に書き込むためにはmtext...
**Windows版 help() の表示方法の変更 [#pbadd398]
>[[X Jr.]] (2007-02-03 (土) 21:31:51)~
~
R-2.4.0以降のWindows版で、インストール時にヘルプの表示形...
//
-Rprofile.siteと言うファイルがR_HOMEのetc以下にあると思い...
- 遅くなりましたが、ありがとうございました。私はヘルプをH...
**wilcox.exact で，P 値が 1 より大 [#ica1e20e]
>ウィリアム (2007-02-02 (金) 07:57:57)~
~
Wilcoxon signed rank test を行おうと思いましたが、同じ値...
> wilcox.exact(c(8,8,-5,-10))
Exact Wilcoxon signed rank test
data: c(8, 8, -5, -10)
V = 5, p-value = 1.125
alternative hypothesis: true mu is not equal to 0
P 値が 1 より大きくなってよいのですか？~
//
-貢献パッケージ中の関数(ですよね)に言及する時はパッケージ...
-すみません。書くのを忘れました。 >library(exactRankTest...
-質問の意図が、バグ報告なのか、どうしても計算したいのに困...
**指定されたベクトルのサイズが長すぎます [#l8bf0ae9]
>[[szk]] (2007-02-01 (木) 09:03:07)~
~
以下のエラーを解決する方法をお尋ね申し上げます。
> dist(rnorm(560327))
以下にエラーvector("double", length) : 指定されたベクト...
以下は正常に計算されます。
> dist(rnorm(560327)[1:10])
//
-dist(1:10) が何を返すかは理解できていますか?~
dist((1:10)[1:3]) とはどこが違うか分かりますか?~
あなたが要求した操作には（実数を格納するには）560327×5603...
そして，あなたのコンピュータはそれだけのメモリを使えます...
-2.3テラバイト! 東工大の誇るグリッドスパコン TSUBAME の総...
-1.2Tぐらいかと. TSUBAMEは共有型クラスタではなかったと思...
終了行:
COLOR(green){SIZE(20){初心者のための R および RjpWiki に...
新規投稿はできません
----
-[[初級Ｑ＆Ａ アーカイブ(6)]] (元記事が 2006-08-24 より 2...
-[[初級Ｑ＆Ａ アーカイブ(5)]] (元記事が 2005-03-13 より 2...
-[[初級Ｑ＆Ａ アーカイブ(4)]] (元記事が 2005-11-09 より 2...
-[[初級Ｑ＆Ａ アーカイブ(3)]] (元記事が 2005-05-02 より 2...
-[[初級Ｑ＆Ａ アーカイブ(2)]] (元記事が 2004-12-13 より 2...
-[[初級Ｑ＆Ａ アーカイブ(1)]] (元記事が 2004-08-03 より 2...
----
#contents
----
**凡例中に数式と数値変数を繰り返し書きたいのですが [#db8d...
>[[TAG]] (2007-06-27 (水) 21:50:40)~
~
使用環境:R version 2.5.0 (2007-04-23), i386-apple-darwin8...
~
#ref(legendmiss.png)~
凡例として~
「1 ミュー g／2 ミュー g／3 ミュー g」~
のように数字と数式（ギリシャ文字）を混ぜて書きたいと思い~
cols <- rainbow(3)
plot(sin, xlim = c(0, pi), ylim = c(-1, 1), ann=F, col=c...
a <- NULL
b <- NULL
a <- expression(paste(1, mu, "g"))
b <- c(b, a)
a <- expression(paste(2, mu, "g"))
b <- c(b, a)
a <- expression(paste(3, mu , "g"))
b <- c(b, a)
legend(0.3, -0.2, b, col = cols, lwd=1)
上のようにコードを書いてみたのですが（プロットの部分はダ...
面倒なので凡例に関する部分を~
a <- NULL
b <- NULL
for(i in 1:3)
{
a <- expression(paste(i, mu, "g"))
b <- c(b, a)
}
legend(0.3, -0.2, b, col = cols, lwd=1)
とすると、「iミューg／iミューg／iミューg」のように変数iが...
paste関数の中なので当然なのかもしれませんが何か良い方法は...
~
RjpWiki内のグラフィックス参考実例集：数式のプロットの例2に~
「数式と数値変数の結合」とあるので参考にしてexpressionの...
a <- NULL
b <- NULL
for(i in 1:3)
{
a <- substitute(paste(x, mu, "g"), list(x=i))
b <- c(b, a)
}
legend(0.3, -0.2, b, col = cols, lwd=1)
このようにすると「paste(1, mu, "g")／paste(2, mu, "g")／p...
数値変数は反映されたのですがその他が駄目になるという結果...
解決策はあるのでしょうか？~
どなたか対策をご存知でしたらご教授願えますでしょうか。~
//
- 今の場合なら，以下のようにしたらいかが？ -- &new{2007-...
cols <- rainbow(3)
plot(sin, xlim = c(0, pi), ylim = c(-1, 1), ann=F, col=c...
b <- paste(1:3, "μg")
legend(0.3, -0.2, b, col = cols, lwd=1)
- なるほど！全角文字を使ってしまえば良いのですね。盲点で...
- legend(0.3, -0.2, sapply(1:3,function(i)(eval(substitut...
- legend(0.3, -0.2, parse(text=sapply(1:3,function(i)past...
- さっそく試してみました。見事に上手くいきました！ありが...
- 一応付け加えておくと、こんなにややこしいのは、legendにv...
- 確かにtext, mtextだと何とかできたのですが、今回は自分の...
**Condor & R [#c0745713]
>[[HPC]] (2007-06-27 (水) 00:31:04)~
~
Condor というPCクラスター上での Rの動作例ってありますか...
//
**NAを含む行列を関数distで計算 [#c99db955]
>[[su]] (2007-06-25 (月) 13:34:43)~
~
使用環境：R version 2.5.0 (2007-04-23) i386-apple-darwin8...
NAを含む行列xyで行xと行yとの間の市街距離を計算しました。~
x = c(1,1,NA)
y = c(1,2,NA)
xy = rbind(x,y)
dist(xy[,-3], method="manhattan") # 3列目を除いた距離は...
dist(xy[,3], method="manhattan") # 3列目だけの距離は、NA
dist(xy, method="manhattan") # 3列目を含めた距離は、1.5
dist(xy, method="manhattan")の値が、1でもなく、NAでもなく...
//
- If some columns are excluded in calculating a Euclidean...
- 意訳すれば、NAを含む列は計算に使われないけど、その分距...
- ありがとうございました -- [[su]] &new{2007-06-25 (月) 2...
**データフレームからデータが取り出せない [#z9b0bbf4]
>[[K]] (2007-06-22 (金) 10:03:40)~
~
度数分布表の値をもとに処理をしようとしているのですが、デ...
> dbh_table <- table(data$DBH.Cls0)
> dbh_df <- data.frame(dbh_table)
> dbh_df
DBH.Cls0 Freq
1 18 5
2 22 10
3 24 35
4 26 45
5 28 30
(中略)
20 58 5
21 60 5
22 72 5
としてデータフレームdbh_dfを作成したのですが、2列目のデー...
> dbh_df[1,2]
[1]5
として欲しい"5"だけが得られるのですが、1列目のデータにつ...
> dbh_df[1,1]
[1] 18
22 Levels: 18 22 24 26 28 30 32 34 36 38 40 42 44 46 48 ...
となってしまい、必要な"18"の他に余分な1行が表示されてしま...
x <- dbh_df[1,1]; y <- dbh_df[1,2]
としても、xに余分な1行が入ってしまい、処理できずに困って...
どなたか、原因と対策をご教授いただけますでしょうか。
よろしくお願いします。
~
//
- 処理が出来ないとは具体的にはどういうことでしょう? 1列目...
- コメント、ありがとうございます。「処理ができない」と書...
- 自己レスです。as.vector(x)とすることで、期待する値を得...
- 18だけを取り出せているんですよ。ただし，18は数値ではな...
data.frame を作った後に， dbh_df[,1] <- as.numeric(as.vec...
- なるほど、Rはfactorのリストを知らせてくれていたわけです...
- ちなみに ?as.factor によれば、因子化された数値ベクトル...
**平行座標プロットの各ラインの選択 [#i0b82e1d]
>[[EDA]] (2007-06-21 (木) 19:01:15)~
~
以前、マウスによるポイントの選択方法が掲載されていまし...
//
- どのようなグラフで，どのように選択して，どのような情報...
parcoord の中を見て，マウスでクリックされた点から一番近い...
結構なプログラム量になりそう。 -- &new{2007-06-21 (木) 1...
**sub()で文字列を置き換えると，特定の文字で文字化けする [...
>[[奥村泰之]] (2007-06-21 (木) 13:38:26)~
~
御世話になっております。
下記のコードに示したように，sub()で文字列を置き換えようと...
文字化けする現象が出て困っております。
対応策を御存じの方，お教え頂けると幸いです。
> sessionInfo()
R version 2.5.0 Patched (2007-04-27 r41346)
i386-pc-mingw32
locale:
LC_COLLATE=Japanese_Japan.932;LC_CTYPE=Japanese_Japan.93...
LC_NUMERIC=C;LC_TIME=Japanese_Japan.932
> sub("AA", "山本", "AA氏")
[1] "山本氏"
> sub("AA", "五十嵐", "AA氏")
[1] "五女虫?x81"
//
- sub("AA", "五十??嵐", "AA氏"); でいいんですが、内部エ...
- ただのバグです. CP932環境で, sub("A","?u8868","A") がう...
- takahashi様，なかま様，ご回答を頂きありがとうございまし...
- 先程パッチ付きレポートを送りました. -- [[なかま]] &new{...
- 2.5.1以降は大丈夫. -- [[なかま]] &new{2007-06-27 (水) 1...
- なかま様。R-develにレポートまでして頂き，誠に感謝してお...
- version 2.5.1 で対応されていることを確認致しました。ま...
**散布図：同じような値を重ならることなく、横並びで表示し...
>[[ちょろ]] (2007-06-20 (水) 14:29:14)~
~
散布図を加工としております。~
しかし、plot関数ではデータが一直線に並んでしまい、同じよ...
一つ一つのドットが重ならないようにしたいのですが、どのよ...
イメージ的には以下のURLの一番上にある画像の左上部のような...
http://www.ne.jp/asahi/room/kuro/ScatterMakerReadMe.html~
~
ご教授いただければ幸いです。~
//
- library(gplots); ?space -- [[takahashi]] &new{2007-06-2...
- 宣伝しておきます。http://aoki2.si.gunma-u.ac.jp/R/dot_p...
- takahashi様、青木先生、アドバイスありがとうございます。...
- 似たようなものとしてはstripchart(Y~X,vert=T,method="sta...
**3次元のデータをX-Yと色分けでplotする方法 [#ga03839a]
>[[ぼう]] (2007-06-19 (火) 17:59:04)~
~
お世話になります。XYZの3次元のデータがあり、Zを色分けして...
~
mat <- matrix(rnorm(30),nrow=10)
plot(mat[,1],mat[,1],col=
ifelse(mat[,3]>0.6,2,
ifelse(mat[,3]>0.3,3,
ifelse(mat[,3]>0,4,
ifelse(mat[,3]>-0.3,5,
ifelse(mat[,3]>-0.6,6,
7
)
)
)
)
)
)
//
- データ数＝色数でよいのなら -- [[g]] &new{2007-06-19 (火...
~
mat <- matrix(rnorm(30),nrow=10)
plot(mat[order(mat[,3]),][,1],mat[order(mat[,3]),][,2],c...
- ありがとうございます。orderでzの順番にするとは妙技です...
- 値を階級に分けて，その階級に色を適用するのでしょうから...
mat <- matrix(rnorm(30), nrow=10)
m = 5 # 階級の個数
# suf は データの値により 1〜m になる
# 例ではいい加減に分類したけど findInterval の第二引数で...
suf <- findInterval(mat[,3], seq(min(mat[,3]), max(mat[...
plot(mat[,1], mat[,2], col=rainbow(m)[suf]) # rainbow で...
- じゃあ，無駄に大げさにして -- [[g]] &new{2007-06-19 (火...
~
plot(mat[order(mat[,3]),][,1],mat[order(mat[,3]),][,2],
col=rep(c("red","orange","yellow","green","blue","violet...
hist(mat[,3],breaks=c(min(mat[,3]),-.6,-.3,0,.3,.6,max(m...
$counts))
- お二方ありがとうございます。分野によっては頻繁に使いそ...
**plotでx軸とy軸を任意の座標で交わらせることはできますか...
>[[くに]] (2007-06-19 (火) 10:24:13)~
~
plotでx軸、y軸を(0,0)などの任意の座標で交わらせることはで...
#ref(kuni.png)
のように。よろしくお願いします。~
//
- axis 関数のヘルプを読むべし -- &new{2007-06-19 (火) 10...
> x <- runif(20, min=40, max=100)
> y <- runif(20, min=10, max=40)
> plot(x, y, axes=FALSE, xlim=c(0,100), ylim=c(0, 80))
> axis(1, pos=0)
> axis(2, pos=0)
- ほんとだ、helpにちゃんと書いてありますね。質問する前に...
**等高線グラフィックのスケールについて [#cc03ca46]
>[[V9]] (2007-06-13 (水) 14:31:14)~
~
たまたま等高線の話題なので、続けさせていただきます。~
contourで描いた等高線に、二次元の点をプロットしたく思って...
単なるcontourではスケールをぴったり合わせられるのですが、...
スケール（軸のメモリ）を合わせるのに、何かうまい方法はあ...
//
- 大変すみませんでした。具体的には、グラフィックス参考実...
- 説明不足で大変申し訳ございませんでした。実は、例３の中...
- 大変すみません。削除の方法を教えて頂けないでしょうか。 ...
- 実は、以下のようなものを書いて実行したら、エラーが出て...
- 上の方の「編集」をクリックして，編集画面の下の方に「テ...
ところで，何回も聞くけど ? filled.contour は読んだの？ -...
-以下のようなものを書いて実行したら、エラーが出てしまいま...
引数の長さが０とは、どういう意味なのでしょうか？
以下にエラーif (axes) { : 引数の長さが0です
filled.contour(
x,y,z,
levels=c(0,10,30,50,70,90,100),
xlim=c(-1,1),ylim=c(-100,100),key.axes=axis(4,seq(0,100,...
plot.title=title(main="main-title",xlab="x",ylab="y"),
key.title=title("ex"),
axes={axis(1);axis(2);points(exdata,pch=16,col="black")},
)
-exdataは、以下のような型式です。
x y
1 -0.2 20
2 0.8 40
3 0.1 30
・・・
- 描画できました。どうもありがとうございました。大変お手...
- ちなみに、プロットした点は、値によって色を変えることは...
- ちなみに、プロットした点は、値によって色を変えることは...
- plot.axis の { } の中を，目的に従って書けば良いだけ。 ...
- ありがとうございました。点の色分けはできたのですが、判...
- 判例の件、以下のようなものでうまく表示されませんでした...
filled.contour(
x,y,z,
levels=c(0,10,30,50,70,90,100),
xlim=c(-1,1),ylim=c(-100,100),key.axes=axis(4,seq(0,100...
plot.axes={axis(1); axis(2); points(exdata, pch=16, col=...
"deeppink", "black", "deepskyblue", "b...
legend(0, 0, paste("value", c("0-20%", "20-40%",
"40-60%", "60-80%", "80-100%")), bg="g...
col=c("blue", "deepskyblue", "black", ...
"red"), pch=16),
plot.title=title(main="main-title",xla...
key.title=title("ex")
)
- そこで、plot.axexの中で修文したところ、うまく表示できま...
**等高線用の外積計算 [#q55775cf]
> (2007-06-13 (水) 13:43:54)~
~
等高線を描くために外積を計算しようとしていますが、length...
何かうまい設定法がございましたら、ご教示頂けると幸いです。
> x<-seq(-4,4,length=20)
> y<-seq(-10000,4000,length=20)
> z<-outer(x,y,function(xi,yj) 100/(1+exp(-(0.279-1.0344...
+ )
以下にエラーouter(x, y, function(xi, yj) 100/(1 + exp(-(...
dim<- : dims [product 400] は object [20] の長さ...
//
- 基本的な間違いで、恥ずかしい限りです、、申し訳ありませ...
**無題 [#x0c3115a]
>[[あ]] (2007-06-10 (日) 20:00:15)~
~
sourceで時系列オブジェクトを読み込もうとしたら、下記のエ...
以下にエラーattributes(.Data) <- c(attributes(.Data), at...
不正な時系列パラメータが指定されました
このエラーはどのような意味なのでしょうか。~
読み込もうとしているデータは以前に読み込みを成功した事が...
//
- 無題とはいやはや・・・
読み込もうとしているデータも，変更の内容もわからないのに...
使用環境，問題の再現可能なデータセット，ロードしているラ...
ないと，誰も答えようがないと思いますが． -- &new{2007-06...
- 理由？不正な時系列パラメータを指定したからです！~
あ さん とは。如何に匿名とはいえ，もう少し考えたら？ -- ...
**ヒストグラムのｙ軸目盛（スケール？）の変更は可能でしょ...
>[[くに]] (2007-06-07 (木) 12:14:30)~
~
ヒストグラムのｙ軸目盛（スケール？）の変更は可能でしょ...
面積が0.3ヘクタールの調査区で樹木の胸高直径を測定しまし...
そのようなことは可能でしょうか。よろしくお願いします。~
//
- 元のデータを 1/0.3 倍して、そのヒストグラムを描けばいい...
- 早速のコメントありがとうございます。元のデータをxとする...
- 余計な御世話かもしれませんが、0.3ha毎の観測値が100個あ...
- ahaha,上のコメントで目が覚めた。私の提案していたのは，...
set.seed(77789)
x <- rnorm(100, mean=20, sd=5)
par(mar=c(4,4,1,4))
hist(x)
axis(4, at=0:4*10, labels=round(0:4*10*(1/0.3),1))
- 色々試行錯誤の末、少し求めているものに近づいてきました。
2cm刻みの直径階別本数を先に求めておいて、そのベクトルを...
#ref(upload.png)
でも、これではもう一息と言うところ。table関数で得られた...
>area=0.3
>layout(matrix(1:2),2)
>DBHclass <- round((data$x+0.05)/2,0)*2
>#15cm以上17cm未満は18cmという、2cm刻みの直径階にするた...
>brks <- seq(14,74,by=2)
>hist(DBHclass, breaks=brks,xlab="胸高直径階(cm)", ylab=...
>
>DBH_table <- table(DBHclass)
>plot(DBH_table/area,type="h",xlab="胸高直径階(cm)", yl...
>
>dbh<-data.frame(DBH_table/area)
>dbh
DBHclass Freq
1 18 5
2 22 10
3 24 35
4 26 45
5 28 30
6 30 55
7 32 45
8 34 60
9 36 55
10 38 90
11 40 75
12 42 65
13 44 85
14 46 30
15 48 10
16 50 10
17 52 10
18 54 15
19 56 5
20 58 5
21 60 5
22 72 5
このデータフレームから、forループか何かで c(18,18,18,18,1...
- x<-rnorm(1000);hist(x,axes=F);axis(2,0:4*50,0:4*500) --...
- x軸も一旦退避しておいてaxis()で書き込む -- &new{2007-0...
set.seed(77789)
x <- rnorm(100, mean=20, sd=5)
temp <- hist(x, breaks=seq(5,40,2), axes=FALSE, xlab="...
axis(2, 0:5*3, round(0:5*3*10/3,0))
axis(1, temp$breaks)
#ref(hist.nag.png)
- parを活用すれば, -- [[なかま]] &new{2007-06-08 (金) 10:...
set.seed(77789)
hist(x,yaxt="n")
axis(2,seq(par()$yaxp[1],par()$yaxp[2],5))
固定値は限りなく少なくなる筈.
- なるほど、
axis(2, 0:5*3, round(0:5*3*10/3,0))
としたところ、欲しいグラフが得られました。皆さん、ありが...
**和や部分リストを得る方法は [#k9c7b21b]
>[[kd]] (2007-06-06 (水) 10:16:48)~
~
リストを与えて，和や部分リストを得たいのですが方法はある...
> dwt<-dwt(x,n.levels=8, boundary="reflection")
> str(dwt@W)
List of 8
$ W1: num [1:300, 1] -0.188 0.166 -0.150 0.170 0.115 ...
$ W2: num [1:300, 1] 0.6416 0.5408 0.0322 -0.3626 -0....
:
:
次に，この dwt@W に対して和や部分リストを得たい思っていま...
dwt@W[1:4]
に対して
dwt@W[[1]]+dwt@W[[2]]+dwt@W[[3]]+dwt@W[[4]]
を得たり，
dwt@W[1:4]
に対して
list(dwt@W[[1]][1:150],dwt@W[[2]][1:150],dwt@W[[3]][1:15...
を得たいと思っています．~
//
- dwt ってのは，どのパッケージに入っているのか？x, n には...
そんなこと教えなくても答えられる人だけ答えてくれというの...
上の方に長々と書いてある注意事項は読んでないんだろうね。...
- 一個目はx<-list(1:3,4:6,9:12,13:15);do.call("+",x[1:2])...
- ただし一個目は"+"が二項関数なので３項以上あると問題だっ...
- そもそも，
dwt@W[[1]]+dwt@W[[2]]+dwt@W[[3]]+dwt@W[[4]]
ってやって正しい答えが出るならそれでいいじゃんか -- &new...
- ご教示いただきありがとうございます．forループで回す関数...
**指定データを除いたプロット [#o1a81cf5]
>[[kd]] (2007-06-03 (日) 15:05:17)~
~
下記のように，指定データを除いてプロットすることは可能で...
~
最初に下記を行いました．
> newx<-modwt(x,n.level=8,boundary="reflection")
> str(newx2)
Formal class 'modwt' [package "wavelets"] with 11 slots
..@ W :List of 8
.. ..$ W1: num [1:3061, 1] -1.82e-06 3.58e-05 -1.12e-...
:
:
..@ V :List of 8
.. ..$ V1: num [1:3061, 1] -2.14e-06 3.91e-05 1.48e-...
:
:
> plot(newx)
このplot(newx) を実行すると，x,newx@W[1],newx@W[2],...,ne...
次にnewx@V[8] だけを除いたプロットを行いたいので，下記の ...
> newx1<-????(newx)
> plot(newx1)
またnewx@W[8] だけをプロットすることも可能でしょうか？
> newx2<-????(newx)
> plot(newx2)
//
- まず str(newx2) は str(newx) の、x,newx@W[1],... は n...
- どうもありがとうございます．計算結果 y1, y2, ... を単に...
- par(mar=c(1, 1, 1, 1)) -- &new{2007-06-05 (火) 15:49:...
**locfit パッケージのロードに失敗 [#b56fb3c5]
>[[kd]] (2007-06-02 (土) 21:23:34)~
~
locfit を v. 2.5 for MacOSX 上でロードしようとすると，下...
要求されたパッケージ akima をロード中です
エラー： パッケージ 'akima' をロードできませんでした
追加情報: Warning message:
'akima' という名前のパッケージはありません in:
library(pkg, character.only = TRUE, logical = TRUE, lib...
locfit は下記からダウンロードしたものです．~
ftp://ftp.u-aizu.ac.jp/pub/lang/R/CRAN/bin/macosx/univers...
ちなみにwavelets等の他のパッケージは正常にロードできて動...
ftp://ftp.u-aizu.ac.jp/pub/lang/R/CRAN/bin/macosx/univers...
//
- 見たところ locfit パッケージは akima パッケージに依...
- この下の方にある「Macでパッケージvcdがロードできません ...
- どうもありがとうございます．「パッケージの依存関係」な...
- ずぼらしないで、まずは akima をインストールしたらど...
- 重ね重ねご教示をいただき本当にありがとうございます．aki...
**カレントディレクトリのファイル名一覧を取得 [#z6e16a93]
>[[mori]] (2007-06-01 (金) 10:41:10)~
~
カレントディレクトリに存在するファイル名の一覧をベクトル...
//
- dir <- dir(path=getwd(), pattern="file") のこと？ -- &...
- ありがとうございます。dir()でできたんですね・・・ -- [[...
- list.files() http://takenaka-akio.cool.ne.jp/doc/r_auto...
- dir と list.files は同じものです。 -- &new{2007-06-04 ...
> identical(dir, list.files)
[1] TRUE
**TIFFファイルをreadLinesで読み込みたい [#yb1adda9]
>[[akira]] (2007-05-31 (木) 20:22:00)~
~
測定器が出力するTiffファイルの1行目に測定条件が記入されて...
> x <- readLines("test.tiff", n=1)
Warning message:
'test.tif' に関する readLines で不完全な最終行が見つか...
画像ファイルをテキストで読み込む方法を教えてください。~
~
OS:WinXPSP2~
R version 2.5.0 Patched (2007-05-04 r41429)
i386-pc-mingw32
です。~
//
- そうですか。MacOS X 10.4 では読めました（何の追加情報も...
こんなファイルを
foo [9] > hexdump -c yellow.tiff
0000000 M M ?0 * ?0 ?0 002 220 301 ?0 301 ?0 ...
0000010 030 001 9 023 ?0 ] ] 037 ?0 253 z ) ...
読んでみました。
R での結果は以下のようになりました。
> x<-readLines("yellow.tiff", n=1)
> x
[1] "MM"
これがあなたの望む結果なんでしょうか。あなたが読んだ tif...
> sessionInfo()
R version 2.5.0 Patched (2007-05-18 r41630)
powerpc-apple-darwin8.9.1
- exifならそんなパッケージがあったような気がしますし, バ...
- ありがとうございます。LINUXでもうまくいきませんでした。~
ubuntu Linux~
R version 2.4.1 (2006-12-18) i486-pc-linux-gnu
本当なら、実際に使用している画像をお見せしたいのですが、...
[[デジカメ画像:http://cweb.canon.jp/camera/ixyd/810is/dat...
rgdalやrtiffのパッケージを用いるとtiffの数値情報をとれま...
> x <- readLines(dir()[1])
Warning message:
'810IS_l0.jpg' に関する readLines で不完全な最終行が見つ...
> x
[1]エラー：マルチバイト文字列が不正です
なお、測定器のTiffファイルにも同じように測定条件の文字列...
なかま様、これはバイナリなのでしょうか？コーディングの問...
- んと, JPEGやTIFFには[[EXIF(ぐぐれば沢山あるでしょう):ht...
- なかまさま、ありがとうございます。バッチ処理をしたいの...
> x <- readBin("test.tiff", what="character", endian="sw...
> x
[1] "II*"
で、scan()を試してみたら、"エラー：マルチバイト文字列が不...
> x <- scan("test.tiff", what="character")
Read 26373 items
> x
[1] "II*"
[2] ""
[3]エラー：マルチバイト文字列が不正です
これならいけそうです。~
これって、「[[ファイルTips大全の特定の文字コードのファイ...
- readBin("test.tiff", what="", n=1) でどうなりますか？わ...
- ありがとうございます。今はOSをubuntuで計算させているの...
- バイナリをバイナリで読むには, -- [[なかま]] &new{2007-0...
con<-file("binary.file",open="rb")
raw<-readBin(con,"raw",8)
raw # Hexで出る
# もし, rawの中がEUCで自分の環境がUTF-8なら"
iconv(rawToChar(raw),"EUC-JP","utf-8")
# めでたし,めでたし
close(con)
などとしないとダメでしょう.
- なかまさまありがとうございます。教えていただいた通り実...
test.tiffの1行目の"II*"から目的とする文字列の間にある制御...
- IIはTiffのフォーマットがLittleEndianと申しておりますか...
- ありがとうございます。生まれがwinなものでfileとかpipeと...
- [[libtiffライブラリ:http://www.libtiff.org/]]の[[tiffin...
libtiffはrtiffパッケージが使用しているのですが、tiffinfo...
- とりあえず、fileとscanだけでもほぼ望む結果が得られるこ...
con <- file(description=i, open="rb", blocking = TRUE)
x <- scan(con, what="logical", n=20, sep=NULL, comment.c...
どの辺が「ほぼ」かというと、^@のコードの直後の文字が取得...
**オンボードRAIDは使ってますか? [#te4364d3]
>[[通りすがりの佐々木]] (2007-05-31 (木) 16:22:29)~
~
R用PC(Athlon64,3200,2G,WXPSP2)をMBオンボードのRAIDを使っ...
//
- 昔家のパソコンで、シミュレーションを夏休み中ぶっつづけ...
- みなさんのRAID環境をしりたかったりする。Rを使うためには...
- 私は，RAID なんていうのは，今の今まで知らなかった。私は...
- RAIDなんてものを信用してはいけない。外付けHDに、こまめ...
- 性能的な物だと, 廉価版RAIDはOS+CPUに処理を丸無げするの...
- やはりオンボードRAIDは負荷がかかるのですね。実はRAIDに...
- その構成ならeSATAが使えればいいですね。NASだとEther自体...
- GentooのHPにSOFTraidについての翻訳がありました。Linux R...
- 実は、先日のクラッシュでHUBもお亡くなりになったので、ギ...
- GigaEtherの欠点はjumbo packetにNIC&HUBが対応していない...
**２変数関数の等高線プロットについて [#n6dc8ef2]
>[[GG]] (2007-05-31 (木) 11:55:06)~
~
２変数（ｘ、ｙ）関数の値（ｚ）を求めることなく、~
式の形から等高線を直接描ける方法はあるのでしょうか？~
例えば、z=1/(1+exp(-(10x-3y)))のようなグラフを直接３次元...
適当に等高線で色分けした後に、~
さらに重ねがきで幾つか点もプロットしたいと思っています。~
//
- 例示のようなきれいな関数なら，z=constant として，y=f(x,...
- 実際には例示よりももう少し複雑な式を考えておりますが、...
- y=もなくて，z も計算しなくて，コンピュータの中にこびと...
**日本語入力してもグラフが表示されない [#u4f683a9]
>[[のり]] (2007-05-27 (日) 20:39:34)~
~
Rインストール後、コマンド（R console）で日本語入力しても...
（例）
x=0:50
y=dbinom（x,50,0.25)
plot（x,y,type=’ｈ’,xlab=’x’,ylab=’ｙ’,main=’二項分布’）
＋
（注、type=ｈはヒストグラム表示、xlab、ylabはｘ、ｙ軸の名...
（使用環境 windowsXP）~
この解決方法をご存じの方がおりましたら、ご教示いただけな...
（尚、この質問は最初、初級者コースが閉鎖されていると認識...
//
- 張り付けられた文字と「日本語を入力しない他のコマンドで...
- type='h' の h までもが，全角じゃ〜〜ないですか。なんで...
x=0:50
y=dbinom(x,50,0.25)
plot(x,y,type='h',xlab='x',ylab='ｙ',main='二項分布')
#ref(baz.png)
- 有り難うございました。日本語に注意を取られすぎ、いろい...
- 本当の正解はこれなのでは？ -- [[okinawa]] &new{2007-05-...
x <- c(0:50)
y <- dbinom(x,50,0.25)
plot(x,y,type='h',xlab='x',ylab='ｙ',main='二項分布')
- あと、sessionInfo( )とコンソールから実行してみてくださ...
- c(0:50) は冗長でしょう。0:50 でよい。= を使って悪いわけ...
- c(0:50)は冗長だとは思いますが、R言語の作法に沿ってコー...
x <- 0:50
y <- dbinom(x,50,0.25)
plot(x,y,type="h",xlab="x",ylab="y",main="二項分布")
- 前にも投稿記事がありましたが、付値演算子として = を使う...
> order(x <- 1:10)
[1] 1 2 3 4 5 6 7 8 9 10
> x
[1] 1 2 3 4 5 6 7 8 9 10
> order(X = 1:10)
[1] 1 2 3 4 5 6 7 8 9 10
> X
エラー： オブジェクト "X" は存在しません
- order(x <- 1:10) がコードの簡略化に役立つとは思いません...
? c で，This is a generic function which combines its arg...
> identical(c(1:10), 1:10)
[1] TRUE
**Macでパッケージvcdがロードできません [#w5c8d728]
>[[ピーター]] (2007-05-25 (金) 18:46:30)~
~
Mac OS X 10.4.9, R 2.5.0 GUI 1.19（4308）を使っています~
「パッケージとデータ」のCRANからvcdをインストールし、「R...
要求されたパッケージ colorspace をロード中です
エラー： パッケージ 'colorspace' をロードできませんでした
追加情報: Warning message:
'colorspace' という名前のパッケージはありません
in: library(pkg, character.only = TRUE, logical = TR...
とあります。また、library(vcd)でも同様でした。~
一体なぜののでしょうか。教えてください。お願いします。~
//
- colorspaceというパッケージをインストールする必要がある...
- パッケージインストーラでインストールするときに，いつも...
- なるほど！colorspaceをインストールしたら、vcdもロードで...
**三次元プロットについて [#g5c55634]
>[[dmc]] (2007-05-23 (水) 13:42:22)~
~
scattaerplot3d関数で三次元プロットした点に、persp関数で二...
点と曲面をうまくひとつのグラフ内に描く方法がございました...
//
- perspだけでやっちゃう。 -- [[takahashi]] &new{2007-05-2...
x <- seq(-10, 10, length= 30)
y <- x
f <- function(x,y) { r <- sqrt(x^2+y^2); 10 * sin(r)/r }
z <- outer(x, y, f)
z[is.na(z)] <- 1
pmat<- persp(x, y, z, theta = 30, phi = 30, expand = 0.5...
ltheta = 120, shade = 0.75, ticktype = "det...
xlab = "X", ylab = "Y", zlab = "Sinc( r )")
p<-list(x=runif(100,-10,10),
y=runif(100,-10,10),
z=runif(100,0,5))
points(trans3d(p$x,p$y,p$z, pmat), col = 2, pch =16)
- ご回答頂き、どうもありがとうございました。無事に描画で...
- 追伸です。曲面を塗らず、透過させるにはどのようなコマン...
- perspの引数colを0とかNULLとか"transparent"とかってこと...
- 引数を０やtransparentにしてみたら、うまく描画できました...
**データがマイナスとなる個数 [#g2824fdf]
>[[port]] (2007-05-23 (水) 12:02:52)~
~
Ｒの超初心者です。~
10000個くらいのデータがあって，そのデータの中でマイナスと...
色々調べたのですが，よくわかりませんでした。~
どのたかご教授ください。~
//
- dat<-rnorm(10^4); sum(dat<0) -- [[takahashi]] &new{2007...
- どうもありがとうございます。助かりました。 -- [[port]] ...
- dat<-rnorm(10^4); length(which(dat<0))、でしょう。 -- [...
**Excelからのコピペが上手くいきません [#cc7998f8]
>[[みかん]] (2007-05-14 (月) 04:58:23)~
~
「Tips紹介」→「R にファイルを読み込む tips 集(暫定版)」→~
「Excel からコピーしたセル範囲を R に読み込む (何でも掲示...
に記載されております方法が上手くいきません。~
記載事項の通り、~
> from.excel <- function(nc)
+ {
+ matrix(scan(""), byrow=TRUE, nc=nc)
+ }
> from.excel(3)
としたのち、
エクセルで
1 F 4
2 F 3
3 M 2
という3×3のセル範囲をコピーし、貼り付けると、~
> from.excel(3)
1: 1F4
1: 2M3
以下にエラーscan("") : scan 関数は '$s' を期待したのに、...
> 3F2
エラー："3F2" に構文エラーがありました
>
というメッセージが出てきます。~
一方同じExcelのコピー範囲を一度テキストエディタに貼り付け...
> x <- read.table(stdin())
の後に貼り付けると、上手くいきます。~
しかしこの方法でも、Excelでコピーしたそのままを貼り付ける...
> x <- read.table(stdin())
0: 1F4
1: 2M3
2: 3F2
3:
>
となってしまいます。~
どうもExcelでコピーした列セル間の区切り子が、R上で認識さ...
どうしてこうなってしまうのか、ご教授をお願いいたします。~
当方WindowsXP上で、ver.2.2.1を使用しております。Excelのバ...
よろしくお願いいたします~
//
- エクセルで範囲選択してコピー後、read.table("clipboard")...
- Windowsね。x をプリントしてみたらどうなっているの？mac ...
> x <- read.table(stdin(),header=FALSE)
0: 1 F 4
1: 2 F 3
2: 3 M 2
3:
> x
V1 V2 V3
1 1 F 4
2 2 F 3
3 3 M 2
- 青木先生の方法ではmatrix関数を使用しているため~
すべての値がすべて数値あるいは文字列になっていないとうま...
~
そこで直接の回答ではありませんが、僕自身は
from.excel2 <- function(h = 0)
{
dat <- (read.delim("clipboard", header=h))
}
のような関数を別に定義してExcelのデータを読み込んでいます...
タイトル行がある場合は header=T とします。~
この場合にはデータフレームとして読み込まれるのでうまくい...
~
また、僕自身はstdin()を知らないです。-- [[あら]] &new{200...
- dat に代入するのではなく（dat って，関数外のオブジェク...
- 教えるつもりが逆に教えていただく型になってしまいました...
**直交実験を計画してくれる関数 [#v29e7c5d]
>[[TH]] (2007-05-13 (日) 09:21:44)~
~
直交実験を計画してくれる関数はどこかのRパッケージに含まれ...
//
- CRAN にあるパッケージ crossdes あたりは使えませんか。 -...
- [[Rseek:http://www.rseek.org/?cx=010923144343702598753%...
- ありがとうございます．でもcrossdesはラテン方格と釣り合...
- oa.design とそれに関する関数を用意して，その段階ではど...
- 移植しようと昨日からいろいろやっているうちに，全体の移...
**確率微分方程式 [#fa2c1ef1]
>[[SDE]] (2007-05-10 (木) 19:02:49)~
~
確率微分方程式で検索してもヒットしなかったのですが、~
解を求めるなど、確率微分方程式を取り扱ったパッケージは無...
Rの守備範囲だと思うのですが…~
//
- どうやって検索したのか？？site: r-project.org "Stochast...
- ありがとうございます。検索が雑でした。CRAN.packagesで一...
**Fizz-Buzz問題 [#x4e95cdc]
>[[X Jr.]] (2007-05-09 (水) 22:20:18)~
~
どのカテゴリに書こうか迷ったのですが、とりあえず初級Q&Aに...
「[[どうしてプログラマに・・・プログラムが書けないのか?:h...
//
- [どうしてプログラマに・・・プログラムが書けないのか?]そ...
> x <- 1:100
> print(ifelse(x%%15==0,"FizzBuzz", ifelse(x%%3==0,"Fizz...
[1] "1" "2" "Fizz" "4"
[5] "Buzz" "Fizz" "7" "8"
[9] "Fizz" "Buzz" "11" "Fizz"
[13] "13" "14" "FizzBuzz" "16"
[17] "17" "Fizz" "19" "Buzz"
:
[93] "Fizz" "94" "Buzz" "Fizz"
[97] "97" "98" "Fizz" "Buzz"
- やり方は幾つもあるけど，どれが一番「プログラマが書いた...
x <- 1:100
x[1:33*3] <- "Fizz"
x[1:20*5] <- "Buzz"
x[1:6*15] <- "FizzBuzz"
print(x)
- せめて，これを拡張して，1000までの素数をエラトステネス...
- cat の方がいいかな。それと ==0 の代わりに ! を使う -- ...
> x<-1:100;cat(ifelse(!x%%15,"FizzBuzz",ifelse(!x%%3,"Fi...
1 2 Fizz 4 Buzz Fizz 7 8 Fizz Buzz 11 Fizz 13 14 FizzBuz...
Fizz 19 Buzz Fizz 22 23 Fizz Buzz 26 Fizz 28 29 FizzBuzz...
Fizz 34 Buzz Fizz 37 38 Fizz Buzz 41 Fizz 43 44 FizzBuzz...
Fizz 49 Buzz Fizz 52 53 Fizz Buzz 56 Fizz 58 59 FizzBuzz...
Fizz 64 Buzz Fizz 67 68 Fizz Buzz 71 Fizz 73 74 FizzBuzz...
Fizz 79 Buzz Fizz 82 83 Fizz Buzz 86 Fizz 88 89 FizzBuzz...
Fizz 94 Buzz Fizz 97 98 Fizz Buzz
2文字だけ短いもの（！？）
x<-1:100;cat(ifelse(x%%15,ifelse(x%%3,ifelse(x%%5,x,"Buz...
- データを考える言語の立場から, 条件分岐はどうでしょうね....
M<-expression(1:100)
x<-eval(M)
x[!eval(M)%%3]<-"Fizz"
x[!eval(M)%%5]<-"Buzz"
x[!eval(M)%%3&!eval(M)%%5]<-"FizzBuzz"
- Rでのベクトルと行列の扱いなどの特徴を活かして(^_^;) -- ...
x<-1:105
dim(x)<-c(3,105/3)
x[3,]<-"Fizz"
dim(x)<-c(5,105/5)
x[5,]<-"Buzz"
dim(x)<-c(15,105/15)
x[15,]<-"FizzBuzz"
cat(x[1:100])
これも，ほんの少しだけ短くなる
a<-function(i,s){dim(x)<<-c(i,105/i);x[i,]<<-s}
x<-1:105
a(3,"Fizz")
a(5,"Buzz")
a(15,"FizzBuzz")
cat(x[1:100])
- いろいろと面白い例をありがとうございました。Rはループを...
- 人に聞く前に，あなたのプログラム例をまな板に載せてから...
- 初心者コーナーとあれば一番Rらしくない（笑）のもあげてお...
> for (i in 1:100) {
+ if(i%%15==0) cat("FizzBuz,")
+ else if (i%%3==0) cat("Fizz,")
+ else if(i%%5==0) cat("Buz,")
+ else cat(i,",")}
- そうですね。順序が逆になりましたが、初心者が書くとこん...
- いや，実際，挙げられた一番単純なものが一番確実に動く，...
- 私が採用係なら直前のコードを書いた人物（=> 実は私本人）...
[1] 2.744 0.000 2.757 0.000 0.000
[1] 2.737 0.000 2.748 0.000 0.000
[1] 2.712 0.000 2.723 0.000 0.000
[1] 2.676 0.004 2.698 0.000 0.000
[1] 2.628 0.000 2.633 0.000 0.000
[1] 2.632 0.000 2.649 0.000 0.000
- 私が採用担当なら, プログラマーになりたいと言う人は, 思...
- > 15,3,5 の順で条件判断すべきではなく、3,5,15の順で~
注意しないと，15 の倍数を3の倍数だと判断して（間違いじゃ...
- 実は，一番起こりやすいのは 3,5いずれの倍数でもないとい...
3,5,15の順で正しく判定するためには，かえって手間が掛かる...
ちょっと計算時間が短くなっても，バグ持ちプログラムでは意...
> a <- function(){
+ for (i in 1:1e7) {
+ if(i%%15==0) 15
+ else if (i%%3==0) 3
+ else if (i%%5==0) 5
+ else 0 }
+ }
> system.time(a())
ユーザ システム 経過
34.522 0.116 34.522
> b <- function() {
+ for (i in 1:1e7) {
+ if (i%%3==0 && i%%5!=0) 3
+ else if(i%%3!=0 && i%%5==0) 5
+ else if(i%%15==0) 15
+ else 0 }
+ }
> system.time(b())
ユーザ システム 経過
46.421 0.133 46.393
> bb <- function() { # こうすれば，無駄な比較をしないで...
+ for (i in 1:1e7) {
+ if(i%%3==0) {
+ if (i%%5==0) 15 else 3
+ }
+ else {
+ if (i%%5==0) 5 else 0}
+ }
+ }
> system.time(bb())
ユーザ システム 経過
29.060 0.138 29.072
- one-liner -- &new{2007-05-12 (土) 19:23:10};
replace(replace(replace(1:100,1:33*3,"Fizz"),1:20*5,"Buz...
**.RData の怪 [#b974396e]
>[[昔の名前で出ています]] (2007-05-09 (水) 16:13:25)~
~
R 2.5.0 をインストールして、使おうとしている途中です。~~
ディレクトリの設定をして、q() で「作業スペースを保存」し...
それでは、と、「作業スペースの保存」をして、できた .RData...
> sessionInfo()
R version 2.5.0 (2007-04-23)
i386-pc-mingw32
locale:
LC_COLLATE=Japanese_Japan.932;LC_CTYPE=Japanese_Japan.93...
;LC_NUMERIC=C;LC_TIME=Japanese_Japan.932
attached base packages:
[1] "stats" "graphics" "grDevices" "utils" "dat...
[7] "base"
> # この後、ワーキングディレクトリを設定
> getwd() # 設定の確認
[1] "C:/Documents and Settings/Owner/デスクトップ/r"
> # 作業スペースの保存
> save.image("C:??Documents and Settings??Owner??デスク...
> #ここで、いったん終了
> q()
デスクトップ??r??.RData をダブルクリックすると、なんと、R...
> sessionInfo()
R version 2.4.0 (2006-10-03)
i386-pc-mingw32
locale:
LC_COLLATE=Japanese_Japan.932;LC_CTYPE=Japanese_Japan.93...
;LC_NUMERIC=C;LC_TIME=Japanese_Japan.932
attached base packages:
[1] "methods" "stats" "graphics" "grDevices" "uti...
[7] "base"
//
- んと. explorerを起動して「ツール：フォルダオプション」...
- お返事ありがとうございました。インストール時の「Ｒを.RD...
- q()で終了したとき、イメージはそのときのワーキングディレ...
- q()で終了したとき イメージは ワーキングディレクトリの ....
単に作業スペースの変更の保存だけでなく，関数の定義なども...
私の Windows が，おかしいのかな。 -- &new{2007-05-10 (木...
- .RDataにはワーキングディレクトリは保存されません。.RDat...
- コメントありがとうございました。以下のようにまとめて宜...
(1) 現在の作業ディレクトリに .RData がない場合には，q() ...
(2) .RData がないときには，明示的に「作業スペースの保存」...
(3) 作業ディレクトリは .RData には保存されない。~
(4) .RData をダブルクリックして R を立ち上げると，立ち上...
(5) 既に .RData がある場合には，q() で「作業スペースを保...
- (1)は違います。(2)は(1)が違うので、別に明示的に作らなく...
- ああ，「一つの例外」に陥っていたのですね。 -- &new{200...
**hist で「'breaks' の数が無効です」とは？ [#waf79fcd]
>[[norion]] (2007-05-06 (日) 21:50:20)~
~
hist(5,x)というカキコミをして、（ｘにはあらかじめ数値列が...
'breaks' の数が無効です というエラーメッセージがでます。
5のところを500とかに変えても同じです。データの数は全部で2...
よろしくお願いします。~
//
- ?hist してみましょう。正しい引数を正しい順で与えていま...
- ?histは見ましたが、breaksというものの意味は（英語が得意...
-- [[norion]] &new{2007-05-06 (日) 22:17:27};
- hist(5,x) というようにしても，仰るようなエラーにはなり...
- 1行目に文字があったことに気づかなかっただけで、これを取...
- breaks とはデータ範囲をどのように分割（簡単にいえば棒と...
**複数ウインドウを開いて、描き出すウインドウを選択したい ...
>[[ななし]] (2007-05-05 (土) 13:35:36)~
~
x11() #ウインドウ一つめ
plot(1:10)
x11() #ウインドウ二つめ
plot(2:11)
このとき一つ目のウインドウにグラフを追記したいのですが可...
ひとつのウインドウを分割する場合、split.screenで分割して...
御存じの方教えてください。~
R version 2.4.0~
os linux-gnu vine/linux4.0~
//
- こんな感じの回答がお望みでしょうか。 -- [[aa]] &new{200...
x11() #ウインドウ一つめ
plot(1:10, main="1つめ")
x11() #ウインドウ二つめ
plot(2:11, main="2つめ")
dev.set(dev.prev()) ##1つめに移動
lines(1:10) #1つめに追記
dev.set(dev.next()) ##2つめに移動
lines(2:11) #2つめに追記
graphics.off()
- まさしくそれです。ありがとうございました。 -- [[ななし]...
**garchFit で「引数 "formula" がありません」というエラー ...
>[[たく]] (2007-05-01 (火) 12:46:50)~
~
gaarchFitで解析したいものがあり、試しにUKgasでやってみた...
どうすれば使えるようになるのでしょうか。
バージョンは2.4.1、OSはXPです。
> data(UKgas)
> library(fSeries)
> UKg.d<-diff(UKgas)
> UKg.m <-garchFit(formula.mean =~arma(2,0),formula.var=...
以下にエラー.modelSeries(fake = FALSE, lhs = TRUE) :
引数 "formula" がありませんし、省略時既定値もあり...
//
- その使い方はどこから探し出してきたのか？UKg.m <-garchFi...
- http://www1.doshisha.ac.jp/~mjin/R/0606_35.pdf
このサイトを参考にしました。ちなみにおっしゃるようなやり...
- どこかのpdfは見るけど，オンラインヘルプは見ていないので...
「おっしゃるようなやり方でやっても同じエラー」おかしいな...
> UKg.d<-diff(UKgas)
> UKg.m <-garchFit( formula = ~arma(2,0)+aparch(1,1),dat...
[1] "arma(2, 0)" "aparch(1, 1)"
Series Initialization:
ARMA model: arma
Formula mean: ~ arma(2, 0)
GARCH model: aparch
Formula var: ~ aparch(1, 1)
ARMA Order: 2 0
Max ARMA Order: 2
途中省略
Hessian Matrix:
mu ar1 ar2 omega
mu 0.4038139 -14.661207 -4.759985 0.005485700
ar1 -14.6612069 3721.675424 502.178082 -0.430102036
ar2 -4.7599849 502.178082 4486.891725 -1.201257840
途中省略
gamma1 -9.5579802 57.61092815 -4.6013946 -8.7055680
beta1 85.2086937 -4.60139464 240.7460409 -101.1148148
delta -22.7605912 -8.70556795 -101.1148148 58.0273266
--- END OF TRACE ---
Warning message:
計算結果が NaN になりました in: sqrt(diag(fit$cvar))
> summary(UKg.m)
Title:
GARCH Modelling
Call:
garchFit(formula = ~arma(2, 0) + aparch(1, 1), data = U...
途中省略
Information Criterion Statistics:
AIC BIC SIC HQIC
-9.893399 -9.693561 -9.903576 -9.812388
Description:
Tue May 1 15:54:32 2007 by user:
- 上のをコピペしてみたらできました・・・。さっきはどこか...
**INDSCAL(個人差MDS，３元データのMDS)の関数 [#v19fdb12]
>[[TH]] (2007-04-29 (日) 08:13:46)~
~
私が探した範囲ではRでは，多次元尺度法を行う関数として，以...
計量的多次元尺度法のための関数
パッケージ'stats'の関数 cmdscale
非計量多次元尺度法のための関数
パッケージ'MASS'の関数 isoMSD ，関数 sammon
パッケージ'vegan'の関数 metaMSD
しかし個人差を扱う，3元データ用のMDSモデルのINDSCALを実行...
//
- SensoMineR に indscal : Construct the Indscal model fo...
- 私の知識不足で，Napping data typeというのがいまひとつ分...
- 上切れました（失礼）私の知識不足で，Napping data typeと...
- 投稿は編集できるんですけどねえ。下のR/qtl も，自分でち...
- [[Rと多次元尺度法:http://www1.doshisha.ac.jp/~mjin/R/27...
- 便乗ですみません。MDSのALSCALって、Rでサポートされてい...
- こういう時は RsiteSearch でキーワード検索します。結果は...
- ALSCAL MDS R でググって，出てくるものを調べれば，http:/...
**R/qtl が load できません [#u95b6b2a]
>[[saoki]] (2007-04-27 (金) 23:29:18)~
~
R/qtlがloadできません。初めてRを使う初心者です。どなたか...
Mac OSX バージョン10.4.7のintel MacにR-2.50をインストール...
> library(qtl)
Error in dyn.load(x, as.logical(local), as.logical(now))...
共有ライブラリ
'/Library/Frameworks/R.framework/Versions/2.5/Resources/...
/qtl/libs/i386/qtl.so' を読み込めません
dlopen(/Library/Frameworks/R.framework/Versions/2.5/Re...
Library not loaded: /Library/Frameworks/R.framework/Vers...
Referenced from: /Library/Frameworks/R.framework/Versi...
Reason: image not found
以下にエラーlibrary(qtl) : .First.lib は 'qtl' に対し...
よろしく御願いいたします。~
//
- intel Mac が今ここにないので確かめられないのですが，iBo...
- ありがととうございます。iBook G4があったので（OSXバージ...
- 以下で解決しました。現在のR/qtl (qtl_1.05-2.tar)はR-2.4...
- おかしいですね。少なくとも iBookG4 の R2.5.0 で qtl_1.0...
> sessionInfo()
R version 2.5.0 Patched (2007-04-26 r41343)
powerpc-apple-darwin8.9.1
locale:
ja_JP.UTF-8/ja_JP.UTF-8/ja_JP.UTF-8/C/ja_JP.UTF-8/ja_JP....
attached base packages:
[1] "stats" "graphics" "grDevices" "utils" "dat...
[6] "methods" "base"
other attached packages:
qtl
"1.05-2"
> library(qtl)
> example(bayesint)
baysnt> data(hyper)
baysnt> ## Don't show:
baysnt> hyper <- subset(hyper, chr=c(1,4))
baysnt> ## End Don't show
baysnt> hyper <- calc.genoprob(hyper, step=0.5)
baysnt> out <- scanone(hyper, method="hk")
baysnt> bayesint(out, chr=1)
chr pos lod
c1.loc33 1 36.3 2.454921
c1.loc44.5 1 47.8 3.562391
c1.loc80.5 1 83.8 2.511073
以下略
- 再度行ったところ前述のエラーがやはり出て参りました。ses...
- おかしいですね. R2.6.0でも試しましたが，動きましたよ。...
R version 2.6.0 Under development (unstable) (2007-04-26...
Copyright (C) 2007 The R Foundation for Statistical Comp...
ISBN 3-900051-07-0
中略
'q()'と入力すればRを終了します。
[Workspace restored from /Users/***/***/***/.RData]
> library(qtl)
> library(help=qtl)
> 要するに，ここまでで何のエラーもないし，この後もちゃん...
- intel Mac でも，ok でしたよ。 -- &new{2007-04-29 (日) ...
- 手動で入れた、[/Library/Frameworks/R.framework/Versions...
**サンプルについて [#d5464f51]
>[[goto]] (2007-04-25 (水) 17:57:25)~
~
どなたか教えてください。~
以下のサンプルでは要素（球体）やノードを大きくしたり加工...
どのようにしたらこのようになるのでしょうか。~
gplot関数のパラメータを変更するのでしょうか。
http://erzuli.ss.uci.edu/R.stuff/sna/screenshots/sna_sam...
#ref(sna_sample3.png)
//
- タイトルが不適切。引用で画像のURLを書いたら表示されると...
- 説明不足の件は大変失礼致しました。snaパッケージを利用し...
- ?gplotとして、ヘルプファイルを熟読する。vertex.cex: exp...
- 引用される図があまりにも大きすぎて不都合なので，ダウン...
- 図の縮小ありがとうございます。gplotのパラメータを弄って...
- 自己レスです。gplot3という関数が用意されていました。解...
**Ｒでお絵かきツール [#e66bd6ae]
>[[ロトリング]] (2007-04-24 (火) 21:24:31)~
~
Ｒで出力したグラフに、ポリゴンやテキストをマウスで描画...
//
- > ポリゴンやテキストをマウスで描画~
その程度のことなら，テキストを指定して，描画位置をマウス...
- ご返答ありがとうございます。返事が遅れて申し訳ございま...
- ポリゴンにもよるけど，ある点がポリゴンの内か外かを判定...
- splancsパッケージにinpip()を使えば楽だと思います。マル...
- みなさん、情報ありがとうございます。ところで、最初の質...
- 単にこういうことを知りたかったとか？ -- &new{2007-04-2...
poly.xy <- locator(type = "l")
**95％信頼楕円 [#y1d9fea6]
>[[suz]] (2007-04-21 (土) 11:23:24)~
~
95％信頼区間楕円を描く方法として以下の操作は正しいでしょ...
x = rnorm(100, sd=1)
y = rnorm(100, sd=3)
plot(x, y)
library(ellipse)
r = cor(x, y)
stdev = c(sd(x), sd(y))
centre = c(mean(x), mean(y))
polygon(ellipse(r, scale=stdev, centre=centre, level=0.9...
//
- （元記事質問者と二番目のコメント投稿者は同じだと仮定し...
> V <- rbind(c(var(x), var(x,y)),c(var(x,y),var(y)))
> a <- ellipse(V, centre=centre, level=0.95)
> b <- ellipse(r, scale=stdev, centre=centre, level=0.95)
> identical(a,b)
[1] TRUE
- 注意：コメントがすでにあった場合、元記事を書き直すのは...
- 再質問します。-- &new{2007-04-22 (日) 12:47:38};
95％信頼区間楕円を描く方法として以下の操作は正しいでしょ...
set.seed(101)
x = rnorm(100, sd=1)
y = rnorm(100, sd=3)
xy = cbind(x,y)
library(ellipse)
stdev = c(sd(x), sd(y))
centre = c(mean(x), mean(y))
plot(x,y)
polygon(ellipse(var(xy), scale=stdev, centre=centre, l...
polygon(ellipse(var(xy), scale=c(1,1), centre=centre, l...
polygon(ellipse(var(xy), centre=centre, l...
polygon(ellipse(cor(x,y), scale=stdev, centre=centre, l...
polygon(ellipse(cor(x,y), scale=c(1,1), centre=centre, l...
polygon(ellipse(cor(x,y), centre=centre, l...
95％信頼区間楕円を描く方法として#1から#6のどれが正しいの...
#3と#4は同じ楕円を描きます。
scaleのデフォルトはc(1,1)ですが、#2と#3は同じ楕円を描きま...
- まあ，基本的には，その関数のオンラインヘルプをちゃんと...
- > ご指摘ありがとうございます。元記事です~
厳しく言えば，Wiki においては編集は無限の自由度があるのだ...
というか，直せば直すほど修正箇所が増えて，どうしようもな...
- ellipseのhelpを読み、exampleを実行してみました。答えの...
- 確率楕円で検索してみたのでしょうか？http://aoki2.si.gun...
- 「統計学関連なんでもあり」の過去ログ //
教えて下さいDensity Ellipse http://aoki2.si.gunma-u.ac.jp...
等確率楕円の長軸 http://aoki2.si.gunma-u.ac.jp/lecture/mb...
これは，確率楕円ですか？ http://aoki2.si.gunma-u.ac.jp/le...
二次元正規分布の確率楕円 http://aoki2.si.gunma-u.ac.jp/le...
散布図（各種の描画機能付き）散布図を描き，棄却楕円（確率...
**mfrowの値とplotmeansの動作 [#b8a84b1f]
>[[suzu]] (2007-04-20 (金) 16:32:43)~
~
グラフィックパラメータpar(mfrow=c(行数, 列数))の行数によ...
library(vegan)
data(BCI)
sp.n = specnumber(BCI)
group = gl(5, 10)
library(gplots)
par(mfrow=c(7,1)) # 行数によってplotmeansの動作が異なる
plotmeans(sp.n ~ group)
において
par(mfrow=c(6,1))
の場合、正常に動作しますが、
par(mfrow=c(7,1))
の場合、group 4のエラーバーが表示されず、以下のWarning me...
Warning messages:
1: 長さゼロの arrow は角度が不定ですので，スキップされま...
2: 長さゼロの arrow は角度が不定ですので，スキップされま...
エラーバーを正常に表示させるためには、どうすればよろしい...
使用環境 R version 2.4.0 (2006-10-03) i386-apple-darwin8....
参考文献 http://cse.naro.affrc.go.jp/takezawa/r-tips/r/68...
//
- 描画領域を大きくする。 -- [[takahashi]] &new{2007-04-20...
- アドバイスありがとうございます。「描画領域を大きくする...
**表形式でないファイルの読み込み [#qbd8961d]
>[[mori]] (2007-04-19 (木) 15:03:13)~
~
ファイルからデータを読み込む方法に関して質問があります。
[test.txt]
milk bread
bread butter
beer
milk bread butter
bread butter
上のようなtest.txtからデータを読み込んで下のように作った...
> x <- list(c("milk","bread"),
c("bread", "butter"),
c("beer"),
c("milk", "bread", "butter"),
c("bread", "butter"))
//
- read.table の fill パラメータを調べましょう。~
スマートではないが，以下のようなものも -- &new{2007-04-1...
x <- as.matrix(read.table("test.txt", header=FALSE, fill...
dimnames(x) <- NULL
apply(x, 1, function(a) a[a != ""])
こんな別解も
a <- sapply(readLines("test.txt"), function(s) unlist(st...
names(a) <- NULL
- もしくは類似解には -- &new{2007-04-19 (木) 16:28:59};
> sapply(as.list(readLines("test.txt")), function(i) str...
[[1]]
[1] "milk" "bread"
[[2]]
[1] "bread" "butter"
[[3]]
[1] "beer"
[[4]]
[1] "milk" "bread" "butter"
[[5]]
[1] "bread" "butter"
Warning message:
'test.txt' に関する readLines で不完全な最終行が見つかり...
注：上の warning は，ファイルの最後に EOF がないために出...
ちゃんとしたファイルなら warning は出ない。
- 類似解ですが，文字数が少ない(^_^;) -- &new{2007-04-19 ...
mapply(strsplit, as.list(readLines("test.txt")), split="...
- readLines()というのを知らなかったので大変勉強になりまし...
**garchFit(パッケージfSeries)のエラーについて [#d83a26c3]
>[[ケンゴ]] (2007-04-14 (土) 16:05:43)~
~
こんにちは、garchFit(パッケージfSeries)関数を使用し時系列...
おこないたいと思い以下のような処理を致しました。~
結果「エラー：!missing(data) is not TRUE」と表示されうま...
どこか間違っているのでしょうか？教えて下さい。~
使っているのはR2.4.1 OSはＸＰです。
> library(fSeries)
> UKgas
Qtr1 Qtr2 Qtr3 Qtr4
1960 160.1 129.7 84.8 120.1
1961 160.1 124.9 84.8 116.9
1962 169.7 140.9 89.7 123.3
1963 187.3 144.1 92.9 120.1
1964 176.1 147.3 89.7 123.3
1965 185.7 155.3 99.3 131.3
1966 200.1 161.7 102.5 136.1
1967 204.9 176.1 112.1 140.9
1968 227.3 195.3 115.3 142.5
1969 244.9 214.5 118.5 153.7
1970 244.9 216.1 188.9 142.5
1971 301.0 196.9 136.1 267.3
1972 317.0 230.5 152.1 336.2
1973 371.4 240.1 158.5 355.4
1974 449.9 286.6 179.3 403.4
1975 491.5 321.8 177.7 409.8
1976 593.9 329.8 176.1 483.5
1977 584.3 395.4 187.3 485.1
1978 669.2 421.0 216.1 509.1
1979 827.7 467.5 209.7 542.7
1980 840.5 414.6 217.7 670.8
1981 848.5 437.0 209.7 701.2
1982 925.3 443.4 214.5 683.6
1983 917.3 515.5 224.1 694.8
1984 989.4 477.1 233.7 730.0
1985 1087.0 534.7 281.8 787.6
1986 1163.9 613.1 347.4 782.8
> UKg.d<-diff(UKgas)
> UKg.m<-garchFit(formula=~arma(2,0)+garch(1,1),series=U...
エラー：!missing(data) is not TRUE
※条件付き平均（mean)=arma(2,0)~
※条件付き分散(variance)=garch(1,1)~
//
- 使ったことがないから適当なコメントですが、引数 formula ...
?garchFit によれば：
garchFit(formula.mean = ~arma(0, 0), formula.var = ~...
series = x, init.rec = c("mci", "uev"), delta = ...
cond.dist = c("dnorm", "dsnorm", "dged", "dsged"...
include.mean = TRUE, include.delta = NULL, inclu...
include.shape = NULL, leverage = NULL, trace = T...
algorithm = c("sqp", "nlminb", "lbfgsb", "nlminb...
control = list(), title = NULL, description = NU...
- アドバイスありがとうございます。 -- [[ケンゴ]] &new{200...
- で、アドバイスは役に立ったんですか、的はずれだったんで...
- どうも，質問の出所というか参考にしているところは同じら...
UKg.m<-garchFit(formula=~arma(2,0)+garch(1,1),data=UKg.d)
のようですよ（？） -- &new{2007-05-01 (火) 18:10:52};
**Ｒにインストールは必要ない？！ [#l0575daa]
>[[豊田秀樹]] (2007-04-14 (土) 07:03:17)~
~
信じられないようなことですが，ちゃんとインストールした~
Ｒのフォルダを単純にＵＳＢフラッシュメモリにコピーして~
おくと，そのＵＳＢフラッシュメモリを持ち歩くだけで，Ｒを...
していない計算機でＲを使えます．（そのフォルダは，全ての~
パッケージがインストールされており1ギガくらいの容量です）~
更にそのＵＳＢフラッシュメモリのＲののフォルダを単純に~
別の計算機にコピーするとその計算機でもＲが使えるようにな...
ここで「コピーで済むなら何のためのインストールだったん~
だろう？」という疑問が沸いてきました．きっとこの方法で~
コピーを繰り返すと「何か困ったことがそのうち起きる」ある...
「何か使えない高度な機能がある」としか思えません．~
文科系の学部なのでインストールはハードルが高い学生もいま...
コピーだけならだれでもでき，各学生に任せられます．~
とても便利なので，つい学生にも薦めて見たくなっていますが，~
あとで変なことが起きると困るので，どなたか「インストール~
しないと，こういう点で後で困るよ」というアドバイス~
がありましたら教えてください．それともインストールなんて~
必要ないのでしょうか？そんなことないですよね．~
//
- すみません環境・バージョンを書き忘れました．OSは Windo...
- 面白いですね。いわれてみればなるほどな、ですね。最近もV...
2.6 RをCDやUSBから起動できますか？
注意すればできます。基本的なRインストールは再配置可能（...
ですから、インストールしたRのイメージをハードディスクに...
フラッシュメモリUSBドライブ等のリムーバブル記憶装置にイ...
できます。（もしパッケージを私的なライブラリ（フォルダー...
いれば、それらの絶対パスがHTMLパッケージリストに記録され...
Rの実行には書き込み可能な一時ディレクトリとホームディレ...
で、もしそうしたものがなければ現在のディレクトリが使われ...
編成されたNTベースのWIndowsのバージョンでは何の問題も起...
もし問題があれば書き込み可能なフォルダーへのショートカッ...
ことはできないかも知れません。
- R は, 全てのOSで再配置可能です. 自分のトップ(R_HOME)が...
- 早速のアドバイスを，どうも有難うございました．了解いた...
- Rguiの引数に, R_LIBS=c:/temp 等を加えれば, パッケージの...
- 貴重な情報をありがとうございます．勉強になりました．計...
- なかまさんが言っているような「R解析キット」が分野別にあ...
- うーん, "「R解析キット」が分野別=>パッケージの追加が面...
- そうそう、私が言いたかったのは「R解析キット分野別作成キ...
- よいことを聞いたと、うれしくなり、とりあえずデスクトッ...
Warning messages:
1: 'lib.loc' 中に如何なるライブラリー木も見つかりません
in: library(package, lib.loc = lib.loc, character.only...
2: package "methods" in options("defaultPackages") was n...
起動準備中です ー 12 件の警告がありました (警告を見るに...
何かが足りないか（レジストリとかをいじらないといけないよ...
- .Rdataをダブルクリックしたまたは, コピーが不完全じゃな...
- お返事ありがとうございます。binディレクトリのRgui.exeを...
- 私の思い付く限りは, せいぜいみょーにフォルダ名の長いと...
- デスクトップの r というディレクトリに移しただけで，prof...
- 日本語ディレクトリはなるべくなら避けてほしいなぁ(希望)....
- 避けて欲しいというか，それが原因のようですね。Cドライブ...
**Windows版コンソールで使えるフォント [#sf37c914]
>[[mori]] (2007-04-12 (木) 14:32:35)~
~
Windows版のR2.4.1なのですが、GUIプリファレンスで指定でき...
//
- たとえば，Terminal とか FixedFont とか，これで，日本語...
- あれはリストボックスじゃなくって, コンボボックスだった...
- コンボボックスに直接指定すればダイジョブ。これは「HGS行...
#ref(ss2.png)
- アドバイスありがとうございます。リストにないフォントで...
**2変量散布図とヒストグラム図 [#uc13cad0]
>[[teru]] (2007-04-11 (水) 19:56:49)~
~
xとyの2変量散布図と，xおよびyのヒストグラムを書きたいので...
Ｒではどのようにすればいいのでしょうか？
**| oo
***| ooo
****| oooo
***| oo
**| oo
*| o
-------------
*********
*******
****
**
//
- 出来合いのものはないと思われるので，自分で関数を書くし...
pairs 関数の example を（実行して）見れば，望むものにちょ...
#ref(pairs.png)
- ヒストグラムではないですが、chplotパッケージがレイアウ...
- [[グラフィックス参考実例集：箱型図]]の「箱型図を軸の装...
- [[ヒストグラムと密度の推定]]の「二変量データのヒストグ...
- par()でmfrowでやるより，layout(matrix(c(1,1,4,2,2,3,2,2...
- 大変助かりました．みなさん，どうもありがとうございます ...
**ハッシュ [#x662ce31]
>[[JR]] (2007-04-11 (水) 14:24:09)~
~
Rにはperlなどのスクリプト言語でいうハッシュに相当するデー...
a => 1
b => 2
…
z => 26
というハッシュと~
v=c("a","b",…"z")
というベクトルがあったとき、vの各要素をハッシュのキーとし...
//
- こういうことでしょうか -- &new{2007-04-11 (水) 14:48:3...
> a <- 1:26
> names(a) <- letters
> a
a b c d e f g h i j k l m n o p q r s...
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19...
> a['b']
b
2
- a[names(a) == "b"] でないと最初に見付かったものしか選ば...
> x <- sample(1:100, 1e6, rep=TRUE)
> y <- sample(letters, 1e6, rep=TRUE)
> xx <- x
> names(xx) <- y
> xxx <- data.frame(x=x, y=y)
> xx["a"]
a
81
> length(xx[names(xx) == 'a'])
[1] 38482
> system.time(xx[names(xx)=="a"])
[1] 0.164 0.000 0.166 0.000 0.000
> system.time(xxx$x[xxx$y =="a"])
[1] 0.248 0.000 0.254 0.000 0.000
- ありがとうございます。ドンピシャなものは無いけど、ベク...
**区切り文字が1文字以上の空白で、#を含むNAに当たる文字が...
>[[TT]] (2007-04-07 (土) 17:30:52)~
~
こんにちは、以下のような、テキストファイルに入ったテーブ...
1 2 3 4 5
11 12 13 -1.#IND 15
これを、
tt<-read.table("hoge.txt")
で読み込もうとすると
'2' 行目には，5 個の要素がありません
というエラーが出ます。どうも2行目の4個目の要素の-1.#INDの...
tt<-read.table("hoge.txt",sep=" ")
とすると、今度は、1行目の1個目と2個目の要素の間の2個の空...
'2' 行目には，6 個の要素がありません
とエラーが出ます。何とか、1文字以上の連続した空白を1個の...
//
- hoge というのはどうも。以下のようにしたらいかが？さらな...
> x <- read.csv("foo.dat", header=FALSE, sep=" ", na.str...
> x
V1 V2 V3 V4 V5
1 1 2 3 4 5
2 11 12 13 NA 15
- read.csv は，read.table にシュガー・コートを着せたもの...
> x <- read.table("foo.dat", sep=" ", na.string="-1.#IND...
- 回答ありがとうございます。read.tableのオンラインマニュ...
> x <- read.csv("foo.dat", header=FALSE, sep=" ", na.str...
> x
V1 V2 V3 V4 V5 V6
1 1 NA 2 3 4 5
2 11 12 13 NA 15 NA
となってしまい、意図した結果になりません。もう少しご教授...
- 見れども見えず。read help. -- &new{2007-04-07 (土) 21:...
sep the field separator character. Values on each line ...
separated by this character. If sep = "" (the d...
the separator is “white space”, that is one or ...
newlines or carriage returns.
ーーーーーーーーーーーーーーーーーーーーーーーーーーーー...
> read.csv("foo.dat", header=FALSE, sep="", na.strings="...
V1 V2 V3 V4 V5
1 1 2 3 4 5
2 11 12 13 NA 15
- もしNA値に置き換えるべき項目いつも"-1.#IND" ではなく、...
> x0 <- scan("temp.txt", what=character(1)) Read 10 items
> x0
[1] "1" "2" "3" "4" "5" "...
[8] "13" "-1.#IND" "15"
> x1 <- gsub(".+#.+", "NA", x0)
> x1
[1] "1" "2" "3" "4" "5" "11" "12" "13" "NA" "15"
> x2 <- as.numeric(x1) # 直後のコメントとからもわかるよ...
> x2
[1] 1 2 3 4 5 11 12 13 NA 15
> x3 <- matrix(x2, ncol=5, by=TRUE)
> x3
[,1] [,2] [,3] [,4] [,5]
[1,] 1 2 3 4 5
[2,] 11 12 13 NA 15
> x4 <- as.data.frame(x3)
> x4
V1 V2 V3 V4 V5
1 1 2 3 4 5
2 11 12 13 NA 15
- そこまでやる必要はないのでは？~
要するに，数値を表す文字列以外を NA にするなら，以下のよ...
> x <- read.csv("foo.dat", header=FALSE, sep="", as.is=T...
> y <- as.data.frame(matrix(as.numeric(unlist(x)), nrow(...
Warning message:
強制変換により NA が生成されました
> y
V1 V2 V3 V4 V5
1 1 2 3 4 5
2 11 12 13 NA 15
無意味な空行はインデントを狂わすのでご注意
- ありがとうございました。オンラインマニュアル、まさしく...
**パッケージのインストール [#u6790519]
>[[tyamada]] (2007-04-06 (金) 18:13:53)~
~
教えてください。Rのパッケージ（たとえば、vegan)をインスト...
使っているのはR2.4.1 OSはビスタです。~
~
chooseCRANmirror()~
と打って、適当なミラーを選択し（日本の3つとも試しました）~
install.packages(c("vegan"))~
と打ちますと、はじめは調子よくいっているように見えるので...
~
URL 'ftp://ftp.u-aizu.ac.jp/pub/lang/R/CRAN/bin/windows/c...
ftp data connection made, file length 1176434 bytes~
開かれた URL~
downloaded 1148Kb~
~
以下にエラーzip.unpack(pkg, tmpDir) : ファイル 'C:/Progra...
~
などというメッセージが返されて、止まってしまいます。~
~
また、ローカルにあるZIPファイルからのインストールでも、同...
どこが間違っているのでしょうか？教えてください。~
//
- Vistaではユーザの権限管理がunix系並みに厳しくなっている...
- ありがとうございます。ご教示していただいた方法で、解決...
**ファイルタイプの識別 [#o10d8165]
>[[cam]] (2007-04-02 (月) 11:33:16)~
~
Rではsave関数でXDR形式のバイナリデータを作ることができま...
データファイルを読み込む時に、まずはファイル名からこのフ...
//
- ?data しかしこれはファイル拡張子でファイルタイプを識別...
- system 関数を使えばよろしいでしょう -- &new{2007-04-02...
> system("file test.XDR")
test.XDR: gzip compressed data, from Unix
> system("file test.data")
test.data: UTF-8 Unicode text
- ありがとうございました。参考になります。 -- [[cam]] &ne...
**作業スペースの日本語が文字化けしてしまいます [#kbd6c17f]
>[[akiba]] (2007-03-31 (土) 21:40:55)~
~
http://cran.md.tsukuba.ac.jp/から、R-2.4.1-win32.exeをイ...
文字化け部分をテキストに貼り付けると日本語としてちゃんと...
どうすれば、文字化けを直すことができるでしょうか。~
ご教示をお願いします。~
//
- まず [[日本語化掲示板]] の過去記事を読みましょう。そ...
- とりあえず直すだけなら「編集」メニューの「GUIプリファレ...
- GUIプリファレンスのFontを変えたら、直りました。ありがと...
**パッケージのインストールができません [#e02b12af]
>[[奥井正俊]] (2007-03-29 (木) 13:02:38)~
~
Rをインストールしたのちに、パッケージのインストールを試み...
//
- 現在、Aizu＆Tukubaにうまくつながっていないようです。CER...
- たぶん質問された方の環境だとプロキシ設定をしないとイン...
**Ｒの2軸プロット [#y1916d05]
>[[nick]] (2007-03-28 (水) 19:52:19)~
~
Ｒのプロットで，y軸の主軸，2軸の設定は可能でしょうか？~
色々調べたのですが，わかりません。~
どなたかご存知の方お願いします。~
//
- 例えば、plot(0,yaxt="n");axis(2,-2:2*0.1);axis(4,-4:4*0...
- お礼が遅くなりました。どうも有難うございます。 -- [[nic...
**プロットしたい行を読みとばしたい [#a23168a6]
>[[内藤]] (2007-03-28 (水) 15:41:16)~
~
Ｒの初心者です。~
プロットしたいデータが10000行あって，行を100行おきとかに...
gnuplotでいうとevery 100みたいな。~
以上，どうぞよろしくお願いします。~
//
- 例えば、 dat<-rnorm(10005); plot(dat[1:100*100]) -- [[t...
- どうもありがとうございます。助かりました。 -- [[内藤]] ...
**行列の1列目の中からある項目だけを取り出す [#z3fdd13e]
>[[norion]] (2007-03-27 (火) 23:56:46)~
~
他で作った表をRに読み込んで作業しようとしています。表全体...
行列の1列目には色々な名前が書いてあるんですが、この中から...
指定した名前に1列目が該当している行の全ての列を取り出した...
どのような関数を書けばそれが可能でしょうか。~
~
不適切な質問だったらすみません。~
//
- なんとなく行列ではなくて、データフレームではないかとい...
> M # 行、列ラベルつき行列
C1 C2 C3
R1 1 4 7
R2 2 5 8
R3 3 6 9
> M[,"C2"] # 列ラベルで第二列を取り出す
R1 R2 R3
4 5 6
> M["R2",] # 行ラベルで第二行を取り出す
C1 C2 C3
2 5 8
> str(M) # 行列なら次のように表示されるはず
int [1:3, 1:3] 1 2 3 4 5 6 7 8 9
- attr(*, "dimnames")=List of 2
..$ : chr [1:3] "R1" "R2" "R3"
..$ : chr [1:3] "C1" "C2" "C3"
> str(M) # データフレームなら次のように表示されるはず
`data.frame': 3 obs. of 3 variables:
$ C1: int 1 2 3
$ C2: int 4 5 6
$ C3: int 7 8 9
- 例えば，
> A
[,1] [,2] [,3]
[1,] "a" "1" "7"
[2,] "b" "2" "8"
[3,] "c" "3" "9"
[4,] "b" "4" "10"
[5,] "c" "5" "11"
[6,] "a" "6" "12"
> A[A[,1]=="b",]
[,1] [,2] [,3]
[1,] "b" "2" "8"
[2,] "b" "4" "10"
こんな感じでいいのでは？ -- [[g]] &new{2007-03-28 (水) 10...
- ありがとうございました。 -- [[norion]] &new{2007-03-28 ...
**関数が効かない [#peb71b5f]
>[[sunny]] (2007-03-26 (月) 23:01:23)~
~
初めて書き込みします。~
Rを使い始めてわずかの、全くの初心者です。~
barplotと入れても、表示をしてくれなくなってしまいました。~
どうやら、何かの弾みにbarplot関数が壊れてしまったようです...
壊れた関数を再設定するにはどうしたらよいのでしょう。~
//
- ls() としてみてください。もし barplot と表示されたら、...
- ls()をやってみましたが、その中にはbarplotはありませんで...
- > 表示をしてくれなくなってしまいました~
とありますが，エラーメッセージも何も表示されないのですか？~
引数なしで，barplot とだけ入力したら何か表示されますか？ ...
- Rをもう一度インストールしてみたらいかがですか。 -- &ne...
- エラーメッセージも何も表示されなかったんですが、インス...
**グラフのX,Yの刻み指定と目盛線の引き方 [#j05cfc33]
>[[のっち]] (2007-03-23 (金) 00:45:59)~
~
目的としては、CPUの負荷傾向を日々のCPU負荷状態(%)およびJO...
Rを入れて色々と調べているのですが、10, 20,30... といったY...
教えて君で申し訳ありませんが、ヒントなどいただけませんで...
//
- 目盛線については、plot(...,tck=1)で引くことが出来ました...
- abline 関数，および線種ならばその引数として渡せる lty ...
- 任意の刻みはaxisで指定しましょう。 -- &new{2007-03-23 ...
- 例として、~
plot(1:10, runif(10)*100, ylim = c(0, 100));abline(h=c(1...
で意図する刻みでの線が引けました。ありがとうございます。...
- axis のヘルプ見ました？第一引数は何だったか分かりません...
c(10,20,30,40,50,60,70,80,90)は1:9*10でよいね。 -- &new{...
- 線を引くのに悩みましたが、axis( 2, 1:9*10, lty="dotted"...
**小数以下が0を明示したいとき [#h7cb08dc]
>[[X Jr.]] (2007-03-21 (水) 21:22:39)~
~
例えば、~
y <- 1:10/2; plot(1:10, y); text(1:10, y, y, pos = 1)
とすると、小数以下が0のものは整数部分のみが表示されますが...
//
- ?format として, nsmall をみませう. -- [[なかま]] &new{2...
-
y <- 1:10/2; plot(1:10, y); text(1:10, y, format(y, nsma...
でできました。ありがとうございました。-- [[X Jr.]] &new{2...
- y <- 1:10/2; plot(1:10, y); text(1:10, y, sprintf("%.1f...
- なるほど。上の例で y <- 1:10/3 とした場合を試してみまし...
- ?formatとして, nsmallを見ると, digitsというのも(略 -- [...
**nlsで出力される回帰係数の読み込みについて [#q7aaab43]
>[[myuhe]] (2007-03-14 (水) 18:55:45)~
~
はじめまして。Rのド素人です。次のことについて、ご教示い...
~
nls()を実行すると、例えば次のような結果が得られます
Nonlinear regression model
model: inverseHeight ~ 1/(A * DBH^B) + 1/C
data: testtrees
A B C
0.2689986 3.1630292 24.6279807
residual sum-of-squares: 0.0003014474
この中の回帰係数（上の例では、0.2689986 3.1630292 24.6...
続けて行いたい処理は、出力された回帰係数で任意の独立変...
~
最初は、この結果をデータフレーム等に変換して、そこから...
~
predict()を使った方法も検討しましたが、任意の独立変数の...
~
初歩的なものかもしれませんが、よろしくお願いします。~
//
- 例えば、nls()が出力したnlsオブジェクトがhogeだとすると...
as.vector(hoge$m$getAllPars())
- coef(hoge) -- [[takahashi]] &new{2007-03-14 (水) 19:25:...
- 何事によらず R オブジェクトの構造を見たければ str() 関...
- coef()を使えば良かったのですね。str()もこれまでデータの...
- coef()を使った方法でできました。多くのアプローチまで教...
**作成したモデル（予測式）の保存・呼び出し方法 [#q17a4b19]
>[[ADMET modeler]] (2007-03-09 (金) 15:37:11)~
~
こんにちは。初めて投稿します。超初心者で非常に基本的な質...
~
仕事の関係で判別予測モデル（手法としてはBayesian, Random ...
~
たとえば、Random Forestの場合、
hlm <- read.csv('C:/work/data.csv',head=TRUE)
trainset <- hlm[1:1952,]
testset <- hlm[1953:2439,]
model <- randomForest(Class ~ ., data = trainset)
というふうにモデルを作成して、~
そのままpredictすると別のデータを予測できますよね。
pred <- predict(model, testset)
しかし、できればここで作成したモデルを保存しておいて、後...
~
この例では、Random Forestですが、一般的な回帰モデルなども...
//
- Rのオブジェクトをそのままディスク上に保存しておくには、...
- あるいは、終了時にそのまま「ワークスペースを保存」して...
- [[ファイルを読み込む tips 集(暫定版)]] で save.image() ...
**連続する整数をまとめる方法 [#s2982f8c]
>[[きゃらっと]] (2007-03-08 (木) 08:31:34)~
~
初めて投稿させていただきます．初心者です．~
(1:3)が(1,2,3)やとなる変換の逆を行わせる方法を教えてくだ...
例~
(1,2,...100, 250,251,...500, 10001, 10002,... 10005)~
変換後~
(1:100, 250:500,10001:10005)~
よろしくお願いします．~
//
- 関数を書くしかないのではないでしょうか？その関数は，以...
vecpack<-function(x) {
S<-x[-1]
E<-x[-length(x)]
TF<-E-S!=-1
s<-c(x[1],S[TF])
e<-c(E[TF],x[length(x)])
myfun<-function(x) {
if(s[x]==e[x])
as.character(s[x])
else
paste(as.character(s[x]),as.char...
}
vlist<-lapply(1:length(s),myfun)
paste("c(",paste(vlist,collapse=","),")")
}
早速付け加えてくれたようですね~
以下のように考えたのですが，ちょっと美しくないけど，せっ...
> set.seed(836478)
> ( x <- unique(sort(sample(100, 90, replace=TRUE))) ) #...
[1] 2 3 5 6 7 8 12 14 15 16 18 19 20 24 28 29 30 ...
[25] 40 42 45 49 50 52 53 54 55 59 61 62 63 65 66 68 69 ...
[49] 82 83 84 86 88 90 91 92 93 97 98 99
> n <- length(x) # ベクトルの長さ
> s <- c(TRUE, x[-n]+1 != x[-1]) # 新しいシークエンスが...
> res <- "" # 結果を得る文字列
> for (i in 1:n) { # いやらしいループで
+ if (s[i]) { # 新しいシークエンスが始まると文字列に「,...
+ res <- paste(res, x[i], sep=", ")
+ last <- x[i] # 最後に加えた整数値を覚えておく
+ }
+ if (i < n) { # シークエンスが終わるときは文字列に「：...
+ if (s[i+1] && x[i] != last) { # ただし，k:k のような...
+ res <- paste(res, x[i], sep=":")
+ }
+ }
+ else if (x[i] != last) { # 最後の要素の処理
+ res <- paste(res, x[i], sep=":")
+ }
+ }
+ res <- paste(res, ")") # 最後に「)」を付けて
> ( res <- sub(",", "c(", res) ) # 先頭が「,」で始まるの...
[1] "c( 2:3, 5:8, 12, 14:16, 18:20, 24, 28:31, 33, 35:40...
59, 61:63, 65:66, 68:74, 77, 80, 82:84, 86, 88, 90:93, 9...
vecpack の vlist lapply を使わなくてもできるので，より分...
vecpack2 <- function(x) {
S <- x[-1]
E <- x[-length(x)]
TF <- E-S != -1
s <- c(x[1], S[TF])
e <- c(E[TF], x[length(x)])
vlist <- ifelse(s==e, s, paste(s, e, sep=":")) # これだ...
paste("c(",paste(vlist,collapse=","),")")
}
巧い!
vecpack3 <- function(x)
{
s <- x[c(TRUE, x[-n]+1 != x[-1])]
e <- x[c(x[-n]+1 != x[-1], TRUE)]
return(paste("c(",paste(ifelse(s == e, s, paste(s, e, s...
collapse=","),")"))
}
マイナーな修正ですが，未定義の n とスペースを除いて
vecpack4 <- function(x)
{
s <- x[c(TRUE, x[-length(x)]+1 != x[-1])]
e <- x[c(x[-length(x)]+1 != x[-1], TRUE)]
return(paste("c(",paste(ifelse(s == e, s, paste(s, e, s...
collapse=","),")", sep=""))
}
- 早速のご返答ありがとうございました．早く皆さんのように...
**実験計画法 [#ia70d912]
>[[ちゃい]] (2007-03-05 (月) 21:56:05)~
~
Rで実験計画法ができるのでしょうか~
//
- 「実験計画法」という言葉で，あなたは何を意味しています...
実験計画を立ててくれというなら，それはたぶん無理なんじゃ...
- CRAN の貢献パッケージリストを最初から読む、または RSit...
- 有難うございました -- [[ちゃい]] &new{2007-03-06 (火) 2...
**integrate()で一様関数を積分できますか？ [#ie859bd7]
>[[藤巻十三]] (2007-02-28 (水) 17:52:22)~
~
こんにちは。~
表題の一様関数の積分は例えで、具体的にやりたいことは他に...
f<-function(x) 1
integrate(f=f,lower=0,upper=1)
とすると答えとして「１」を期待したいところでしたが、
以下にエラーintegrate(f = f, lower = 0, upper = 1) :
evaluation of function gave a result of wrong len...
となります。この原因と回避する方法を教えて頂けませんか？~
尚、「f<-function(x) ceiling(x)」なら答えは１が返ってきま...
「f<-function(x) max(x,1)」では駄目です。~
//
- すみません。わかりました。 -- [[藤巻十三]] &new{2007-02...
- なぜエラーになるかはおわかり（ヘルプ文章に書いてあるよ...
- そうじゃないです。fが返すベクトルの長さの問題です。f<-f...
- 失礼しました．なるほど、離散点での関数値を同時にサンプ...
**リスト内の各ベクトルから特定の要素のみを抽出する方法 [#...
>[[生物系大学院生]] (2007-02-26 (月) 06:11:57)~
~
こんにちは、初めて質問させていただきます。~
下記のようにして作成したブートストラップデータのリスト(z)...
> x<-c(1:100)
> y<-function(){sample(x,100,replace=TRUE)}
> z<-lapply(1:100,function(i)try(y()))
塩基配列データに基づくパラメータの推定値の分散をブートス...
~
各ベクトルについては、
> z[[1]][z[[1]]<11]
として11未満の数字のみを抜き出す方法や、
~
> sort(z[[1]])
として先頭にinformative siteを集める方法にはたどり着いた...
どなたかいい方法をご存知ないでしょうか？
//
- もう一歩。初心者ということならまずは for ループを使う...
> z <- lapply(1:100, function(i) sample(1:100, 100, rep=...
> A <- list(rep(NULL,100)) # 空リスト
> for (i in 1:100) A[[i]] <- z[[i]][z[[i]] < 11]
> B <- lapply(1:100, function(i) z[[i]][z[[i]] < 11])
> identical(A, B)
[1] TRUE
- 早速解法を教えていただきどうもありがとうございました。
新しいリストを構築するということまで考えがおよびませんで...
教えていただいた方法で思い描いていた結果が得られました。
どうもありがとうございました。 -- [[生物系大学院生]] &new...
**陰関数の表示について [#a6580f67]
>[[もも]] (2007-02-24 (土) 17:27:16)~
~
単発の質問で申し訳ないのですが~
陰関数を表示する関数、あるいはパッケージをご存知の方教え...
//
- 「陰関数を表示する」とはどういう意味ですか。 -- &new{2...
- すみません。自己解決しました。等高線＝０に気づいていま...
- そのノウハウを[[知っているといつか役にたつ(?)関数達(2)]...
**他言語で作成したプログラムのコンパイル方法について。 [#...
>[[Y.Takenaka]] (2007-02-18 (日) 16:52:55)~
~
mt19937の乱数の改良版が発表されていましたので、ラッパーを...
~
altivecやsse2と行ったものに対応していないものでは使えるよ...
http://sun.s167.xrea.com/Archives/HomePage.html (r_sfmt19...
//
- もう少し状況が分かるように説明した方がよろしいのでは？ ...
- コンパイルオプションを付けるなら,
R CMD COMPILE CFLAGS=-msse2 a.c b.c c.c
R CMD SHLIB -o abc.so a.o b.o c.o
こんな感じになります.
R CMD COMPILE -h
R CMD SHLIB -h
を見てください.-- [[なかま]] &new{2007-02-18 (日) 19:05:2...
- 情報ありがとうございます。2つめの方法でうまくいったです...
- Rで使えるようにと言う意味ならば, help(Random.user)する...
- なるほど、おっしゃるように関数のオーバヘッドは大きい事...
**自己回帰係数の算出方法 [#fb0da746]
>[[千成]] (2007-02-16 (金) 13:17:45)~
~
ARモデルの自己回帰係数の算出方法として、~
Yule-Walker法、最小二乗法、最尤法、Burg法などがありますが...
~
Yule-Walker法が他の方法に比べてどの点で秀でているのでしょ...
助言をお願い致します。~
//
- ここで尋ねるような話では無いですね。時系列の本を（何冊...
- やはりそうですか…。ありがとうございます。 -- [[千成]] &...
**整数の割り算について [#qa6ada38]
>[[廣瀬_敏之]] (2007-02-15 (木) 02:16:56)~
次のような結果をどう考えればいいのでしょう？
98-BASICの整数の割り算の結果とあまりに違うのでとまどって...
どちらが正しいか、という問題でない事はわかっていますし、...
各々のround(x,0)をとってから割り算をしていると考えても少...
この計算の定義までマニュアルに書いてあったでしょうか？
いわゆる『Ｒ言語のマニュアル』と、過去の『初心者向けＱ＆...
使用環境はＤＯＳ／Ｖ系のパソコンで、『ＸＰ』と『Ｖｉｓｔ...
R version 2.2.0, 2005-10-06, i386-pc-mingw32
attached base packages:
[1] "methods" "stats" "graphics" "grDevices" "uti...
[7] "base"
です。本当は複数のバージョンでチェックすべきなのでしょう...
なお『Ｒ』自体のソースコードをどうぞ、と言われましても、...
> -2.3%/%1
[1] -3
> test<-c(-5.99999,-5.00001,5.00001,5.99999)
> test%/%1
[1] -6 -6 5 5
> test%/%2
[1] -3 -3 2 2
> test%/%2.5
[1] -3 -3 2 2
> test%/%2.2
[1] -3 -3 2 2
> test%/%2.7
[1] -3 -2 1 2
> test%/%-1
[1] 5 5 -6 -6
> test%/%-2
[1] 2 2 -3 -3
> test%/%-2.5
[1] 2 2 -3 -3
> test%/%-2.2
[1] 2 2 -3 -3
> test%/%-2.7
[1] 2 1 -2 -2
- もうしわけありませんでした。helpで分かるそうです。 -- [...
- 演算記号に対してまでhelpが使えるとは思いませんでした。...
- %/% はいわゆる整数商ですが、98-basic （懐かしい響き！）...
**plot()のy軸ラベルと目盛りの表示について [#ff2f6052]
>[[kokko]] (2007-02-14 (水) 18:22:11)~
~
x <- 1:5
y <- c(100,1000,10000,100000,1000000)
plot(x,y, log="y", xlab="X_num", ylab="Y_num", las=1)
こんにちは。このようにy軸を対数目盛にしたグラフを作製した...
また目盛りの値が "1 e+05" のように1とeの間にスペースが入...
どなたかよろしくお願いします。~
//
- ＞ y軸のラベル "Y_num"を時計回りに45°回転させてx軸と平行~
時計回りに 90 度回転ではないですか？それはともかく，text ...
＞ "1 e+05" のように1とeの間にスペースが~
お使いの R のバージョンと，OS の種類およびバージョンは？~
私は R 2.5.0 で Mac ですが，スペースが入っているようには...
もし，どうしても気になるということなら，axis 関数で，文字...
図や表の体裁も，自分ならこのように描く（書く）のにどうし...
- ご回答ありがとうございます。余白に書き込むためにはmtext...
**Windows版 help() の表示方法の変更 [#pbadd398]
>[[X Jr.]] (2007-02-03 (土) 21:31:51)~
~
R-2.4.0以降のWindows版で、インストール時にヘルプの表示形...
//
-Rprofile.siteと言うファイルがR_HOMEのetc以下にあると思い...
- 遅くなりましたが、ありがとうございました。私はヘルプをH...
**wilcox.exact で，P 値が 1 より大 [#ica1e20e]
>ウィリアム (2007-02-02 (金) 07:57:57)~
~
Wilcoxon signed rank test を行おうと思いましたが、同じ値...
> wilcox.exact(c(8,8,-5,-10))
Exact Wilcoxon signed rank test
data: c(8, 8, -5, -10)
V = 5, p-value = 1.125
alternative hypothesis: true mu is not equal to 0
P 値が 1 より大きくなってよいのですか？~
//
-貢献パッケージ中の関数(ですよね)に言及する時はパッケージ...
-すみません。書くのを忘れました。 >library(exactRankTest...
-質問の意図が、バグ報告なのか、どうしても計算したいのに困...
**指定されたベクトルのサイズが長すぎます [#l8bf0ae9]
>[[szk]] (2007-02-01 (木) 09:03:07)~
~
以下のエラーを解決する方法をお尋ね申し上げます。
> dist(rnorm(560327))
以下にエラーvector("double", length) : 指定されたベクト...
以下は正常に計算されます。
> dist(rnorm(560327)[1:10])
//
-dist(1:10) が何を返すかは理解できていますか?~
dist((1:10)[1:3]) とはどこが違うか分かりますか?~
あなたが要求した操作には（実数を格納するには）560327×5603...
そして，あなたのコンピュータはそれだけのメモリを使えます...
-2.3テラバイト! 東工大の誇るグリッドスパコン TSUBAME の総...
-1.2Tぐらいかと. TSUBAMEは共有型クラスタではなかったと思...
ページ名:
![[RjpWiki]](image/../Rlogo.png "[RjpWiki]")
