初級Ｑ＆Ａ アーカイブ(1)
をテンプレートにして作成
[
トップ
] [
新規
|
一覧
|
検索
|
最終更新
|
ヘルプ
]
開始行:
COLOR(green){SIZE(20){初心者のための R および RjpWiki に...
新規投稿はできません
----
#contents
----
**丸め方向の指定 [#tf1abe4a]
>[[abc]] (2004-11-30 (火) 19:25:46)~
~
windows版のRを使用しているのですが、丸めモードの変換を行...
どのようにすればよろしいのでしょうか？教えてください。~
よろしくお願いします。~
//
-かなり意味不明。 -- &new{2004-11-30 (火) 20:22:47};
-「丸めモードの変換」ってのがよく分かりませんが、丸めた結...
-浮動小数点数を_controlfp(_RC_DOWN,_MCW_RC); -- &new{200...
-もしかして、round, floor, ceil 関数の比較という意味です...
-何を意図しているかはわかりませんね。誠意のないところを見...
-確かに何を意図している質問かわかりませんが，あまりにも初...
-R(LAPACKも)はIEEE演算の丸めを基に開発されていますので,CP...
**エクセルへ値コピー [#qc29c6e2]
>[[ABA]] (2004-11-29 (月) 23:47:56)~
~
Ｒで数値データを一列に出力させたとき、値の横に一行目には[...
２行目には[2,]・・・と表示されるのですが、出力されたデー...
エクセルにコピー＆ペーストするには、どうすればよいのでし...
色々と試行錯誤したのですが、分からないので教えてください。~
宜しくお願いします。~
//
-Windows環境だとすれば，Rで数値データがベクトルXに入って...
-Macな人は，apply(X, 1, function(i) cat(paste(i, "?t"), "...
d.list <- function(d)
{
invisible(apply(d, 1, function(i) cat(paste(i, c(rep...
}
質問を読み返してみたら，数値データはベクトルに入っている...
-ありがとうございます。中澤先生と青木先生の両方の手法で解...
-ご丁寧な対応，ありがとうございます。わかっている人には，...
**標準化回帰係数 [#g5920031]
>[[B]] (2004-11-26 (金) 13:58:42)~
~
重回帰をした際に標準化回帰係数を知りたいのですが、summary...
//
-[[ここ:http://phi.med.gunma-u.ac.jp/swtips/R.html#MISC]]...
-発想の転換で，scale を使って標準化したデータセットで同じ...
d がデータフレームだとして，
ds <- data.frame(scale(d)) # 標準化する
lm(y ~ x1+x2+x5, d)
lm(y ~ x1+x2+x5, ds) # 第二引数だけが違う
# 厳密に言えば，後者のモデルの最後に -1 を付けておくのが...
前者は回帰係数，後者は標準化回帰係数が得られる(^_^) -- [[...
-ありがとうございます。両方ともやってみます -- [[B]] &new...
**rimage インストール時のエラー [#i91a88cf]
>[[ささかま]] (2004-11-26 (金) 00:10:45)~
~
R2.0.1 mac版を使っています。OS10.3.6 Xcode 1.1を入れてあ...
rimageをインストールするため、Finkを使ってfftwとlibjpegを...
checking jpeglib.h usability... no
checking jpeglib.h presence... no
checking for jpeglib.h... no
configure: error: Sorry, can't find jpeglib header
と出てしまい、インストールできませんでした。fftwをcheckす...
この掲示板の使い方がよくわからず上手く書き込めませんでし...
~
//
-BTW(by the way)ですが，[[Fink:http://fink.sf.net/]]のペ...
-レスありがとうございます。Xcodeを1.5にしてgccのパッチも...
-こんにちは。ただいま、ささかまさんと同じく、Mac版R2.1にr...
**統計的有意の記号 [#z4e59028]
>[[y.a.]] (2004-11-20 (土) 22:37:54)~
~
例えば，回帰分析の結果に対して summary.lm()を signif.star...
オプション付きで実行するとt値の各有意水準に応じて星をつけ...
この星の水準を変更するにはどうすればよいでしょうか？~
具体的には，10%, 5%, 1% 水準でそれぞれ，*, **, ***~
をつけたいのです（経済学ではこれがデフォルトです）．~
~
summary.lm をみると，係数の行列は ans$coefficients ですが...
これに星をつけている箇所はよく分かりませんでした．~
//
-私も不思議(経済学の標準も不思議ですが、経済データでは 1%...
For example, calling summary(lmObject) dispatches on met...
hwich creates an object of class "summary.lm". The latte...
method print.summary.lm() which calls printCoefmat().
The stars are hard-coded there, and I don't think anybod...
that. I suggest to turn of the printing of siginificant ...
print(summary(.....), signif.stars = FALSE) or by settin...
option().
Uwe Ligges
It would be possible to re-define 'printCoefmat' private...
the lines
cutpoints = c(0, 0.001, 0.01, 0.05, 0.1, 1),
symbols = c("***", "**", "*", ".", " "))
towards the end of its code into whatever you prefer, e.g.
cutpoints = c(0, 0.01, 0.05, 0.1, 1),
symbols = c("***", "**", "*", " "))
or
cutpoints = c(0, 0.001, 0.01, 0.05, 0.1, 1),
symbols = c("****", "***", "**", "*", " "))
(both of which are compatible with your description of w...
The most straightforward way of redefining it is to copy...
'printCoefmat' into a file, e.g.
sink("printCoefmat.R")
printCoefmat
sink()
and then edit that file.
NOTE that the code written to the file does not include
the name of the function, i.e. it starts
function (x, digits = max(3, getOption("digits") - 2),...
so the first modification has to be
printCoefmat<-function(x, digits = .... )
Then, when you want your private version, simply do
source("printCoefmat.R")
and it will overlay the original version. (Experts will ...
whether this clashes with any "namespace" issues.
On my reading of the code, it doesn't seem to; but I'm n...
If your friend wants to use this new definition all the ...
then one way to arrange this is to put the revised funct...
definition (as in the edited file) into his .Rprofile,
or put the command
source("printCoefmat")
into that file.
Best wishes,
Ted.
-早速のお返事ありがとうございます． -- [[y.a.]] &new{2004...
経済学（や他の社会科学？）では，5%水準で有意であればだい...
ところで，上記の上書きする方法を試したのですが，どうもう...
以下の例では，printCoefmat が上書きされているにも関わらず...
> tail(printCoefmat, 15)
74 Signif <- symnum(pv, corr = FALSE, na...
75 cutpoints = c(0, 0.001, 0.01, 0.05,...
76 symbols = c("***", "**", "*", ".", ...
（以下略）
> data(swiss)
> summary(lm(swiss$Fertility ~ swiss$Education + swiss$E...
(前略)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 85.2533 3.0855 27.630 <2e-16 ***
swiss$Education -0.5395 0.1924 -2.803 0.0075 **
swiss$Examination -0.5572 0.2319 -2.402 0.0206 *
---
Signif. codes: 0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1...
(略)
> # 上記の方法でprintCoefmatを上書きする
> source(file="printCoefmat.R")
> tail(printCoefmat, 15) # 上書きできているか確認
74 Signif <- symnum(pv, corr = FALSE, na...
75 cutpoints = c(0, 0.01, 0.05, 0.1), ...
76 "**", "*")) ...
(略)
> summary(lm(swiss$Fertility ~ swiss$Education + swiss$E...
(前略)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 85.2533 3.0855 27.630 <2e-16 ***
swiss$Education -0.5395 0.1924 -2.803 0.0075 **
swiss$Examination -0.5572 0.2319 -2.402 0.0206 *
---
Signif. codes: 0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1...
(略)
-私の理解を越えていますので単なる推測ですが第二回答者が書...
-今回のパズルは前にも首をかしげた R のメカニズムに関係す...
-ありがとうございます．どうもprint.summary.lmが呼び出して...
-r-help への回答記事の続きです。-- [[間瀬茂]] &new{2004-1...
第二回答者のやりかたでは駄目との U. Ligges 氏(第一回答者)...
Ted, it "clashes"! Functions in the namespace are looked...
そしてこうすれば良いという G. Grothendieck 氏(r-help の常...
True, but one can still get the effect by using assignIn...
example, run these two lines (the body(...) <- line is j...
here. You want to ultimately replace that line with you...
printCoefmat: printCoefmat <- function... as discussed ...
body(printCoefmat) <- parse(text = "cat('Greetings from ...
## 次の命令が namespace を操作する魔法の呪文らしい？！
assignInNamespace("printCoefmat", printCoefmat, "stats")
Now running summary.lm as shown below displays the desir...
R> example(lm)...snip...
R> summary(lm.D90)
Call:lm(formula = weight ~ group - 1)
Residuals: Min 1Q Median 3Q Max -1.071...
Coefficients:Greetings from printCoefmat!!!
Residual standard error: 0.6964 on 18 degrees of freedom
Multiple R-Squared: 0.9818, Adjusted R-squared: 0.97...
-例えば print.summary.lm 等の関数コードがそのままでは表示...
-他にも, stats:::print.summary.lm とすると内部の名前空間...
-G. Gro氏の提案でうまくいきますね。勉強になりました。 -- ...
T. H氏のやりかたで plotCoefmat 関数（の一部）を次のように...
cutpoints = c(0, 0.01, 0.05, 0.1, 1),
symbols = c("TottemoYuui!", "KanariYuui!", "MaamaaYuui!"...
次に
> assignInNamespace("printCoefmat", printCoefmat, "stats...
> example(lm)
--- snip ---
> summary(lm.D9)
Call:
lm(formula = weight ~ group)
Residuals:
Min 1Q Median 3Q Max
-1.0710 -0.4938 0.0685 0.2462 1.3690
Coefficients:
Estimate Std. Error t value Pr(>|t|) ...
(Intercept) 5.0320 0.2202 22.850 9.55e-15 Tottemo...
groupTrt -0.3710 0.3114 -1.191 0.249 Yuuijya...
---
Signif. codes: 0 `TottemoYuui!' 0.01 `KanariYuui!' 0.05...
Residual standard error: 0.6964 on 18 degrees of freedom
Multiple R-Squared: 0.07308, Adjusted R-squared: 0.02...
F-statistic: 1.419 on 1 and 18 DF, p-value: 0.249
注意：こうして printCoefmat 関数を変更すると rm(printCoef...
- getAnywhere(print.summary.lm) 等としてもコードが得られ...
-namespace メカニズムはアドオンパッケージ等で、他のパッケ...
-Signif. codes: 0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1...
という行は，分野の違いにかかわらず，共通認識だと思います。~
記号の付け方の慣習が違うというのは私はみたことがありませ...
そもそも，記号で有意性を示すのは前近代的であり，P値を明示...
-namespace，奥が深いですね… どうも，ありがとうございまし...
**ラベルの縦書き [#b3829d75]
>[[学生]] (2004-11-19 (金) 19:30:47)~
~
日本語版R 2.00(Win版)を使用しています．~
~
グラフを描くに，x軸やy軸のラベルを縦書きにしたいです．~
以前，[[掲示板:http://www.okada.jp.org/RWiki/index.php?cm...
どなたか，ご存知の方，ご助言頂ければ幸いです．~
//
-こんなのじゃ駄目ですよねぇ・・・ -- &new{2004-11-20 (土...
> plot(1:10, type="n", ylab="")
> text(3,3, "A?nB?nC?nD?nE")
-基本的にはこうした日本 specific な工夫は労多くして益少な...
まず R で
x <- rnorm(10)
pictex("test.tex")
plot(x)
dev.off()
これで作られたファイル test.tex の最初に次を加える
?documentclass[a4paper]{jarticle}
?usepackage{pictex}
?begin{document}
当然最後に次を加える
?end{document}
次に test.tex からx,y軸ラベルに相当するところを見つけ編集...
?put {Index} [rB] <0.00pt,0.00pt> at 169.59 5.17
?put {x} [lB] <0.00pt,0.00pt> at 9.72 154.73
をたとえば
?put {インデックス} [lB] <0.00pt,0.00pt> at 169.59 5.17
?put {?shortstack[r]{ス??[-7mm] ク??[-7mm] ッ??[-7mm] エ...
<0.00pt,0.00pt> at 2 154.73
に変更する。at 以下はテキストをおく位置の x,y 座標ですか...
#ref(test.png)
-丁寧なご解説ありがとうございました．MS Wordしか使ったこ...
-TeXの勉強から始めるくらいなら，PowerpointとかOpenOffice....
-おかげさまで，emf形式で保存して，MS Power pointで，「挿...
-文字化けのWORDのバージョンにもよると思いますが,ツール+オ...
-Word2000を使用しています．ツール＋オプション＋編集と日本...
-「Drawとかに貼り付け」とか「emf形式で保存」とは MSW だけ...
-emf自体はOpenOffice for Linuxでも編集可能です. -- [[なか...
**Decisiontreeはどうすればできるのでしょうか？ [#y8d65551]
>[[MI]] (2004-11-19 (金) 17:51:39)~
~
ＲでDecisiontreeができると人づてに聞きました。~
調べるとtreeのような感じですが、よくわかりません。~
どなたか、ご存知ないでしょか？~
//
-使ったことはありませんが，rpartというパッケージでできる...
-[[Jin's Page:http://su10.sgu.ac.jp/~jin/2003/DA/tree/ind...
-Jin's Pageっていいですね。 -- &new{2004-11-19 (金) 22:4...
-読みやすいですよね． -- &new{2004-11-19 (金) 23:03:17};
-忘れてましたが，英語に抵抗がなければ，http://www.liacc.u...
-ありがとうございます。Jin'sPageは丁寧にまとめてあり初級...
**Scatterplot3Dについて [#s559f803]
>[[みち]] (2004-11-15 (月) 16:12:57)~
~
色々調べてみても分からなかったので、こちらに投稿させて頂...
~
１．Scatterlot3D で Type=h として、プロットした点から底...
~
２．Scatterlot3Dで軸ラベルを回転させたり移動させたりする...
~
宜しくお願いします。~
//
-type=h でプロットした後、par(new=TRUE) として、今度は t...
-ありがとうございます。やってみましたら、それらしくなりま...
**複数条件下でのベクトルの取り出し [#h1016980]
>[[P]] (2004-11-13 (土) 18:14:19)~
~
初歩的な質問と思いますが，あれこれと調べてどうしても分か...
~
例えば，dataという配列もしくはデータフレームがあったとし...
data[,3][data$A==1][data$B==2]~
では上手く行かないのですが，どのようにすればよいのでしょ...
ご教示いただければ幸いです．~
//
- たとえば行列なら -- [[QDU]] &new{2004-11-13 (土) 18:28:...
> x = matrix(c(2,1,2,2,2, 3,4,3,4,3, 5,6,7,8,9), 5,3)
> x
[,1] [,2] [,3]
[1,] 2 3 5
[2,] 1 4 6
[3,] 2 3 7
[4,] 2 4 8
[5,] 2 3 9
> x[x[,1]==2 & x[,2]==3, ]
[,1] [,2] [,3]
[1,] 2 3 5
[2,] 2 3 7
[3,] 2 3 9
> x[x[,1]==2 & x[,2]==3, 3]
[1] 5 7 9
データフレームなら
> x = data.frame(C1=c(2,1,2,2,2), C2=c(3,4,3,4,3), C3=c(...
> x
C1 C2 C3
1 2 3 5
2 1 4 6
3 2 3 7
4 2 4 8
5 2 3 9
> x[x$C1==2 & x$C2==3,] # 結果はデータフレームになる
C1 C2 C3
1 2 3 5
3 2 3 7
5 2 3 9
> x[x$C1==2 & x$C2==3,]$C3
[1] 5 7 9
-大変参考になりました． -- [[P]] &new{2004-11-13 (土) 20:...
**macでソースファイルからパッケージをインストールする方法...
>[[ささかま]] (2004-11-12 (金) 20:33:39)~
~
R2.0.0 mac版を使用しています。OSは 10.3.6です。xcode1.0を...
lpSolveというパッケージを導入したいのですが、CRANにはmac...
In file included from commonlib.c:17:~
commonlib.h:6:20: malloc.h: No such file or directory~
make: *** [commonlib.o] Error 1~
ERROR: compilation failed for package 'lpSolve'~
というエラーが出てしまい、出来ませんでした。~
ソースを手動でダウンロードし、ターミナルから~
R CMD INSTALL lpsolve~
と実行しても同様のエラーが出てしまいました。ルートになっ...
Rを使い始めて数週間ですし、UNIXは素人です。色々と検索して...
//
-/usr/include/malloc/malloc.h を /usr/local/include にコ...
-ありがとうございます。教えて頂いたようにしましたら無事イ...
**RExcelのインストール [#c657f290]
>[[NK]] (2004-11-10 (水) 14:44:56)~
~
R(D)COM Server V1.2を今まで使っていて、V1.2をアンインスト...
//
-ロード済のDLLを置き換えられないと言っているので該当ファ...
-お忙しいところ有り難うございました。お蔭様でインストール...
**出力画像のサイズを変更したい [#nadfffa4]
>[[fujim]] (2004-11-04 (木) 22:15:36)~
~
初めまして。近頃Rをはじめようと思い、R-tips,データサイエ...
などを読みつついろいろと試しているものです。~
環境：Windows2000,R2.0.0jp 実行方法：test.Rファイルを作成...
分からない点：例えばグラフィックス参考実例集barplotの項の...
調べた点：parパラメータを以下のように変更してみた。~
1) par(pin=c(5,10))~
2) par(fin=c(5,10))~
いずれも.outファイルにエラーが出て画像が出力できなくなり...
エラーは以下の通りです。~
Error in plot.new() : Figure region too large~
Execution halted~
Warning message: ~
calling par(new=) with no plot~
宜しくお願いします。~
//
-見つけました。済みません。png("outfilename.png",width=60...
-便乗質問ですが、pngの場合はwidthやheightの単位がpxになっ...
-こんな感じでしょうか．
| 関数 | 単位 | default width & height | defaul...
| bitmap() | inch | 6 | ‐（適当）|
| bmp() | pixel | 480 | 12 |
| dev2bitmap() | inch | 6 | ‐ |
| jpg() | pixel | 480 | 12 |
| pdf() | inch | 6 | ‐ |
| pictex() | inch | w=5, h=4 | なし |
| png() | pixel | 480 | 12 |
| postscript() | inch | ‐ | ‐ |
| quartz() | inch | 5 | 12 |
| windows() | inch | 7 | 12 |
| win.graph() | inch | 7 | 12 |
| win.metafile()| inch | 7 | 12 |
| win.print() | inch | 7 | 12 |
| x11() | inch | 7 | 12 |
| X11() | inch | 7 | 12 |
| xfig() | inch | ‐ | ‐ |
本に載せるいいネタを戴きました．ありがとうございます！ --...
**スクリプト実行 [#u6e4a442]
>[[taku]] (2004-11-04 (木) 14:19:56)~
~
perlやrubyの様に、一行目に次のようなRのパスを書いて、~
一連のコマンドをスクリプトとして実行したいのですが、(ファ...
#!/usr/bin/R~
~
ARGUMENT './hello' __ignored__~
というメッセージが出て、普通にRのプロンプトが出てしまいま...
~
Rプログラムをスクリプトとして実行するにはどうしたらいいの...
//
-スクリプトとして実行なんかできるんですか？~
まずもって，R --help で，オプション確認 -- &new{2004-11-...
-この話は r-help 等で繰返し議論されています。一度 batch ...
-いや，だから。バッチと「Rプログラムをスクリプトとして実...
taku さんは，どういうスクリプトファイルを用意したのですか...
#!/usr/bin/R
x <- rnorm(2000)
print(mean(x)
などとして，それを chmod 755 hello なんかしても，./hello ...
例えば prog というファイルにやりたいことを書いておいて
$ cat prog
x <- rnorm(2000)
print(mean(x))
そこから入力して R を実行する
$ R --vanilla --quiet < prog
> x <- rnorm(2000)
> print(mean(x))
[1] -0.007029891
こういうことをやりたいのでしょう？~
quiet を slave にすれば，入力ファイルからの入力行は出力さ...
-[[PerlスクリプトからRを起動する]]などがヒントになるでし...
-いろいろありがとうございます。
どうやら、#!/usr/bin/Rというのはできないみたいですね。
perlやrubyを使う身としては、実行のたびにいちいち'R --vani...
このあたりは何とかして欲しいものですが…
結局、今のところはshを使うのがよさそうです。
#!/usr/bin/sh
R --vanilla --slave <<EOF
x <- rnorm(50)
y <- rnorm(x)
plot(x,y)
EOF
-面倒ならば，
alias rb 'R --vanila --slave ?!*'~
しておいて，
rb < ファイル名
で動くようにしておいたらいかが。~
sh のなかでいちいち
#!/usr/bin/sh
R --vanilla --slave <<EOF
なんたら
EOF
なんて書く方が面倒だ。 -- &new{2004-11-05 (金) 11:12:45};
-ああ、こういう人もいるんじゃないかと思ってはいましたが。...
-どうぞ。おすきなように。Cプログラムを記述してコンパイル...
内容が
R --vanilla --slave < $1
な実行ファイルを/usr/bin/Rsとでもしておけば，
#!/usr/bin/Rs
で動く． -- &new{2004-12-01 (水) 11:15:44};
- 今は Rscript があります。#!/usr/bin/Rscript -- &new{20...
**関数の中身を見る関数 [#e8ccf8da]
>[[akosuke]] (2004-10-29 (金) 10:45:01)~
~
自分で作った関数で名前はわかるけど中身を忘れてしまった、...
ここでも調べてみたのですが調べ方が悪いのか見つかりません...
くだらない質問で申し訳ありませんがどうぞよろしくお願いし...
//
-コンソールにその関数の名前だけを入力してみる.... これは...
-おわかりになりましたか？ -- &new{2004-11-04 (木) 14:36:...
**AICの計算 [#t0d6c154]
>[[AIC]] (2004-10-26 (火) 21:50:37)~
~
The R Bookの第5章を読んでいる者です。p.109の上から2行目で...
また、AICのパラメータ数はこの場合npram<-2と定義しています...
場合、式（5.14）にあるようにω、α、βの3と定義するべきでは...
それから（長くてすみません）GARCHの派生モデルを"fSeries"...
gjr(data,order(p,q,γ))？？)と思ってやってみたのですができ...
ご存知でしたら教えていただけないでしょうか。~
//
-本来ならば「The R bookご意見」でいただくべきご質問ですが...
-まず1番目のご質問は以下のようにしてください。
for (p in 1:3){
for ( q in 1:3){
tmp <- garch(dnikkei, order=c(p,q))
cat("p=", p, "q=",q, ", AIC=", (tmp$n.likeli+p+q),"?...
}
}
warningなどが出て少し見づらいですが、とりあえず。-- [[矢...
-2番目のご質問に対しては、AICはモデルの当てはまりの良さに...
-3番目のご質問のfSeriesに関しては詳しくないのですが、apar...
-ありがとうございます。上のプログラムで実行した結果、1番...
-３番目の質問に関しては、fSeriesをインストールしたあと、...
-fSeriesを用いたGJR推定は以下のようになります。
The R book第5章のデータを使った例(GJR(1,1))を示します。
library(fSeries)
tmp <- aparchFit(x=dnikkei,
order = list(alpha.lags = 1, beta.lags = 1, delta...
opt = list(gamma = TRUE, delta = FALSE,
disparm = FALSE), doprint = FALSE)
tmp
ポイントは上でgamma=TRUEとなっている部分です。この
部分がGJR定義式でのγを示します。gamma=FALSEとすると
ただのGARCHになります。GJRで日経225を推定すると非対称性が...
分かります。aparchFitについてより詳しくは
help(aparchFit)
としてhelpの内容を参照してください。-- [[矢野]] &new{2004...
-懇切丁寧な説明、ありがとうございます。頑張ってみます。 -...
**0を埋める [#z99a4268]
> (2004-10-25 (月) 18:52:25)~
~
とってもFAQだと思うのですが…~
1 -> 001，42 -> 042 ~
のように有効桁数（この例だと3ケタ）を決めて，~
それに満たない場合は0で埋めるにはどうすればいいでしょうか...
//
-sprintfでこのサイトを検索してみて下さい。 -- [[谷村]] &n...
-ありがとうございました．formatC という関数がありました：
x <- c(1:20)
formatC(x,width=2,flag="0")
-関数の紹介ありがとうございます。
> x <- c(1:20)
> formatC(x,width=2,flag="*")
flag にいろいろな文字を指定するとおもしろい結果が出てきま...
本来は，無効な文字ははじかなくてはならないのでしょうが。~
ベクトル化されているという点では sprintf より使いでがある...
**変数名に別の変数の内容を含めたい [#wc8d2de9]
>[[質問者]] (2004-10-23 (土) 15:43:00)~
~
大変初歩的な質問と思いますが、うまくいかないのでよろしく...
month <- c("Jan","Feb","Mar","Apr","May","Jun",
"Jul","Aug","Sep","Oct","Nov","Dec")
for (i in 1:12) paste("Hoge.",month[i],".total",sep="") ...
として、例えばHoge.Jan.totalに0が格納されるようにしたいの...
//
-[[RプログラミングTips大全]] に定石構文が紹介されています...
> month <- c("Jan","Feb","Mar","Apr","May","Jun",
"Jul","Aug","Sep","Oct","Nov","Dec")
> for (i in 1:12) eval(parse(text=paste("Hoge.",month[i]...
> Hoge.Apr.total
[1] 0
ついでですが次のようにしても良いわけです。
> for (i in 1:12) eval(parse(text=paste("Hoge",month[i],...
また R は月名を表す組み込みの文字列ベクトルを最初から持っ...
>month.abb
[1] "Jan" "Feb" "Mar" "Apr" "May" "Jun" "Jul" "Aug" "Sep...
> month.name
[1] "January" "February" "March" "April" "May...
[7] "July" "August" "September" "October" "Nov...
-ありがとうございました。evalも使って試行錯誤もやってので...
**データフレームの分割処理について [#xddecfdc]
>[[R-Fan]] (2004-10-18 (月) 14:38:21)~
~
このサイトでRを知りました。いつも有益な情報をありがとうご...
表題についてですが、大きなデータをキー項目（下記例ではdat...
(案1)
for (key in sort(unique(maindata$date))){
subdata <- subset(maindata, date == key)
#subdataの処理
}
自己レスになりますが、subset()を使用するかわりにキー項目...
(案2)
idx <- tapply(1:nrow(maindata), maindata$date,
function(subidx){return (subidx)})
for (key in sort(unique(maindata$date))){
subdata <- maindata[idx[[toString(key)]],]
#subdataの処理
}
もっと根本的に改善する方法がありましたらご教授いただけれ...
どうぞよろしくお願いいたします。
//
-ケースにもよるでしょうけど,unique(maindata$date)をlevels...
-ご回答いただきましてありがとうございます。後者のご提案に...
> d <- data.frame(x = factor(c(1,1,1,2,2,3)), y = 10:15)
> d[levels(d$x) == 1,]
x y
1 1 10
4 2 13
-次のようになるのでは？
> d[levels(d$x)[1]==d[,1],]
x y
1 1 10
2 1 11
3 1 12
>
> d[levels(d$x)[2]==d[,1],]
x y
4 2 13
5 2 14
> d[levels(d$x)[3]==d[,1],]
x y
6 3 15
データファイルが分析のたびに付け加わっていくのでなければ...
そうそう，思い出した。by という関数，速いか遅いか知りませ...
- levels()関数の例は大変よく分かりました。ご紹介いただい...
(案3)
by(maindata, maindata$date,
function(subdata){
#subdataの処理
}
)
**関数の読み込み [#x20303d3]
>[[恐縮です]] (2004-10-17 (日) 00:24:05)~
~
Rを使い始めて間もない者です。~
Rで色々な関数をinstall.packages()でインストールするのです...
一度インストールした関数をそのまま使い続けるにはどう保存...
//
-library(パッケージ名） という使い方ができないということ...
-すみません。library(パッケージ名)を実行して使えるように...
-library は毎回，必要なときにやらないといけないでしょう。...
-OSを書かれていませんので、Linuxを前提に説明するとinstall...
-library()を毎回実行するのが面倒なら、青木先生の助言のよ...
-普通，ライブラリのインストール時には管理者権限でやるので...
-.First()関数に初期化処理を定義しておく方法はだめでしょう...
-参考URL http://www.okada.jp.org/RWebRef/RManualJP/Startu...
-R 2.0 からは、たしかパッケージは実際に使う際になって初め...
**積集合(INTERSECT)と和集合(UNION)の中間的なものはありま...
>[[撃墜王]] (2004-10-14 (木) 15:39:55)~
~
R言語には、INTERSECTとUNIONといったコマンドが存在します。~
ただ、入力されるリストが非常に多く存在した場合に、「それ...
また、「それら全てのUNION」を取ると、当然ながらアウトプッ...
そこで、質問させて頂きたいのは、INTERSECTとUNIONの中間的...
~
なお、自分の中のイメージとしては~
~
入力されるリスト 100個~
1リストあたりに含まれる変数名(Aomori,Miyagi,Hokkaido...et...
全てのINTERSECTを取った場合 最終的に残る変数名 3個(少な...
全てのUNIONを取った場合、 最終的に残る変数名 1024個(多...
~
そこで、単純なUNIONではなくて、1つのリストにしか含まれな...
具体的な例をさらにあげますと~
~
list1 Kanagawa,Tokyo,Saitama~
list2 Chiba,Kanagawa,Tokyo~
list3 Tokyo,Irabaki,Gunma~
~
Intersect = Tokyo,Union=Kanagawa,Tokyo,Saitama,Chiba,Ibar...
となってしまいますが、上のようなリストの場合Kanagawaはlis...
//
-ある要素が，いくつのリストに含まれているかを勘定し，リス...
> top.n <- function(lst, n)
+ {
+ result <- table(unlist(lst))
+ result[result >= n]
+ }
> lst <- list(no1=c("aa", "bb", "cc"),
no2 =c("bb", "dd", "ff", "hh"),
no3=c("cc", "ff", "xx", "yy"),
no4=c("aa", "cc"))
> top.n(lst, 3)
cc
3
> top.n(lst, 2)
aa bb cc ff
2 2 3 2
探すのがめんどいから作っちゃった。関数にするまでもないが...
-青木先生有難うございます。先生の手に掛かるとあっというま...
**CSVファイルのインポートで errorが。 [#pa3d6bf6]
>[[竹内]] (2004-10-10 (日) 19:42:29)~
~
こんにちは。Ｒを使い始めてまだ5日目の初心者です。使用環境...
Error in file(file, "r") : unable to open connection
In addition: Warning message:
cannot open file `test.csv'
R-Bookや googleで一応は調べた（つもり）なのですが、 error...
//
-ファイルが開けないといっているのです。どこにそのファイル...
Tools --> Change Working Directory で，そのファイルがある...
-さっそくのご回答ありがとうございました。教えていただいた...
**frequencyについて [#w6a1889e]
>[[R初心者]] (2004-10-07 (木) 23:49:37)~
~
Rで株価の時系列分析をしようとしている者ですが、frequency...
//
-土日が無いなら freq = 5 とする他無いのでは。祭日は NA で...
> ts(1:10, frequency = 5)
Time Series:
Start = c(1, 1) # 始点は第一週の第１日と解釈
End = c(2, 5) # 終点は第二週の第５日と解釈
Frequency = 5
[1] 1 2 3 4 5 6 7 8 9 10
-ご質問は株価データが日次データ（土日祝日データなし）が何...
-御返事有難うございます。矢野先生のおっしゃる通り株価が何...
-えーっと、まず「先生」というのはやめてください。先生では...
-それでfrequencyを250にすると良い理由ですが、それは簡単で...
-どうもありがとうございます。納得しました。ただ何年か前ま...
**GARCHモデルのパラメータ推定 [#wd87da8d]
>[[ピエロ]] (2004-10-07 (木) 22:18:58)~
~
RでGARCHモデル分析する際、パラメータ推定は何の方法を用い...
擬似最尤法とかいうやつでしょうか？~
~
あと、Rでの時系列分析の細かな手法を扱っている本あるいはサ...
宜しくお願い致します。~
//
-ソースを御覧になるより仕方がないのでは。こうした質問にす...
-help(garch)とすると"maximum-likelihood estimates of the ...
-それと時系列解析の大まかな手法でよければ、矢野が書かせて...
~
時系列解析法:http://web.sfc.keio.ac.jp/~kogure/courses/20...
-それと上のコメントの方がおっしゃるとおり、Rの日本でのユ...
-それとRjpの以下のページは非常に便利ですばらしいです（作...
~
[[Rの基本パッケージ中の時系列オブジェクト一覧:http://www....
-↑を作成されているのは間瀬先生ですね。ご尽力ありがとうご...
-間瀬先生でしたか。改めてありがとうございます。 -- [[矢野...
-実は私の先生は時系列の専門家（藤井光昭先生）なのですが、...
-矢野先生の書かれた本の第5章、読ませて頂きました。5.7.2節...
-こちらも同じで、「先生」というのは勘弁してください -- [[...
-plot(predict(ng))というのはパラメータを決定したgarchモデ...
-御返事有難うございます。矢野さん。上記のplot(predict(ng)...
-ちょっとは，自分で調べましょう。 -- &new{2004-10-15 (金...
-私の勉強不足で初歩的な質問になってしまい、申し訳ありませ...
-縦軸が2003年1月から12月までになってしまうのは昔からある...
**pairsのプロット範囲 [#a4f0ce6f]
>[[ぷーすか]] (2004-10-07 (木) 20:17:03)~
~
散布図行列で、プロット点をドットに変え、y=xの直線を書きも...
~
pairs(x,panel=function(x,y){points(x,y,pch=".");abline(0,...
~
でうまくできました。~
ここでさらに各図のプロット範囲を同じにしようとxlim、ylim...
~
pairs(x,xlim=c(-3,3), ylim=c(-3,3), panel=function(x,y){p...
~
ではエラーが出てしまいます。~
プロット範囲以外のカスタマイズはあきらめ、~
~
pairs(x,xlim=c(-3,3), ylim=c(-3,3))~
~
とすれば、warningは出ますが、一応プロット範囲は制御できま...
~
プロット内容とプロット範囲のカスタマイズを両立させる方法...
//
-エラーメッセージはパネル関数は追加パラメータを受け付けな...
data(USJudgeRatings)
pairs(USJudgeRatings[1:3],
panel=function(x,y, ...) {points(x, y, pch="."); a...
xlim=c(5, 10), ylim=c(5, 10) )
#ref(pairspanel.png)
-ありがとうございます。panelに引数の追加を指定しつつも、x...
**eclipseの使い方 [#i5f83d3a]
> (2004-10-04 (月) 19:18:27)~
~
eclipse (eclipse-SDK-3.0-linux-gtk.zip)をダウンロードし、...
//
-自己フォローです。File->New->R Script fileでエディタが現...
**エクセルファイルのデータをRに取り込む方法は？ [#h1ace67d]
>[[初心者]] (2004-09-27 (月) 11:51:31)~
~
Rの初心者ですが、どうやってエクセルファイルのデータをRに...
//
-このページの下の方に，「CSVファイルのインポート」と「Rcm...
-出来たら初心者の為のデータインポートとか初心者が躓きやす...
-R の公式マニュアルの 「R Data Import/Export」に他システ...
-最近の [[なんでも掲示板]] の青木さんの関数も使えるかも知...
-みなさんありがとうございました。何とかなりそうです。でも...
-なぜか繰り返される質問と言うのは現状のドキュメントでは不...
-「よろしくお願いいしマース〜〜」などと甘まったれた文句を...
-「2万5000くらいのレコード」ということは，2万5000個の数値...
質問者の常だが，「できるのでしょうか」なんて聞くより，や...
初心者には優しくしましょうね(^_^)（発言タイトルはなおして...
-ごめんなさい！いろいろ試してみたんだけど、だめだったので...
-いつもこのようなメッセジーが出て、困っています。Warning ...
-初心者の一人ですが，なかまさんのアドバイスに従って，「エ...
-なるほど,MsysのシェルにPATH追加してRとか打ってるようでは...
----
~
**ヒストグラムＸ軸目盛について。 [#o01284d8]
>[[tofu]] (2004-09-21 (火) 11:11:26)~
~
度数分布をヒストグラムにしたときののＸ軸の目盛の数字、目...
//
-axis 関数で，任意の場所に任意の文字列を書くことができま...
-せっかく [[グラフィックス参考実例集]] があるのですから、...
-すみません…いろいろ勉強不足なもので。アドバイスありがと...
-実例集のページもずいぶん大きくなったのですね。見るのも大...
簡単に一例を示すと，以下のようになるか（縦軸は自分で考え...
x <- rnorm(2000)
hist(x, axes=F)
axis(1, at=c(-1.96,0,1.96), labels=c("-1.96 SD", "Zero",...
#ref(hist3.png)
-丁寧な解説ありがとうございました。Ｒとともにこのサイトの...
**plotでドットだけを打ちたい [#g180f43c]
>[[KabuTaro]] (2004-09-20 (月) 04:02:13)~
~
plot(c(1,2,3),c(4,5,6),type="p")~
とすると(1,4),(2,5),(3,6)に丸印でプロットされますが~
これを単にドットだけにするにはどうしたらよいでしょうか？~
pchを弄くっても形が変わるだけでうまくいきません。~
よろしくお願いいたします。~
//
-ドットだけの意味が分かりませんが，http://cran.r-project....
-pch="."ですね。できました。ありがとうございました。これ...
-plot(1,pch="A",cex=34) と plot(1,pch=".",cex=34) を比べ...
-勝手注： 中間さんのコメントの通り pch="." は真正のドット...
#ref(pchtest.png)
-なるほど、pch="."には特別な意味があるんですね。てっきり...
**package(gregmisc)のheatmap.2の引数設定について [#da53a4...
>[[ちょろべぇ]] (2004-09-14 (火) 18:57:31)~
~
gregmiscパッケージのheatmap.2を用いてヒートマップを作図し...
heatmap.2は様々なパラメータ設定が可能であり、自分としては...
(希望)~
?distfun→距離関数だと思われるので、ユークリッド距離を計算...
?hclustfun→average linkageで結合したい~
?できればdivisiveではなくagglomerativeでクラスタしたい~
~
ただ、
heatmap.2(data,distfun=dist(method="euclidian"),
hclustfun=hclust(dist(method="euclidian"),meth...
dendrogram="both")
のように設定しても、hclustfun部分でエラーが帰ってきてしま...
かなり愚かな質問だとは思いますが、アドバイスありましたら...
//
-こうした質問にはエラーメッセージをそえるのが礼儀作法です...
Hclust <- function(d) hclust(d, method = "average", memb...
heatmap.2(data, hclustfun=Hclust, dendrogram="both")
# もしくは一度に済ませるには
heatmap.2(data, hclustfun=function(d) hclust(d, method =...
dendrogram="both")
-こんな風なエラーメッセージが出てしまいますね --[[ちょろ...
Error in as.matrix(x) : Argument "x" is missing, with no...
-euclidian?euclidean? どちらが正しいのでしょうか？~
ちょっと違いますがこんな感じでどうでしょう-- [[○田五郎]]&...
my.dist <- function(x) dist(x, method="binary")
my.hclust <- function(d) hclust(d, method="ward")
hm <- heatmap(blablabla, distfun=my.dist, hclustfun=my.h...
-確認しないでコメントしました(最初のコメントのバグを修正...
library(gregmisc)
data(mtcars)
x <- as.matrix(mtcars)
rc <- rainbow(nrow(x), start=0, end=.3)
cc <- rainbow(ncol(x), start=0, end=.3)
my.hclust <- function(d) hclust(d, method = "average")
hv <- heatmap.2(x, col = cm.colors(256), scale="column",
hclustfun=my.hclust,
RowSideColors = rc, ColSideColors = cc, ...
xlab = "specification variables", ylab= ...
main = "heatmap(<Mtcars data>, ..., scal...
tracecol="green", density="density", den...
# 形式 hclustfun=function(d) hclust(d, method = "average...
#ref(hm2.png, left)
-諸先輩方、親切なアドバイスありがとうございました。ご教授...
**エラーバー [#a79042cd]
>[[Kei]] (2004-09-11 (土) 16:54:52)~
~
実験のデータをエラーバーで表現したいですが、R言語でどのふ...
mybar <- function(X) {
n <- nrow(X)
df <- n-1
X.mean <- apply(X,2,mean)
X.var <- apply(X,2,var)
d.error<-qt(0.975,df)*sqrt(X.var/n)
d.error.x <- (0.2 + 1) * 1:5 - 1/2 # エラーバーを表...
barplot(X.mean, names.arg=colnames(X),
ylim=c(0,ceiling(max(X.mean+(d.error)))),
col=c("blue", "green", "red", "yellow"))
arrows(d.error.x, X.mean, d.error.x, X.mean+ d.error, ...
}
mybar(分析するデータ[,列：列])
軸の調整方についてちょっと質問したいですが、教えてくださ...
//
-自分で作れるというのは素晴らしいことですが、こうした問題...
-縦軸をどのように調整したいのでしょうか？~
描画される範囲を，必要・十分に確保するというとことでしょ...
僭越ながら，[[私の書いた関数:http://aoki2.si.gunma-u.ac.j...
要点は，描画する前にある軸（今回は縦軸）において，描画対...
なお，必要なら，描画範囲を決めるだけで，軸も実際には描か...
plot(c(x軸の最小値, 同最大値), c(y軸の最小値, 同最大値),
type="n", # 描画しない
xaxt="n", # x 軸を描かない（必要なら）
yaxt="n", # y 軸を描かない（必要なら）
その他のパラメータ)
:
points(...)
lines(...)
axis(...)
その他諸々の描画関数
健康のために使いすぎに注意。 -- [[青木繁伸]] &new{2004-09...
**棒グラフの重ね合わせ [#l6008258]
>[[TETSU]] (2004-09-10 (金) 11:51:50)~
~
２つの棒グラフをbarplot関数を用いて描画したのですが、これ...
//
-barplot の beside という引数について，よく調べてください...
それでもよくわからなかったら，[[グラフィックス参考実例集...
このサイトで何かを調べようと思ったら，[[単語検索:http://w...
-ありがとうございます！ただ、besideをTrueにするとグラフは...
-「色々調べてみた」-> まず ?barplot でたくさんある引数を...
barplot(VADeaths, border = "dark blue", space=c(-0.5))
#ref(bp.png)
-青木先生、QDUさん、ありがとうございます。全ては私の日本...
-polygon 関数を使って，自分でヒストグラムを描く関数を書け...
**層別の相関値 [#ge316df8]
>[[HOJO]] (2004-09-10 (金) 11:39:45)~
~
欠損値を含む二つの変数について，相関値をcor(x,y,use=pai...
//
-上位30件というのはどのように決めたのでしょうか。~
それはともかく，データの散布図を描けば問題解決の手がかり...
下の図の右上30個ほどの点は相関係数0.910，次の70個ほどは0....
#ref(plot.png)
R はグラフを描くのがとっても簡単なので，数値解析する前に...
-青木先生いつも有り難うございます。散布図の特徴は，(1)原...
-＞一般に全体の相関係数は，層別した相関係数より大きくはな...
そんなことありませんよ。例えば，原点付近と(10,10)を中心と...
> gendat2 # 任意の相関係数を持つデータを生成する関数
function(nc, r)
{
z <- matrix(rnorm(2*nc), ncol=2)
res <- eigen(r2 <- cor(z))
coeff <- solve(r2) %*% (sqrt(matrix(res$values, 2, 2,
byrow=TRUE))*res$vectors)
z <- t((t(z)-apply(z, 2, mean))/sqrt(apply(z, 2, var)*(...
z %*% chol(matrix(c(1, r, r, 1), ncol=2))
}
> d1 <- gendat2(10,0)
> d2 <- gendat2(10,0)+10
> cor(rbind(d1,d2))
[,1] [,2]
[1,] 1.0000000 0.9615385 # プールすると相関は0.962
[2,] 0.9615385 1.0000000
> cor(d1)
[,1] [,2]
[1,] 1.000000e+00 -1.010639e-16 # プールする前は無相関
[2,] -1.010639e-16 1.000000e+00
> cor(d2)
[,1] [,2]
[1,] 1.00000e+00 3.27047e-16 # プールする前は無相関
[2,] 3.27047e-16 1.00000e+00
極端な場合，層別では負の相関，全体では正の相関ということ...
＞（10000,600)付近に他と離れたデータが１件~
というのは，すさまじいような気がしますが。その1つのデータ...
ともかく，HOJOさんが今扱っているデータではそのように計算...
散布図をここにペーストしてくれると，後学のためにもなるか...
-＞そんなことありませんよ。例えば，原点付近と(10,10)を --...
-＞そんなことありませんよ。例えば，原点付近と(10,10)を．...
-#ref(D:?My Documents?40joho?scatterplot1.bmp) -- &new{...
-この wiki に画像を張り付けるには、jpg か png 形式で、ペ...
-すこし極端な例を。全体では相関 0, 層別すると相関が 1, -1...
#ref(cor.png)
-#ref(scplot1.jpeg) -- [[HOJO]] &new{2004-09-11 (土) 20:4...
#ref(scplot1.jpeg)
-「編修」を使ってようやくできました。 アドバイスありがと...
-この図から見るにやはり「外れ値」とみなすべき例外的データ...
-青木先生，ＱＤＵさん他の皆様方アドバイス有り難うございま...
**単純主効果 [#i4379192]
> (2004-09-10 (金) 02:07:41)~
~
2元配置の分散分析（完全無作為）をする場合に~
交互作用がある場合，~
１．「単純主効果の検定」，~
２．「単純主効果の多重比較」~
をしたいのですが，Rではどのよなプログラムを書けば良いので...
~
参考までに，
SASでの単純主効果についての情報と，そのデータを以下にリン...
例：http://koko15.hus.osaka-u.ac.jp/~kano/lecture/faq/int...
データ：http://mat.isc.chubu.ac.jp/R/tech.html#CRFpq~
ご存知の方，教えて頂ければ幸いです~
~
//
-単純主効果という言葉の意味が一寸わからないのですが、R の...
-aov関数では得られないようです．単純主効果の説明は，http:...
-CRAN の検索では simple main effect + anova では一件も引...
**Rcmdrでグラフの変数選択欄 [#wb1d1e2a]
>[[HOJO]] (2004-08-31 (火) 16:06:53)~
~
たびたびですみません．RcmdrでIndexPlotなどを用いてグラフ...
//
-変数選択欄に出るのは数値として認識出来た列のみになるかと...
-なかまさん，何度もお手を煩わして申し訳ありません．因子か...
-変数毎にas.vectorを行うと，文字データになってしまいます...
-こうした質問はもとのデータ（長過ぎるなら問題になりそうな...
-どうしても意図した通によみこみたいならscanを使ってくださ...
-数値そのものに，あまり意味があるとは思えないので掲載しな...
-scanを使用して解決できました。有り難うございました。 -- ...
-もとのデータがただの数値で、単純に因子化されているならば...
-詳細に記載すると，以下の処理をしてからscanを施し解決して...
-詳細に記載すると，csvファイルを作成する際に事前に以下の...
-(1)は，カンマ区切りの数値は引用符で囲まれるから（そうし...
(2)は文字として % が実際にファイルに書き出されているため...
(3)空セルは文字通り，なんにも出力されないので，ケースあた...
全て，Excel に罪はないとも言える（制御のためのオプション...
いちいち手作業で(1),(2),(3)の前処理を行うのも面倒なので（...
-解説有り難うございました．相変わらず，やみなべ（Ｒ）をつ...
**CSVファイルのインポート [#y06b905e]
>[[HOJO]] (2004-08-31 (火) 11:26:52)~
~
エクセルで作成した２００行，２０列程度の数値データのみの...
１． Error: unable to open connection~
２． data is not a data frame and cannot be attached~
これはどのような点に気をつければよいのでしょうか。~
//
-多分ですが,日本語化してないRで日本語ファイル名のファイル...
-なかまさん，いつも有り難うございます．ご指摘のとおりでし...
**非心パラメータ [#xe80473c]
> (2004-08-24 (火) 19:12:06)~
~
正確な用語がわからないのですが~
SASではTNONCT関数で算出できる，~
「t分布の非心パラメータ（引数は，x, df, prob）を求める関...
Rにないでしょうか？~
http://support.sas.com/91doc/getDoc/lrdict.hlp/a000245960...
~
例えば，~
x <- 3.1; df<- 58; prob <- .975;~
?(x,df,prob)~
> 1.04844~
と求めたいのです．~
//
-そうした関数はないですが、要するに pt(3.1, 58, ncp) = 0....
> TNONCT <- function(x, df, prob, interval = c(-30, 30)) {
temp <- function(ncp) pt(x, df, ncp...
uniroot(temp, interval)[[1]] }
> TNONCT(3.1, 58, 0.975)
[1] 1.048444
-同じようなものになりましたが，参考までに。 -- [[青木繁伸...
> TNONCT <- function(x, df, prob, interval=c(0, 37.62))
+ {
+ uniroot(function(ncp) pt(x, df, ncp)-prob, interva...
+ }
> TNONCT(3.1, 58, 0.975)
[1] 1.048444
-非心分布など普段意識して使ったことがないので、いまひとつ...
-ありがとうございます．おかげさまで，目的の計算ができまし...
-私も，非心度を求めるようなことはやったことがありませんで...
> TNONCT(3.1, 58, 0.975)
[1] 1.048444
> TNONCT(-3.1, 58, 0.025)
[1] -1.048444
という関係があるようです。interval=c(-37.62, 37.62) とし...
質問者が示した URL での解説から，SAS でも同じような関数を...
google で site:r-project.org ncp pt で，もろ回答がありま...
[[R Help 2001:http://www.r-project.org/nocvs/mail/r-help/...
これでは，c(-1e7,1e7) としていますが，QDU さんの言うよう...
**行列のグラフィカル表示image()の色に任意の値を割り当てる...
>[[川村健介]] (2004-08-18 (水) 13:14:01)~
~
お世話になります。~
行列をグラフィカルに表示するimage()関数の色調を任意に割り...
例えば配列の値が0〜0.88の間にあった場合，0.0〜0.5に白を割...
//
- 色見本付きのイメージ図は example(filled.contour) 参照。...
-ありがとうございます。col=grey(c(1,1,1,1,1,0.4,0.3,0.2,0...
-数値と濃淡が正しく対応しているか慎重に確認して下さい。な...
**plot()で軸の数値をなくしたい [#uba2f2b7]
> (2004-08-18 (水) 00:22:26)~
~
10ほど並べた散布図で、最初の１つだけ軸の数値を表示し、残...
//
par(mfrow=c(2,5))
plot(sin)
plot(sin, xaxt = "n", yaxt = "n", xlab="", ylab="")
plot(sin, xaxt = "n", yaxt = "n", xlab="", ylab="")
# あと 7 回繰り返し
上のような感じで如何でしょうか？ご参考になれば幸いです．-...
-ありがとうございます。目的の図を得ることができました。 -...
**Ｒで扱える統計手法 [#k9c2485f]
>[[石川誠]] (2004-08-16 (月) 01:26:14)~
~
お世話になります。~
Ｒで扱える統計手法を知りたいのですが、何処かに日本語で一...
//
-まさにこのサイトでしょう。[[stats(R 統計)パッケージ中の...
ほかにどのようなページがあるかは，[[ページの一覧:http://w...
-R 本体のオブジェクト数は千以上、アドオンパッケージは少な...
-Ｒ に含まれる統計手法全部を知りたいというわけでもないで...
- Rで扱える統計手法でしたら，まず見なければいけないのは [...
-え〜〜，本人から申し上げますが，あそこはまんまり最初っか...
-またまた〜、ご謙遜されてますねぇ〜。「R」＋「統計」＝「R...
**行列の値上位5つの要素の位置を知りたい [#p6a1c4f8]
>[[川村健介]] (2004-08-14 (土) 09:05:51)~
~
初めて投稿させていただきます。~
630*630の配列の値の中から，上位5位までの要素の位置(x,y)を...
which.maxだと最も高い値の位置，例えば39726と帰ってきます。~
できれば，上位から[326,52],[555,156]...といった感じで位置...
//
- 欠損値 NA がなければもう少し簡単になるのですが。(PS. も...
> y [,1] [,2] [,3]
[1,] NA 0.7122784 0.7761505
[2,] 0.47935820 0.4700606 0.4547836
[3,] 0.02593524 0.8849352 0.3914027
> yy = matrix(length(y)+1-rank(y, na.last=F), dim(y), by...
> y [,1] [,2] [,3]
[1,] 9 3 2
[2,] 4 5 6
[3,] 8 1 7
> lapply(1:5, function(i) which(yy==i, arr.ind=T))
[[1]]
row col
[1,] 3 2
[[2]]
row col
[1,] 1 3
[[3]]
row col
[1,] 1 2
[[4]]
row col
[1,] 2 1
[[5]]
row col
[1,] 2 2
-すばやいレスありがとうございます。うまくいきました。ちな...
-ここはやはりベクトルの n 番目に大きな値を求める関数を定...
# ベクトル(行列、配列も可)の n 番目に大きな値を求める関数
# タイ や NA があっても使える(はず？)
> nthmax
function(x, n) {
if ( !(n %in% 1:length(x)) )
stop( "n is not in 1:length(x)" )
if (n == 1) return(max(x, na.rm=TRUE))
else
for ( i in 1:(n-1) ) {
x <- x[x != max(x, na.rm=TRUE)]
if ( length(x) == 0 )
stop("no 'n'-th max element in x")
}
return(max(x, na.rm=TRUE))
}
> x = matrix(runif(9), 3,3)
> for (i in 1:length(x)) print( which(x==nthmax(x,i), a...
row col
[1,] 1 1
row col
[1,] 3 1
row col
(途中略)
row col
[1,] 1 2
> x[1,2] <- x[1,1] # タイを入れる
> for (i in 1:length(x)) print( which(x==nthmax(x,i), ar...
row col
[1,] 1 1 # <- 最大値が二つある
[2,] 1 2
row col
(途中略)
row col
[1,] 1 3
Error in nthmax(x, i) : no 'n'-th max element in x
> x[1,2] <- NA # 欠損値を入れる
> for (i in 1:length(x)) print( which(x==nthmax(x,i), ar...
row col
[1,] 1 1
row col
(途中略)
row col
[1,] 1 3
numeric(0) # 最後がエレガントでない(欠損値のせい)
Warning message:
no finite arguments to max; returning -Inf
-初等的なやりかたは見つけた最大値を -Inf に置き換えて繰り...
> x = matrix(runif(9), 3,3)
> which(x==max(x), arr.ind=T)
row col
[1,] 2 3
> x[which(x==max(x), arr.ind=T)] <- -Inf
> which(x==max(x), arr.ind=T)
row col
[1,] 1 3
> x[which(x==max(x), arr.ind=T)] <- -Inf
> which(x==max(x), arr.ind=T)
row col
[1,] 3 2
-----
> x[ print(which(x==max(x), arr.ind=T)) ] <- -Inf # 一行...
-わかりやすい解説ありがとうございます。岐阜大で講義がある...
-そんなに実行時間の短いやり方が必要なんですか。スローライ...
作りつけの速い関数を求めるのもいいが，それを問い合わせる...
作りつけの関数でできないことも，いろいろ付け加えられます...
こういう考え方を許容できない人は無視すればいいだけです(^_...
> x <- matrix(runif(630*630),630)
>
> max.at <- function (x, top=10)
+ {
+ result <- matrix(0, top, 3)
+ for (loop in 1:top) {
+ mx <- x[1,1]
+ for(i in 1:630) {
+ for (j in 1:630) {
+ if (x[i,j] > mx) {
+ mx <- x[i,j]
+ mxi <- i
+ mxj <- j
+ }
+ }
+ }
+ result[loop,] <- c(mxi, mxj, mx)
+ x[mxi, mxj] <- -Inf
+ }
+ result
+ }
>
> system.time(max.at(x))
[ 1] 78.30 4.31 89.11 0.00 0.00
> max.at(x)
[,1] [,2] [,3]
[1,] 466 95 0.9999976
[2,] 239 31 0.9999970
[3,] 205 620 0.9999952
[4,] 595 183 0.9999935
[5,] 500 609 0.9999923
[6,] 37 556 0.9999901
[7,] 292 3 0.9999889
[8,] 515 278 0.9999818
[9,] 539 237 0.9999802
[10,] 483 331 0.9999786
速くしたいなら，C ででもプログラムすればいいだけです。プ...
-こんなのはいかが？
max.at2 <- function (x, top=10)
{
n <- nrow(x)*ncol(x)
col <- col(x)
row <- row(x)
o <- order(x)
cbind(row[o], col[o], x[o])[n:(n-top+1),]
}
> system.time(max.at2(x))
[1] 2.52 0.18 2.80 0.00 0.00 ＃ すんばらしく，速いわ！
> max.at3(x)
[,1] [,2] [,3]
[1,] 2 3 NA # NA はこの位置に入る
[2,] 18 178 99.99992
[3,] 549 252 99.99988
[4,] 208 115 99.99891
[5,] 182 140 99.99866
[6,] 308 597 99.99823
[7,] 207 156 99.99820
[8,] 616 112 99.99811
[9,] 499 249 99.99801
[10,] 355 79 99.99794
本当は，以下のようにしたいのだが，decleasing=TRUE が効か...
max.at3 <- function (x, top=10)
{
n <- nrow(x)*ncol(x)
col <- col(x)
row <- row(x)
o <- order(x, decreasing = TRUE)
cbind(row[o], col[o], x[o])[1:top,]
}
いずれにせよ，
max.at4 <- function (x)
{
o <- order(x)
cbind(row(x)[o], col(x)[o], x[o])
}
とすれば，小さい順の結果は全部得られるので，必要なところ...
-青木先生，ありがとうございました。大変，勉強になりました...
**statsライブラリ内のheatmap関数について [#r429e061]
>[[○田五郎]] (2004-08-10 (火) 01:45:09)~
~
こんばんわ。○田です。~
この度はstasライブラリ内のheatmap関数について相談させてい...
statsライブラリをロードして、heatmapがかけることはわかっ...
?heatmap内のデモスクリプトを書くなどしてみましたが、~
heatmap関数の引数であるcol引数の種類が良くわかりません。...
自分としては正の値が赤系、マイナスの値を青系で表示したい...
デモのcol = cm.colors(256)はレインボー系、、、~
col = topo.colors(16)は緑と薄いこげ茶系と、自分の望んでい...
どうぞよろしくお願いしますm(_ _)m。~
//
-こんなのありますよー http://gdds.pharm.kyoto-u.ac.jp/le...
-アドオンパッケージ RColorBrewer には様々な配色のパレット...
rainbow(n, s = 1, v = 1, start = 0, end = max(1,n - 1)/n...
heat.colors(n)
terrain.colors(n)
topo.colors(n)
cm.colors(n)
と同じように使えます。一般的には個人で配色を決めるよりは...
## help(brewer.pal) より
mypalette<-brewer.pal(7,"Greens")
image(1:7,1,as.matrix(1:7),col=mypalette,xlab="Greens (s...
ylab="",xaxt="n",yaxt="n",bty="n")
-[[グラフィックス参考実例集：カラーパレット]] を見て下さ...
-書きたいことは，ちゃんと書きましょう。主張したいことがよ...
- 元記事投稿者が topo.colors では青い色がでないといってい...
data(volcano)
x <- 10*(1:nrow(volcano))
y <- 10*(1:ncol(volcano))
image(x, y, volcano, col = c(topo.colors(12)[1:4],
heat.colors(8)[c(1,3,5,8)]), axes = FALSE)
#ref(topoheat.png)
**mvtnorm内で帰ってくる値が小数点4桁以下がその都度異なる...
>[[ちょろちゃん]] (2004-08-08 (日) 15:34:16)~
~
お世話になっております。~
Rのmultcompライブラリのsimtest、method=dunnettを~
利用しております。~
同じデータセットを入力した際に、simtestのダネットの結果中...
T値は常に一定の値が帰ってくるのですが、P値の小数点4桁より...
計算するたびに異なっております。~
どこで値が異なってしまったか追跡したところ、mulcompから呼...
mvtnorm内でP値を計算するときにどうやら異なった値が帰って...
これは仕様なのでしょうか、それともmvtnormのバグなのでしょ...
もしご存知の方がいらっしゃいましたらアドバイスいただけれ...
どうぞよろしくお願いいたします。~
//
-どこかで乱数でも使っていない限り、仕様にせよバグにせよ、...
-mvtnorm って，multcomp から呼ばれるんですか。よくわかり...
R の関数には，uniroot だっけを呼び出して解を求めるなどと...
あなたが利用した関数の定義を参照すれば，毎回異なった解が...
P値などというものは，0.05や0.01と比べて大きいのか小さいの...
いろいろな統計量に関するP値も，近似値なので（統計ソフト間...
いずれにしろ，近似値だと言うことで。そんなに目くじらたて...
-自分がテストに使ったデータは青木先生のページに記載してあ...
-灯台もと暗しでしたか。mvt 関数が疑惑の張本人ですね。主要...
? pmvt を読みますと（めんどくさいので斜め読み。従って精度...
Randomized quasi-Monte Carlo methods are used for
the computations.
とか
Value:
error: estimated absolute error
とか読めますが。~
また，pmvt の引数には，abseps = 0.001, releps = 0 という...
-なるほど。青木先生ありがとうございますm(_ _)m。感謝感激...
-青木さんのおっしゃる通り、乱数が絡んでいるようですね。試...
**繰り返しのｔ検定の結果を表示させるには [#bc1ff9ad]
>[[はまちゃん]] (2004-08-03 (火) 18:12:58)~
~
お世話になっております。
a<-matrix(rnorm(96),c(16,6))
for (i in 2:6) t.test(a[,1],a[,i])
のように、繰り返しでｔ検定を行う場合に、結果が表示されま...
なにか方法があったら、教えてお願いできますか？~
ほかの検定でも同じことが起こるのでしょうか？宜しくお願い...
//
-オブジェクトを書き出すには print を使うのが基本であると...
a<-matrix(rnorm(96),c(16,6))
for (i in 2:6) print(t.test(a[,1],a[,i]))
これで，うまくいくでしょう。~
統計学的な観点から言うと，このようなt検定の繰り返しが正し...
- -- [[MKR]] &new{2004-08-03 (火) 20:28:05};
x = sapply(2:6, function(i) {a<-matrix(rnorm(96),c(16,6)...
t.test(a[,1],a[,i])} )
-sapply は実行結果の正確な出力が得られないのでお勧めでは...
結果は
[,1] [,2] ...
statistic -0.3054952 0.3751029
終了行:
COLOR(green){SIZE(20){初心者のための R および RjpWiki に...
新規投稿はできません
----
#contents
----
**丸め方向の指定 [#tf1abe4a]
>[[abc]] (2004-11-30 (火) 19:25:46)~
~
windows版のRを使用しているのですが、丸めモードの変換を行...
どのようにすればよろしいのでしょうか？教えてください。~
よろしくお願いします。~
//
-かなり意味不明。 -- &new{2004-11-30 (火) 20:22:47};
-「丸めモードの変換」ってのがよく分かりませんが、丸めた結...
-浮動小数点数を_controlfp(_RC_DOWN,_MCW_RC); -- &new{200...
-もしかして、round, floor, ceil 関数の比較という意味です...
-何を意図しているかはわかりませんね。誠意のないところを見...
-確かに何を意図している質問かわかりませんが，あまりにも初...
-R(LAPACKも)はIEEE演算の丸めを基に開発されていますので,CP...
**エクセルへ値コピー [#qc29c6e2]
>[[ABA]] (2004-11-29 (月) 23:47:56)~
~
Ｒで数値データを一列に出力させたとき、値の横に一行目には[...
２行目には[2,]・・・と表示されるのですが、出力されたデー...
エクセルにコピー＆ペーストするには、どうすればよいのでし...
色々と試行錯誤したのですが、分からないので教えてください。~
宜しくお願いします。~
//
-Windows環境だとすれば，Rで数値データがベクトルXに入って...
-Macな人は，apply(X, 1, function(i) cat(paste(i, "?t"), "...
d.list <- function(d)
{
invisible(apply(d, 1, function(i) cat(paste(i, c(rep...
}
質問を読み返してみたら，数値データはベクトルに入っている...
-ありがとうございます。中澤先生と青木先生の両方の手法で解...
-ご丁寧な対応，ありがとうございます。わかっている人には，...
**標準化回帰係数 [#g5920031]
>[[B]] (2004-11-26 (金) 13:58:42)~
~
重回帰をした際に標準化回帰係数を知りたいのですが、summary...
//
-[[ここ:http://phi.med.gunma-u.ac.jp/swtips/R.html#MISC]]...
-発想の転換で，scale を使って標準化したデータセットで同じ...
d がデータフレームだとして，
ds <- data.frame(scale(d)) # 標準化する
lm(y ~ x1+x2+x5, d)
lm(y ~ x1+x2+x5, ds) # 第二引数だけが違う
# 厳密に言えば，後者のモデルの最後に -1 を付けておくのが...
前者は回帰係数，後者は標準化回帰係数が得られる(^_^) -- [[...
-ありがとうございます。両方ともやってみます -- [[B]] &new...
**rimage インストール時のエラー [#i91a88cf]
>[[ささかま]] (2004-11-26 (金) 00:10:45)~
~
R2.0.1 mac版を使っています。OS10.3.6 Xcode 1.1を入れてあ...
rimageをインストールするため、Finkを使ってfftwとlibjpegを...
checking jpeglib.h usability... no
checking jpeglib.h presence... no
checking for jpeglib.h... no
configure: error: Sorry, can't find jpeglib header
と出てしまい、インストールできませんでした。fftwをcheckす...
この掲示板の使い方がよくわからず上手く書き込めませんでし...
~
//
-BTW(by the way)ですが，[[Fink:http://fink.sf.net/]]のペ...
-レスありがとうございます。Xcodeを1.5にしてgccのパッチも...
-こんにちは。ただいま、ささかまさんと同じく、Mac版R2.1にr...
**統計的有意の記号 [#z4e59028]
>[[y.a.]] (2004-11-20 (土) 22:37:54)~
~
例えば，回帰分析の結果に対して summary.lm()を signif.star...
オプション付きで実行するとt値の各有意水準に応じて星をつけ...
この星の水準を変更するにはどうすればよいでしょうか？~
具体的には，10%, 5%, 1% 水準でそれぞれ，*, **, ***~
をつけたいのです（経済学ではこれがデフォルトです）．~
~
summary.lm をみると，係数の行列は ans$coefficients ですが...
これに星をつけている箇所はよく分かりませんでした．~
//
-私も不思議(経済学の標準も不思議ですが、経済データでは 1%...
For example, calling summary(lmObject) dispatches on met...
hwich creates an object of class "summary.lm". The latte...
method print.summary.lm() which calls printCoefmat().
The stars are hard-coded there, and I don't think anybod...
that. I suggest to turn of the printing of siginificant ...
print(summary(.....), signif.stars = FALSE) or by settin...
option().
Uwe Ligges
It would be possible to re-define 'printCoefmat' private...
the lines
cutpoints = c(0, 0.001, 0.01, 0.05, 0.1, 1),
symbols = c("***", "**", "*", ".", " "))
towards the end of its code into whatever you prefer, e.g.
cutpoints = c(0, 0.01, 0.05, 0.1, 1),
symbols = c("***", "**", "*", " "))
or
cutpoints = c(0, 0.001, 0.01, 0.05, 0.1, 1),
symbols = c("****", "***", "**", "*", " "))
(both of which are compatible with your description of w...
The most straightforward way of redefining it is to copy...
'printCoefmat' into a file, e.g.
sink("printCoefmat.R")
printCoefmat
sink()
and then edit that file.
NOTE that the code written to the file does not include
the name of the function, i.e. it starts
function (x, digits = max(3, getOption("digits") - 2),...
so the first modification has to be
printCoefmat<-function(x, digits = .... )
Then, when you want your private version, simply do
source("printCoefmat.R")
and it will overlay the original version. (Experts will ...
whether this clashes with any "namespace" issues.
On my reading of the code, it doesn't seem to; but I'm n...
If your friend wants to use this new definition all the ...
then one way to arrange this is to put the revised funct...
definition (as in the edited file) into his .Rprofile,
or put the command
source("printCoefmat")
into that file.
Best wishes,
Ted.
-早速のお返事ありがとうございます． -- [[y.a.]] &new{2004...
経済学（や他の社会科学？）では，5%水準で有意であればだい...
ところで，上記の上書きする方法を試したのですが，どうもう...
以下の例では，printCoefmat が上書きされているにも関わらず...
> tail(printCoefmat, 15)
74 Signif <- symnum(pv, corr = FALSE, na...
75 cutpoints = c(0, 0.001, 0.01, 0.05,...
76 symbols = c("***", "**", "*", ".", ...
（以下略）
> data(swiss)
> summary(lm(swiss$Fertility ~ swiss$Education + swiss$E...
(前略)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 85.2533 3.0855 27.630 <2e-16 ***
swiss$Education -0.5395 0.1924 -2.803 0.0075 **
swiss$Examination -0.5572 0.2319 -2.402 0.0206 *
---
Signif. codes: 0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1...
(略)
> # 上記の方法でprintCoefmatを上書きする
> source(file="printCoefmat.R")
> tail(printCoefmat, 15) # 上書きできているか確認
74 Signif <- symnum(pv, corr = FALSE, na...
75 cutpoints = c(0, 0.01, 0.05, 0.1), ...
76 "**", "*")) ...
(略)
> summary(lm(swiss$Fertility ~ swiss$Education + swiss$E...
(前略)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 85.2533 3.0855 27.630 <2e-16 ***
swiss$Education -0.5395 0.1924 -2.803 0.0075 **
swiss$Examination -0.5572 0.2319 -2.402 0.0206 *
---
Signif. codes: 0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1...
(略)
-私の理解を越えていますので単なる推測ですが第二回答者が書...
-今回のパズルは前にも首をかしげた R のメカニズムに関係す...
-ありがとうございます．どうもprint.summary.lmが呼び出して...
-r-help への回答記事の続きです。-- [[間瀬茂]] &new{2004-1...
第二回答者のやりかたでは駄目との U. Ligges 氏(第一回答者)...
Ted, it "clashes"! Functions in the namespace are looked...
そしてこうすれば良いという G. Grothendieck 氏(r-help の常...
True, but one can still get the effect by using assignIn...
example, run these two lines (the body(...) <- line is j...
here. You want to ultimately replace that line with you...
printCoefmat: printCoefmat <- function... as discussed ...
body(printCoefmat) <- parse(text = "cat('Greetings from ...
## 次の命令が namespace を操作する魔法の呪文らしい？！
assignInNamespace("printCoefmat", printCoefmat, "stats")
Now running summary.lm as shown below displays the desir...
R> example(lm)...snip...
R> summary(lm.D90)
Call:lm(formula = weight ~ group - 1)
Residuals: Min 1Q Median 3Q Max -1.071...
Coefficients:Greetings from printCoefmat!!!
Residual standard error: 0.6964 on 18 degrees of freedom
Multiple R-Squared: 0.9818, Adjusted R-squared: 0.97...
-例えば print.summary.lm 等の関数コードがそのままでは表示...
-他にも, stats:::print.summary.lm とすると内部の名前空間...
-G. Gro氏の提案でうまくいきますね。勉強になりました。 -- ...
T. H氏のやりかたで plotCoefmat 関数（の一部）を次のように...
cutpoints = c(0, 0.01, 0.05, 0.1, 1),
symbols = c("TottemoYuui!", "KanariYuui!", "MaamaaYuui!"...
次に
> assignInNamespace("printCoefmat", printCoefmat, "stats...
> example(lm)
--- snip ---
> summary(lm.D9)
Call:
lm(formula = weight ~ group)
Residuals:
Min 1Q Median 3Q Max
-1.0710 -0.4938 0.0685 0.2462 1.3690
Coefficients:
Estimate Std. Error t value Pr(>|t|) ...
(Intercept) 5.0320 0.2202 22.850 9.55e-15 Tottemo...
groupTrt -0.3710 0.3114 -1.191 0.249 Yuuijya...
---
Signif. codes: 0 `TottemoYuui!' 0.01 `KanariYuui!' 0.05...
Residual standard error: 0.6964 on 18 degrees of freedom
Multiple R-Squared: 0.07308, Adjusted R-squared: 0.02...
F-statistic: 1.419 on 1 and 18 DF, p-value: 0.249
注意：こうして printCoefmat 関数を変更すると rm(printCoef...
- getAnywhere(print.summary.lm) 等としてもコードが得られ...
-namespace メカニズムはアドオンパッケージ等で、他のパッケ...
-Signif. codes: 0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1...
という行は，分野の違いにかかわらず，共通認識だと思います。~
記号の付け方の慣習が違うというのは私はみたことがありませ...
そもそも，記号で有意性を示すのは前近代的であり，P値を明示...
-namespace，奥が深いですね… どうも，ありがとうございまし...
**ラベルの縦書き [#b3829d75]
>[[学生]] (2004-11-19 (金) 19:30:47)~
~
日本語版R 2.00(Win版)を使用しています．~
~
グラフを描くに，x軸やy軸のラベルを縦書きにしたいです．~
以前，[[掲示板:http://www.okada.jp.org/RWiki/index.php?cm...
どなたか，ご存知の方，ご助言頂ければ幸いです．~
//
-こんなのじゃ駄目ですよねぇ・・・ -- &new{2004-11-20 (土...
> plot(1:10, type="n", ylab="")
> text(3,3, "A?nB?nC?nD?nE")
-基本的にはこうした日本 specific な工夫は労多くして益少な...
まず R で
x <- rnorm(10)
pictex("test.tex")
plot(x)
dev.off()
これで作られたファイル test.tex の最初に次を加える
?documentclass[a4paper]{jarticle}
?usepackage{pictex}
?begin{document}
当然最後に次を加える
?end{document}
次に test.tex からx,y軸ラベルに相当するところを見つけ編集...
?put {Index} [rB] <0.00pt,0.00pt> at 169.59 5.17
?put {x} [lB] <0.00pt,0.00pt> at 9.72 154.73
をたとえば
?put {インデックス} [lB] <0.00pt,0.00pt> at 169.59 5.17
?put {?shortstack[r]{ス??[-7mm] ク??[-7mm] ッ??[-7mm] エ...
<0.00pt,0.00pt> at 2 154.73
に変更する。at 以下はテキストをおく位置の x,y 座標ですか...
#ref(test.png)
-丁寧なご解説ありがとうございました．MS Wordしか使ったこ...
-TeXの勉強から始めるくらいなら，PowerpointとかOpenOffice....
-おかげさまで，emf形式で保存して，MS Power pointで，「挿...
-文字化けのWORDのバージョンにもよると思いますが,ツール+オ...
-Word2000を使用しています．ツール＋オプション＋編集と日本...
-「Drawとかに貼り付け」とか「emf形式で保存」とは MSW だけ...
-emf自体はOpenOffice for Linuxでも編集可能です. -- [[なか...
**Decisiontreeはどうすればできるのでしょうか？ [#y8d65551]
>[[MI]] (2004-11-19 (金) 17:51:39)~
~
ＲでDecisiontreeができると人づてに聞きました。~
調べるとtreeのような感じですが、よくわかりません。~
どなたか、ご存知ないでしょか？~
//
-使ったことはありませんが，rpartというパッケージでできる...
-[[Jin's Page:http://su10.sgu.ac.jp/~jin/2003/DA/tree/ind...
-Jin's Pageっていいですね。 -- &new{2004-11-19 (金) 22:4...
-読みやすいですよね． -- &new{2004-11-19 (金) 23:03:17};
-忘れてましたが，英語に抵抗がなければ，http://www.liacc.u...
-ありがとうございます。Jin'sPageは丁寧にまとめてあり初級...
**Scatterplot3Dについて [#s559f803]
>[[みち]] (2004-11-15 (月) 16:12:57)~
~
色々調べてみても分からなかったので、こちらに投稿させて頂...
~
１．Scatterlot3D で Type=h として、プロットした点から底...
~
２．Scatterlot3Dで軸ラベルを回転させたり移動させたりする...
~
宜しくお願いします。~
//
-type=h でプロットした後、par(new=TRUE) として、今度は t...
-ありがとうございます。やってみましたら、それらしくなりま...
**複数条件下でのベクトルの取り出し [#h1016980]
>[[P]] (2004-11-13 (土) 18:14:19)~
~
初歩的な質問と思いますが，あれこれと調べてどうしても分か...
~
例えば，dataという配列もしくはデータフレームがあったとし...
data[,3][data$A==1][data$B==2]~
では上手く行かないのですが，どのようにすればよいのでしょ...
ご教示いただければ幸いです．~
//
- たとえば行列なら -- [[QDU]] &new{2004-11-13 (土) 18:28:...
> x = matrix(c(2,1,2,2,2, 3,4,3,4,3, 5,6,7,8,9), 5,3)
> x
[,1] [,2] [,3]
[1,] 2 3 5
[2,] 1 4 6
[3,] 2 3 7
[4,] 2 4 8
[5,] 2 3 9
> x[x[,1]==2 & x[,2]==3, ]
[,1] [,2] [,3]
[1,] 2 3 5
[2,] 2 3 7
[3,] 2 3 9
> x[x[,1]==2 & x[,2]==3, 3]
[1] 5 7 9
データフレームなら
> x = data.frame(C1=c(2,1,2,2,2), C2=c(3,4,3,4,3), C3=c(...
> x
C1 C2 C3
1 2 3 5
2 1 4 6
3 2 3 7
4 2 4 8
5 2 3 9
> x[x$C1==2 & x$C2==3,] # 結果はデータフレームになる
C1 C2 C3
1 2 3 5
3 2 3 7
5 2 3 9
> x[x$C1==2 & x$C2==3,]$C3
[1] 5 7 9
-大変参考になりました． -- [[P]] &new{2004-11-13 (土) 20:...
**macでソースファイルからパッケージをインストールする方法...
>[[ささかま]] (2004-11-12 (金) 20:33:39)~
~
R2.0.0 mac版を使用しています。OSは 10.3.6です。xcode1.0を...
lpSolveというパッケージを導入したいのですが、CRANにはmac...
In file included from commonlib.c:17:~
commonlib.h:6:20: malloc.h: No such file or directory~
make: *** [commonlib.o] Error 1~
ERROR: compilation failed for package 'lpSolve'~
というエラーが出てしまい、出来ませんでした。~
ソースを手動でダウンロードし、ターミナルから~
R CMD INSTALL lpsolve~
と実行しても同様のエラーが出てしまいました。ルートになっ...
Rを使い始めて数週間ですし、UNIXは素人です。色々と検索して...
//
-/usr/include/malloc/malloc.h を /usr/local/include にコ...
-ありがとうございます。教えて頂いたようにしましたら無事イ...
**RExcelのインストール [#c657f290]
>[[NK]] (2004-11-10 (水) 14:44:56)~
~
R(D)COM Server V1.2を今まで使っていて、V1.2をアンインスト...
//
-ロード済のDLLを置き換えられないと言っているので該当ファ...
-お忙しいところ有り難うございました。お蔭様でインストール...
**出力画像のサイズを変更したい [#nadfffa4]
>[[fujim]] (2004-11-04 (木) 22:15:36)~
~
初めまして。近頃Rをはじめようと思い、R-tips,データサイエ...
などを読みつついろいろと試しているものです。~
環境：Windows2000,R2.0.0jp 実行方法：test.Rファイルを作成...
分からない点：例えばグラフィックス参考実例集barplotの項の...
調べた点：parパラメータを以下のように変更してみた。~
1) par(pin=c(5,10))~
2) par(fin=c(5,10))~
いずれも.outファイルにエラーが出て画像が出力できなくなり...
エラーは以下の通りです。~
Error in plot.new() : Figure region too large~
Execution halted~
Warning message: ~
calling par(new=) with no plot~
宜しくお願いします。~
//
-見つけました。済みません。png("outfilename.png",width=60...
-便乗質問ですが、pngの場合はwidthやheightの単位がpxになっ...
-こんな感じでしょうか．
| 関数 | 単位 | default width & height | defaul...
| bitmap() | inch | 6 | ‐（適当）|
| bmp() | pixel | 480 | 12 |
| dev2bitmap() | inch | 6 | ‐ |
| jpg() | pixel | 480 | 12 |
| pdf() | inch | 6 | ‐ |
| pictex() | inch | w=5, h=4 | なし |
| png() | pixel | 480 | 12 |
| postscript() | inch | ‐ | ‐ |
| quartz() | inch | 5 | 12 |
| windows() | inch | 7 | 12 |
| win.graph() | inch | 7 | 12 |
| win.metafile()| inch | 7 | 12 |
| win.print() | inch | 7 | 12 |
| x11() | inch | 7 | 12 |
| X11() | inch | 7 | 12 |
| xfig() | inch | ‐ | ‐ |
本に載せるいいネタを戴きました．ありがとうございます！ --...
**スクリプト実行 [#u6e4a442]
>[[taku]] (2004-11-04 (木) 14:19:56)~
~
perlやrubyの様に、一行目に次のようなRのパスを書いて、~
一連のコマンドをスクリプトとして実行したいのですが、(ファ...
#!/usr/bin/R~
~
ARGUMENT './hello' __ignored__~
というメッセージが出て、普通にRのプロンプトが出てしまいま...
~
Rプログラムをスクリプトとして実行するにはどうしたらいいの...
//
-スクリプトとして実行なんかできるんですか？~
まずもって，R --help で，オプション確認 -- &new{2004-11-...
-この話は r-help 等で繰返し議論されています。一度 batch ...
-いや，だから。バッチと「Rプログラムをスクリプトとして実...
taku さんは，どういうスクリプトファイルを用意したのですか...
#!/usr/bin/R
x <- rnorm(2000)
print(mean(x)
などとして，それを chmod 755 hello なんかしても，./hello ...
例えば prog というファイルにやりたいことを書いておいて
$ cat prog
x <- rnorm(2000)
print(mean(x))
そこから入力して R を実行する
$ R --vanilla --quiet < prog
> x <- rnorm(2000)
> print(mean(x))
[1] -0.007029891
こういうことをやりたいのでしょう？~
quiet を slave にすれば，入力ファイルからの入力行は出力さ...
-[[PerlスクリプトからRを起動する]]などがヒントになるでし...
-いろいろありがとうございます。
どうやら、#!/usr/bin/Rというのはできないみたいですね。
perlやrubyを使う身としては、実行のたびにいちいち'R --vani...
このあたりは何とかして欲しいものですが…
結局、今のところはshを使うのがよさそうです。
#!/usr/bin/sh
R --vanilla --slave <<EOF
x <- rnorm(50)
y <- rnorm(x)
plot(x,y)
EOF
-面倒ならば，
alias rb 'R --vanila --slave ?!*'~
しておいて，
rb < ファイル名
で動くようにしておいたらいかが。~
sh のなかでいちいち
#!/usr/bin/sh
R --vanilla --slave <<EOF
なんたら
EOF
なんて書く方が面倒だ。 -- &new{2004-11-05 (金) 11:12:45};
-ああ、こういう人もいるんじゃないかと思ってはいましたが。...
-どうぞ。おすきなように。Cプログラムを記述してコンパイル...
内容が
R --vanilla --slave < $1
な実行ファイルを/usr/bin/Rsとでもしておけば，
#!/usr/bin/Rs
で動く． -- &new{2004-12-01 (水) 11:15:44};
- 今は Rscript があります。#!/usr/bin/Rscript -- &new{20...
**関数の中身を見る関数 [#e8ccf8da]
>[[akosuke]] (2004-10-29 (金) 10:45:01)~
~
自分で作った関数で名前はわかるけど中身を忘れてしまった、...
ここでも調べてみたのですが調べ方が悪いのか見つかりません...
くだらない質問で申し訳ありませんがどうぞよろしくお願いし...
//
-コンソールにその関数の名前だけを入力してみる.... これは...
-おわかりになりましたか？ -- &new{2004-11-04 (木) 14:36:...
**AICの計算 [#t0d6c154]
>[[AIC]] (2004-10-26 (火) 21:50:37)~
~
The R Bookの第5章を読んでいる者です。p.109の上から2行目で...
また、AICのパラメータ数はこの場合npram<-2と定義しています...
場合、式（5.14）にあるようにω、α、βの3と定義するべきでは...
それから（長くてすみません）GARCHの派生モデルを"fSeries"...
gjr(data,order(p,q,γ))？？)と思ってやってみたのですができ...
ご存知でしたら教えていただけないでしょうか。~
//
-本来ならば「The R bookご意見」でいただくべきご質問ですが...
-まず1番目のご質問は以下のようにしてください。
for (p in 1:3){
for ( q in 1:3){
tmp <- garch(dnikkei, order=c(p,q))
cat("p=", p, "q=",q, ", AIC=", (tmp$n.likeli+p+q),"?...
}
}
warningなどが出て少し見づらいですが、とりあえず。-- [[矢...
-2番目のご質問に対しては、AICはモデルの当てはまりの良さに...
-3番目のご質問のfSeriesに関しては詳しくないのですが、apar...
-ありがとうございます。上のプログラムで実行した結果、1番...
-３番目の質問に関しては、fSeriesをインストールしたあと、...
-fSeriesを用いたGJR推定は以下のようになります。
The R book第5章のデータを使った例(GJR(1,1))を示します。
library(fSeries)
tmp <- aparchFit(x=dnikkei,
order = list(alpha.lags = 1, beta.lags = 1, delta...
opt = list(gamma = TRUE, delta = FALSE,
disparm = FALSE), doprint = FALSE)
tmp
ポイントは上でgamma=TRUEとなっている部分です。この
部分がGJR定義式でのγを示します。gamma=FALSEとすると
ただのGARCHになります。GJRで日経225を推定すると非対称性が...
分かります。aparchFitについてより詳しくは
help(aparchFit)
としてhelpの内容を参照してください。-- [[矢野]] &new{2004...
-懇切丁寧な説明、ありがとうございます。頑張ってみます。 -...
**0を埋める [#z99a4268]
> (2004-10-25 (月) 18:52:25)~
~
とってもFAQだと思うのですが…~
1 -> 001，42 -> 042 ~
のように有効桁数（この例だと3ケタ）を決めて，~
それに満たない場合は0で埋めるにはどうすればいいでしょうか...
//
-sprintfでこのサイトを検索してみて下さい。 -- [[谷村]] &n...
-ありがとうございました．formatC という関数がありました：
x <- c(1:20)
formatC(x,width=2,flag="0")
-関数の紹介ありがとうございます。
> x <- c(1:20)
> formatC(x,width=2,flag="*")
flag にいろいろな文字を指定するとおもしろい結果が出てきま...
本来は，無効な文字ははじかなくてはならないのでしょうが。~
ベクトル化されているという点では sprintf より使いでがある...
**変数名に別の変数の内容を含めたい [#wc8d2de9]
>[[質問者]] (2004-10-23 (土) 15:43:00)~
~
大変初歩的な質問と思いますが、うまくいかないのでよろしく...
month <- c("Jan","Feb","Mar","Apr","May","Jun",
"Jul","Aug","Sep","Oct","Nov","Dec")
for (i in 1:12) paste("Hoge.",month[i],".total",sep="") ...
として、例えばHoge.Jan.totalに0が格納されるようにしたいの...
//
-[[RプログラミングTips大全]] に定石構文が紹介されています...
> month <- c("Jan","Feb","Mar","Apr","May","Jun",
"Jul","Aug","Sep","Oct","Nov","Dec")
> for (i in 1:12) eval(parse(text=paste("Hoge.",month[i]...
> Hoge.Apr.total
[1] 0
ついでですが次のようにしても良いわけです。
> for (i in 1:12) eval(parse(text=paste("Hoge",month[i],...
また R は月名を表す組み込みの文字列ベクトルを最初から持っ...
>month.abb
[1] "Jan" "Feb" "Mar" "Apr" "May" "Jun" "Jul" "Aug" "Sep...
> month.name
[1] "January" "February" "March" "April" "May...
[7] "July" "August" "September" "October" "Nov...
-ありがとうございました。evalも使って試行錯誤もやってので...
**データフレームの分割処理について [#xddecfdc]
>[[R-Fan]] (2004-10-18 (月) 14:38:21)~
~
このサイトでRを知りました。いつも有益な情報をありがとうご...
表題についてですが、大きなデータをキー項目（下記例ではdat...
(案1)
for (key in sort(unique(maindata$date))){
subdata <- subset(maindata, date == key)
#subdataの処理
}
自己レスになりますが、subset()を使用するかわりにキー項目...
(案2)
idx <- tapply(1:nrow(maindata), maindata$date,
function(subidx){return (subidx)})
for (key in sort(unique(maindata$date))){
subdata <- maindata[idx[[toString(key)]],]
#subdataの処理
}
もっと根本的に改善する方法がありましたらご教授いただけれ...
どうぞよろしくお願いいたします。
//
-ケースにもよるでしょうけど,unique(maindata$date)をlevels...
-ご回答いただきましてありがとうございます。後者のご提案に...
> d <- data.frame(x = factor(c(1,1,1,2,2,3)), y = 10:15)
> d[levels(d$x) == 1,]
x y
1 1 10
4 2 13
-次のようになるのでは？
> d[levels(d$x)[1]==d[,1],]
x y
1 1 10
2 1 11
3 1 12
>
> d[levels(d$x)[2]==d[,1],]
x y
4 2 13
5 2 14
> d[levels(d$x)[3]==d[,1],]
x y
6 3 15
データファイルが分析のたびに付け加わっていくのでなければ...
そうそう，思い出した。by という関数，速いか遅いか知りませ...
- levels()関数の例は大変よく分かりました。ご紹介いただい...
(案3)
by(maindata, maindata$date,
function(subdata){
#subdataの処理
}
)
**関数の読み込み [#x20303d3]
>[[恐縮です]] (2004-10-17 (日) 00:24:05)~
~
Rを使い始めて間もない者です。~
Rで色々な関数をinstall.packages()でインストールするのです...
一度インストールした関数をそのまま使い続けるにはどう保存...
//
-library(パッケージ名） という使い方ができないということ...
-すみません。library(パッケージ名)を実行して使えるように...
-library は毎回，必要なときにやらないといけないでしょう。...
-OSを書かれていませんので、Linuxを前提に説明するとinstall...
-library()を毎回実行するのが面倒なら、青木先生の助言のよ...
-普通，ライブラリのインストール時には管理者権限でやるので...
-.First()関数に初期化処理を定義しておく方法はだめでしょう...
-参考URL http://www.okada.jp.org/RWebRef/RManualJP/Startu...
-R 2.0 からは、たしかパッケージは実際に使う際になって初め...
**積集合(INTERSECT)と和集合(UNION)の中間的なものはありま...
>[[撃墜王]] (2004-10-14 (木) 15:39:55)~
~
R言語には、INTERSECTとUNIONといったコマンドが存在します。~
ただ、入力されるリストが非常に多く存在した場合に、「それ...
また、「それら全てのUNION」を取ると、当然ながらアウトプッ...
そこで、質問させて頂きたいのは、INTERSECTとUNIONの中間的...
~
なお、自分の中のイメージとしては~
~
入力されるリスト 100個~
1リストあたりに含まれる変数名(Aomori,Miyagi,Hokkaido...et...
全てのINTERSECTを取った場合 最終的に残る変数名 3個(少な...
全てのUNIONを取った場合、 最終的に残る変数名 1024個(多...
~
そこで、単純なUNIONではなくて、1つのリストにしか含まれな...
具体的な例をさらにあげますと~
~
list1 Kanagawa,Tokyo,Saitama~
list2 Chiba,Kanagawa,Tokyo~
list3 Tokyo,Irabaki,Gunma~
~
Intersect = Tokyo,Union=Kanagawa,Tokyo,Saitama,Chiba,Ibar...
となってしまいますが、上のようなリストの場合Kanagawaはlis...
//
-ある要素が，いくつのリストに含まれているかを勘定し，リス...
> top.n <- function(lst, n)
+ {
+ result <- table(unlist(lst))
+ result[result >= n]
+ }
> lst <- list(no1=c("aa", "bb", "cc"),
no2 =c("bb", "dd", "ff", "hh"),
no3=c("cc", "ff", "xx", "yy"),
no4=c("aa", "cc"))
> top.n(lst, 3)
cc
3
> top.n(lst, 2)
aa bb cc ff
2 2 3 2
探すのがめんどいから作っちゃった。関数にするまでもないが...
-青木先生有難うございます。先生の手に掛かるとあっというま...
**CSVファイルのインポートで errorが。 [#pa3d6bf6]
>[[竹内]] (2004-10-10 (日) 19:42:29)~
~
こんにちは。Ｒを使い始めてまだ5日目の初心者です。使用環境...
Error in file(file, "r") : unable to open connection
In addition: Warning message:
cannot open file `test.csv'
R-Bookや googleで一応は調べた（つもり）なのですが、 error...
//
-ファイルが開けないといっているのです。どこにそのファイル...
Tools --> Change Working Directory で，そのファイルがある...
-さっそくのご回答ありがとうございました。教えていただいた...
**frequencyについて [#w6a1889e]
>[[R初心者]] (2004-10-07 (木) 23:49:37)~
~
Rで株価の時系列分析をしようとしている者ですが、frequency...
//
-土日が無いなら freq = 5 とする他無いのでは。祭日は NA で...
> ts(1:10, frequency = 5)
Time Series:
Start = c(1, 1) # 始点は第一週の第１日と解釈
End = c(2, 5) # 終点は第二週の第５日と解釈
Frequency = 5
[1] 1 2 3 4 5 6 7 8 9 10
-ご質問は株価データが日次データ（土日祝日データなし）が何...
-御返事有難うございます。矢野先生のおっしゃる通り株価が何...
-えーっと、まず「先生」というのはやめてください。先生では...
-それでfrequencyを250にすると良い理由ですが、それは簡単で...
-どうもありがとうございます。納得しました。ただ何年か前ま...
**GARCHモデルのパラメータ推定 [#wd87da8d]
>[[ピエロ]] (2004-10-07 (木) 22:18:58)~
~
RでGARCHモデル分析する際、パラメータ推定は何の方法を用い...
擬似最尤法とかいうやつでしょうか？~
~
あと、Rでの時系列分析の細かな手法を扱っている本あるいはサ...
宜しくお願い致します。~
//
-ソースを御覧になるより仕方がないのでは。こうした質問にす...
-help(garch)とすると"maximum-likelihood estimates of the ...
-それと時系列解析の大まかな手法でよければ、矢野が書かせて...
~
時系列解析法:http://web.sfc.keio.ac.jp/~kogure/courses/20...
-それと上のコメントの方がおっしゃるとおり、Rの日本でのユ...
-それとRjpの以下のページは非常に便利ですばらしいです（作...
~
[[Rの基本パッケージ中の時系列オブジェクト一覧:http://www....
-↑を作成されているのは間瀬先生ですね。ご尽力ありがとうご...
-間瀬先生でしたか。改めてありがとうございます。 -- [[矢野...
-実は私の先生は時系列の専門家（藤井光昭先生）なのですが、...
-矢野先生の書かれた本の第5章、読ませて頂きました。5.7.2節...
-こちらも同じで、「先生」というのは勘弁してください -- [[...
-plot(predict(ng))というのはパラメータを決定したgarchモデ...
-御返事有難うございます。矢野さん。上記のplot(predict(ng)...
-ちょっとは，自分で調べましょう。 -- &new{2004-10-15 (金...
-私の勉強不足で初歩的な質問になってしまい、申し訳ありませ...
-縦軸が2003年1月から12月までになってしまうのは昔からある...
**pairsのプロット範囲 [#a4f0ce6f]
>[[ぷーすか]] (2004-10-07 (木) 20:17:03)~
~
散布図行列で、プロット点をドットに変え、y=xの直線を書きも...
~
pairs(x,panel=function(x,y){points(x,y,pch=".");abline(0,...
~
でうまくできました。~
ここでさらに各図のプロット範囲を同じにしようとxlim、ylim...
~
pairs(x,xlim=c(-3,3), ylim=c(-3,3), panel=function(x,y){p...
~
ではエラーが出てしまいます。~
プロット範囲以外のカスタマイズはあきらめ、~
~
pairs(x,xlim=c(-3,3), ylim=c(-3,3))~
~
とすれば、warningは出ますが、一応プロット範囲は制御できま...
~
プロット内容とプロット範囲のカスタマイズを両立させる方法...
//
-エラーメッセージはパネル関数は追加パラメータを受け付けな...
data(USJudgeRatings)
pairs(USJudgeRatings[1:3],
panel=function(x,y, ...) {points(x, y, pch="."); a...
xlim=c(5, 10), ylim=c(5, 10) )
#ref(pairspanel.png)
-ありがとうございます。panelに引数の追加を指定しつつも、x...
**eclipseの使い方 [#i5f83d3a]
> (2004-10-04 (月) 19:18:27)~
~
eclipse (eclipse-SDK-3.0-linux-gtk.zip)をダウンロードし、...
//
-自己フォローです。File->New->R Script fileでエディタが現...
**エクセルファイルのデータをRに取り込む方法は？ [#h1ace67d]
>[[初心者]] (2004-09-27 (月) 11:51:31)~
~
Rの初心者ですが、どうやってエクセルファイルのデータをRに...
//
-このページの下の方に，「CSVファイルのインポート」と「Rcm...
-出来たら初心者の為のデータインポートとか初心者が躓きやす...
-R の公式マニュアルの 「R Data Import/Export」に他システ...
-最近の [[なんでも掲示板]] の青木さんの関数も使えるかも知...
-みなさんありがとうございました。何とかなりそうです。でも...
-なぜか繰り返される質問と言うのは現状のドキュメントでは不...
-「よろしくお願いいしマース〜〜」などと甘まったれた文句を...
-「2万5000くらいのレコード」ということは，2万5000個の数値...
質問者の常だが，「できるのでしょうか」なんて聞くより，や...
初心者には優しくしましょうね(^_^)（発言タイトルはなおして...
-ごめんなさい！いろいろ試してみたんだけど、だめだったので...
-いつもこのようなメッセジーが出て、困っています。Warning ...
-初心者の一人ですが，なかまさんのアドバイスに従って，「エ...
-なるほど,MsysのシェルにPATH追加してRとか打ってるようでは...
----
~
**ヒストグラムＸ軸目盛について。 [#o01284d8]
>[[tofu]] (2004-09-21 (火) 11:11:26)~
~
度数分布をヒストグラムにしたときののＸ軸の目盛の数字、目...
//
-axis 関数で，任意の場所に任意の文字列を書くことができま...
-せっかく [[グラフィックス参考実例集]] があるのですから、...
-すみません…いろいろ勉強不足なもので。アドバイスありがと...
-実例集のページもずいぶん大きくなったのですね。見るのも大...
簡単に一例を示すと，以下のようになるか（縦軸は自分で考え...
x <- rnorm(2000)
hist(x, axes=F)
axis(1, at=c(-1.96,0,1.96), labels=c("-1.96 SD", "Zero",...
#ref(hist3.png)
-丁寧な解説ありがとうございました。Ｒとともにこのサイトの...
**plotでドットだけを打ちたい [#g180f43c]
>[[KabuTaro]] (2004-09-20 (月) 04:02:13)~
~
plot(c(1,2,3),c(4,5,6),type="p")~
とすると(1,4),(2,5),(3,6)に丸印でプロットされますが~
これを単にドットだけにするにはどうしたらよいでしょうか？~
pchを弄くっても形が変わるだけでうまくいきません。~
よろしくお願いいたします。~
//
-ドットだけの意味が分かりませんが，http://cran.r-project....
-pch="."ですね。できました。ありがとうございました。これ...
-plot(1,pch="A",cex=34) と plot(1,pch=".",cex=34) を比べ...
-勝手注： 中間さんのコメントの通り pch="." は真正のドット...
#ref(pchtest.png)
-なるほど、pch="."には特別な意味があるんですね。てっきり...
**package(gregmisc)のheatmap.2の引数設定について [#da53a4...
>[[ちょろべぇ]] (2004-09-14 (火) 18:57:31)~
~
gregmiscパッケージのheatmap.2を用いてヒートマップを作図し...
heatmap.2は様々なパラメータ設定が可能であり、自分としては...
(希望)~
?distfun→距離関数だと思われるので、ユークリッド距離を計算...
?hclustfun→average linkageで結合したい~
?できればdivisiveではなくagglomerativeでクラスタしたい~
~
ただ、
heatmap.2(data,distfun=dist(method="euclidian"),
hclustfun=hclust(dist(method="euclidian"),meth...
dendrogram="both")
のように設定しても、hclustfun部分でエラーが帰ってきてしま...
かなり愚かな質問だとは思いますが、アドバイスありましたら...
//
-こうした質問にはエラーメッセージをそえるのが礼儀作法です...
Hclust <- function(d) hclust(d, method = "average", memb...
heatmap.2(data, hclustfun=Hclust, dendrogram="both")
# もしくは一度に済ませるには
heatmap.2(data, hclustfun=function(d) hclust(d, method =...
dendrogram="both")
-こんな風なエラーメッセージが出てしまいますね --[[ちょろ...
Error in as.matrix(x) : Argument "x" is missing, with no...
-euclidian?euclidean? どちらが正しいのでしょうか？~
ちょっと違いますがこんな感じでどうでしょう-- [[○田五郎]]&...
my.dist <- function(x) dist(x, method="binary")
my.hclust <- function(d) hclust(d, method="ward")
hm <- heatmap(blablabla, distfun=my.dist, hclustfun=my.h...
-確認しないでコメントしました(最初のコメントのバグを修正...
library(gregmisc)
data(mtcars)
x <- as.matrix(mtcars)
rc <- rainbow(nrow(x), start=0, end=.3)
cc <- rainbow(ncol(x), start=0, end=.3)
my.hclust <- function(d) hclust(d, method = "average")
hv <- heatmap.2(x, col = cm.colors(256), scale="column",
hclustfun=my.hclust,
RowSideColors = rc, ColSideColors = cc, ...
xlab = "specification variables", ylab= ...
main = "heatmap(<Mtcars data>, ..., scal...
tracecol="green", density="density", den...
# 形式 hclustfun=function(d) hclust(d, method = "average...
#ref(hm2.png, left)
-諸先輩方、親切なアドバイスありがとうございました。ご教授...
**エラーバー [#a79042cd]
>[[Kei]] (2004-09-11 (土) 16:54:52)~
~
実験のデータをエラーバーで表現したいですが、R言語でどのふ...
mybar <- function(X) {
n <- nrow(X)
df <- n-1
X.mean <- apply(X,2,mean)
X.var <- apply(X,2,var)
d.error<-qt(0.975,df)*sqrt(X.var/n)
d.error.x <- (0.2 + 1) * 1:5 - 1/2 # エラーバーを表...
barplot(X.mean, names.arg=colnames(X),
ylim=c(0,ceiling(max(X.mean+(d.error)))),
col=c("blue", "green", "red", "yellow"))
arrows(d.error.x, X.mean, d.error.x, X.mean+ d.error, ...
}
mybar(分析するデータ[,列：列])
軸の調整方についてちょっと質問したいですが、教えてくださ...
//
-自分で作れるというのは素晴らしいことですが、こうした問題...
-縦軸をどのように調整したいのでしょうか？~
描画される範囲を，必要・十分に確保するというとことでしょ...
僭越ながら，[[私の書いた関数:http://aoki2.si.gunma-u.ac.j...
要点は，描画する前にある軸（今回は縦軸）において，描画対...
なお，必要なら，描画範囲を決めるだけで，軸も実際には描か...
plot(c(x軸の最小値, 同最大値), c(y軸の最小値, 同最大値),
type="n", # 描画しない
xaxt="n", # x 軸を描かない（必要なら）
yaxt="n", # y 軸を描かない（必要なら）
その他のパラメータ)
:
points(...)
lines(...)
axis(...)
その他諸々の描画関数
健康のために使いすぎに注意。 -- [[青木繁伸]] &new{2004-09...
**棒グラフの重ね合わせ [#l6008258]
>[[TETSU]] (2004-09-10 (金) 11:51:50)~
~
２つの棒グラフをbarplot関数を用いて描画したのですが、これ...
//
-barplot の beside という引数について，よく調べてください...
それでもよくわからなかったら，[[グラフィックス参考実例集...
このサイトで何かを調べようと思ったら，[[単語検索:http://w...
-ありがとうございます！ただ、besideをTrueにするとグラフは...
-「色々調べてみた」-> まず ?barplot でたくさんある引数を...
barplot(VADeaths, border = "dark blue", space=c(-0.5))
#ref(bp.png)
-青木先生、QDUさん、ありがとうございます。全ては私の日本...
-polygon 関数を使って，自分でヒストグラムを描く関数を書け...
**層別の相関値 [#ge316df8]
>[[HOJO]] (2004-09-10 (金) 11:39:45)~
~
欠損値を含む二つの変数について，相関値をcor(x,y,use=pai...
//
-上位30件というのはどのように決めたのでしょうか。~
それはともかく，データの散布図を描けば問題解決の手がかり...
下の図の右上30個ほどの点は相関係数0.910，次の70個ほどは0....
#ref(plot.png)
R はグラフを描くのがとっても簡単なので，数値解析する前に...
-青木先生いつも有り難うございます。散布図の特徴は，(1)原...
-＞一般に全体の相関係数は，層別した相関係数より大きくはな...
そんなことありませんよ。例えば，原点付近と(10,10)を中心と...
> gendat2 # 任意の相関係数を持つデータを生成する関数
function(nc, r)
{
z <- matrix(rnorm(2*nc), ncol=2)
res <- eigen(r2 <- cor(z))
coeff <- solve(r2) %*% (sqrt(matrix(res$values, 2, 2,
byrow=TRUE))*res$vectors)
z <- t((t(z)-apply(z, 2, mean))/sqrt(apply(z, 2, var)*(...
z %*% chol(matrix(c(1, r, r, 1), ncol=2))
}
> d1 <- gendat2(10,0)
> d2 <- gendat2(10,0)+10
> cor(rbind(d1,d2))
[,1] [,2]
[1,] 1.0000000 0.9615385 # プールすると相関は0.962
[2,] 0.9615385 1.0000000
> cor(d1)
[,1] [,2]
[1,] 1.000000e+00 -1.010639e-16 # プールする前は無相関
[2,] -1.010639e-16 1.000000e+00
> cor(d2)
[,1] [,2]
[1,] 1.00000e+00 3.27047e-16 # プールする前は無相関
[2,] 3.27047e-16 1.00000e+00
極端な場合，層別では負の相関，全体では正の相関ということ...
＞（10000,600)付近に他と離れたデータが１件~
というのは，すさまじいような気がしますが。その1つのデータ...
ともかく，HOJOさんが今扱っているデータではそのように計算...
散布図をここにペーストしてくれると，後学のためにもなるか...
-＞そんなことありませんよ。例えば，原点付近と(10,10)を --...
-＞そんなことありませんよ。例えば，原点付近と(10,10)を．...
-#ref(D:?My Documents?40joho?scatterplot1.bmp) -- &new{...
-この wiki に画像を張り付けるには、jpg か png 形式で、ペ...
-すこし極端な例を。全体では相関 0, 層別すると相関が 1, -1...
#ref(cor.png)
-#ref(scplot1.jpeg) -- [[HOJO]] &new{2004-09-11 (土) 20:4...
#ref(scplot1.jpeg)
-「編修」を使ってようやくできました。 アドバイスありがと...
-この図から見るにやはり「外れ値」とみなすべき例外的データ...
-青木先生，ＱＤＵさん他の皆様方アドバイス有り難うございま...
**単純主効果 [#i4379192]
> (2004-09-10 (金) 02:07:41)~
~
2元配置の分散分析（完全無作為）をする場合に~
交互作用がある場合，~
１．「単純主効果の検定」，~
２．「単純主効果の多重比較」~
をしたいのですが，Rではどのよなプログラムを書けば良いので...
~
参考までに，
SASでの単純主効果についての情報と，そのデータを以下にリン...
例：http://koko15.hus.osaka-u.ac.jp/~kano/lecture/faq/int...
データ：http://mat.isc.chubu.ac.jp/R/tech.html#CRFpq~
ご存知の方，教えて頂ければ幸いです~
~
//
-単純主効果という言葉の意味が一寸わからないのですが、R の...
-aov関数では得られないようです．単純主効果の説明は，http:...
-CRAN の検索では simple main effect + anova では一件も引...
**Rcmdrでグラフの変数選択欄 [#wb1d1e2a]
>[[HOJO]] (2004-08-31 (火) 16:06:53)~
~
たびたびですみません．RcmdrでIndexPlotなどを用いてグラフ...
//
-変数選択欄に出るのは数値として認識出来た列のみになるかと...
-なかまさん，何度もお手を煩わして申し訳ありません．因子か...
-変数毎にas.vectorを行うと，文字データになってしまいます...
-こうした質問はもとのデータ（長過ぎるなら問題になりそうな...
-どうしても意図した通によみこみたいならscanを使ってくださ...
-数値そのものに，あまり意味があるとは思えないので掲載しな...
-scanを使用して解決できました。有り難うございました。 -- ...
-もとのデータがただの数値で、単純に因子化されているならば...
-詳細に記載すると，以下の処理をしてからscanを施し解決して...
-詳細に記載すると，csvファイルを作成する際に事前に以下の...
-(1)は，カンマ区切りの数値は引用符で囲まれるから（そうし...
(2)は文字として % が実際にファイルに書き出されているため...
(3)空セルは文字通り，なんにも出力されないので，ケースあた...
全て，Excel に罪はないとも言える（制御のためのオプション...
いちいち手作業で(1),(2),(3)の前処理を行うのも面倒なので（...
-解説有り難うございました．相変わらず，やみなべ（Ｒ）をつ...
**CSVファイルのインポート [#y06b905e]
>[[HOJO]] (2004-08-31 (火) 11:26:52)~
~
エクセルで作成した２００行，２０列程度の数値データのみの...
１． Error: unable to open connection~
２． data is not a data frame and cannot be attached~
これはどのような点に気をつければよいのでしょうか。~
//
-多分ですが,日本語化してないRで日本語ファイル名のファイル...
-なかまさん，いつも有り難うございます．ご指摘のとおりでし...
**非心パラメータ [#xe80473c]
> (2004-08-24 (火) 19:12:06)~
~
正確な用語がわからないのですが~
SASではTNONCT関数で算出できる，~
「t分布の非心パラメータ（引数は，x, df, prob）を求める関...
Rにないでしょうか？~
http://support.sas.com/91doc/getDoc/lrdict.hlp/a000245960...
~
例えば，~
x <- 3.1; df<- 58; prob <- .975;~
?(x,df,prob)~
> 1.04844~
と求めたいのです．~
//
-そうした関数はないですが、要するに pt(3.1, 58, ncp) = 0....
> TNONCT <- function(x, df, prob, interval = c(-30, 30)) {
temp <- function(ncp) pt(x, df, ncp...
uniroot(temp, interval)[[1]] }
> TNONCT(3.1, 58, 0.975)
[1] 1.048444
-同じようなものになりましたが，参考までに。 -- [[青木繁伸...
> TNONCT <- function(x, df, prob, interval=c(0, 37.62))
+ {
+ uniroot(function(ncp) pt(x, df, ncp)-prob, interva...
+ }
> TNONCT(3.1, 58, 0.975)
[1] 1.048444
-非心分布など普段意識して使ったことがないので、いまひとつ...
-ありがとうございます．おかげさまで，目的の計算ができまし...
-私も，非心度を求めるようなことはやったことがありませんで...
> TNONCT(3.1, 58, 0.975)
[1] 1.048444
> TNONCT(-3.1, 58, 0.025)
[1] -1.048444
という関係があるようです。interval=c(-37.62, 37.62) とし...
質問者が示した URL での解説から，SAS でも同じような関数を...
google で site:r-project.org ncp pt で，もろ回答がありま...
[[R Help 2001:http://www.r-project.org/nocvs/mail/r-help/...
これでは，c(-1e7,1e7) としていますが，QDU さんの言うよう...
**行列のグラフィカル表示image()の色に任意の値を割り当てる...
>[[川村健介]] (2004-08-18 (水) 13:14:01)~
~
お世話になります。~
行列をグラフィカルに表示するimage()関数の色調を任意に割り...
例えば配列の値が0〜0.88の間にあった場合，0.0〜0.5に白を割...
//
- 色見本付きのイメージ図は example(filled.contour) 参照。...
-ありがとうございます。col=grey(c(1,1,1,1,1,0.4,0.3,0.2,0...
-数値と濃淡が正しく対応しているか慎重に確認して下さい。な...
**plot()で軸の数値をなくしたい [#uba2f2b7]
> (2004-08-18 (水) 00:22:26)~
~
10ほど並べた散布図で、最初の１つだけ軸の数値を表示し、残...
//
par(mfrow=c(2,5))
plot(sin)
plot(sin, xaxt = "n", yaxt = "n", xlab="", ylab="")
plot(sin, xaxt = "n", yaxt = "n", xlab="", ylab="")
# あと 7 回繰り返し
上のような感じで如何でしょうか？ご参考になれば幸いです．-...
-ありがとうございます。目的の図を得ることができました。 -...
**Ｒで扱える統計手法 [#k9c2485f]
>[[石川誠]] (2004-08-16 (月) 01:26:14)~
~
お世話になります。~
Ｒで扱える統計手法を知りたいのですが、何処かに日本語で一...
//
-まさにこのサイトでしょう。[[stats(R 統計)パッケージ中の...
ほかにどのようなページがあるかは，[[ページの一覧:http://w...
-R 本体のオブジェクト数は千以上、アドオンパッケージは少な...
-Ｒ に含まれる統計手法全部を知りたいというわけでもないで...
- Rで扱える統計手法でしたら，まず見なければいけないのは [...
-え〜〜，本人から申し上げますが，あそこはまんまり最初っか...
-またまた〜、ご謙遜されてますねぇ〜。「R」＋「統計」＝「R...
**行列の値上位5つの要素の位置を知りたい [#p6a1c4f8]
>[[川村健介]] (2004-08-14 (土) 09:05:51)~
~
初めて投稿させていただきます。~
630*630の配列の値の中から，上位5位までの要素の位置(x,y)を...
which.maxだと最も高い値の位置，例えば39726と帰ってきます。~
できれば，上位から[326,52],[555,156]...といった感じで位置...
//
- 欠損値 NA がなければもう少し簡単になるのですが。(PS. も...
> y [,1] [,2] [,3]
[1,] NA 0.7122784 0.7761505
[2,] 0.47935820 0.4700606 0.4547836
[3,] 0.02593524 0.8849352 0.3914027
> yy = matrix(length(y)+1-rank(y, na.last=F), dim(y), by...
> y [,1] [,2] [,3]
[1,] 9 3 2
[2,] 4 5 6
[3,] 8 1 7
> lapply(1:5, function(i) which(yy==i, arr.ind=T))
[[1]]
row col
[1,] 3 2
[[2]]
row col
[1,] 1 3
[[3]]
row col
[1,] 1 2
[[4]]
row col
[1,] 2 1
[[5]]
row col
[1,] 2 2
-すばやいレスありがとうございます。うまくいきました。ちな...
-ここはやはりベクトルの n 番目に大きな値を求める関数を定...
# ベクトル(行列、配列も可)の n 番目に大きな値を求める関数
# タイ や NA があっても使える(はず？)
> nthmax
function(x, n) {
if ( !(n %in% 1:length(x)) )
stop( "n is not in 1:length(x)" )
if (n == 1) return(max(x, na.rm=TRUE))
else
for ( i in 1:(n-1) ) {
x <- x[x != max(x, na.rm=TRUE)]
if ( length(x) == 0 )
stop("no 'n'-th max element in x")
}
return(max(x, na.rm=TRUE))
}
> x = matrix(runif(9), 3,3)
> for (i in 1:length(x)) print( which(x==nthmax(x,i), a...
row col
[1,] 1 1
row col
[1,] 3 1
row col
(途中略)
row col
[1,] 1 2
> x[1,2] <- x[1,1] # タイを入れる
> for (i in 1:length(x)) print( which(x==nthmax(x,i), ar...
row col
[1,] 1 1 # <- 最大値が二つある
[2,] 1 2
row col
(途中略)
row col
[1,] 1 3
Error in nthmax(x, i) : no 'n'-th max element in x
> x[1,2] <- NA # 欠損値を入れる
> for (i in 1:length(x)) print( which(x==nthmax(x,i), ar...
row col
[1,] 1 1
row col
(途中略)
row col
[1,] 1 3
numeric(0) # 最後がエレガントでない(欠損値のせい)
Warning message:
no finite arguments to max; returning -Inf
-初等的なやりかたは見つけた最大値を -Inf に置き換えて繰り...
> x = matrix(runif(9), 3,3)
> which(x==max(x), arr.ind=T)
row col
[1,] 2 3
> x[which(x==max(x), arr.ind=T)] <- -Inf
> which(x==max(x), arr.ind=T)
row col
[1,] 1 3
> x[which(x==max(x), arr.ind=T)] <- -Inf
> which(x==max(x), arr.ind=T)
row col
[1,] 3 2
-----
> x[ print(which(x==max(x), arr.ind=T)) ] <- -Inf # 一行...
-わかりやすい解説ありがとうございます。岐阜大で講義がある...
-そんなに実行時間の短いやり方が必要なんですか。スローライ...
作りつけの速い関数を求めるのもいいが，それを問い合わせる...
作りつけの関数でできないことも，いろいろ付け加えられます...
こういう考え方を許容できない人は無視すればいいだけです(^_...
> x <- matrix(runif(630*630),630)
>
> max.at <- function (x, top=10)
+ {
+ result <- matrix(0, top, 3)
+ for (loop in 1:top) {
+ mx <- x[1,1]
+ for(i in 1:630) {
+ for (j in 1:630) {
+ if (x[i,j] > mx) {
+ mx <- x[i,j]
+ mxi <- i
+ mxj <- j
+ }
+ }
+ }
+ result[loop,] <- c(mxi, mxj, mx)
+ x[mxi, mxj] <- -Inf
+ }
+ result
+ }
>
> system.time(max.at(x))
[ 1] 78.30 4.31 89.11 0.00 0.00
> max.at(x)
[,1] [,2] [,3]
[1,] 466 95 0.9999976
[2,] 239 31 0.9999970
[3,] 205 620 0.9999952
[4,] 595 183 0.9999935
[5,] 500 609 0.9999923
[6,] 37 556 0.9999901
[7,] 292 3 0.9999889
[8,] 515 278 0.9999818
[9,] 539 237 0.9999802
[10,] 483 331 0.9999786
速くしたいなら，C ででもプログラムすればいいだけです。プ...
-こんなのはいかが？
max.at2 <- function (x, top=10)
{
n <- nrow(x)*ncol(x)
col <- col(x)
row <- row(x)
o <- order(x)
cbind(row[o], col[o], x[o])[n:(n-top+1),]
}
> system.time(max.at2(x))
[1] 2.52 0.18 2.80 0.00 0.00 ＃ すんばらしく，速いわ！
> max.at3(x)
[,1] [,2] [,3]
[1,] 2 3 NA # NA はこの位置に入る
[2,] 18 178 99.99992
[3,] 549 252 99.99988
[4,] 208 115 99.99891
[5,] 182 140 99.99866
[6,] 308 597 99.99823
[7,] 207 156 99.99820
[8,] 616 112 99.99811
[9,] 499 249 99.99801
[10,] 355 79 99.99794
本当は，以下のようにしたいのだが，decleasing=TRUE が効か...
max.at3 <- function (x, top=10)
{
n <- nrow(x)*ncol(x)
col <- col(x)
row <- row(x)
o <- order(x, decreasing = TRUE)
cbind(row[o], col[o], x[o])[1:top,]
}
いずれにせよ，
max.at4 <- function (x)
{
o <- order(x)
cbind(row(x)[o], col(x)[o], x[o])
}
とすれば，小さい順の結果は全部得られるので，必要なところ...
-青木先生，ありがとうございました。大変，勉強になりました...
**statsライブラリ内のheatmap関数について [#r429e061]
>[[○田五郎]] (2004-08-10 (火) 01:45:09)~
~
こんばんわ。○田です。~
この度はstasライブラリ内のheatmap関数について相談させてい...
statsライブラリをロードして、heatmapがかけることはわかっ...
?heatmap内のデモスクリプトを書くなどしてみましたが、~
heatmap関数の引数であるcol引数の種類が良くわかりません。...
自分としては正の値が赤系、マイナスの値を青系で表示したい...
デモのcol = cm.colors(256)はレインボー系、、、~
col = topo.colors(16)は緑と薄いこげ茶系と、自分の望んでい...
どうぞよろしくお願いしますm(_ _)m。~
//
-こんなのありますよー http://gdds.pharm.kyoto-u.ac.jp/le...
-アドオンパッケージ RColorBrewer には様々な配色のパレット...
rainbow(n, s = 1, v = 1, start = 0, end = max(1,n - 1)/n...
heat.colors(n)
terrain.colors(n)
topo.colors(n)
cm.colors(n)
と同じように使えます。一般的には個人で配色を決めるよりは...
## help(brewer.pal) より
mypalette<-brewer.pal(7,"Greens")
image(1:7,1,as.matrix(1:7),col=mypalette,xlab="Greens (s...
ylab="",xaxt="n",yaxt="n",bty="n")
-[[グラフィックス参考実例集：カラーパレット]] を見て下さ...
-書きたいことは，ちゃんと書きましょう。主張したいことがよ...
- 元記事投稿者が topo.colors では青い色がでないといってい...
data(volcano)
x <- 10*(1:nrow(volcano))
y <- 10*(1:ncol(volcano))
image(x, y, volcano, col = c(topo.colors(12)[1:4],
heat.colors(8)[c(1,3,5,8)]), axes = FALSE)
#ref(topoheat.png)
**mvtnorm内で帰ってくる値が小数点4桁以下がその都度異なる...
>[[ちょろちゃん]] (2004-08-08 (日) 15:34:16)~
~
お世話になっております。~
Rのmultcompライブラリのsimtest、method=dunnettを~
利用しております。~
同じデータセットを入力した際に、simtestのダネットの結果中...
T値は常に一定の値が帰ってくるのですが、P値の小数点4桁より...
計算するたびに異なっております。~
どこで値が異なってしまったか追跡したところ、mulcompから呼...
mvtnorm内でP値を計算するときにどうやら異なった値が帰って...
これは仕様なのでしょうか、それともmvtnormのバグなのでしょ...
もしご存知の方がいらっしゃいましたらアドバイスいただけれ...
どうぞよろしくお願いいたします。~
//
-どこかで乱数でも使っていない限り、仕様にせよバグにせよ、...
-mvtnorm って，multcomp から呼ばれるんですか。よくわかり...
R の関数には，uniroot だっけを呼び出して解を求めるなどと...
あなたが利用した関数の定義を参照すれば，毎回異なった解が...
P値などというものは，0.05や0.01と比べて大きいのか小さいの...
いろいろな統計量に関するP値も，近似値なので（統計ソフト間...
いずれにしろ，近似値だと言うことで。そんなに目くじらたて...
-自分がテストに使ったデータは青木先生のページに記載してあ...
-灯台もと暗しでしたか。mvt 関数が疑惑の張本人ですね。主要...
? pmvt を読みますと（めんどくさいので斜め読み。従って精度...
Randomized quasi-Monte Carlo methods are used for
the computations.
とか
Value:
error: estimated absolute error
とか読めますが。~
また，pmvt の引数には，abseps = 0.001, releps = 0 という...
-なるほど。青木先生ありがとうございますm(_ _)m。感謝感激...
-青木さんのおっしゃる通り、乱数が絡んでいるようですね。試...
**繰り返しのｔ検定の結果を表示させるには [#bc1ff9ad]
>[[はまちゃん]] (2004-08-03 (火) 18:12:58)~
~
お世話になっております。
a<-matrix(rnorm(96),c(16,6))
for (i in 2:6) t.test(a[,1],a[,i])
のように、繰り返しでｔ検定を行う場合に、結果が表示されま...
なにか方法があったら、教えてお願いできますか？~
ほかの検定でも同じことが起こるのでしょうか？宜しくお願い...
//
-オブジェクトを書き出すには print を使うのが基本であると...
a<-matrix(rnorm(96),c(16,6))
for (i in 2:6) print(t.test(a[,1],a[,i]))
これで，うまくいくでしょう。~
統計学的な観点から言うと，このようなt検定の繰り返しが正し...
- -- [[MKR]] &new{2004-08-03 (火) 20:28:05};
x = sapply(2:6, function(i) {a<-matrix(rnorm(96),c(16,6)...
t.test(a[,1],a[,i])} )
-sapply は実行結果の正確な出力が得られないのでお勧めでは...
結果は
[,1] [,2] ...
statistic -0.3054952 0.3751029
ページ名:
![[RjpWiki]](image/../Rlogo.png "[RjpWiki]")
