Ｑ＆Ａ(旧1)
をテンプレートにして作成
[
トップ
] [
新規
|
一覧
|
検索
|
最終更新
|
ヘルプ
]
開始行:
COLOR(blue){SIZE(20){R および RjpWiki に関する質問コーナ...
COLOR(red){注意}：新規記事用の入力欄は以下の目次の直後に...
もし新規に質問をされる場合は、効率的にコメントを得るため...
- 広い意味で R に関係する話題に限定!
- 学校の宿題は自分で考えるべきです
- 具体的な状況がわかる様に背景説明をけちらない
- [[CRAN の検索エンジン:http://finzi.psych.upenn.edu/sear...
- 自分には馴染みでも、他の人には馴染みの無い手法・概念も...
- もし問題が解決した(結局解決しない場合を含め)場合はその...
~
~
COLOR(red){SIZE(18){注意：このコーナーは重量オーバーで表...
~
----
#contents
//#article
~
~
**群散布図（Group Scatter）を作成するパッケージはあります...
> (2004-05-13 (木) 22:15:48)~
~
[[ここ:http://www.osaka-amt.or.jp/lecture/stat/index2.htm...
//
- パッケージ ade4 にそれらしき関数が多数用意されています...
> library(ade4)
> data(lascaux)
> scatter(dudi.acm(lascaux$ornem, sca = FALSE), csub = 3)
#ref(plot-00.png,center)
例えばこんな関数がありました．
> x <- c(0.5,0.2,-0.5,-0.2) ; y <- c(0.2,0.5,-0.2,-0.5)
> eti <- c("toto", "kjbk", "gdgiglgl", "sdfg")
> plot(x, y, xlim = c(-1,1), ylim = c(-1,1))
> scatterutil.eti.circ(x, y, eti, 2.5)
> abline(0, 1, lty = 2) ; abline(0, -1, lty = 2)
> x <- c(0.5,0.2,-0.5,-0.2) ; y <- c(0.2,0.5,-0.2,-0.5)
> eti <- c("toto", "kjbk", "gdgiglgl", "sdfg")
> plot(x, y, xlim = c(-1,1), ylim = c(-1,1))
> scatterutil.eti(x, y, eti, 1.5)
> plot(runif(10,-3,5), runif(10,-1,1), asp = 1)
> scatterutil.grid(2)
> abline(h = 0, v = 0, lwd = 3)
> x <- runif(10,0,1) ; y <- rnorm(10) ; z <- rep(1,10)
> plot(x,y) ; scatterutil.star(x, y, z, 0.5)
> plot(x,y) ; scatterutil.star(x, y, z, 1)
> x <- c(runif(10,0,0.5), runif(10,0.5,1))
> y <- runif(20)
> plot(x, y, asp = 1) # asp=1 is essential to have perpe...
> scatterutil.ellipse(x, y, rep(c(1,0), c(10,10)), cell ...
> scatterutil.ellipse(x, y, rep(c(0,1), c(10,10)), cell ...
> x <- c(runif(100,0,0.75), runif(100,0.25,1))
> y <- c(runif(100,0,0.75), runif(100,0.25,1))
> z <- factor(rep(c(1,2), c(100,100)))
> plot(x, y, pch = rep(c(1,20), c(100,100)))
> scatterutil.chull(x, y, z, opt = c(0.25,0.50,0.75,1))
> par(mfrow = c(1,1))
#ref(plot-01.png,center)
関数が多数あるので全部は調べ切れませんでした．すみません...
-いろいろ調べるのが面倒なので，作る方が楽。 -- [[青木繁伸...
# テストデータの作成
mean <- c(3,5,6,8,7)
sd <-c(0.3, 0.4, 0.5, 0.2, 0.4)
x <- rep(1:5, each=20) # 5群（20ケースずつ）のデータ，平...
y <- NULL
for (i in 1:5) y <- as.integer(100*c(y, rnorm(20, mean=m...
#accu 測定値（縦軸）の精度（階級幅）
#stp 横軸方向のずらし量
# 以上二つは，試行錯誤で決定(^_^;)
graph <- function(x, y, accu, stp)
{
y <- round(y/accu)*accu
x1 <- unique(x)
for (i in 1:length(x1)) {
freq <- table(y[x==x1[i]])
for (j in 1:length(freq)) {
if (freq[j] >= 2) {
for (k in 1:length(y)) {
if (abs(y[k]-as.numeric(names(freq)[j])) < 1e-...
x[k] <- x[k]+(freq[j] <- freq[j]-1)*stp
}
}
}
}
}
plot(x, y)
}
graph(x, y, 0.2, 0.05)
#ref(graph.png,center)
-質問した者です。ありがとうございました。青木先生の関数を...
-凄いっ！いつもながら鮮やかなソースコードですね．簡潔さと...
**巡回セールスマン問題 [#icaee477]
> (2004-05-06 (木) 17:15:50)~
~
巡回セールスマン問題用の関数がRのあるパッケージにあったよ...
//
-?optim です。なんでも掲示板で以前教えていただきました。 ...
-example(optim) を実行すると、組み込みデータセット eurodi...
**matplotの軸 [#g1ba8468]
>[[さち]] (2004-05-03 (月) 00:00:58)~
~
matplotで軸の正の方向を左向きにしたいと思い、~
~
matplot(c(30,-80),c(50,-50))~
~
のようにしたのですが、グラフにしてみると軸は正の方向が右...
matplotで軸の方向を変える方法を知っている方がいたらアドバ...
//
-x 軸の正の方向を左向きにする場合は matplot() に引数 xlim...
xlim で調節してください．この方法は他の高水準作図関数 (例...
でも使えます．
sines <- outer(1:20, 1:4, function(x, y) sin(x / 20 * pi...
matplot(sines, xlim=c(20,1), type="o", pch=1:4, col = ra...
#ref(matplot-00.png,center)
y 軸側も同様にして 引数 ylim で調節してください．-- [[舟...
sines <- outer(1:20, 1:4, function(x, y) sin(x / 20 * pi...
matplot(sines, pch = 1:4, xlim=c(20,1), ylim=c(1,-1),
type = "o", col = rainbow(ncol(sines)))
-xlimを求めたい場合
matplot(c(30,-80),c(50,-50))
matplot(c(30,-80),c(50,-50),xlim=par("xaxp")[2:1]
とか、2度描画させるのが一番楽なのでしょうか? -- &new{200...
-matplot() で描いた直後は座標系の値がグラフィックスパラメ...
par(xaxp=c(0,1,5))
par("xaxp") # xlim は [0 1]
[1] 0 1 5
plot(1:10) # xlim ( c(xaxp[1],xaxp[2]) ) がこの...
par("xaxp") # xlim は [2 10] に
[1] 2 10 4
par(xaxp=c(0,1,5)) # 自分好みの座標系にする場合はxaxpを...
plot(1:10, xlim=c(xaxp[1],xaxp[2]))
par("xaxp")
[1] 0 1 5
余談ですが，par() で設定して (永続的に) 自分好みの座標を...
のでしょうか？"usr" や "xaxp" は値がコロコロ変わってしま...
仕方なく xlim , ylim で一時的に変更しているのですが．どう...
せいか，"usr" や "xaxp" の適当な使い方が思いつきません．-...
# 原点を通る座標軸を描く
plot(rnorm(30), rnorm(30), axes = F, xlab = "", ylab = "")
axis(1, pos = 0, at = pretty(par("usr")[1:2]), adj = 0)
axis(2, pos = 0, at = pretty(par("usr")[3:4]), las = 2, ...
box()
**[質問]ベクトル [#vd5e3dc9]
>[[さち]] (2004-04-30 (金) 05:19:49)~
~
ある変数に値を代入していき、その変数に値をベクトルのように~
記憶させていくプログラムを作りたいです。例えば、~
x<-1~
x<-3~
x<-9~
x~
1 3 9~
のようにしたいです。このとき、~
x[3]~
のように長さを与えれればよいのですが、~
いま代入する変数の数がいくつになるのかわからない状態です。~
このような時はどうしたらよいのでしょうか？~
初歩的な質問ですが、アドバイスなどお願いします。~
//
- R のベクトルは任意に長さを延長できますから気にせずに代...
> x <- numeric(0) # x <- NULL でもOK （この時点では x ...
> x
numeric(0) # この時点では x は空の変数
> x[1] <- 1
> x[2] <- 2
> x[3] <- 3
> x
[1] 1 2 3
> x[5] <- 5 # 途中飛ばしてもOK
> x
[1] 1 2 3 NA 5 # 未定義の x[4] には自動的に NA 値が...
-append という関数を使うと吉-- [[青木繁伸]] &new{2004-04-...
> x <- NULL
> x <- append(x, 1)
> x <- append(x, 3)
> x
[1] 1 3
> x <- append(x, 2, after=1)
> x
[1] 1 2 3
-ありがとうございます。 -- [[さち]] &new{2004-04-30 (金) ...
-ありがとうございました。早速ためしてみます！ -- [[さち]]...
**princompのloadingsは固有ベクトルそのもの？ [#a1152f15]
>[[石原茂和]] (2004-04-27 (火) 07:42:09)~
~
（すみません，同内容を旧Rwikiの掲示板にまちがえてかいてし...
~
こんにちは．いま気がついたのですが，princompのloadingsは~
固有ベクトルそのものですね．~
helpでも~
loadings: the matrix of variable loadings (i.e., a matri...
contain the eigenvectors).
~
と，そのように書いてあります．~
しかし，主成分負荷量は，加工していないもとの変数での値と...
のはずですが，~
確かに固有ベクトルとたいがいの場合，同じような傾向になり...
固有ベクトルそのものを負荷量というのは誤解のもとでは？~
//
-example(Harman74.cor) を実行してみてください。おっしゃる...
-出ますけど，原質問者の質問の答えになっているのかな？~
それと，必ず回転された答えが出るんでしょうか。何回転にな...
二変数の分析してみると，教科書と違った解になるようですが。~
> x
[,1] [,2]
[1,] -0.24451455 0.5090637
[2,] 0.04563039 -0.5012150
[3,] 0.97775313 -1.1832550
[4,] 1.26989942 1.4814664
[5,] -0.14870401 -1.3809846
> loadings(princomp(x))
Loadings:
Comp.1 Comp.2
[1,] 0.273 0.962
[2,] 0.962 -0.273
Comp.1 Comp.2
SS loadings 1.0 1.0
Proportion Var 0.5 0.5
Cumulative Var 0.5 1.0
-- &new{2004-04-29 (木) 19:22:43};
-失礼、Harman74.cor の example は princomp ではなく、fact...
**数量化１類で消費されるメモリ量について。 [#m90801fb]
>[[koichi]] (2004-04-19 (月) 23:48:32)~
~
はじめまして。Rの初心者です。~
Rで数量化1類を行おうとしているのですが、メモリの問題で手...
具体的には、~
１列目(V1)を『0から1までの８桁の乱数』、２列目(V2)を『1か...
~
d0$V2 <- factor(d0$V2)
result <- summary(lm(V1 ~ V2,data = d0))
~
すると~
~
Error:cannot allocate vector of size 390625 Kb
Reached total allocation of 1022Mb
~
と出て、タスクマネージャーを見るとRは「応答なし」で止まっ...
~
このデータ量で390625 Kbもの行列を生成してしまうのでしょう...
~
メモリ関連の操作や、計算法の工夫を通してなんとか計算させ...
~
当方の環境はWindowsXP,Pentium4 2.8GHz,1GBRAMです。~
よろしくお願いいたします。~
//
-変数の長さが n ならば、回帰分析で登場するハット行列は nx...
-一つのやりかたは、必要な情報だけを公式通りに自前で計算す...
# 参考 10055007 は文字列、数値の項目数ですから、
# メモリー量はざっとその 8 倍程度になると考えるべきでし...
> n=100; x=rnorm(n); y=factor(floor(1000*runif(n))); z=l...
[1] 10902
> n=1000; x=rnorm(n); y=factor(floor(1000*runif(n))); z=...
[1] 661272
> n=10000; x=rnorm(n); y=factor(floor(1000*runif(n))); z...
[1] 10055007 # z に含まれる項目総数
~
-提示された条件に従ってデータを生成して数量化I類を（定義...
そもそも，このような設定の数量化I類が必要なのか適切なのか...
-解析対象のデータがわかりませんから、野次馬コメントですが...
-みなさん、早速のご回答ありがとうございます。また、つたな...
参考のために提示された式、とてもためになりました。Rに習熟...
提示された条件に従ってデータを生成して数量化I類を（定義に...
-数量化I類のプログラムは，[[ここ:http://aoki2.si.gunma-u....
for を使わないとかえって複雑怪奇になるところもあるので（...
実行したのは Mac G5 です。~
必要メモリー数ですが，原データが50000×2×8バイト，原データ...
cov 関数で分散共分散行列を計算するようにしたので，展開し...
プログラミング言語で書くときには，こんなにメモリを贅沢に...
**.Call/.Externalで使用するDLLの作成方法について [#p9d0c3...
>どうむ (2004-04-19 (月) 15:20:39)~
~
現在R（Windows版）から呼び出せるCで書いたDLLを作成中です...
■R1.8.1~
■コンパイラー：「Borland C++ Compiler 5.5」~
■OS：WindowsXP~
です。いままでRやSPlusでCで書いたルーティンを作ってきたの...
まず、簡単なCルーティンですが
extern "C" __declspec(dllexport) void R_TEST001(int *x, ...
__declspec(dllexport) void R_TEST001(int x[], int *ret)
{
*ret = x[0]+x[1]+x[2]+x[3];
}
整数型のベクトルを引数にとってその合計を「ret」で返してい...
DLLを作成するのにMakefileを作っておりコンパイルとリンクは...
bcc32 -u- -I"***" -O2 /c ***.cpp
ilink32 /w /m /Tpd /L"***?Lib"
"C:?borland?bcc55?Lib?C0D32.OBJ" ***.obj,
***.dll,
CW32.LIB IMPORT32.LIB
（mapファイルの指定が欠如していますが気になさらぬよう。...
「***.cpp」が該当のCルーティンです。~
これを呼び出すRのプログラムは
> dyn.load("***.dll")
> is.loaded("R_TEST001") #確認
[1] TRUE
と、なります。結果は
> ret <- .C("R_TEST001",as.integer(c(1,2,3,4)),ret=integ...
> ret$ret
[1] 10
です。~
ここまでは、良いのですがRの「.Call」関数を使用できるよう...
#include <R.h>
#include <Rinternals.h>
extern "C" __declspec(dllexport) void R_TEST001(int *x, ...
__declspec(dllexport) void R_TEST001(int x[], int *ret)
{
*ret = x[0]+x[1]+x[2]+x[3];
}
たった、これだけで上のコンパイルが通らなくなりました。~
が以下のように「-A」オプション（ANSI 準拠の予約語のみを使...
bcc32 -u- -I"***" -O2 -A /c ***.cpp
ilink32 /w /m /Tpd /L"***?Lib"
"C:?borland?bcc55?Lib?C0D32.OBJ" ***.obj, ・・・
今度は、「Writing R Extensions」の第４章にあるように関数...
#include <R.h>
#include <Rinternals.h>
extern "C" __declspec(dllexport) void R_TEST001(int *x, ...
extern "C" __declspec(dllexport) SEXP R_TEST003(SEXP x, ...
__declspec(dllexport) void R_TEST001(int x[], int *ret)
{
*ret = x[0]+x[1]+x[2]+x[3];
}
__declspec(dllexport) SEXP R_TEST003(SEXP x, SEXP ret)
{
int i, j, nx, ny;
SEXP ans;
nx = length(x);
PROTECT(ans = allocMatrix(REALSXP, nx, 1));
UNPROTECT(1);
return(ans);
}
このようにすると、コンパイルは通るのですがリンクができま...
ilink32 /w /m /Tpd /L"***?Lib" "C:?borland?bcc55?Lib?C0D...
***.dll, CW32.LIB IMPORT32.LIB
Turbo Incremental Link 5.00 Copyright (c) 1997, 2000 Bor...
Error: 外部シンボル 'Rf_length(SEXPREC *)' が未解決（***...
Error: 外部シンボル 'Rf_allocMatrix(unsigned int, int, i...
Error: 外部シンボル 'Rf_protect(SEXPREC *)' が未解決（**...
Error: 外部シンボル 'Rf_unprotect(int)' が未解決（***.OB...
** error 2 ** deleting ***.dll
R.DLLから「implib.exe」でR.libを作成してそれを指定しても...
このようなエラーが出るということは、致命的な間違いを、犯...
解決につながるような情報がございましたら、よろしくお願い...
//
-関数が「Rprintf("Hello!?n");」のみだったら通りますね。 -...
> ret <- .Call("R_TEST003",as.double(c(1,2,3,4)),ret=dou...
Hello!
> ret
[1] 1 2 3 4
ということは、
length
allocMatrix
PROTECT
UNPROTECT
が、いけないことになりますね。適切なヘッダーを指定してい...
なにか、大きな過ちを犯しているような。。。~
何を見落としているのでしょう？？？
//
-自己レスばかりですがなんとなくわかりました。 -- どうむ &...
~Rprintfが定義してある「Print.h」には
#ifdef __cplusplus
extern "C" {
#endif
とあって、C++でコンパイルされることを想定しているようです...
ですので、C++ではなくCだと思わせてコンパイルをかけること...
「Rinternals.h」に「extern "C"」を書けばC++でもコンパイル...
同一ファイルでC++のクラスが使えないのは不便かなぁ。
-こうした使い方ができることは噂には聞いていましたが、実際...
-これにあたっては[[Rから他言語利用]]を、参考にさせていた...
**Mac OS X の作図デバイス [#c2c090f8]
>[[舟尾]] (2004-04-18 (日) 15:49:17)~
~
R を起動しますと，UNIX では x11 が，Windows では windows ...
//
-Mac OS X では、x11 が起動することも、Aqua (Mac OS X のグ...
-回答ありがとうございます．なるほど，Aqua を起動させるこ...
-Linux系の場合、Rの起動コマンドオプションに --gui=X11(デ...
-回答ありがとうございます．勉強になりました．基本はどのOS...
-バイナリのレベルでは、RAqua.dmg (1.9.0ではR.dmg)からイン...
**GLMM()でできるオブジェクトの参照方法 [#qe73ec2c]
>林 (2004-04-14 (水) 13:38:20)~
~
lme4パッケージのGLMM()でできるオブジェクトは~
is.list()で調べるとリストなのですが、names()では名前はな...
list()
となります。通常のリストでは
List of 5
などといくつ要素があるとか返ってくるので、GLMM()でできる...
//
**1.9.0がmakeできない [#j38a01fc]
>[[谷村]] (2004-04-13 (火) 00:42:16)~
~
VineSeedでR-1.9.0がmakeできません。~
making dataentry.d from dataentry.c
making devX11.d from devX11.c
making rotated.d from rotated.c
making rbitmap.d from rbitmap.c
make[4]: ディレクトリ `/home/hoge/tmp/R-1.9.0/src/module...
make[4]: ディレクトリ `/home/hoge/tmp/R-1.9.0/src/module...
gcc -I. -I../../../src/include -I../../../src/include -...
-I/usr/local/include -DHAVE_CONFIG_H -D__NO_MATH_INLI...
-g -O2 -c dataentry.c -o dataentry.lo
In file included from dataentry.c:31:
/usr/X11R6/include/X11/Xlib.h:1400: error: 文法エラー be...
/usr/X11R6/include/X11/Xlib.h:1488: error: 文法エラー be...
/usr/X11R6/include/X11/Xlib.h:1516: error: 文法エラー be...
/usr/X11R6/include/X11/Xlib.h:1520: error: 文法エラー be...
/usr/X11R6/include/X11/Xlib.h:1542: error: 文法エラー be...
[snip]
誰かVineSeedでmakeに成功している方はいらっしゃいますか？~
XはXOrg-6.7.0-0vl2です。
//
-コンパイルが通りました。dataentry.cの#include <X11/X.h> ...
**RjpWiki 頁のアーカイブのうまい作り方？ [#r1012a9e]
>[[間瀬茂]] (2004-04-04 (日) 23:06:46)~
~
例えばこのＱ＆Ａコーナーがかなり長くなってきたので、アー...
//
**SJava [#u4f6bbf5]
>[[高階知巳]] (2004-03-31 (水) 10:56:07)~
~
S/RとJavaのインタフェース SJava を、1年半ほど前に試したと...
//
-SJavaうまくいきました。omegahatは今は大丈夫なようで、し...
import org.omegahat.R.Java.*;
public class REvalSample {
public static void main(String [] args) {
String [] rargs = {"--slave", "--vanilla"};
System.out.println("Sample program to call R engi...
ROmegahatInterpreter interp
= new ROmegahatInterpreter(ROmegahatInterpret...
false);
REvaluator e = new REvaluator();
Object val = e.eval("x <- sin(seq(0, 2*pi, length...
val = e.eval("x * 2.0");
if (val != null) {
double[] objects = (double[])val;
for (int i=0; i<objects.length; i++) {
System.err.println("("+i+") " + objects[i...
}
}
}
}
-R-1.9.0+SJava-0.67.3でコンパイルは成功しました．が，実行...
**pairs(data.frm, upper.panel = points, lower.panel = 残...
>[[みやむら]] (2004-03-26 (金) 13:42:27)~
~
pairs() の描画のようにして，上三角部分に通常の散布図を，...
//
-[[pairs用 回帰直線つきpanel]]を参考に対応する lower.pane...
-単回帰の残差であればよいのですが，この場合重回帰ですので...
**help.startができなくなりました。 [#mb2d87a7]
> (2004-03-26 (金) 13:15:11)~
~
help.startでブラウザが起動しなくなりました。どこをチェッ...
ご存知の方は教えてもらえないでしょうか？~
こんなコメントが出ます。Making links in per-session dir ....
使用ソフトRAqua ver1.8です。~
//
-Macintosh HD/private/tmp/Rtmp988/.R/doc/html/index.html ...
-GUI で見えないのでどうでもいいのですが，Rtmp988 みたいな...
-私の場合だと，その index.html ファイルを jedit で開くよ...
-コメントどうもです。アドバイスどおりやってみましたが、、...
-open index.htmlだとちゃんとひらくのですが、、この前OpenO...
**2変量のヒストグラム [#y279fd4d]
>[[sak]] (2004-03-24 (水) 23:53:51)~
~
２変量のヒストグラムを描く事は可能でしょうか？~
//
-何種類かパッケージがありますが，その一例は[[ここ:http://...
-うーん、なんでもできそうですね。すごい。これからRの世界...
**[質問]DOSの |more や /p にあたるコマンドは？ [#cfdd0...
>[[inagaki00]] (2004-03-23 (火) 22:27:03)~
~
長ーいデータを見るとき、DOSでは~
>type hogehoge.txt | more とか~
>type hogehoge.txt /p のように打ちます（UNIXだとcatですね...
これと同じスイッチがRにはありますでしょうか。~
~
ありそうで見つからないので困っています。~
ご存知の方、どうぞよろしくお願い致します。~
//
-普段は x[1:100] のように一部だけ表示しますのでおっしゃる...
page(base) Invoke a Pager on an R o...
ObjecttaskCallbackManager(base) Create an R-level task c...
managergrid.newpage(grid) Move to a New Page on a ...
Devicetkpager(tcltk) Page file using Tk text ...
-なるほど！Pagerというんでしたか。 早速のご回答有り難う...
**ページのヘッダやフッタ描画、またはページの４隅のユーザ...
>[[Я・・R]] (2004-03-23 (火) 18:38:15)~
~
レイアウト、グラフ、余白設定などに関係なく､ページの特定位...
ヘッダやフッタを text(ユーザ座標) で描こうと思ったのです...
ところが~
~
box(which="outer", lwd=4, lty=2, col="blue")~
~
とやると必ず正しく４隅を囲むので、どうやって４隅の座標を...
~
?box~
~
としましたがソースコードが見れませんでした。~
ヘッダやフッタ描画、またはページの４隅のユーザ座標を得る...
//
-identify、locator 関数でインタラクティブに指定位置の座標...
-はい、ページ数が非常に多いので、マウス入力はちょっと無理...
-フリーソフトですから、ソースは(その気にさえなれば)見るの...
-一応解決しました。１ページ分を作図した後に par(new=T, ce...
**ESSの中の日本語 [#o2747e03]
> (2004-03-23 (火) 11:33:40)~
~
gnome-terminal上では~
> a <- "日本語"
> a
[1] "日本語"
と日本語が使えますが、ESS上では~
#ref(ess-20040323.png)~
という感じになります。gnome-ternimal上では問題がないのでE...
何とかなればと思っています。ご助言などがございましたらよ...
[環境]~
VineLinux2.6~
ess-emacs-5.2.0beta1-1vl1~
R-1.8.1-4vl1(本家版)~
//
-[[なかま様のご助言:http://www.okada.jp.org/RWiki/index.p...
**ファイル書き出し時の文字コードの指定方法 [#n01ba267]
>[[matsu]] (2004-03-18 (木) 19:06:00)~
~
はじめまして、matsuと申します。~
~
早速ですが、write()やwrite.table()でオブジェクトをファイル~
に書き出す場合、ファイルの文字コードを制御することはでき...
でしょうか。~
~
私の環境( R ver1.6.1 Win2k )では sjis(改行コードCRLF)で~
書き出されるのですが R側の設定でeuc(改行コードLF)に変更~
はできるものなのでしょうか。~
~
ポインタ等でもお教えいただければ幸いです。~
よろしくお願い致します。~
//
-オリジナルの R は当然日本語コードなど知るはずもないから...
-中からperl（WindowsならActivePerl）を呼び出すことにしてj...
-中澤コメントに補足。例えば write 関数には file="|cmd" の...
-みなさまありがとうございます。ソースをEUCで記述したり、w...
-以前、save(file="hoge.txt",ascii=T,list=hoge)ですと、asc...
**グラフィックウインドウの名前 [#w3dd2d3f]
> (2004-03-18 (木) 11:20:43)~
~
windows()などで表示されたウインドウ~
に自動的にふられる名称~
「R Graphics: Device x」~
を変更することは可能でしょうか？~
//
-Rソースコードのsrc/modules/X11/devX11.c の中で、Xライブ...
-そうですか。作成グラフ数が多いときにできれば便利です。上...
**非線形回帰の初期値 [#c9768ff8]
>[[ぶんしょ]] (2004-03-17 (水) 17:40:55)~
~
非線形回帰の初期値がすぐ入らないのですが、いい方法があれ...
今やっているのが~
tpot<-nls(HAV~(a*SI+b)*(1-exp(k*AGE))^m,start=list(a=1,b...
という式です。
見当をつけて初期値を入れれば大体はすぐ当てはまるのですが...
Error in numericDeriv(form[[3]], names(ind), env) :
Missing value or an Infinity produced when evaluating th...
や
Error in nls(HAV ~ (a * SI + b) * (1 - exp(k * AGE))^m, ...
singular gradient
や
Error in nls(HAV ~ (a * SI + b) * (1 - exp(k * AGE))^m, ...
step factor 0.000488281 reduced below `minFactor' of 0.0...
といったエラー文が出ます。これらの意味もできれば教えてく...
//
-最適化の初期値の選び方についてはチップス [[関数の最大・...
-四つの変数を適当な範囲と、刻目で変えながら、非線形誤差自...
-ありがとうございます。やはり神頼み的な思考錯誤しかないよ...
-誤差の入り方を無視し log(HAV) ~ log(a*SI+b)+m*log(1-exp(...
-すみません勉強不足なんですが、これは対数で回帰させるので...
-そうです。単なる思いつきです。目的変数を log(HAV) にとり...
-何度もすみません。対数で試したら~
Error in numericDeriv(form[[3]], names(ind), env) :
Missing value or an Infinity produced when evalua...
In addition: Warning messages:
1: NaNs produced in: log(x)
2: NaNs produced in: log(x)
といったエラー文が出るのですが何がいけないのでしょうか？ ...
-これは当然 1-exp(k*AGE)) が試行中に非正値になったという...
-もし差し支えなければ、そしてデータがあまり大量でなければ...
-
AGE HAV SI
18 197.98 21.48
20 266.02 23.78
24 695.26 39.88
26 429.21 25.14
27 312.04 22.79
33 484.55 24.15
35 708.11 30.79
25 743.29 35.23
28 714.40 32.50
29 377.93 19.36
29 684.48 29.50
42 557.06 20.37
47 483.70 16.62
50 725.47 20.88
24 400.58 33.11
25 458.79 32.69
25 296.57 24.12
27 382.04 27.50
32 343.38 21.45
33 450.72 26.33
33 375.25 22.31
34 420.39 21.82
36 384.44 20.51
が当てはまらない一例なんですが、初期値a,ｍ=正、ｂ,ｋ=負の...
-- [[ぶんしょ]] &new{2004-03-19 (金) 09:06:58};
-少し解析をしてみましたが、やはり最初に注釈があったように...
> optim(c(0.1, 0.1, -0.1, 0.1), SF)
$par
[1] 10.907251 -1.975936 -0.057869 1.013191
$value
[1] 403297
$counts
function gradient
313 NA
$convergence
[1] 0
$message
NULL
-わざわざありがとうございます。確かに粗悪なデータです。自...
**ｌｍ（）関数の不思議 [#y49c88e7]
>[[杉本]] (2004-03-17 (水) 16:31:34)~
x<-c(8,3,4,2,1,4,8,4,9)
y<-c(4,1,5,1,1,2,5,2,3)
として x=a*y+bの形で線形回帰したいときは~
lm(x~y)
ですが x=a*y^2+bの形で線形回帰するときはどうすればよいの...
lm(x~y^2)
とすると上と同じ結果になってしまいます。もちろん
x=a*y+b*y^2+cの形で線形回帰するときも~
lm(x~y+y^2)
Call:
lm(formula = x ~ y + y^2)
Coefficients:
(Intercept) y
1.586 1.197
となり、x=a*y+bの形で線形回帰してしまっています。新たにデ...
//
-[[Rの統計解析関数Tips]] を見て下さい。R のモデル公式では...
-ありがとうございます。うまくいきました。 -- [[杉本]] &ne...
**無題 [#r9128296]
>[[大熊 大郎]] (2004-03-12 (金) 23:58:36)~
~
関数 width.SJ はなんですか？おしえてください。~
//
-悪い質問の典型ですね。標準パッケージ以外のパッケージ中の...
**ヒストグラムを描くhist() [#eb449b55]
>[[Rin]] (2004-03-04 (木) 11:57:28)~
~
hist()を使いヒストグラムを描くとき、matrixデータを与える...
よろしくおねがいします。
//
-意味がいまひとつ理解できません。具体例をあげていただけま...
-上の方が仰るとおり，hist()では多次元データはベクトルに強...
> x <- matrix(1:24,6,4)
> barplot(x)
> barplot(x,col=rainbow(24))
> barplot(x,col=rainbow(12))
> barplot(x,col=rainbow(6))
> barplot(x,col=rainbow(3))
-理解しにくい質問の仕方をして申し訳ありません。凡例はどの...
-[[ここ:http://www.okada.jp.org/RWiki/?cmd=read&page=%A5%...
-はい。そうです。早速、legend()を使い試してみましたが、や...
- histogram につける気の利いた凡例が思いつかなかったのが...
> data(islands)
> hist(islands)
> legend( 5000,y=30,paste("data",c(1:3)), col = c(1:3), ...
> legend(10000,y=30,paste("data",c(1:3)), col = c(1:3), ...
> legend(10000,y=40,paste("data",c(1:3)), col = c(1:3), ...
#ref(hist.png,center)
//-凡例をつけることができました！やっと解決できて嬉しいで...
//&new{2004-03-09 (火) 17:26:53};
-凡例をつけることができました！やっと解決できて嬉しいです...
-いえいえ．?legend や help(legend) でさらなる情報が得られ...
-初めまして。matlab初心者の者です。棒グラフで値を表示した...
-続き）matlabでは棒１本１本の色をかえることができずその方...
**帯グラフ [#mc629b65]
> (2004-03-02 (火) 12:57:28)~
~
グラフィックス参考実例集には、なぜか帯グラフがありません...
//
-リンク集、グラフィックス参考実例集にもリンクが張られてい...
> data(HairEyeColor)
> x <- apply(HairEyeColor, 2, sum)
> par(fin=c(6,2)) # グラフィックス画面を 6x4 インチに指定
> barplot(as.matrix(x), horiz=TRUE, col=rainbow(length(x...
#ref(obigraph.jpg,left)
-ありがとうございます。私のイメージとはちょっと異なるので...
> data(HairEyeColor)
> a <- as.table( apply(HairEyeColor, c(1,2), sum) )
> b <- a / apply(a, 1, sum)
> barplot(t(b))
#ref(obi.png,left)
> d.seg.xs <- rep((0.2 + 1) * 1:(ncol(t(b)) - 1), each=n...
> d.seg.xe <- rep((0.2 + 1) * 2:ncol(t(b)) - 1, each=nro...
> d.seg.ys <- apply(t(b)[,1:(ncol(t(b)) - 1)], 2, cumsum)
> d.seg.ye <- apply(t(b)[,2:ncol(t(b))], 2, cumsum)
> segments(d.seg.xs, d.seg.ys, d.seg.xe, d.seg.ye)
#ref(obi2.png,left)
-おぉ、horiz=TRUEなんてあるのですか。早速、使ってみました...
> barplot(t(b),horiz=TRUE)
> d.seg.ys <- rep((0.2 + 1) * 1:(ncol(t(b)) - 1), each=n...
> d.seg.ye <- rep((0.2 + 1) * 2:ncol(t(b)) - 1, each=nro...
> d.seg.xs <- apply(t(b)[,1:(ncol(t(b)) - 1)], 2, cumsum)
> d.seg.xe <- apply(t(b)[,2:ncol(t(b))], 2, cumsum)
> segments(d.seg.xs, d.seg.ys, d.seg.xe, d.seg.ye)
#ref(obi3.png,left)
-あとは、0.0-1.0ではなくてパーセント表示にしたいのですが...
-% 表示はやはり手作業技でしょうね。ところで元記事質問者の...
-名乗らなくてすみません。元記事質問者です。上の3つの帯グ...
-barplot()の中でaxes=Fとしておいて，barplot()後にaxis(1,s...
-ありがとうございます！なるほど。%が欲しかったので、axis(...
-質問者の分際でなんなのですが、しかも自力ではないのですが...
**Rが動作するOS [#lcec4b52]
> (2004-03-01 (月) 18:58:11)~
~
Rの紹介記事にRが動作するOS一覧を載せたいのですがどこかに...
現在の把握状況は下記の通りです。~
- MacOS X、Windows、Linux、FreeBSD、NetBSD、OpenBSD、Sola...
- 古いMacOSは古いRなら動く？~
- Linux版のザウルスは、制限がありながらも動作する。~
~
上記以外にRが動作したOSがありますか？~
LindowsやBeOS上でRは動作しますか？超漢字はどうですか？~
//
- 2002年2月の r-help 記事に WindowsCE ではビルドできない...
-Mac OS 9 以下では，R 1.7.1 までしか動きません。 -- &new...
**ESSでのコメント行の色 [#m589a9ba]
>[[さとう]] (2004-02-28 (土) 09:46:19)~
~
Emacs+ESS+Rを利用してます。emacsでRファイルを編集するとき...
とりあえずの対処療法として、~
M-x hi-lock-mode
M-x highlight-regexp [RET] ### [RET] hi-red-b
M-x highlight-regexp [RET] ## [RET] hi-green-b
などとしてます。(順序が大事)~
拡張子がR(またはr)のファイルを開いたときには、いつもこの...
//
**Rで並列処理計算（グリッド・コンピューティング？） [#i53...
> (2004-02-23 (月) 15:18:03)~
~
Rで利用できる並列処理のチュートリアルのようなものはある...
//
-例えば R news 2002/2 号に "Rmpi: Parallel Statistical Co...
**６４ビット環境でのRの動作 [#w32f2127]
> (2004-02-23 (月) 15:16:40)~
~
最近、６４ビットのPCが盛んに話題になっていますが、Rは３...
//
-C言語のdoubleは現在64bitであることが多いですから、double...
~
**階級間隔が異なるヒストグラム [#tf77c6ee]
>[[さかもと]] (2004-02-20 (金) 11:30:56)~
~
Rにて階級区間が異なるヒストグラムを作成したいのですが，ど...
(度数)*d(基準とする標準間隔とでも言いますか･･･)/f~
を間隔ｆの高さにすればいいのですが，どうすればRでそれが可...
ちなみに，初心者です･･･．~
//
//-タイトルの"感覚"=>"間隔"です(-_-;) -- [[さかもと]] &ne...
-引数 breaks (vector型)に、binの区切りとなる値を与えれば...
-ありがとうございます。さっそく試してみます。小生が試した...
-確か引数 width は S の頃に(関数density()などで)使われて...
> x <- rnorm(50) # データ
> h <- 0.8
> bins <- seq(min(x)-0.1, max(x)+0.1+h, by=h)
> pointlist <- hist(x, breaks=bins, prob=T, plot=F)$counts
> hist(x, breaks=bins) # 普通にヒストグラムを描く場合(区...
リストpointlistには，長さ h ごとに区切られた場合の左から...
**免疫アルゴリズム [#h86a6ecd]
> (2004-02-17 (火) 18:01:23)~
~
R や S-plus で免疫アルゴリズム（IA）を扱った例はあるので...
//
**ライブラリのバージョン表示 [#c441c04d]
> (2004-02-17 (火) 13:16:00)~
~
ライブラリパッケージが最新バージョンかどうかを確認するた...
しているライブラリのバージョンを表示させたいのですが、方...
//
-update.packages() で解決しました。 -- &new{2004-02-17 (...
-read.dcf(file=system.file("DESCRIPTION",package="tree"))...
**Rmapやmaptoolsでの縮尺バーや方位記号の表示 [#w57d3f4c]
> (2004-02-17 (火) 00:17:51)~
~
Rmapやmaptoolsで表示した地図には縮尺バーや方位記号を入れ...
//
-arrows などを使って自作するしかないのでは -- &new{2004-...
-凡例はデモにハワイの火山表示のコードが参考になるかも -- ...
**plotで縦軸と横軸の単位幅を同じ長さにする方法 [#ia0b6a17]
> (2004-02-16 (月) 19:10:03)~
~
FAQでしたら、削除しください。~
x軸とy軸がどちらも単位がm(メートル)であるplotで、x軸の1m...
長さにしたいのですが、勝手に調整されてそろいません。x軸と...
同じにする方法はありませんか。~
//
-xlim と ylim で座標軸範囲を同じにすれば良い？-- &new{20...
plot(runif(10), 2*runif(10),xlim=c(0,2),ylim=c(0,2))
-ありがとうございました。確かに、座標範囲を指定すると軸と...
-サンプルコードを紹介して頂かないと、コメントのしようがあ...
-例えば、下のような場合では、上に大きな余白ができます。 -...
plot(2*runif(10), runif(10,min=0,max=.1),xlim=c(0,2),yli...
-x,y 座標を同じスケールにする以上これは避け難い。グラフ領...
-申し訳ありません。最初にいうべきでしたが、x,y 座標が同じ...
-par(mai=c(1,2,6,2))ってしてみるとか，postscriptデバイス...
-dev.copy(x11,whdth=値,height=値)で書き直させる? -- &new...
-widthやheightを指定する方法は、作図領域の指定ではなく、...
-library(grid)を使えばできそうな気がします -- &new{2004-...
-正確な指定の仕方はわかりませんが、次のようにするとそれら...
> par(fin=c(6,3))
> plot(6*runif(100),1.5*runif(100),xlim=c(0,6),ylim=c(0...
比例定数 3/6 というのは本質的ではなく、実際の画面を見て適...
oldpar <- par(no.readonly = TRUE) # 現在のグラフィックス...
... 作業 ...
par(oldpar) # パラメータ復帰
を実行しないと後で困るかも知れません。
#ref(par.fin.jpg, left)
-追加。これでも実際は余計な余白がたっぷり付いている可能性...
-大体様子がわかったような気がします。グラフィックス画面の...
par()$fin[1]*(par()$plt[2]-par()$plt[1]) # x 方向のサイ...
par()$fin[2]*(par()$plt[4]-par()$plt[3]) # y 方向のサイ...
と決まります。従って x,y 座標のスケールを同じにするには x...
> par(fin=c(6,4)) # グラフィックス画面を 6x4 インチに...
> par(plt=c(0.1,0.9,0.1,0.9)) # 作図領域をグラフィックス...
> Xlim <- 10*par()$fin[1]*c(0, par()$plt[2]-par()$plt[1])
> Ylim <- 10*par()$fin[2]*c(0, par()$plt[4]-par()$plt[3])
> plot(1:20,1:20, xlim=Xlim, ylim=Ylim)
#ref(finplt.jpg, left)
参考までにメモすれば plt, fin の既定値は次のようになって...
> par()$plt
[1] 0.1172537 0.9399432 0.1457143 0.8828571
> par()$fin
[1] 6.990803 6.997416
従って既定値での作図画面のサイズはインチ単位でそれぞれ 5....
-お陰様でうまくいきました。ありがとうございました。 -- [[...
-その後、ただ単にasp=1というオプションを付ければよいこと...
例）
> x <- 1:10
> y <- x * 10
> plot(x,y,asp=1)
結果)
#ref(asp.png,center)~ -- &new{2004-05-14 (金) 23:14:18};
**正準相関分析での第 i 正準変量の求め方 [#r7a06ebf]
> (2004-02-14 (土) 20:55:31)~
~
library(mva)の関数cancor()を用いれば正準相関係数行列が得...
第 i 正準変量を求める場合はどうすればよいのでしょう？例え...
[[このような分析:http://bstat.f7.ems.okayama-u.ac.jp/stat...
（exampleでは x %*% xcoef と計算していますが何か違う気が...
//
-一般化固有値問題を解く短い関数を書いて，[[正準相関分析を...
~
**latex形式のヘルプファイルの見方 [#c96fa4c7]
> (2004-02-06 (金) 23:25:05)~
~
libraryごとにhtmlやhelpなどヘルプ文書がありますが、latex...
//
-ソースから R システムを構築するとき make install-dvi と...
-もしかすると [[森さんのWeb「R の日本語文章 (pdf 版)」:ht...
-全く無知で申し訳ないのですが、R CMD Rd2dvi パッケージ名....
-例. library/nlme/latex/nlme.tex だけを make する手順 (R ...
-なお, 上記の手順は, bin/Rd2dvi の内部で実行されているこ...
-わざわざ latex 形式のヘルプファイルを一揃い用意している...
-謎が解けました。やはり help 関数には tex 形式のヘルプフ...
> options(latexcmd="/usr/bin/latex")
> options(dvipscmd="/usr/bin/dvips")
> help(optim, offline=TRUE)
Saving help page to 'optim.ps' # optim.tex から optim.p...
思いついて次のようにしたら dvi ファイルのオンライン表示が...
> options(dvipscmd="/usr/bin/xdvi")
-おぉ、すごいです。無事にxdviでコマンドのhelpが読めました...
-さらにdvips を dvipdfm に変えてみました。すると無事にpdf...
-パッケージのhelpと各々のコマンドのhelpをxdviで表示する方...
-課題としては面白いですが、まず気になるのは、例えば base ...
-help() に offline という引数があったとは．勉強させてもら...
-dvipscmd をだましてやると、dvi ファイルも保存できます。...
cp ./*.dvi /tmp
というスクリプトファイル(名前を仮に makedvi )に実行属性を...
> options(dvipscmd="./makedvi")
としてやる。もう少しスクリプトを工夫する余地？
**文字列のオブジェクトへの変換方法 (なんでも掲示版より移...
>[[sakamoto]] (2004-01-30 (金) 19:35:15)~
~
まずは次のプログラムを見てください。
money<-1:10
>a<-readline()
money
>money[a>5]
[1] 1 2 3 4 5 6 7 8 9 10
> a<-as.name(a)
> a
money
> money[a>5]
[1] 1 2 3 4 5 6 7 8 9 10
ここでしたい事は、以下の事である。
>money[a>5]
[1] 6 7 8 9 10
つまり、文字列として取り込んだものをオブジェクトとして扱...
-ご質問の趣旨はこういうことですね。すでに a という名前の...
> a <- runif(10) # a という名前の数値ベクトルを作る
> a
[1] 0.6313190 0.4538556 0.3934363 0.4643964 0.3375608 0...
> aa <- readline() # その名前をコンマンドラインから入力
a # 名前は a
> aa
[1] "a" # aa は文字列 "a"
> get(aa) # aa という文字列を名前に持つオブジ...
[1] 0.6313190 0.4538556 0.3934363 0.4643964 0.3375608 0...
まとめてみると
> a <- runif(10)
> b <- runif(10)
> x <- 1:10
> x[get(readline())<0.5]
a # コンマンドラインから a を入...
[1] 3 4 5 6 10 # x[a<0.5] が得られた
> x[get(readline())<0.5]
b # コマンドラインから b を入力
[1] 1 3 4 5 6 7 8 # x[b<0.5] が得られた
-money[eval(parse(text=a))>5] としたいと言ってるのかな？ ...
-中間さんの解を注釈すると -- &new{2004-01-31 (土) 00:05:...
> a
[1] 0.557755694 0.756125389 0.062838319 0.116221907 0.28...
> y <- readline()
a
> y
[1] "a"
> eval(parse(text=y)) # 文字列 y をオブジェクト名 a に変...
[1] 0.557755694 0.756125389 0.062838319 0.116221907 0.28...
-ご回答ありがとうございます。おかげさまで問題が解決しまし...
**hist()の引数breaksの決め方 [#hc0fa05f]
> (2004-01-29 (木) 17:57:01)~
~
ヒストグラムを描くhist()の引数のbreaks=seq(min,max,by=bin...
//
-データ数が大きい場合(50以上)はあまり気にしなくても良いは...
-代替案的なものですが，ASH という方法があります．ただ，こ...
**ベクトルと行列 [#c6999be7]
>[[Soichi]] (2004-01-29 (木) 05:39:56)~
~
Matlabから乗り換えようかとRを試し始めました。~
ベクトルと行列のどちらも引数として受け取ることができ、~
行列を受け取ったときは各列に処理を、~
ベクトルを受け取ったときは(n,1)行列に変換した上で、その後...
Rでそのようなことは可能でしょうか。~
//
- as.matrix()をかませておけばいいのでは？ -- &new{2004-0...
- 次のようにすれば良い。問題によりますがいちいち行列に変...
if(is.vector(x)) x <- as.matrix(x) #もしくは
if(!is.matrix(x)) x <- as.matrix(x)
ちなみに縦ベクトルにしたければ
if(is.vector(x)) x <- t(as.matrix(x))
判断する必要がないのなら単に
x <- as.matrix(x) # x がすでに行列なら何もしない
**R の sleep() 機能 [#l2cf0b29]
>[[φ(oo-)o]] (2004-01-28 (水) 09:11:25)~
~
このサイトの「R プログラミング Tips 大全」のページで~
> 指定秒数プログラムをアイドリングする関数~
は標準関数の Sys.sleep(sec) ではどうでしょうか。ヘルプには~
> This function may not be implemented on all systems.~
と書かれていますが、Windows2000では動くようです。~
(R 1.8.1を使用しています)~
//
-本当ですね。知りませんでした。「R プログラミング Tips 大...
**作ったページの削除は？ [#ld06b03b]
>[[mnyu]] (2004-01-23 (金) 01:08:50)~
~
作ったページの削除は どうしますか~
//
-編集モードで内容をすべて削除し、更新すると一人でになくな...
**非線形には対応していますか [#jc5b6ba7]
>[[tia]] (2004-01-22 (木) 17:22:56)~
~
Rに関してまったくの初心者です．~
Rで非線形問題を解くことは可能でしょうか．~
ご教授いただければ幸いです．~
//
-問題によりますが、汎用最適化関数 optim や nlm を使えばた...
-早速のお返事ありがとうございます．いろいろと試してみます...
-非線形最小自乗法用関数 nls については [[Rの統計解析関数T...
**CCDFについて。 [#m5bd524c]
>[[sato]] (2004-01-22 (木) 14:24:52)~
~
Rを使って、1‐累積密度分布であるCCDFを描きたいのですが、パ...
またもしない場合は、パッケージstepfunの中のecdfを使って描...
//
-stepfun パッケージの ecdf では困る、もしくは避けたい理由...
-ecdfで困るということではないのですが、それ用に用意されて...
-無いようですね。累積和関数 cumsum を使えば自作は十分可能...
-わかりました。ありがとうございました。 -- [[sato]] &new{...
**Rでビット演算 [#je04e68d]
> (2004-01-19 (月) 16:29:43)~
~
Rでビット演算をするには、どうすれば、良いのでしょう~
//
-R にはそのものずばりのパッケージはないようです。具体的に...
-Colorオブジェクト同士のビット演算をしたいのですが。 -- ...
-
> heat.colors(9)
[1] "#FF0000" "#FF2A00" "#FF5500" "#FF8000" "#FFAA00" "#...
[8] "#FFFF40" "#FFFFBF"
のようなカラーを与える16進数(を表す文字列)間の and, or, x...
**直線回帰の信頼区間について [#z762ae44]
>[[ぶんしょ]] (2004-01-18 (日) 20:38:47)~
~
直線回帰の95％信頼区間をRで計算したいのですがどうしたらよ...
//
-predict()関数を使うのが便利です。直線回帰でいいなら，exa...
-ありがとうございます。さっそく試してみます。 -- [[ぶんし...
-プロットまでやる例をhttp://phi.ypu.jp/swtips/sumfunc.htm...
-久しぶりに見ました。丁寧にありがとうございます -- [[ぶん...
**スレッド表示 [#b83a13c4]
> (2004-01-13 (火) 21:19:06)~
~
Q & Aではできないのでしょうか？~
//
**[質問] lm() 関数の返り値のF統計量のp-value <= なんでも...
>[[舟尾]] (2004-01-10 (土) 21:12:33)~
~
lm()関数で重回帰分析を行った際，~
~
result <- summary(lm(y~x))~
~
として result［［k］］ や result$names (k:整数，names:名...
をすることで推定値やt値などが取り出せますが，F-statistics~
のp-valueだけがどうやっても取り出せません．names()で調べ...
F-statisticsのp-valueにあたるものが見当たらないのですが，~
この件に関して情報をお持ちの方，レス頂ければ幸いです．~
//
> [[間瀬]] &new{2004-01-10 (土) 22:52:45};
~
たしかに summary.lm 関数は F 統計量の p 値自体は計算しな...
> x=1:10
> y=c(rnorm(5),rnorm(5,mean=1))
> z=summary(lm(y~x))
> w=z$fstatistic
> pf(w[1], w[2], w[3], lower.tail=FALSE)
value
0.1643454
//
> [[中澤]] &new{2004-1-12（月）11:42:00};
~
print.summary.lmの中のようです。プロンプトでprint.summary...
if (!is.null(x$fstatistic)) {
cat("Multiple R-Squared:", formatC(x$r.squared, d...
cat(", Adjusted R-squared:", formatC(x$adj.r.squ...
digits = digits), "?nF-statistic:", formatC(x...
digits = digits), "on", x$fstatistic[2], "and...
"DF, p-value:", format.pval(pf(x$fstatistic[...
x$fstatistic[2], x$fstatistic[3], lower.t...
digits = digits), "?n")
}
となっているのがわかります。
-なるほど。気づきませんでした。R は部品化されているのでこ...
-間瀬さん、中澤さん、非常に丁寧な回答をありがとうございま...
**ファイルが読み込めません [#ca5b1182]
>[[haru]] (2004-01-09 (金) 14:52:09)~
~
初心者なのですが、R1.8.1をダウンロードして早速Excelでデー...
> customernumber <- read.table("costomernumber.txt")~
としましたが~
Error in file(file, "r") : unable to open connection~
In addition: Warning message: ~
cannot open file `costomertable.txt'~
というエラーメッセージがでてしまい読み込むことができませ...
read.delimなども使ってみましたが同じ結果でした。~
http://phi.ypu.jp/swtips/Data2R.htmlに書いてある説明通り...
Rに関する質問なのかもよくわかりませんが、よろしくお願いし...
//
-原因は，ファイル名の指定が正しくないことだと思われます。...
-できました！ご説明本当にありがとうございました！ -- [[ha...
**R初心者 [#y11bdb10]
>[[文字の読み込み]] (2003-12-16 (火) 01:33:52)~
~
文字列（例えば3030A6-0010など）の入ったファイルから各々の...
stnid <- cdataとしても、Error in scan("worldstnid.dat") :...
//
-scan("worldstnid.dat", what=character())とwhat=を引数で...
-scan にデータがどういうタイプか教える必要があります。既...
**macでRのproxy設定 [#kece459e]
>[[mac]] (2003-12-14 (日) 05:55:22)~
~
OSXのパンサーを使っています。~
Rでproxyを利用したいのですが、どうすればいいのでしょうか。~
教えてください。~
//
-すでに試されたかもしれませんが、googleでproxy site:r-pro...
-切れてしました。続き。proxyの設定方法が書いてありました...
-ありがとうございました。出来ました。他にも困っている人が...
-まず、インストールした直後に、Sys.putenv("http_proxy"="h...
-ftpも同様に入れてください。Rをインストール直後にこれを行...
-環境変数はプログラムの起動時ごとに設定されるので、インス...
-また、毎回起動時に環境変数を設定する必要をなくすには、/A...
-あと、インストールしたパッケージとかは ~/Libaray/RAqua/ ...
**同時実行 [#n213d61b]
> (2003-12-12 (金) 16:58:27)~
~
Windows用のRを2つ立ち上げて、~
それぞれを操作した場合に、~
それぞれの実行結果は完全に独立だと思ってよいのでしょうか？~
//
-独立とは統計的な独立（乱数を使用した場合に限り意味のある...
**グラフがかけない X11 is not available [#od5d3d06]
>[[TAKU]] (2003-12-12 (金) 11:33:19)~
~
超初心者で、はじめてLinuxをインストールし、Rをemacs上で立...
とこまでできたのですが、postscript以外のグラフィックコマ...
使えないのです。png("hoge")やx11()どを入力すると~
X11 is not available~
とエラーメッセージがでてしまいます。Rの質問ではないのかも...
助け舟をよろしくお願いします。~
ディストリビューションはdebianです。~
//
-X Window Systemは起動しているんですね？それならば X11(di...
-Error in x11(display = ";0") : X11 is not available と...
-私もdebianです。同じことがおこったことがあります。もしか...
-で、debianのX11クライアント開発環境のインストールですが...
-もしインストールされていなければ、apt-get -s install xli...
-開発環境がインストールできたら、Rのsourceに戻って、./con...
-私は試していませんが、make が面倒なら、www.r-project.org...
-detourさんほんとうにいろいろありがとうございます！（;_;)...
-ありがとうございました！解決しました。このことでずっと迷...
-いえいえ、Knoppix-jp CDROM から HD へインストールすれば...
-[[テスト中CRANミラー:http://cran.okada.jp.org/bin/linux/...
-ほんとにいろいろありがとうございました。まだまだ初心者で...
**Rプログラミングで得られた結果を別のファイルへ出力するに...
>[[ぽいと]] (2003-12-11 (木) 18:21:16)~
~
現在、R言語でプログラムを書いたのはいいのですが、その結果...
//
-プログラムの冒頭でsink("c:/temp.txt")のようにすれば、tem...
-Tips 集の [[R出力の記録]] が参考になります。 -- &new{20...
-早速の回答をありがとうございます。投稿後、色々調べてみた...
**Rパッケージのインストール＆ロード方法 [#i0c435fe]
>[[R初心者]] (2003-12-11 (木) 16:03:04)~
~
現在Linux上でパッケージ「evd」をinstall.package("evd")で...
//
-該当パッケージをダウンロードした上、ルートになって R CMD...
**大規模データからランダムにデータをサンプリングする [#ma...
>[[msasalo]] (2003-12-10 (水) 21:48:07)~
~
はじめまして。最近Rを使いはじめるようになった者です。~
200万件38Mbのデータを処理しようと思っていたのですが、読み...
そこで、データのサブセットを作って36万件6Mbのデータを作っ...
マシンのスペックは、P4-1.6GHz,512Mで、ページングメモリー...
データファイルを元にランダムにとる方法がRにあるとありがた...
//
- sample 関数が探しているものでは？-- &new{2003-12-10 (...
> x=rnorm(2000000) # 200万個の正規乱数発生
> length(x)
[1] 2000000
> 2000000*8/1024/1024
[1] 15.25879 # 15 Mbyte
> xx = sample(x, 10000) # その中から一万個を無作為抽出
> length(xx)
[1] 10000
//
-Rに読み込んであるデータならばsample()という関数が使えま...
-そうなのです、以前はPostgreSQLを立ち上げてデータをODBCと...
-ありがとうございました!! なんとか解決しました!! Rで216万...
**[質問]2変量t乱数の作り方 [#r417af38]
> (2003-12-10 (水) 16:59:27)~
~
2変量乱数の発生でt乱数の作り方を教えて下さい。
//
-CRAN にあるアドオンパッケージ mvtnorm は多変量の正規分布...
アルゴリズムを書く代わりに命題みたいなものを示します。
- R：自由度mのカイ二乗分布に従う確率変数
- Z：p変量正規分布N(0,Ip)に従う確率ベクトル（独立な標準正...
- V：正定値対称行列（ちらばりを表す行列で分散共分散行列と...
- V＝C%*%t(C)：コレスキー分解
-- X = sqrt(m/R)*Z はp変量楕円t分布met(p,m,0,Ip)に従います
-- Y=μ+C%*%X はp変量楕円t分布met(p,m,μ,V)に従います
拡張性があるように2変量ではなく3変量楕円t乱数に従う乱数を...
met3 <- function(m, mu, V, n) {
# m : 自由度
# mu : 平均ベクトル
# V : 散らばり行列
# n : 乱数の個数
U <- svd(V)$u
V1 <- svd(V)$v
D <- diag(sqrt(svd(V)$d))
B <- U %*% D %% t(V1)
w <- c()
for (i in 1:n) {
R <- 0
for (j in 1:m) {
R <- R + rnorm(1)^2
}
w <- append(w, list(mu + B %*% (cbind(rnorm(3))*sqrt...
}
return(w)
}
# 実行例
mu <- cbind(c(1,1,1))
V <- array(c(2,1,1,1,2,1,1,1,2), dim=c(3,3))
n <- 1000
w <- met3(5, mu, V, n)
sm <- 0
for (i in 1:n) { sm <- sm + w[[i]] }
sm <- sm/n
sv <- 0
for (i in 1:n) { sv <- sv + (w[[i]]-sm) %*% t(w[[i]]-sm) }
sv <- sv/n
この分布の平均ベクトルはμ、分散共分散行列はm/(m-2)*V（m>2...
多変量楕円t乱数がうまく出来ているかどうかはこれで確かめる...
**[質問]hist()関数とnclass()関数の関係 [#bdaab073]
> (2003-12-06 (土) 11:09:38)~
~
どちらの関数でもSturgesの方法によって階級が求められるよう...
//
-ご質問の意味が今一不明確ですが、?hist に書いてあるように...
-- (1) 単に hist(x) は hist.default(x, breaks = "Sturges...
-- (2) 明示的に分割点ベクトルを与えても良い、
-- (3) 単に整数 n を一つを与えれば、区間が n 等分される、
-- (4) 分割点を計算するユーザー定義関数を与える
ということでしょうか。 -- &new{2003-12-06 (土) 17:00:57};
//
-具体的には(1)に関してのことになりますが、以下ではsは常に...
x<-rnorm(80,mean=0,sd=1)
s<-nclass.Sturges(x)
b<-length(hist.default(x,breaks = "Sturges")$breaks)
s;b
-たしかに変ですね。関係する関数は内部関数なので詳細はソー...
> test <- function(x) {c(nclass.Sturges(x),nclass.scott...
+ length(hist(x)$break...
> test(rnorm(1000))
[1] 11 21 26 16 16
> test(runif(1000))
[1] 11 10 10 11 11
> test(rnorm(1000))
[1] 11 19 25 15 15
> test(rexp(1000))
[1] 11 20 32 15 15
> test(rcauchy(1000))
[1] 11 109 3057 13 13 # ずいぶんと違うもんですね...
-R の ML に Sturges 方式を hist の既定にするほうが良いか...
- おなじことは他の方式でもおきているようです。一つのあり...
> nclass.scott(x)
[1] 9
> length(hist(x,breaks="scott")$breaks)
[1] 8
> nclass.FD(x)
[1] 18
> length(hist(x,breaks="FD")$breaks)
[1] 14
> nclass.Sturges(x)
[1] 9
> length(hist(x)$breaks)
[1] 8
**R(D)COMの使い方 [#l18c0057]
>[[f]] (2003-12-06 (土) 00:37:19)~
~
R(D)COMを使って、VC++で作ったプログラムにRを組み込みたい...
R(D)COMをお使いの方がおられましたら、使用法等、ご助言お願...
//
**[質問]条件分岐 [#l0e4043d]
>[[さかい]] (2003-12-04 (木) 18:59:26)~
~
i<500のときx=1、500<=i<1000のときy=1、~
それ以外のときz=1とするプログラムを作成したかったので、~
以下のようなプログラムを作成したのですが、~
iが1000以上のときx=0,y=1,Z=0となり、うまく値を返してくれ...
どのようにしたら解決できるのか分かる方がいたらアドバイス...
~
以下のプログラムはRプログラミングTips大全集の条件実行 if...
#プログラム
x <- 0
y <- 0
z <- 0
i <- 2000
if (i < 500) {x <- 1} else if (500 <= i < 1000) {y <- 1...
#結果
> x
[1] 0
> y
[1] 1
> z
[1] 0
//
-disc <- function(i) { x <- 0; y <- 0; z <- 0; if (i<500)...
-２番目のifは１番目のifに対するelseの後なので，500<=iの部...
-ありがとうございます。やってみます。 -- [[さかい]] &new{...
-eval(500 <= 2000 < 1000)はなぜTRUEを返すのですか？ -- [[...
-if (i < 500) {x <- 1} else if (500 <= i && i < 1000) {y...
x <- ifelse(i < 500, 1, 0)
y <- ifelse(i < 1000, 1, 0)
z <- ifelse(i >= 1000, 1, 0)
ついでですが、今チェックをしていて気づきましたが R 1.8 か...
> test <- function(x,y) return(x,y)
> test(1,2)
$x
[1] 1
$y
[1] 2
Warning message:
multi-argument returns are deprecated in: return(x, y)
-- &new{2003-12-04 (木) 21:38:49};
-もっと簡単に(+FALSE は論理値を整数に強制変換するトリック...
> x = (i < 500)+FALSE; y = (i < 1000)+FALSE; z=(i >=100...
> x;y;z
[1] 0
[1] 0
[1] 1
-- &new{2003-12-04 (木) 22:02:08};
**Rdbi RdbiPgSQL [#r0441a96]
>[[なかま]] (2003-12-04 (木) 10:42:45)~
~
http://www.bioconductor.org/ のが 1.0.2 で cran.r-project...
bioconductor のが最新?~
//
-ところで、RdbiPgsql はCRAN にはバイナリがありませんが、...
-[[PostgreSQL Native Win32:http://momjian.postgresql.org/...
-やってみると、面倒くさかったのでVine2.6(PostgresSQL-7.2....
-Rdbiはbioconductorの方にバイナリがあります。 -- [[なかま...
-やっぱりbioconductorのが最新のようですね(他にもあったり...
-なかまさん、RdbiPgSQLのバイナリありがとうございます。わ...
-libpqとpostmasterのバージョンが異なると通信出来なかった...
**［質問］棒グラフで対数軸 [#k47e8d60]
> (2003-11-29 (土) 17:20:27)~
~
棒グラフを作図するのには barplot(1:10) のようにすれば良い...
そこで、Y 軸を対数軸として作図しようと log="Y" を指定した...
> barplot(1:10, log="y")~
Error in plot.window(xlim, ylim, log = "", ...) :~
formal argument "log" matched by multiple actual arguments~
対数軸の棒グラフはどのように指示すれば描くことが出来るの...
よろしくお願いします。~
//
-残念ながら base パッケージの barplot 関数は対数軸をサポ...
-ありがとうございます。アドオンパッケージを調べてみます。...
**subsetで値を取り除く [#p4d4e573]
>[[detour]] (2003-11-28 (金) 14:49:55)~
~
subsetで名義変数の特定の値を除いたサブセットを作成したと...
a<-data.frame(choice=c("ab","ab","ab","c","c","d","d","d"...
summary(a)~
summary(subset(a,choice!="c"))~
"c"というのが一応残っているようです。これはどうやったら消...
//
-次のようにすれば一応ラベル c は消えますが、おそらくこれ...
> b=data.frame(choice=a[a!="c"])
> summary(b)
choice
ab:3
d :3
-使用の目的にもよるのでしょうが、気にしないのもありかと。...
-消せないせいでうまく行かなかったことがありましたが、勘違...
-おまけ:subsetを見付ける前はbquoteを使ってみようと思った...
-気にしない事にしましたがやはりうまく行かない事もありまし...
-で、できました。b<-subset(a,choice!="c")だとすると、leve...
-おめでとうございます。苦労して得られたノウハウは是非然る...
**Rのインストール方法（UNIX） [#ja1a0abf]
>[[R初心者]] (2003-11-28 (金) 14:21:53)~
~
ただいまRをUNIXにインストールしようとしています。R-1.8.1...
> make~
Make: Cannot open /share/make/vars.mk. Stop.~
vars.mkファイルは確かに/share/make/に存在しましたし、原因...
//
-Debian GNU/Linux 3.0ですが、問題なくmakeできました。Unix...
-使っているUNIXマシンはCompaq Tru64です。 -- [[R初心者]] ...
終了行:
COLOR(blue){SIZE(20){R および RjpWiki に関する質問コーナ...
COLOR(red){注意}：新規記事用の入力欄は以下の目次の直後に...
もし新規に質問をされる場合は、効率的にコメントを得るため...
- 広い意味で R に関係する話題に限定!
- 学校の宿題は自分で考えるべきです
- 具体的な状況がわかる様に背景説明をけちらない
- [[CRAN の検索エンジン:http://finzi.psych.upenn.edu/sear...
- 自分には馴染みでも、他の人には馴染みの無い手法・概念も...
- もし問題が解決した(結局解決しない場合を含め)場合はその...
~
~
COLOR(red){SIZE(18){注意：このコーナーは重量オーバーで表...
~
----
#contents
//#article
~
~
**群散布図（Group Scatter）を作成するパッケージはあります...
> (2004-05-13 (木) 22:15:48)~
~
[[ここ:http://www.osaka-amt.or.jp/lecture/stat/index2.htm...
//
- パッケージ ade4 にそれらしき関数が多数用意されています...
> library(ade4)
> data(lascaux)
> scatter(dudi.acm(lascaux$ornem, sca = FALSE), csub = 3)
#ref(plot-00.png,center)
例えばこんな関数がありました．
> x <- c(0.5,0.2,-0.5,-0.2) ; y <- c(0.2,0.5,-0.2,-0.5)
> eti <- c("toto", "kjbk", "gdgiglgl", "sdfg")
> plot(x, y, xlim = c(-1,1), ylim = c(-1,1))
> scatterutil.eti.circ(x, y, eti, 2.5)
> abline(0, 1, lty = 2) ; abline(0, -1, lty = 2)
> x <- c(0.5,0.2,-0.5,-0.2) ; y <- c(0.2,0.5,-0.2,-0.5)
> eti <- c("toto", "kjbk", "gdgiglgl", "sdfg")
> plot(x, y, xlim = c(-1,1), ylim = c(-1,1))
> scatterutil.eti(x, y, eti, 1.5)
> plot(runif(10,-3,5), runif(10,-1,1), asp = 1)
> scatterutil.grid(2)
> abline(h = 0, v = 0, lwd = 3)
> x <- runif(10,0,1) ; y <- rnorm(10) ; z <- rep(1,10)
> plot(x,y) ; scatterutil.star(x, y, z, 0.5)
> plot(x,y) ; scatterutil.star(x, y, z, 1)
> x <- c(runif(10,0,0.5), runif(10,0.5,1))
> y <- runif(20)
> plot(x, y, asp = 1) # asp=1 is essential to have perpe...
> scatterutil.ellipse(x, y, rep(c(1,0), c(10,10)), cell ...
> scatterutil.ellipse(x, y, rep(c(0,1), c(10,10)), cell ...
> x <- c(runif(100,0,0.75), runif(100,0.25,1))
> y <- c(runif(100,0,0.75), runif(100,0.25,1))
> z <- factor(rep(c(1,2), c(100,100)))
> plot(x, y, pch = rep(c(1,20), c(100,100)))
> scatterutil.chull(x, y, z, opt = c(0.25,0.50,0.75,1))
> par(mfrow = c(1,1))
#ref(plot-01.png,center)
関数が多数あるので全部は調べ切れませんでした．すみません...
-いろいろ調べるのが面倒なので，作る方が楽。 -- [[青木繁伸...
# テストデータの作成
mean <- c(3,5,6,8,7)
sd <-c(0.3, 0.4, 0.5, 0.2, 0.4)
x <- rep(1:5, each=20) # 5群（20ケースずつ）のデータ，平...
y <- NULL
for (i in 1:5) y <- as.integer(100*c(y, rnorm(20, mean=m...
#accu 測定値（縦軸）の精度（階級幅）
#stp 横軸方向のずらし量
# 以上二つは，試行錯誤で決定(^_^;)
graph <- function(x, y, accu, stp)
{
y <- round(y/accu)*accu
x1 <- unique(x)
for (i in 1:length(x1)) {
freq <- table(y[x==x1[i]])
for (j in 1:length(freq)) {
if (freq[j] >= 2) {
for (k in 1:length(y)) {
if (abs(y[k]-as.numeric(names(freq)[j])) < 1e-...
x[k] <- x[k]+(freq[j] <- freq[j]-1)*stp
}
}
}
}
}
plot(x, y)
}
graph(x, y, 0.2, 0.05)
#ref(graph.png,center)
-質問した者です。ありがとうございました。青木先生の関数を...
-凄いっ！いつもながら鮮やかなソースコードですね．簡潔さと...
**巡回セールスマン問題 [#icaee477]
> (2004-05-06 (木) 17:15:50)~
~
巡回セールスマン問題用の関数がRのあるパッケージにあったよ...
//
-?optim です。なんでも掲示板で以前教えていただきました。 ...
-example(optim) を実行すると、組み込みデータセット eurodi...
**matplotの軸 [#g1ba8468]
>[[さち]] (2004-05-03 (月) 00:00:58)~
~
matplotで軸の正の方向を左向きにしたいと思い、~
~
matplot(c(30,-80),c(50,-50))~
~
のようにしたのですが、グラフにしてみると軸は正の方向が右...
matplotで軸の方向を変える方法を知っている方がいたらアドバ...
//
-x 軸の正の方向を左向きにする場合は matplot() に引数 xlim...
xlim で調節してください．この方法は他の高水準作図関数 (例...
でも使えます．
sines <- outer(1:20, 1:4, function(x, y) sin(x / 20 * pi...
matplot(sines, xlim=c(20,1), type="o", pch=1:4, col = ra...
#ref(matplot-00.png,center)
y 軸側も同様にして 引数 ylim で調節してください．-- [[舟...
sines <- outer(1:20, 1:4, function(x, y) sin(x / 20 * pi...
matplot(sines, pch = 1:4, xlim=c(20,1), ylim=c(1,-1),
type = "o", col = rainbow(ncol(sines)))
-xlimを求めたい場合
matplot(c(30,-80),c(50,-50))
matplot(c(30,-80),c(50,-50),xlim=par("xaxp")[2:1]
とか、2度描画させるのが一番楽なのでしょうか? -- &new{200...
-matplot() で描いた直後は座標系の値がグラフィックスパラメ...
par(xaxp=c(0,1,5))
par("xaxp") # xlim は [0 1]
[1] 0 1 5
plot(1:10) # xlim ( c(xaxp[1],xaxp[2]) ) がこの...
par("xaxp") # xlim は [2 10] に
[1] 2 10 4
par(xaxp=c(0,1,5)) # 自分好みの座標系にする場合はxaxpを...
plot(1:10, xlim=c(xaxp[1],xaxp[2]))
par("xaxp")
[1] 0 1 5
余談ですが，par() で設定して (永続的に) 自分好みの座標を...
のでしょうか？"usr" や "xaxp" は値がコロコロ変わってしま...
仕方なく xlim , ylim で一時的に変更しているのですが．どう...
せいか，"usr" や "xaxp" の適当な使い方が思いつきません．-...
# 原点を通る座標軸を描く
plot(rnorm(30), rnorm(30), axes = F, xlab = "", ylab = "")
axis(1, pos = 0, at = pretty(par("usr")[1:2]), adj = 0)
axis(2, pos = 0, at = pretty(par("usr")[3:4]), las = 2, ...
box()
**[質問]ベクトル [#vd5e3dc9]
>[[さち]] (2004-04-30 (金) 05:19:49)~
~
ある変数に値を代入していき、その変数に値をベクトルのように~
記憶させていくプログラムを作りたいです。例えば、~
x<-1~
x<-3~
x<-9~
x~
1 3 9~
のようにしたいです。このとき、~
x[3]~
のように長さを与えれればよいのですが、~
いま代入する変数の数がいくつになるのかわからない状態です。~
このような時はどうしたらよいのでしょうか？~
初歩的な質問ですが、アドバイスなどお願いします。~
//
- R のベクトルは任意に長さを延長できますから気にせずに代...
> x <- numeric(0) # x <- NULL でもOK （この時点では x ...
> x
numeric(0) # この時点では x は空の変数
> x[1] <- 1
> x[2] <- 2
> x[3] <- 3
> x
[1] 1 2 3
> x[5] <- 5 # 途中飛ばしてもOK
> x
[1] 1 2 3 NA 5 # 未定義の x[4] には自動的に NA 値が...
-append という関数を使うと吉-- [[青木繁伸]] &new{2004-04-...
> x <- NULL
> x <- append(x, 1)
> x <- append(x, 3)
> x
[1] 1 3
> x <- append(x, 2, after=1)
> x
[1] 1 2 3
-ありがとうございます。 -- [[さち]] &new{2004-04-30 (金) ...
-ありがとうございました。早速ためしてみます！ -- [[さち]]...
**princompのloadingsは固有ベクトルそのもの？ [#a1152f15]
>[[石原茂和]] (2004-04-27 (火) 07:42:09)~
~
（すみません，同内容を旧Rwikiの掲示板にまちがえてかいてし...
~
こんにちは．いま気がついたのですが，princompのloadingsは~
固有ベクトルそのものですね．~
helpでも~
loadings: the matrix of variable loadings (i.e., a matri...
contain the eigenvectors).
~
と，そのように書いてあります．~
しかし，主成分負荷量は，加工していないもとの変数での値と...
のはずですが，~
確かに固有ベクトルとたいがいの場合，同じような傾向になり...
固有ベクトルそのものを負荷量というのは誤解のもとでは？~
//
-example(Harman74.cor) を実行してみてください。おっしゃる...
-出ますけど，原質問者の質問の答えになっているのかな？~
それと，必ず回転された答えが出るんでしょうか。何回転にな...
二変数の分析してみると，教科書と違った解になるようですが。~
> x
[,1] [,2]
[1,] -0.24451455 0.5090637
[2,] 0.04563039 -0.5012150
[3,] 0.97775313 -1.1832550
[4,] 1.26989942 1.4814664
[5,] -0.14870401 -1.3809846
> loadings(princomp(x))
Loadings:
Comp.1 Comp.2
[1,] 0.273 0.962
[2,] 0.962 -0.273
Comp.1 Comp.2
SS loadings 1.0 1.0
Proportion Var 0.5 0.5
Cumulative Var 0.5 1.0
-- &new{2004-04-29 (木) 19:22:43};
-失礼、Harman74.cor の example は princomp ではなく、fact...
**数量化１類で消費されるメモリ量について。 [#m90801fb]
>[[koichi]] (2004-04-19 (月) 23:48:32)~
~
はじめまして。Rの初心者です。~
Rで数量化1類を行おうとしているのですが、メモリの問題で手...
具体的には、~
１列目(V1)を『0から1までの８桁の乱数』、２列目(V2)を『1か...
~
d0$V2 <- factor(d0$V2)
result <- summary(lm(V1 ~ V2,data = d0))
~
すると~
~
Error:cannot allocate vector of size 390625 Kb
Reached total allocation of 1022Mb
~
と出て、タスクマネージャーを見るとRは「応答なし」で止まっ...
~
このデータ量で390625 Kbもの行列を生成してしまうのでしょう...
~
メモリ関連の操作や、計算法の工夫を通してなんとか計算させ...
~
当方の環境はWindowsXP,Pentium4 2.8GHz,1GBRAMです。~
よろしくお願いいたします。~
//
-変数の長さが n ならば、回帰分析で登場するハット行列は nx...
-一つのやりかたは、必要な情報だけを公式通りに自前で計算す...
# 参考 10055007 は文字列、数値の項目数ですから、
# メモリー量はざっとその 8 倍程度になると考えるべきでし...
> n=100; x=rnorm(n); y=factor(floor(1000*runif(n))); z=l...
[1] 10902
> n=1000; x=rnorm(n); y=factor(floor(1000*runif(n))); z=...
[1] 661272
> n=10000; x=rnorm(n); y=factor(floor(1000*runif(n))); z...
[1] 10055007 # z に含まれる項目総数
~
-提示された条件に従ってデータを生成して数量化I類を（定義...
そもそも，このような設定の数量化I類が必要なのか適切なのか...
-解析対象のデータがわかりませんから、野次馬コメントですが...
-みなさん、早速のご回答ありがとうございます。また、つたな...
参考のために提示された式、とてもためになりました。Rに習熟...
提示された条件に従ってデータを生成して数量化I類を（定義に...
-数量化I類のプログラムは，[[ここ:http://aoki2.si.gunma-u....
for を使わないとかえって複雑怪奇になるところもあるので（...
実行したのは Mac G5 です。~
必要メモリー数ですが，原データが50000×2×8バイト，原データ...
cov 関数で分散共分散行列を計算するようにしたので，展開し...
プログラミング言語で書くときには，こんなにメモリを贅沢に...
**.Call/.Externalで使用するDLLの作成方法について [#p9d0c3...
>どうむ (2004-04-19 (月) 15:20:39)~
~
現在R（Windows版）から呼び出せるCで書いたDLLを作成中です...
■R1.8.1~
■コンパイラー：「Borland C++ Compiler 5.5」~
■OS：WindowsXP~
です。いままでRやSPlusでCで書いたルーティンを作ってきたの...
まず、簡単なCルーティンですが
extern "C" __declspec(dllexport) void R_TEST001(int *x, ...
__declspec(dllexport) void R_TEST001(int x[], int *ret)
{
*ret = x[0]+x[1]+x[2]+x[3];
}
整数型のベクトルを引数にとってその合計を「ret」で返してい...
DLLを作成するのにMakefileを作っておりコンパイルとリンクは...
bcc32 -u- -I"***" -O2 /c ***.cpp
ilink32 /w /m /Tpd /L"***?Lib"
"C:?borland?bcc55?Lib?C0D32.OBJ" ***.obj,
***.dll,
CW32.LIB IMPORT32.LIB
（mapファイルの指定が欠如していますが気になさらぬよう。...
「***.cpp」が該当のCルーティンです。~
これを呼び出すRのプログラムは
> dyn.load("***.dll")
> is.loaded("R_TEST001") #確認
[1] TRUE
と、なります。結果は
> ret <- .C("R_TEST001",as.integer(c(1,2,3,4)),ret=integ...
> ret$ret
[1] 10
です。~
ここまでは、良いのですがRの「.Call」関数を使用できるよう...
#include <R.h>
#include <Rinternals.h>
extern "C" __declspec(dllexport) void R_TEST001(int *x, ...
__declspec(dllexport) void R_TEST001(int x[], int *ret)
{
*ret = x[0]+x[1]+x[2]+x[3];
}
たった、これだけで上のコンパイルが通らなくなりました。~
が以下のように「-A」オプション（ANSI 準拠の予約語のみを使...
bcc32 -u- -I"***" -O2 -A /c ***.cpp
ilink32 /w /m /Tpd /L"***?Lib"
"C:?borland?bcc55?Lib?C0D32.OBJ" ***.obj, ・・・
今度は、「Writing R Extensions」の第４章にあるように関数...
#include <R.h>
#include <Rinternals.h>
extern "C" __declspec(dllexport) void R_TEST001(int *x, ...
extern "C" __declspec(dllexport) SEXP R_TEST003(SEXP x, ...
__declspec(dllexport) void R_TEST001(int x[], int *ret)
{
*ret = x[0]+x[1]+x[2]+x[3];
}
__declspec(dllexport) SEXP R_TEST003(SEXP x, SEXP ret)
{
int i, j, nx, ny;
SEXP ans;
nx = length(x);
PROTECT(ans = allocMatrix(REALSXP, nx, 1));
UNPROTECT(1);
return(ans);
}
このようにすると、コンパイルは通るのですがリンクができま...
ilink32 /w /m /Tpd /L"***?Lib" "C:?borland?bcc55?Lib?C0D...
***.dll, CW32.LIB IMPORT32.LIB
Turbo Incremental Link 5.00 Copyright (c) 1997, 2000 Bor...
Error: 外部シンボル 'Rf_length(SEXPREC *)' が未解決（***...
Error: 外部シンボル 'Rf_allocMatrix(unsigned int, int, i...
Error: 外部シンボル 'Rf_protect(SEXPREC *)' が未解決（**...
Error: 外部シンボル 'Rf_unprotect(int)' が未解決（***.OB...
** error 2 ** deleting ***.dll
R.DLLから「implib.exe」でR.libを作成してそれを指定しても...
このようなエラーが出るということは、致命的な間違いを、犯...
解決につながるような情報がございましたら、よろしくお願い...
//
-関数が「Rprintf("Hello!?n");」のみだったら通りますね。 -...
> ret <- .Call("R_TEST003",as.double(c(1,2,3,4)),ret=dou...
Hello!
> ret
[1] 1 2 3 4
ということは、
length
allocMatrix
PROTECT
UNPROTECT
が、いけないことになりますね。適切なヘッダーを指定してい...
なにか、大きな過ちを犯しているような。。。~
何を見落としているのでしょう？？？
//
-自己レスばかりですがなんとなくわかりました。 -- どうむ &...
~Rprintfが定義してある「Print.h」には
#ifdef __cplusplus
extern "C" {
#endif
とあって、C++でコンパイルされることを想定しているようです...
ですので、C++ではなくCだと思わせてコンパイルをかけること...
「Rinternals.h」に「extern "C"」を書けばC++でもコンパイル...
同一ファイルでC++のクラスが使えないのは不便かなぁ。
-こうした使い方ができることは噂には聞いていましたが、実際...
-これにあたっては[[Rから他言語利用]]を、参考にさせていた...
**Mac OS X の作図デバイス [#c2c090f8]
>[[舟尾]] (2004-04-18 (日) 15:49:17)~
~
R を起動しますと，UNIX では x11 が，Windows では windows ...
//
-Mac OS X では、x11 が起動することも、Aqua (Mac OS X のグ...
-回答ありがとうございます．なるほど，Aqua を起動させるこ...
-Linux系の場合、Rの起動コマンドオプションに --gui=X11(デ...
-回答ありがとうございます．勉強になりました．基本はどのOS...
-バイナリのレベルでは、RAqua.dmg (1.9.0ではR.dmg)からイン...
**GLMM()でできるオブジェクトの参照方法 [#qe73ec2c]
>林 (2004-04-14 (水) 13:38:20)~
~
lme4パッケージのGLMM()でできるオブジェクトは~
is.list()で調べるとリストなのですが、names()では名前はな...
list()
となります。通常のリストでは
List of 5
などといくつ要素があるとか返ってくるので、GLMM()でできる...
//
**1.9.0がmakeできない [#j38a01fc]
>[[谷村]] (2004-04-13 (火) 00:42:16)~
~
VineSeedでR-1.9.0がmakeできません。~
making dataentry.d from dataentry.c
making devX11.d from devX11.c
making rotated.d from rotated.c
making rbitmap.d from rbitmap.c
make[4]: ディレクトリ `/home/hoge/tmp/R-1.9.0/src/module...
make[4]: ディレクトリ `/home/hoge/tmp/R-1.9.0/src/module...
gcc -I. -I../../../src/include -I../../../src/include -...
-I/usr/local/include -DHAVE_CONFIG_H -D__NO_MATH_INLI...
-g -O2 -c dataentry.c -o dataentry.lo
In file included from dataentry.c:31:
/usr/X11R6/include/X11/Xlib.h:1400: error: 文法エラー be...
/usr/X11R6/include/X11/Xlib.h:1488: error: 文法エラー be...
/usr/X11R6/include/X11/Xlib.h:1516: error: 文法エラー be...
/usr/X11R6/include/X11/Xlib.h:1520: error: 文法エラー be...
/usr/X11R6/include/X11/Xlib.h:1542: error: 文法エラー be...
[snip]
誰かVineSeedでmakeに成功している方はいらっしゃいますか？~
XはXOrg-6.7.0-0vl2です。
//
-コンパイルが通りました。dataentry.cの#include <X11/X.h> ...
**RjpWiki 頁のアーカイブのうまい作り方？ [#r1012a9e]
>[[間瀬茂]] (2004-04-04 (日) 23:06:46)~
~
例えばこのＱ＆Ａコーナーがかなり長くなってきたので、アー...
//
**SJava [#u4f6bbf5]
>[[高階知巳]] (2004-03-31 (水) 10:56:07)~
~
S/RとJavaのインタフェース SJava を、1年半ほど前に試したと...
//
-SJavaうまくいきました。omegahatは今は大丈夫なようで、し...
import org.omegahat.R.Java.*;
public class REvalSample {
public static void main(String [] args) {
String [] rargs = {"--slave", "--vanilla"};
System.out.println("Sample program to call R engi...
ROmegahatInterpreter interp
= new ROmegahatInterpreter(ROmegahatInterpret...
false);
REvaluator e = new REvaluator();
Object val = e.eval("x <- sin(seq(0, 2*pi, length...
val = e.eval("x * 2.0");
if (val != null) {
double[] objects = (double[])val;
for (int i=0; i<objects.length; i++) {
System.err.println("("+i+") " + objects[i...
}
}
}
}
-R-1.9.0+SJava-0.67.3でコンパイルは成功しました．が，実行...
**pairs(data.frm, upper.panel = points, lower.panel = 残...
>[[みやむら]] (2004-03-26 (金) 13:42:27)~
~
pairs() の描画のようにして，上三角部分に通常の散布図を，...
//
-[[pairs用 回帰直線つきpanel]]を参考に対応する lower.pane...
-単回帰の残差であればよいのですが，この場合重回帰ですので...
**help.startができなくなりました。 [#mb2d87a7]
> (2004-03-26 (金) 13:15:11)~
~
help.startでブラウザが起動しなくなりました。どこをチェッ...
ご存知の方は教えてもらえないでしょうか？~
こんなコメントが出ます。Making links in per-session dir ....
使用ソフトRAqua ver1.8です。~
//
-Macintosh HD/private/tmp/Rtmp988/.R/doc/html/index.html ...
-GUI で見えないのでどうでもいいのですが，Rtmp988 みたいな...
-私の場合だと，その index.html ファイルを jedit で開くよ...
-コメントどうもです。アドバイスどおりやってみましたが、、...
-open index.htmlだとちゃんとひらくのですが、、この前OpenO...
**2変量のヒストグラム [#y279fd4d]
>[[sak]] (2004-03-24 (水) 23:53:51)~
~
２変量のヒストグラムを描く事は可能でしょうか？~
//
-何種類かパッケージがありますが，その一例は[[ここ:http://...
-うーん、なんでもできそうですね。すごい。これからRの世界...
**[質問]DOSの |more や /p にあたるコマンドは？ [#cfdd0...
>[[inagaki00]] (2004-03-23 (火) 22:27:03)~
~
長ーいデータを見るとき、DOSでは~
>type hogehoge.txt | more とか~
>type hogehoge.txt /p のように打ちます（UNIXだとcatですね...
これと同じスイッチがRにはありますでしょうか。~
~
ありそうで見つからないので困っています。~
ご存知の方、どうぞよろしくお願い致します。~
//
-普段は x[1:100] のように一部だけ表示しますのでおっしゃる...
page(base) Invoke a Pager on an R o...
ObjecttaskCallbackManager(base) Create an R-level task c...
managergrid.newpage(grid) Move to a New Page on a ...
Devicetkpager(tcltk) Page file using Tk text ...
-なるほど！Pagerというんでしたか。 早速のご回答有り難う...
**ページのヘッダやフッタ描画、またはページの４隅のユーザ...
>[[Я・・R]] (2004-03-23 (火) 18:38:15)~
~
レイアウト、グラフ、余白設定などに関係なく､ページの特定位...
ヘッダやフッタを text(ユーザ座標) で描こうと思ったのです...
ところが~
~
box(which="outer", lwd=4, lty=2, col="blue")~
~
とやると必ず正しく４隅を囲むので、どうやって４隅の座標を...
~
?box~
~
としましたがソースコードが見れませんでした。~
ヘッダやフッタ描画、またはページの４隅のユーザ座標を得る...
//
-identify、locator 関数でインタラクティブに指定位置の座標...
-はい、ページ数が非常に多いので、マウス入力はちょっと無理...
-フリーソフトですから、ソースは(その気にさえなれば)見るの...
-一応解決しました。１ページ分を作図した後に par(new=T, ce...
**ESSの中の日本語 [#o2747e03]
> (2004-03-23 (火) 11:33:40)~
~
gnome-terminal上では~
> a <- "日本語"
> a
[1] "日本語"
と日本語が使えますが、ESS上では~
#ref(ess-20040323.png)~
という感じになります。gnome-ternimal上では問題がないのでE...
何とかなればと思っています。ご助言などがございましたらよ...
[環境]~
VineLinux2.6~
ess-emacs-5.2.0beta1-1vl1~
R-1.8.1-4vl1(本家版)~
//
-[[なかま様のご助言:http://www.okada.jp.org/RWiki/index.p...
**ファイル書き出し時の文字コードの指定方法 [#n01ba267]
>[[matsu]] (2004-03-18 (木) 19:06:00)~
~
はじめまして、matsuと申します。~
~
早速ですが、write()やwrite.table()でオブジェクトをファイル~
に書き出す場合、ファイルの文字コードを制御することはでき...
でしょうか。~
~
私の環境( R ver1.6.1 Win2k )では sjis(改行コードCRLF)で~
書き出されるのですが R側の設定でeuc(改行コードLF)に変更~
はできるものなのでしょうか。~
~
ポインタ等でもお教えいただければ幸いです。~
よろしくお願い致します。~
//
-オリジナルの R は当然日本語コードなど知るはずもないから...
-中からperl（WindowsならActivePerl）を呼び出すことにしてj...
-中澤コメントに補足。例えば write 関数には file="|cmd" の...
-みなさまありがとうございます。ソースをEUCで記述したり、w...
-以前、save(file="hoge.txt",ascii=T,list=hoge)ですと、asc...
**グラフィックウインドウの名前 [#w3dd2d3f]
> (2004-03-18 (木) 11:20:43)~
~
windows()などで表示されたウインドウ~
に自動的にふられる名称~
「R Graphics: Device x」~
を変更することは可能でしょうか？~
//
-Rソースコードのsrc/modules/X11/devX11.c の中で、Xライブ...
-そうですか。作成グラフ数が多いときにできれば便利です。上...
**非線形回帰の初期値 [#c9768ff8]
>[[ぶんしょ]] (2004-03-17 (水) 17:40:55)~
~
非線形回帰の初期値がすぐ入らないのですが、いい方法があれ...
今やっているのが~
tpot<-nls(HAV~(a*SI+b)*(1-exp(k*AGE))^m,start=list(a=1,b...
という式です。
見当をつけて初期値を入れれば大体はすぐ当てはまるのですが...
Error in numericDeriv(form[[3]], names(ind), env) :
Missing value or an Infinity produced when evaluating th...
や
Error in nls(HAV ~ (a * SI + b) * (1 - exp(k * AGE))^m, ...
singular gradient
や
Error in nls(HAV ~ (a * SI + b) * (1 - exp(k * AGE))^m, ...
step factor 0.000488281 reduced below `minFactor' of 0.0...
といったエラー文が出ます。これらの意味もできれば教えてく...
//
-最適化の初期値の選び方についてはチップス [[関数の最大・...
-四つの変数を適当な範囲と、刻目で変えながら、非線形誤差自...
-ありがとうございます。やはり神頼み的な思考錯誤しかないよ...
-誤差の入り方を無視し log(HAV) ~ log(a*SI+b)+m*log(1-exp(...
-すみません勉強不足なんですが、これは対数で回帰させるので...
-そうです。単なる思いつきです。目的変数を log(HAV) にとり...
-何度もすみません。対数で試したら~
Error in numericDeriv(form[[3]], names(ind), env) :
Missing value or an Infinity produced when evalua...
In addition: Warning messages:
1: NaNs produced in: log(x)
2: NaNs produced in: log(x)
といったエラー文が出るのですが何がいけないのでしょうか？ ...
-これは当然 1-exp(k*AGE)) が試行中に非正値になったという...
-もし差し支えなければ、そしてデータがあまり大量でなければ...
-
AGE HAV SI
18 197.98 21.48
20 266.02 23.78
24 695.26 39.88
26 429.21 25.14
27 312.04 22.79
33 484.55 24.15
35 708.11 30.79
25 743.29 35.23
28 714.40 32.50
29 377.93 19.36
29 684.48 29.50
42 557.06 20.37
47 483.70 16.62
50 725.47 20.88
24 400.58 33.11
25 458.79 32.69
25 296.57 24.12
27 382.04 27.50
32 343.38 21.45
33 450.72 26.33
33 375.25 22.31
34 420.39 21.82
36 384.44 20.51
が当てはまらない一例なんですが、初期値a,ｍ=正、ｂ,ｋ=負の...
-- [[ぶんしょ]] &new{2004-03-19 (金) 09:06:58};
-少し解析をしてみましたが、やはり最初に注釈があったように...
> optim(c(0.1, 0.1, -0.1, 0.1), SF)
$par
[1] 10.907251 -1.975936 -0.057869 1.013191
$value
[1] 403297
$counts
function gradient
313 NA
$convergence
[1] 0
$message
NULL
-わざわざありがとうございます。確かに粗悪なデータです。自...
**ｌｍ（）関数の不思議 [#y49c88e7]
>[[杉本]] (2004-03-17 (水) 16:31:34)~
x<-c(8,3,4,2,1,4,8,4,9)
y<-c(4,1,5,1,1,2,5,2,3)
として x=a*y+bの形で線形回帰したいときは~
lm(x~y)
ですが x=a*y^2+bの形で線形回帰するときはどうすればよいの...
lm(x~y^2)
とすると上と同じ結果になってしまいます。もちろん
x=a*y+b*y^2+cの形で線形回帰するときも~
lm(x~y+y^2)
Call:
lm(formula = x ~ y + y^2)
Coefficients:
(Intercept) y
1.586 1.197
となり、x=a*y+bの形で線形回帰してしまっています。新たにデ...
//
-[[Rの統計解析関数Tips]] を見て下さい。R のモデル公式では...
-ありがとうございます。うまくいきました。 -- [[杉本]] &ne...
**無題 [#r9128296]
>[[大熊 大郎]] (2004-03-12 (金) 23:58:36)~
~
関数 width.SJ はなんですか？おしえてください。~
//
-悪い質問の典型ですね。標準パッケージ以外のパッケージ中の...
**ヒストグラムを描くhist() [#eb449b55]
>[[Rin]] (2004-03-04 (木) 11:57:28)~
~
hist()を使いヒストグラムを描くとき、matrixデータを与える...
よろしくおねがいします。
//
-意味がいまひとつ理解できません。具体例をあげていただけま...
-上の方が仰るとおり，hist()では多次元データはベクトルに強...
> x <- matrix(1:24,6,4)
> barplot(x)
> barplot(x,col=rainbow(24))
> barplot(x,col=rainbow(12))
> barplot(x,col=rainbow(6))
> barplot(x,col=rainbow(3))
-理解しにくい質問の仕方をして申し訳ありません。凡例はどの...
-[[ここ:http://www.okada.jp.org/RWiki/?cmd=read&page=%A5%...
-はい。そうです。早速、legend()を使い試してみましたが、や...
- histogram につける気の利いた凡例が思いつかなかったのが...
> data(islands)
> hist(islands)
> legend( 5000,y=30,paste("data",c(1:3)), col = c(1:3), ...
> legend(10000,y=30,paste("data",c(1:3)), col = c(1:3), ...
> legend(10000,y=40,paste("data",c(1:3)), col = c(1:3), ...
#ref(hist.png,center)
//-凡例をつけることができました！やっと解決できて嬉しいで...
//&new{2004-03-09 (火) 17:26:53};
-凡例をつけることができました！やっと解決できて嬉しいです...
-いえいえ．?legend や help(legend) でさらなる情報が得られ...
-初めまして。matlab初心者の者です。棒グラフで値を表示した...
-続き）matlabでは棒１本１本の色をかえることができずその方...
**帯グラフ [#mc629b65]
> (2004-03-02 (火) 12:57:28)~
~
グラフィックス参考実例集には、なぜか帯グラフがありません...
//
-リンク集、グラフィックス参考実例集にもリンクが張られてい...
> data(HairEyeColor)
> x <- apply(HairEyeColor, 2, sum)
> par(fin=c(6,2)) # グラフィックス画面を 6x4 インチに指定
> barplot(as.matrix(x), horiz=TRUE, col=rainbow(length(x...
#ref(obigraph.jpg,left)
-ありがとうございます。私のイメージとはちょっと異なるので...
> data(HairEyeColor)
> a <- as.table( apply(HairEyeColor, c(1,2), sum) )
> b <- a / apply(a, 1, sum)
> barplot(t(b))
#ref(obi.png,left)
> d.seg.xs <- rep((0.2 + 1) * 1:(ncol(t(b)) - 1), each=n...
> d.seg.xe <- rep((0.2 + 1) * 2:ncol(t(b)) - 1, each=nro...
> d.seg.ys <- apply(t(b)[,1:(ncol(t(b)) - 1)], 2, cumsum)
> d.seg.ye <- apply(t(b)[,2:ncol(t(b))], 2, cumsum)
> segments(d.seg.xs, d.seg.ys, d.seg.xe, d.seg.ye)
#ref(obi2.png,left)
-おぉ、horiz=TRUEなんてあるのですか。早速、使ってみました...
> barplot(t(b),horiz=TRUE)
> d.seg.ys <- rep((0.2 + 1) * 1:(ncol(t(b)) - 1), each=n...
> d.seg.ye <- rep((0.2 + 1) * 2:ncol(t(b)) - 1, each=nro...
> d.seg.xs <- apply(t(b)[,1:(ncol(t(b)) - 1)], 2, cumsum)
> d.seg.xe <- apply(t(b)[,2:ncol(t(b))], 2, cumsum)
> segments(d.seg.xs, d.seg.ys, d.seg.xe, d.seg.ye)
#ref(obi3.png,left)
-あとは、0.0-1.0ではなくてパーセント表示にしたいのですが...
-% 表示はやはり手作業技でしょうね。ところで元記事質問者の...
-名乗らなくてすみません。元記事質問者です。上の3つの帯グ...
-barplot()の中でaxes=Fとしておいて，barplot()後にaxis(1,s...
-ありがとうございます！なるほど。%が欲しかったので、axis(...
-質問者の分際でなんなのですが、しかも自力ではないのですが...
**Rが動作するOS [#lcec4b52]
> (2004-03-01 (月) 18:58:11)~
~
Rの紹介記事にRが動作するOS一覧を載せたいのですがどこかに...
現在の把握状況は下記の通りです。~
- MacOS X、Windows、Linux、FreeBSD、NetBSD、OpenBSD、Sola...
- 古いMacOSは古いRなら動く？~
- Linux版のザウルスは、制限がありながらも動作する。~
~
上記以外にRが動作したOSがありますか？~
LindowsやBeOS上でRは動作しますか？超漢字はどうですか？~
//
- 2002年2月の r-help 記事に WindowsCE ではビルドできない...
-Mac OS 9 以下では，R 1.7.1 までしか動きません。 -- &new...
**ESSでのコメント行の色 [#m589a9ba]
>[[さとう]] (2004-02-28 (土) 09:46:19)~
~
Emacs+ESS+Rを利用してます。emacsでRファイルを編集するとき...
とりあえずの対処療法として、~
M-x hi-lock-mode
M-x highlight-regexp [RET] ### [RET] hi-red-b
M-x highlight-regexp [RET] ## [RET] hi-green-b
などとしてます。(順序が大事)~
拡張子がR(またはr)のファイルを開いたときには、いつもこの...
//
**Rで並列処理計算（グリッド・コンピューティング？） [#i53...
> (2004-02-23 (月) 15:18:03)~
~
Rで利用できる並列処理のチュートリアルのようなものはある...
//
-例えば R news 2002/2 号に "Rmpi: Parallel Statistical Co...
**６４ビット環境でのRの動作 [#w32f2127]
> (2004-02-23 (月) 15:16:40)~
~
最近、６４ビットのPCが盛んに話題になっていますが、Rは３...
//
-C言語のdoubleは現在64bitであることが多いですから、double...
~
**階級間隔が異なるヒストグラム [#tf77c6ee]
>[[さかもと]] (2004-02-20 (金) 11:30:56)~
~
Rにて階級区間が異なるヒストグラムを作成したいのですが，ど...
(度数)*d(基準とする標準間隔とでも言いますか･･･)/f~
を間隔ｆの高さにすればいいのですが，どうすればRでそれが可...
ちなみに，初心者です･･･．~
//
//-タイトルの"感覚"=>"間隔"です(-_-;) -- [[さかもと]] &ne...
-引数 breaks (vector型)に、binの区切りとなる値を与えれば...
-ありがとうございます。さっそく試してみます。小生が試した...
-確か引数 width は S の頃に(関数density()などで)使われて...
> x <- rnorm(50) # データ
> h <- 0.8
> bins <- seq(min(x)-0.1, max(x)+0.1+h, by=h)
> pointlist <- hist(x, breaks=bins, prob=T, plot=F)$counts
> hist(x, breaks=bins) # 普通にヒストグラムを描く場合(区...
リストpointlistには，長さ h ごとに区切られた場合の左から...
**免疫アルゴリズム [#h86a6ecd]
> (2004-02-17 (火) 18:01:23)~
~
R や S-plus で免疫アルゴリズム（IA）を扱った例はあるので...
//
**ライブラリのバージョン表示 [#c441c04d]
> (2004-02-17 (火) 13:16:00)~
~
ライブラリパッケージが最新バージョンかどうかを確認するた...
しているライブラリのバージョンを表示させたいのですが、方...
//
-update.packages() で解決しました。 -- &new{2004-02-17 (...
-read.dcf(file=system.file("DESCRIPTION",package="tree"))...
**Rmapやmaptoolsでの縮尺バーや方位記号の表示 [#w57d3f4c]
> (2004-02-17 (火) 00:17:51)~
~
Rmapやmaptoolsで表示した地図には縮尺バーや方位記号を入れ...
//
-arrows などを使って自作するしかないのでは -- &new{2004-...
-凡例はデモにハワイの火山表示のコードが参考になるかも -- ...
**plotで縦軸と横軸の単位幅を同じ長さにする方法 [#ia0b6a17]
> (2004-02-16 (月) 19:10:03)~
~
FAQでしたら、削除しください。~
x軸とy軸がどちらも単位がm(メートル)であるplotで、x軸の1m...
長さにしたいのですが、勝手に調整されてそろいません。x軸と...
同じにする方法はありませんか。~
//
-xlim と ylim で座標軸範囲を同じにすれば良い？-- &new{20...
plot(runif(10), 2*runif(10),xlim=c(0,2),ylim=c(0,2))
-ありがとうございました。確かに、座標範囲を指定すると軸と...
-サンプルコードを紹介して頂かないと、コメントのしようがあ...
-例えば、下のような場合では、上に大きな余白ができます。 -...
plot(2*runif(10), runif(10,min=0,max=.1),xlim=c(0,2),yli...
-x,y 座標を同じスケールにする以上これは避け難い。グラフ領...
-申し訳ありません。最初にいうべきでしたが、x,y 座標が同じ...
-par(mai=c(1,2,6,2))ってしてみるとか，postscriptデバイス...
-dev.copy(x11,whdth=値,height=値)で書き直させる? -- &new...
-widthやheightを指定する方法は、作図領域の指定ではなく、...
-library(grid)を使えばできそうな気がします -- &new{2004-...
-正確な指定の仕方はわかりませんが、次のようにするとそれら...
> par(fin=c(6,3))
> plot(6*runif(100),1.5*runif(100),xlim=c(0,6),ylim=c(0...
比例定数 3/6 というのは本質的ではなく、実際の画面を見て適...
oldpar <- par(no.readonly = TRUE) # 現在のグラフィックス...
... 作業 ...
par(oldpar) # パラメータ復帰
を実行しないと後で困るかも知れません。
#ref(par.fin.jpg, left)
-追加。これでも実際は余計な余白がたっぷり付いている可能性...
-大体様子がわかったような気がします。グラフィックス画面の...
par()$fin[1]*(par()$plt[2]-par()$plt[1]) # x 方向のサイ...
par()$fin[2]*(par()$plt[4]-par()$plt[3]) # y 方向のサイ...
と決まります。従って x,y 座標のスケールを同じにするには x...
> par(fin=c(6,4)) # グラフィックス画面を 6x4 インチに...
> par(plt=c(0.1,0.9,0.1,0.9)) # 作図領域をグラフィックス...
> Xlim <- 10*par()$fin[1]*c(0, par()$plt[2]-par()$plt[1])
> Ylim <- 10*par()$fin[2]*c(0, par()$plt[4]-par()$plt[3])
> plot(1:20,1:20, xlim=Xlim, ylim=Ylim)
#ref(finplt.jpg, left)
参考までにメモすれば plt, fin の既定値は次のようになって...
> par()$plt
[1] 0.1172537 0.9399432 0.1457143 0.8828571
> par()$fin
[1] 6.990803 6.997416
従って既定値での作図画面のサイズはインチ単位でそれぞれ 5....
-お陰様でうまくいきました。ありがとうございました。 -- [[...
-その後、ただ単にasp=1というオプションを付ければよいこと...
例）
> x <- 1:10
> y <- x * 10
> plot(x,y,asp=1)
結果)
#ref(asp.png,center)~ -- &new{2004-05-14 (金) 23:14:18};
**正準相関分析での第 i 正準変量の求め方 [#r7a06ebf]
> (2004-02-14 (土) 20:55:31)~
~
library(mva)の関数cancor()を用いれば正準相関係数行列が得...
第 i 正準変量を求める場合はどうすればよいのでしょう？例え...
[[このような分析:http://bstat.f7.ems.okayama-u.ac.jp/stat...
（exampleでは x %*% xcoef と計算していますが何か違う気が...
//
-一般化固有値問題を解く短い関数を書いて，[[正準相関分析を...
~
**latex形式のヘルプファイルの見方 [#c96fa4c7]
> (2004-02-06 (金) 23:25:05)~
~
libraryごとにhtmlやhelpなどヘルプ文書がありますが、latex...
//
-ソースから R システムを構築するとき make install-dvi と...
-もしかすると [[森さんのWeb「R の日本語文章 (pdf 版)」:ht...
-全く無知で申し訳ないのですが、R CMD Rd2dvi パッケージ名....
-例. library/nlme/latex/nlme.tex だけを make する手順 (R ...
-なお, 上記の手順は, bin/Rd2dvi の内部で実行されているこ...
-わざわざ latex 形式のヘルプファイルを一揃い用意している...
-謎が解けました。やはり help 関数には tex 形式のヘルプフ...
> options(latexcmd="/usr/bin/latex")
> options(dvipscmd="/usr/bin/dvips")
> help(optim, offline=TRUE)
Saving help page to 'optim.ps' # optim.tex から optim.p...
思いついて次のようにしたら dvi ファイルのオンライン表示が...
> options(dvipscmd="/usr/bin/xdvi")
-おぉ、すごいです。無事にxdviでコマンドのhelpが読めました...
-さらにdvips を dvipdfm に変えてみました。すると無事にpdf...
-パッケージのhelpと各々のコマンドのhelpをxdviで表示する方...
-課題としては面白いですが、まず気になるのは、例えば base ...
-help() に offline という引数があったとは．勉強させてもら...
-dvipscmd をだましてやると、dvi ファイルも保存できます。...
cp ./*.dvi /tmp
というスクリプトファイル(名前を仮に makedvi )に実行属性を...
> options(dvipscmd="./makedvi")
としてやる。もう少しスクリプトを工夫する余地？
**文字列のオブジェクトへの変換方法 (なんでも掲示版より移...
>[[sakamoto]] (2004-01-30 (金) 19:35:15)~
~
まずは次のプログラムを見てください。
money<-1:10
>a<-readline()
money
>money[a>5]
[1] 1 2 3 4 5 6 7 8 9 10
> a<-as.name(a)
> a
money
> money[a>5]
[1] 1 2 3 4 5 6 7 8 9 10
ここでしたい事は、以下の事である。
>money[a>5]
[1] 6 7 8 9 10
つまり、文字列として取り込んだものをオブジェクトとして扱...
-ご質問の趣旨はこういうことですね。すでに a という名前の...
> a <- runif(10) # a という名前の数値ベクトルを作る
> a
[1] 0.6313190 0.4538556 0.3934363 0.4643964 0.3375608 0...
> aa <- readline() # その名前をコンマンドラインから入力
a # 名前は a
> aa
[1] "a" # aa は文字列 "a"
> get(aa) # aa という文字列を名前に持つオブジ...
[1] 0.6313190 0.4538556 0.3934363 0.4643964 0.3375608 0...
まとめてみると
> a <- runif(10)
> b <- runif(10)
> x <- 1:10
> x[get(readline())<0.5]
a # コンマンドラインから a を入...
[1] 3 4 5 6 10 # x[a<0.5] が得られた
> x[get(readline())<0.5]
b # コマンドラインから b を入力
[1] 1 3 4 5 6 7 8 # x[b<0.5] が得られた
-money[eval(parse(text=a))>5] としたいと言ってるのかな？ ...
-中間さんの解を注釈すると -- &new{2004-01-31 (土) 00:05:...
> a
[1] 0.557755694 0.756125389 0.062838319 0.116221907 0.28...
> y <- readline()
a
> y
[1] "a"
> eval(parse(text=y)) # 文字列 y をオブジェクト名 a に変...
[1] 0.557755694 0.756125389 0.062838319 0.116221907 0.28...
-ご回答ありがとうございます。おかげさまで問題が解決しまし...
**hist()の引数breaksの決め方 [#hc0fa05f]
> (2004-01-29 (木) 17:57:01)~
~
ヒストグラムを描くhist()の引数のbreaks=seq(min,max,by=bin...
//
-データ数が大きい場合(50以上)はあまり気にしなくても良いは...
-代替案的なものですが，ASH という方法があります．ただ，こ...
**ベクトルと行列 [#c6999be7]
>[[Soichi]] (2004-01-29 (木) 05:39:56)~
~
Matlabから乗り換えようかとRを試し始めました。~
ベクトルと行列のどちらも引数として受け取ることができ、~
行列を受け取ったときは各列に処理を、~
ベクトルを受け取ったときは(n,1)行列に変換した上で、その後...
Rでそのようなことは可能でしょうか。~
//
- as.matrix()をかませておけばいいのでは？ -- &new{2004-0...
- 次のようにすれば良い。問題によりますがいちいち行列に変...
if(is.vector(x)) x <- as.matrix(x) #もしくは
if(!is.matrix(x)) x <- as.matrix(x)
ちなみに縦ベクトルにしたければ
if(is.vector(x)) x <- t(as.matrix(x))
判断する必要がないのなら単に
x <- as.matrix(x) # x がすでに行列なら何もしない
**R の sleep() 機能 [#l2cf0b29]
>[[φ(oo-)o]] (2004-01-28 (水) 09:11:25)~
~
このサイトの「R プログラミング Tips 大全」のページで~
> 指定秒数プログラムをアイドリングする関数~
は標準関数の Sys.sleep(sec) ではどうでしょうか。ヘルプには~
> This function may not be implemented on all systems.~
と書かれていますが、Windows2000では動くようです。~
(R 1.8.1を使用しています)~
//
-本当ですね。知りませんでした。「R プログラミング Tips 大...
**作ったページの削除は？ [#ld06b03b]
>[[mnyu]] (2004-01-23 (金) 01:08:50)~
~
作ったページの削除は どうしますか~
//
-編集モードで内容をすべて削除し、更新すると一人でになくな...
**非線形には対応していますか [#jc5b6ba7]
>[[tia]] (2004-01-22 (木) 17:22:56)~
~
Rに関してまったくの初心者です．~
Rで非線形問題を解くことは可能でしょうか．~
ご教授いただければ幸いです．~
//
-問題によりますが、汎用最適化関数 optim や nlm を使えばた...
-早速のお返事ありがとうございます．いろいろと試してみます...
-非線形最小自乗法用関数 nls については [[Rの統計解析関数T...
**CCDFについて。 [#m5bd524c]
>[[sato]] (2004-01-22 (木) 14:24:52)~
~
Rを使って、1‐累積密度分布であるCCDFを描きたいのですが、パ...
またもしない場合は、パッケージstepfunの中のecdfを使って描...
//
-stepfun パッケージの ecdf では困る、もしくは避けたい理由...
-ecdfで困るということではないのですが、それ用に用意されて...
-無いようですね。累積和関数 cumsum を使えば自作は十分可能...
-わかりました。ありがとうございました。 -- [[sato]] &new{...
**Rでビット演算 [#je04e68d]
> (2004-01-19 (月) 16:29:43)~
~
Rでビット演算をするには、どうすれば、良いのでしょう~
//
-R にはそのものずばりのパッケージはないようです。具体的に...
-Colorオブジェクト同士のビット演算をしたいのですが。 -- ...
-
> heat.colors(9)
[1] "#FF0000" "#FF2A00" "#FF5500" "#FF8000" "#FFAA00" "#...
[8] "#FFFF40" "#FFFFBF"
のようなカラーを与える16進数(を表す文字列)間の and, or, x...
**直線回帰の信頼区間について [#z762ae44]
>[[ぶんしょ]] (2004-01-18 (日) 20:38:47)~
~
直線回帰の95％信頼区間をRで計算したいのですがどうしたらよ...
//
-predict()関数を使うのが便利です。直線回帰でいいなら，exa...
-ありがとうございます。さっそく試してみます。 -- [[ぶんし...
-プロットまでやる例をhttp://phi.ypu.jp/swtips/sumfunc.htm...
-久しぶりに見ました。丁寧にありがとうございます -- [[ぶん...
**スレッド表示 [#b83a13c4]
> (2004-01-13 (火) 21:19:06)~
~
Q & Aではできないのでしょうか？~
//
**[質問] lm() 関数の返り値のF統計量のp-value <= なんでも...
>[[舟尾]] (2004-01-10 (土) 21:12:33)~
~
lm()関数で重回帰分析を行った際，~
~
result <- summary(lm(y~x))~
~
として result［［k］］ や result$names (k:整数，names:名...
をすることで推定値やt値などが取り出せますが，F-statistics~
のp-valueだけがどうやっても取り出せません．names()で調べ...
F-statisticsのp-valueにあたるものが見当たらないのですが，~
この件に関して情報をお持ちの方，レス頂ければ幸いです．~
//
> [[間瀬]] &new{2004-01-10 (土) 22:52:45};
~
たしかに summary.lm 関数は F 統計量の p 値自体は計算しな...
> x=1:10
> y=c(rnorm(5),rnorm(5,mean=1))
> z=summary(lm(y~x))
> w=z$fstatistic
> pf(w[1], w[2], w[3], lower.tail=FALSE)
value
0.1643454
//
> [[中澤]] &new{2004-1-12（月）11:42:00};
~
print.summary.lmの中のようです。プロンプトでprint.summary...
if (!is.null(x$fstatistic)) {
cat("Multiple R-Squared:", formatC(x$r.squared, d...
cat(", Adjusted R-squared:", formatC(x$adj.r.squ...
digits = digits), "?nF-statistic:", formatC(x...
digits = digits), "on", x$fstatistic[2], "and...
"DF, p-value:", format.pval(pf(x$fstatistic[...
x$fstatistic[2], x$fstatistic[3], lower.t...
digits = digits), "?n")
}
となっているのがわかります。
-なるほど。気づきませんでした。R は部品化されているのでこ...
-間瀬さん、中澤さん、非常に丁寧な回答をありがとうございま...
**ファイルが読み込めません [#ca5b1182]
>[[haru]] (2004-01-09 (金) 14:52:09)~
~
初心者なのですが、R1.8.1をダウンロードして早速Excelでデー...
> customernumber <- read.table("costomernumber.txt")~
としましたが~
Error in file(file, "r") : unable to open connection~
In addition: Warning message: ~
cannot open file `costomertable.txt'~
というエラーメッセージがでてしまい読み込むことができませ...
read.delimなども使ってみましたが同じ結果でした。~
http://phi.ypu.jp/swtips/Data2R.htmlに書いてある説明通り...
Rに関する質問なのかもよくわかりませんが、よろしくお願いし...
//
-原因は，ファイル名の指定が正しくないことだと思われます。...
-できました！ご説明本当にありがとうございました！ -- [[ha...
**R初心者 [#y11bdb10]
>[[文字の読み込み]] (2003-12-16 (火) 01:33:52)~
~
文字列（例えば3030A6-0010など）の入ったファイルから各々の...
stnid <- cdataとしても、Error in scan("worldstnid.dat") :...
//
-scan("worldstnid.dat", what=character())とwhat=を引数で...
-scan にデータがどういうタイプか教える必要があります。既...
**macでRのproxy設定 [#kece459e]
>[[mac]] (2003-12-14 (日) 05:55:22)~
~
OSXのパンサーを使っています。~
Rでproxyを利用したいのですが、どうすればいいのでしょうか。~
教えてください。~
//
-すでに試されたかもしれませんが、googleでproxy site:r-pro...
-切れてしました。続き。proxyの設定方法が書いてありました...
-ありがとうございました。出来ました。他にも困っている人が...
-まず、インストールした直後に、Sys.putenv("http_proxy"="h...
-ftpも同様に入れてください。Rをインストール直後にこれを行...
-環境変数はプログラムの起動時ごとに設定されるので、インス...
-また、毎回起動時に環境変数を設定する必要をなくすには、/A...
-あと、インストールしたパッケージとかは ~/Libaray/RAqua/ ...
**同時実行 [#n213d61b]
> (2003-12-12 (金) 16:58:27)~
~
Windows用のRを2つ立ち上げて、~
それぞれを操作した場合に、~
それぞれの実行結果は完全に独立だと思ってよいのでしょうか？~
//
-独立とは統計的な独立（乱数を使用した場合に限り意味のある...
**グラフがかけない X11 is not available [#od5d3d06]
>[[TAKU]] (2003-12-12 (金) 11:33:19)~
~
超初心者で、はじめてLinuxをインストールし、Rをemacs上で立...
とこまでできたのですが、postscript以外のグラフィックコマ...
使えないのです。png("hoge")やx11()どを入力すると~
X11 is not available~
とエラーメッセージがでてしまいます。Rの質問ではないのかも...
助け舟をよろしくお願いします。~
ディストリビューションはdebianです。~
//
-X Window Systemは起動しているんですね？それならば X11(di...
-Error in x11(display = ";0") : X11 is not available と...
-私もdebianです。同じことがおこったことがあります。もしか...
-で、debianのX11クライアント開発環境のインストールですが...
-もしインストールされていなければ、apt-get -s install xli...
-開発環境がインストールできたら、Rのsourceに戻って、./con...
-私は試していませんが、make が面倒なら、www.r-project.org...
-detourさんほんとうにいろいろありがとうございます！（;_;)...
-ありがとうございました！解決しました。このことでずっと迷...
-いえいえ、Knoppix-jp CDROM から HD へインストールすれば...
-[[テスト中CRANミラー:http://cran.okada.jp.org/bin/linux/...
-ほんとにいろいろありがとうございました。まだまだ初心者で...
**Rプログラミングで得られた結果を別のファイルへ出力するに...
>[[ぽいと]] (2003-12-11 (木) 18:21:16)~
~
現在、R言語でプログラムを書いたのはいいのですが、その結果...
//
-プログラムの冒頭でsink("c:/temp.txt")のようにすれば、tem...
-Tips 集の [[R出力の記録]] が参考になります。 -- &new{20...
-早速の回答をありがとうございます。投稿後、色々調べてみた...
**Rパッケージのインストール＆ロード方法 [#i0c435fe]
>[[R初心者]] (2003-12-11 (木) 16:03:04)~
~
現在Linux上でパッケージ「evd」をinstall.package("evd")で...
//
-該当パッケージをダウンロードした上、ルートになって R CMD...
**大規模データからランダムにデータをサンプリングする [#ma...
>[[msasalo]] (2003-12-10 (水) 21:48:07)~
~
はじめまして。最近Rを使いはじめるようになった者です。~
200万件38Mbのデータを処理しようと思っていたのですが、読み...
そこで、データのサブセットを作って36万件6Mbのデータを作っ...
マシンのスペックは、P4-1.6GHz,512Mで、ページングメモリー...
データファイルを元にランダムにとる方法がRにあるとありがた...
//
- sample 関数が探しているものでは？-- &new{2003-12-10 (...
> x=rnorm(2000000) # 200万個の正規乱数発生
> length(x)
[1] 2000000
> 2000000*8/1024/1024
[1] 15.25879 # 15 Mbyte
> xx = sample(x, 10000) # その中から一万個を無作為抽出
> length(xx)
[1] 10000
//
-Rに読み込んであるデータならばsample()という関数が使えま...
-そうなのです、以前はPostgreSQLを立ち上げてデータをODBCと...
-ありがとうございました!! なんとか解決しました!! Rで216万...
**[質問]2変量t乱数の作り方 [#r417af38]
> (2003-12-10 (水) 16:59:27)~
~
2変量乱数の発生でt乱数の作り方を教えて下さい。
//
-CRAN にあるアドオンパッケージ mvtnorm は多変量の正規分布...
アルゴリズムを書く代わりに命題みたいなものを示します。
- R：自由度mのカイ二乗分布に従う確率変数
- Z：p変量正規分布N(0,Ip)に従う確率ベクトル（独立な標準正...
- V：正定値対称行列（ちらばりを表す行列で分散共分散行列と...
- V＝C%*%t(C)：コレスキー分解
-- X = sqrt(m/R)*Z はp変量楕円t分布met(p,m,0,Ip)に従います
-- Y=μ+C%*%X はp変量楕円t分布met(p,m,μ,V)に従います
拡張性があるように2変量ではなく3変量楕円t乱数に従う乱数を...
met3 <- function(m, mu, V, n) {
# m : 自由度
# mu : 平均ベクトル
# V : 散らばり行列
# n : 乱数の個数
U <- svd(V)$u
V1 <- svd(V)$v
D <- diag(sqrt(svd(V)$d))
B <- U %*% D %% t(V1)
w <- c()
for (i in 1:n) {
R <- 0
for (j in 1:m) {
R <- R + rnorm(1)^2
}
w <- append(w, list(mu + B %*% (cbind(rnorm(3))*sqrt...
}
return(w)
}
# 実行例
mu <- cbind(c(1,1,1))
V <- array(c(2,1,1,1,2,1,1,1,2), dim=c(3,3))
n <- 1000
w <- met3(5, mu, V, n)
sm <- 0
for (i in 1:n) { sm <- sm + w[[i]] }
sm <- sm/n
sv <- 0
for (i in 1:n) { sv <- sv + (w[[i]]-sm) %*% t(w[[i]]-sm) }
sv <- sv/n
この分布の平均ベクトルはμ、分散共分散行列はm/(m-2)*V（m>2...
多変量楕円t乱数がうまく出来ているかどうかはこれで確かめる...
**[質問]hist()関数とnclass()関数の関係 [#bdaab073]
> (2003-12-06 (土) 11:09:38)~
~
どちらの関数でもSturgesの方法によって階級が求められるよう...
//
-ご質問の意味が今一不明確ですが、?hist に書いてあるように...
-- (1) 単に hist(x) は hist.default(x, breaks = "Sturges...
-- (2) 明示的に分割点ベクトルを与えても良い、
-- (3) 単に整数 n を一つを与えれば、区間が n 等分される、
-- (4) 分割点を計算するユーザー定義関数を与える
ということでしょうか。 -- &new{2003-12-06 (土) 17:00:57};
//
-具体的には(1)に関してのことになりますが、以下ではsは常に...
x<-rnorm(80,mean=0,sd=1)
s<-nclass.Sturges(x)
b<-length(hist.default(x,breaks = "Sturges")$breaks)
s;b
-たしかに変ですね。関係する関数は内部関数なので詳細はソー...
> test <- function(x) {c(nclass.Sturges(x),nclass.scott...
+ length(hist(x)$break...
> test(rnorm(1000))
[1] 11 21 26 16 16
> test(runif(1000))
[1] 11 10 10 11 11
> test(rnorm(1000))
[1] 11 19 25 15 15
> test(rexp(1000))
[1] 11 20 32 15 15
> test(rcauchy(1000))
[1] 11 109 3057 13 13 # ずいぶんと違うもんですね...
-R の ML に Sturges 方式を hist の既定にするほうが良いか...
- おなじことは他の方式でもおきているようです。一つのあり...
> nclass.scott(x)
[1] 9
> length(hist(x,breaks="scott")$breaks)
[1] 8
> nclass.FD(x)
[1] 18
> length(hist(x,breaks="FD")$breaks)
[1] 14
> nclass.Sturges(x)
[1] 9
> length(hist(x)$breaks)
[1] 8
**R(D)COMの使い方 [#l18c0057]
>[[f]] (2003-12-06 (土) 00:37:19)~
~
R(D)COMを使って、VC++で作ったプログラムにRを組み込みたい...
R(D)COMをお使いの方がおられましたら、使用法等、ご助言お願...
//
**[質問]条件分岐 [#l0e4043d]
>[[さかい]] (2003-12-04 (木) 18:59:26)~
~
i<500のときx=1、500<=i<1000のときy=1、~
それ以外のときz=1とするプログラムを作成したかったので、~
以下のようなプログラムを作成したのですが、~
iが1000以上のときx=0,y=1,Z=0となり、うまく値を返してくれ...
どのようにしたら解決できるのか分かる方がいたらアドバイス...
~
以下のプログラムはRプログラミングTips大全集の条件実行 if...
#プログラム
x <- 0
y <- 0
z <- 0
i <- 2000
if (i < 500) {x <- 1} else if (500 <= i < 1000) {y <- 1...
#結果
> x
[1] 0
> y
[1] 1
> z
[1] 0
//
-disc <- function(i) { x <- 0; y <- 0; z <- 0; if (i<500)...
-２番目のifは１番目のifに対するelseの後なので，500<=iの部...
-ありがとうございます。やってみます。 -- [[さかい]] &new{...
-eval(500 <= 2000 < 1000)はなぜTRUEを返すのですか？ -- [[...
-if (i < 500) {x <- 1} else if (500 <= i && i < 1000) {y...
x <- ifelse(i < 500, 1, 0)
y <- ifelse(i < 1000, 1, 0)
z <- ifelse(i >= 1000, 1, 0)
ついでですが、今チェックをしていて気づきましたが R 1.8 か...
> test <- function(x,y) return(x,y)
> test(1,2)
$x
[1] 1
$y
[1] 2
Warning message:
multi-argument returns are deprecated in: return(x, y)
-- &new{2003-12-04 (木) 21:38:49};
-もっと簡単に(+FALSE は論理値を整数に強制変換するトリック...
> x = (i < 500)+FALSE; y = (i < 1000)+FALSE; z=(i >=100...
> x;y;z
[1] 0
[1] 0
[1] 1
-- &new{2003-12-04 (木) 22:02:08};
**Rdbi RdbiPgSQL [#r0441a96]
>[[なかま]] (2003-12-04 (木) 10:42:45)~
~
http://www.bioconductor.org/ のが 1.0.2 で cran.r-project...
bioconductor のが最新?~
//
-ところで、RdbiPgsql はCRAN にはバイナリがありませんが、...
-[[PostgreSQL Native Win32:http://momjian.postgresql.org/...
-やってみると、面倒くさかったのでVine2.6(PostgresSQL-7.2....
-Rdbiはbioconductorの方にバイナリがあります。 -- [[なかま...
-やっぱりbioconductorのが最新のようですね(他にもあったり...
-なかまさん、RdbiPgSQLのバイナリありがとうございます。わ...
-libpqとpostmasterのバージョンが異なると通信出来なかった...
**［質問］棒グラフで対数軸 [#k47e8d60]
> (2003-11-29 (土) 17:20:27)~
~
棒グラフを作図するのには barplot(1:10) のようにすれば良い...
そこで、Y 軸を対数軸として作図しようと log="Y" を指定した...
> barplot(1:10, log="y")~
Error in plot.window(xlim, ylim, log = "", ...) :~
formal argument "log" matched by multiple actual arguments~
対数軸の棒グラフはどのように指示すれば描くことが出来るの...
よろしくお願いします。~
//
-残念ながら base パッケージの barplot 関数は対数軸をサポ...
-ありがとうございます。アドオンパッケージを調べてみます。...
**subsetで値を取り除く [#p4d4e573]
>[[detour]] (2003-11-28 (金) 14:49:55)~
~
subsetで名義変数の特定の値を除いたサブセットを作成したと...
a<-data.frame(choice=c("ab","ab","ab","c","c","d","d","d"...
summary(a)~
summary(subset(a,choice!="c"))~
"c"というのが一応残っているようです。これはどうやったら消...
//
-次のようにすれば一応ラベル c は消えますが、おそらくこれ...
> b=data.frame(choice=a[a!="c"])
> summary(b)
choice
ab:3
d :3
-使用の目的にもよるのでしょうが、気にしないのもありかと。...
-消せないせいでうまく行かなかったことがありましたが、勘違...
-おまけ:subsetを見付ける前はbquoteを使ってみようと思った...
-気にしない事にしましたがやはりうまく行かない事もありまし...
-で、できました。b<-subset(a,choice!="c")だとすると、leve...
-おめでとうございます。苦労して得られたノウハウは是非然る...
**Rのインストール方法（UNIX） [#ja1a0abf]
>[[R初心者]] (2003-11-28 (金) 14:21:53)~
~
ただいまRをUNIXにインストールしようとしています。R-1.8.1...
> make~
Make: Cannot open /share/make/vars.mk. Stop.~
vars.mkファイルは確かに/share/make/に存在しましたし、原因...
//
-Debian GNU/Linux 3.0ですが、問題なくmakeできました。Unix...
-使っているUNIXマシンはCompaq Tru64です。 -- [[R初心者]] ...
ページ名:
![[RjpWiki]](image/../Rlogo.png "[RjpWiki]")
