なんでも掲示板アーカイブ（３）
をテンプレートにして作成
[
トップ
] [
新規
|
一覧
|
検索
|
最終更新
|
ヘルプ
]
開始行:
COLOR(blue){SIZE(20){RjpWiki自体へのコメントでも、雑談で...
COLOR(red){注意}：新規記事用の入力欄は以下の目次の直後に...
----
COLOR(red){過去の記事のアーカイブ} 。アーカイブ中の記事へ...
-[[なんでも掲示板アーカイブ（１）]] (元記事が 2003-07-17 ...
-[[なんでも掲示板アーカイブ（２）]] (元記事が 2003-11-27 ...
----
#contents
----
回答を期待する記事は原則 COLOR(red){質問コーナー [[Ｑ＆Ａ...
~
#article
**なるほど！ [#ecceab77]
> (2004-06-30 (水) 23:58:28)~
~
r-help 記事でまた嬉しくなる発見をしました。言われてみれば...
test1 <- function(x) 1+(x-2)*2*(x-2<0)+(x-2)*(-2)*(x-2>=...
plot(test, -1, 5)
こうした論理値の使い方を知らなければ次のように定義したで...
test2 <- function(x) ifelse(x < 2, 1+2*(x-2), 1-2*(x-2))
ついでに奇妙なことに気づきました。test1 は二回の論理判断...
> set.seed(222); system.time(test1(runif(1e6)))
[1] 0.38 0.30 0.69 0.00 0.00
> set.seed(222); system.time(test2(runif(1e6)))
[1] 0.72 0.33 1.06 0.00 0.00
//
-計算速度も速くなるんですねぇ！勉強になりました．例として...
f(x) = 0.6φ(x)+0.4ψ(x)
となる密度 f(x) に従う乱数を生成する関数を定義してみます．
generator <- function(n) {
data1 <- rnorm(n)-1 # φ(x) に従う乱数
data2 <- rnorm(n)+2 # ψ(x) に従う乱数
data3 <- (runif(n) <= 0.6) # 確率0.6でdata1を採択し，
data1*data3+data2*(1-data3) # 確率0.4でdata2を採択する...
}
mydata <- generator(1000)
([[ヒストグラムと密度の推定:http://www.okada.jp.org/RWiki...
-R は lisp で動いているのでしたっけ。知能的な最適化が組み...
上の関数は，以前に拝見して，次のような書き換えをしてみた...
> generator <- function(n) {
+ data1 <- rnorm(n)-1 # φ(x) に従う乱数
+ data2 <- rnorm(n)+2 # ψ(x) に従う乱数
+ data3 <- (runif(n) <= 0.6) # 確率0.6でdata1を採択...
+ data1*data3+data2*(1-data3) # 確率0.4でdata2を採択...
+ }
> system.time(mydata <- generator(100000))
[1] 0.24 0.02 0.28 0.00 0.00
> generator2 <- function(n) {
+ (rnorm(n)-1)+3*(runif(n) > 0.6)
+ }
> system.time(mydata2 <- generator2(100000))
[1] 0.15 0.00 0.15 0.00 0.00
ついでに，先の関数も
Rweb:> test1 <- function(x) 1+(x-2)*2*(x-2<0)+(x-2)*(-2)...
Rweb:> system.time(test1(rnorm(1000000)+2))
[1] 2.79 0.13 5.94 0.00 0.00
Rweb:> test2 <- function(x) 1+2*(x-2)*((x<2)-(x>=2))
Rweb:> system.time(test2(rnorm(1000000)+2))
[1] 1.36 0.09 1.45 0.00 0.00
人手で最適化しても，たいした効果はないようにも見えます。 ...
-R は Scheme を手本にしているようですが、Lisp インタプリ...
-す、素晴らしい！鳥肌が立ちました！関数 generator2() は速...
-上の書き換え (rnorm(n)-1)+3*(runif(n) > 0.6) は元のコー...
-私も最初「間違いなのでは？」と思ったのですが、
--確率0.6でφ(x)に従う乱数を、確率0.4でψ(x)に従う乱数を生...
という手順は、
--まず平均-1，分散1の正規乱数（φ(x)に従う乱数）を生成する
--この乱数を、確率0.4で平均2，分散1の正規乱数（ψ(x)に従う...
ことと同じになると思います。変更手順は単に「φ(x)に従う乱...
-「元のコードの二つの rnorm(n) が同じベクトルという間違い...
-なるほど！ 勘違い失礼。お詫びに（？）平均、分散がともに...
# 2つの正規分布の混合分布乱数
# n: 発生する混合正規乱数の総数
# p1: 最初の正規分布の比率
# N1: 最初の正規分布の平均、標準偏差のベクトル c(mu1, sd1)
# N2: 2番目の正規分布の平均、標準偏差のベクトル c(mu2, s...
rmixnorm <- function(n, p1, N1, N2) {
P <- (runif(n) < p1)
(P*N1[2] + (1-P)*N2[2])*rnorm(n) + (P*N1[1] + (1-P)*N...
}
-[[Rコードの最適化例：混合正規乱数の発生コード]] にありが...
-問題の解析
Rweb:> rmixnorm <- function(n, p1, N1, N2) {
+ P <- (runif(n) < p1)
+ (P*N1[2] + (1-P)*N2[2])*rnorm(n) + (P*N1[1] + (1-P...
+ }
Rweb:> mean(sapply(1:100, function(i) system.time(rmixno...
[1] 0.10932
100回やって平均を求めた。~
単に混合分布のデータを作るだけなら，乱数の発生順までラン...
Rweb:> rmixnorm2 <- function(n, p1, N1, N2)
+ {
+ n1 <- floor(n*p1)
+ c(rnorm(n1)*N1[2]+N1[1], rnorm(n-n1)*N2[2]+N2[1])
+ }
Rweb:> mean(sapply(1:100, function(i) system.time(rmixno...
[1] 0.04744
これじゃ，それぞれの分布の個数が決まってしまうじゃないか...
Rweb:> rmixnorm3 <- function(n, p1, N1, N2)
+ {
+ n1 <- sum(runif(n) < p1)
+ c(rnorm(n1)*N1[2]+N1[1], rnorm(n-n1)*N2[2]+N2[1])
+ }
Rweb:> mean(sapply(1:100, function(i) system.time(rmixno...
[1] 0.06088
みたいな。 -- &new{2004-07-03 (土) 18:28:01};
-(なるほど)^3、感服！ -- &new{2004-07-03 (土) 19:22:16};
-ああ，わたしってほんとにうっかりもの。平均値も標準偏差も...
Rweb:> rmixnorm4 <- function(n, p1, N1, N2)
+ {
+ n1 <- sum(runif(n) < p1)
+ c(rnorm(n1, mean=N1[1], sd=N1[2]), rnorm(n-n1, mean...
+ }
Rweb:> mean(sapply(1:100, function(i) system.time(rmixno...
[1] 0.05924
時間は大して変わらないけど，可読性はずいぶん違う。 -- &n...
-コメント -- &new{2004-07-04 (日)00:01:00};
> mean(sapply(1:100, function(i) system.time(rmixnorm4(1...
[1] 0.05618
ではなく(これは 100x5 行列の平均になってしまう)
> mean( sapply( 1:100, function(i) system.time( rmixnorm...
[1] 0.1298
としないとまずい。
-最適化された関数をいろいろ紹介して頂きましたが，「本当に...
generator <- function(n) {
set.seed(7777); data1 <- rnorm(n)-1
set.seed(7777); data2 <- rnorm(n)+2
set.seed(7777); data3 <- (runif(n) <= 0.6)
data1*data3+data2*(1-data3)
}
mydata <- generator(1000) # ...
plot(density(mydata)) # ...
generator2 <- function(n) {
set.seed(7777); data1 <- rnorm(n)-1
set.seed(7777); data1 <- data1 + 3*(runif(n) > 0.6)
}
mydata <- generator2(1000) # ...
plot(density(mydata)) # ...
rmixnorm <- function(n, p1, N1, N2) {
set.seed(7777); P <- (runif(n) < p1)
set.seed(7777)
(P*N1[2] + (1-P)*N2[2])*rnorm(n) + (P*N1[1] + (1-P)*N2...
}
mydata <- rmixnorm(1000,0.6,c(-1,1),c(2,1)) # ...
plot(density(mydata)) # ...
#ref(density.gif, center)
- 上と同じ動作をするような，rmixnorm2() 〜 rmixnorm4() の...
-こういうコーディング例の検証を考えると，「[[RUnit を使う...
-set.seed(7777); rnorm(20) と set.seed(7777); rnorm(15); ...
set.seed(7777); data1 <- rnorm(n)-1
set.seed(7777); data2 <- rnorm(n)+2
は，
set.seed(7777); data1 <- rnorm(n)-1; data2 <- data1+3
と全く同じということではないでしょうか。実際にやってみる...
> n <- 1000
> set.seed(7777); data1 <- rnorm(n)-1
> set.seed(7777); data2 <- rnorm(n)+2
> set.seed(7777); data11 <- rnorm(n)-1; data12 <- data11+3
> sum(data12==data2)
[1] 905
差はどれくらいかというと
[1] -4.440892e-16 -2.220446e-16 -2.220446e-16 2.220446e...
程度 -- &new{2004-07-04 (日) 16:52:40};
-乱数の不思議~
同種類の乱数だと以下のようになるが
> set.seed(1)
> rnorm(10)
[1] -0.6264538 0.1836433 -0.8356286 1.5952808 0.3295...
[9] 0.5757814 -0.3053884
> rnorm(5, mean=5, sd=10)
[1] 20.117812 8.898432 -1.212406 -17.146999 16.249309
この乱数列の最後の5つは
> set.seed(1)
> rnorm(15, mean=5, sd=10)
[1] -1.264538 6.836433 -3.356286 20.952808 8.295...
[9] 10.757814 1.946116 20.117812 8.898432 -1.212...
の，最後の5つと同じになっている。
しかし，
> set.seed(1); rnorm(5); runif(5); rnorm(5, mean=5, sd=10)
[1] -0.6264538 0.1836433 -0.8356286 1.5952808 0.3295078
[1] 0.2059746 0.1765568 0.6870228 0.3841037 0.7698414
[1] 4.942328 29.046534 12.635935 -2.990092 -6.476570
のように途中で，異なった種類の乱数発生を挟むと，違ったも...
-set.seed(n) は指定した正整数 n から乱数の初期値を適当に...
-「というか何度も同じ系列を使うことになりむしろ有害」と書...
-私の言葉足らずでした．すみません…．私がやりたかったこと...
よく考えると，rmixnorm2() 〜 rmixnorm4() は，生成すべき片...
よく考えると，rmixnorm2() 〜 rmixnorm4() は，生成すべき片...
（追記）generator() の data1 と data2 には同じ乱数列を与...
#comment
**日本人による新パッケージ [#vf5e7758]
> (2004-06-30 (水) 08:21:13)~
~
We would like to announce the release of our software, w...
available through CRAN.
MNP: R Package for Fitting the Multinomial Probit Models
Abstract:
MNP is a publicly available R package that fits the Baye...
probit models via Markov chain Monte Carlo. Along with t...
multinomial probit model, it can also fit models with di...
sets for each observation, and complete or partial order...
available alternatives. The computation is based on the ...
data augmentation algorithm that is developed by Imai an...
``A Bayesian Analysis of the Multinomial Probit Model Us...
Augmentation,'' Journal of Econometrics, Forthcoming.
Kosuke Imai, Department of Politics, Princeton University
Jordan R. Vance, Department of Computer Science, Princet...
David A. van Dyk, Department of Statistics, University o...
//
#comment
**R for Bioinformatics [#s76a2aea]
>[[nami]] (2004-06-28 (月) 07:59:24)~
~
オライリー本家から「R for Bioinformatics」が、11月に出る...
[[amazon.co.jp:http://www.amazon.co.jp/exec/obidos/ASIN/0...
日本語化されないかな〜~
//
-オライリージャパンの方に直接聞いたのですが、現時点では和...
-どうもありがとうございました。なんとかがんばっていただき...
#comment
**R の一様乱数 runif() [#ie70b59f]
>[[間瀬茂]] (2004-06-27 (日) 00:35:43)~
~
r-help にオヤッと思う記事がありました。~
> set.seed(222)
> sum(duplicated(runif(1e5)))
[1] 4
つまり、10万個の一様乱数中にタイが４個(実際は２組のタイ)...
//
-32ビットマシンで，8バイト実数だと，表現できる正の実数は ...
-MT法は内部的にはもっとたくさんのビットを使っています。さ...
-周期と相異なる乱数の個数は無関係でしょう。 -- &new{2004...
-Ripley が注意しているのは
> qbirthday(classes=2^(53))
[1] 111743588
> qbirthday(classes=2^(32))
[1] 77163
つまり３２ビット実数一様乱数ではタイが生じる確率が 0.5 以...
-- &new{2004-06-27 (日) 10:59:34};
#comment
**package [#g9c79caf]
> (2004-06-26 (土) 23:33:46)~
~
パッケージquadprogの関数solve.QP()を利用してQPを解く際、~
そのパッケージ例に当てはまらないような目的関数を解くこと...
//
-これは独り言ですか、質問ですか？質問ならば場所が不適当で...
#comment
**関数 getAnywhere() [#g6b38f3c]
>[[間瀬]] (2004-06-24 (木) 21:09:49)~
~
また不思議な関数を発見しました。getAnywhere()。引数で与え...
//
#comment
**R fortunes [#lebbe91d]
> (2004-06-20 (日) 08:25:24)~
~
パッケージ更新案内で見つけ、R fortunes というパッケージを...
//
-SAS や「お間抜けな」r-help 投稿記事が格好の的になってい...
Some of us feel that type III sum of squares and so-call...
statistical nonsense which should have been left in SAS.
-- Brian D. Ripley
s-news (May 1999)
To paraphrase provocatively, `machine learning is statis...
of models and assumptions'.
-- Brian D. Ripley (about the difference between machi...
statistics)
useR! 2004, Vienna (May 2004)
Can one be a good data analyst without being a half-good...
answer to that is, `No.' The long answer to that is, `No.'
-- Frank Harrell
1999 S-PLUS User Conference, New Orleans (October 1...
#comment
**RCC (R compiler) ! [#s8315127]
>[[間瀬茂]] (2004-06-19 (土) 23:51:47)~
~
r-help 記事に驚くべきものを発見しました。R コンパイラー！...
[[参考１:http://www.cs.rice.edu/Colloquia/garvin-04-21.sh...
[[参考２:http://hipersoft.cs.rice.edu/rcc/]]~
//
-gcc, bcc に続いて rcc ですか！また凄まじい題材で修士論文...
-参考ホームページの内容だけからは良くわからないが、R コー...
#comment
**R は S(-plus) 互換？ [#a1655aa3]
>[[間瀬]] (2004-06-13 (日) 07:16:11)~
~
最近、J. Chambers の「データによるプログラミング」森北出...
//
#comment
**ｚｉｐのパスワード [#ja6562ba]
>[[初心者]] (2004-06-10 (木) 23:21:48)~
~
unix上でzip圧縮した際にパスワードを設定したいのですが、~
どなたかご存知な方いらっしゃいいますか？~
windowsだと基本的にzip圧縮の際にパスワードを設定できるよ...
一応manコマンドで調べると-eという暗号化のオプションが存在...
環境はsolarisです。よろしくお願いします。~
//
-せめて R との関連を申し訳でも書かないと顰蹙をかううこと...
$ zip -e abc abc.tex # abc.zip という名...
Enter password: # パスワード入力
Verify password: # もう一度
updating: abc.tex (deflated 67%)
$ ls
./ abc.tex abc.zip
$ unzip abc.zip
Archive: abc.zip
[abc.zip] abc.tex password: # 設定パスワードを...
replace abc.tex? [y]es, [n]o, [A]ll, [N]one, [r]ename: y
inflating: abc.tex
-少し強引にRの話題に振りますが、圧縮アルゴリズムって公開...
-?zip.file.extract もちろん R で書いてあるわけではなく、...
-なるほど、内部にunzipのコードは持っていて、すべてのプラ...
#comment
**R-jp, RjpWiki の来歴について [#u6fe8b45]
>[[間瀬茂]] (2004-06-10 (木) 21:19:54)~
~
明日は RjpWiki の開設一周年記念日です。管理者の岡田さん、...
//
-Rの日本での普及の歴史をそろそろまとめるといいかもしれま...
#comment
**予想したとおりになるでしょうか [#w264e6ae]
> (2004-06-10 (木) 13:27:27)~
~
以下のコード，出力はどうなると思いますか。~
やってみて，予想したとおりになったでしょうか。~
私は，驚きました。~
for (i in 1:10) {
if (i == 3) i <- i+5
print(i)
}
//
-スコープについて学ぶときにとてもよい教材ですね（print(i)...
-if 文の { } の有無は，私が意外と思ったことには関係ない...
先に挙げたコードは，以下のコードと同じことを表しているの...
for (i in 1:10) {
if (i == 3) {
i <- i+5
}
print(i)
}
ついでに，試してみたらこんなことも。
n <- 10
for (i in 1:n) {
if (i == 5) {
n <- 20
}
print(i)
}
どうしても for を使わなくてはならなかった局面で（笑），出...
~
help を見れば，書いてあることなんですけどね。~
The index 'seq' in a 'for' loop is evaluated at the start...
loop; changing it subsequently does not affect the loop.
-なるほど、forループではforループの最初のところに戻ってこ...
-駄目といわれると、掟破りをしたくなるのが悲しい人間のさが...
-次も不思議ではありませんか。ループ変数が最終値で残ってい...
> rm(i); for (i in 1:10) {i^2}; i
[1] 10
-私のループ変数の理解がぐらついてしまって非常に混乱してお...
-あくまでも推測ですが，for の最初に制御変数の集合を作って...
> x <- c(3,2,3,8,5)
> for (i in x) { if (i==2) x<-c(4:6); print(i)}
[1] 3
[1] 2
[1] 3
[1] 8
[1] 5
> x
[1] 4 5 6
> i
[1] 5
-i^2 は副作用なし（i に付値してませんから）です。 -- &ne...
> rm(i); for (i in 1:10) {i<-i^2}; i
[1] 100
こうしてみる？
> rm(i); sum<-0; for (i in 1:10) {sum<-sum+i^2}; print(i...
[1] 10
[1] 385
-C では確か、ループ変数をループ中で変更できませんでしたっ...
-C でも，BASIC でも，FORTRAN でも，AWK でも可能です。とい...
- 少し R をからかってみました。文句をいうが、へたれないと...
> rm(i); for (i in 1:3) {print(i); rm(i)}
Warning message: )}remove: variable "i" was not found
NULL
[1] 1
[1] 2
[1] 3
//-記事の衝突を調整中舟尾さんの記事をうっかり消してしまっ...
//-キャッシュに残っていたので，復帰させてみました。 -- &...
//-三人が集中してたのね -- &new{2004-06-10 (木) 22:39:11};
-in 以降の式を評価してスタックして、ループの間だ 指定した...
# 呪われた言語と言われるゆえん?
@ary = (1,2,3);
for $i ( (@ary,4) ){
push @ary,5 if ( $i == 2 );
print $i,"?n";
}
print join(",",@ary),"?n";
for $i ( "a","b","c" )
{
$i="c";
print $i,"?n";
}
1
2
3
4
1,2,3,5
Modification of a read-only value attempted at a.pl line...
-確かに以下のようなへんなことが出来る以上 「制御変数」と...
> r <- list(1, "a", matrix(1:4,nc=2))
> for (i in r) print(i)
[1] 1
[1] "a"
[,1] [,2]
[1,] 1 3
[2,] 2 4
#comment
**CRAN更新パッケージ一覧の表示方法 [#e2c27c40]
>[[岡田]] (2004-06-04 (金) 17:37:34)~
~
とりあえず、MenuBarの最新30件に混ざって出るようにしてみま...
//
-別ページ表示に一票 -- &new{2004-06-04 (金) 21:40:46};
#comment
**無相関は独立性を意味しない、例 [#l8e28809]
> (2004-06-01 (火) 22:58:00)~
~
R とは直接関係ない話題で済みません。また学識深い R ユーザ...
//
-正規分布のような裾の広い分布で示すより，一様分布のような...
-一様分布でも原点対称（例えば [-0.5, 0.5] 上の一様分布）...
#comment
**「R-Tips」のPDF版とその原稿 (04/06/07に修正版をup) [#y0...
>[[舟尾]] (2004-02-03 (火) 22:00:00)~
~
竹澤先生の御厚意により、竹澤さんのウェブスペースに[[R-Tip...
//
-素晴らしい。２８０ページの力作！もともとこんなに量があり...
-竹澤（中央農業総合研究センター・農業情報研究部）です。 ...
-改めて全体を見ました。感想は「すごい」。内容の豊富さも素...
-非常に恐縮であります…。「Introduction to R」の日本語訳版...
このPDF(+html)をネットに上げるにあたって、管理等を全てお...
-リンク切れでしょうか？PDF を見ることができないのですが・...
-復帰したようですね。 -- &new{2004-02-08 (日) 20:57:41};
-ですね。リンクの張り間違いではなかったようです。 -- &ne...
-ありがとうございます。これはものすごいPDFが出ましたね。 ...
- 2/6 に竹澤先生にお願いして up させてもらった PDF は，煩...
-内容についてではなく、つまらないことで恐縮です。印刷して...
-こちらのレーザープリンタ・インクジェットプリンタでは印刷...
-私も１２ページだけ印刷できませんでした。軽い版でも駄目で...
-↑この書き込みの者です。原因は分かりませんでしたが，解決...
-大変申し訳ありません．お手数をおかけしております…．上の...
-すみません。最初に12頁が印刷できないと書き込みした者です...
-お役に立てなくて申し訳ありませんでした…．また原因が分か...
#comment
**統計グラフで特許出願！ [#t93e7aa8]
>[[間瀬茂]] (2004-05-18 (火) 10:15:03)~
~
r-help 記事によると、三次元棒グラフの代替グラフである dia...
グラフが特許出願中だそうです。~
~
[[実物はここを参照:http://www.jhsph.edu/Press_Room/Press_...
//
-引用されたグラフ中の数字が読みとれないので，ネットサーフ...
https://stat.ethz.ch/pipermail/r-help/2003-August/036697....
そのほか，http://maths.newcastle.edu.au/~rking/R/help/03b...
-統計グラフで特許申請というアイデアそのものは「目から鱗」...
-ついでに。特許は有用性ではなく、新規性に対して与えられる...
-このグラフの特許問題に関する R-help 記事に Ross Ihaka が...
（R をオープンソースとする決定をしたこと）は私には考えて...
わかりました。私の所属する機関は知的な財産、特にソフトウ...
(diabolical ママ)な方針を持っています。我々がソフトを密...
漏らした時、私は私の作品の自由な公開について承認を求める...
要とし、辞職から給料のカットまでいくつもの脅しを受けまし...
#comment
**Bioconductor 1.4 released [#s38ae704]
> (2004-05-18 (火) 07:08:01)~
~
r-help でBioconductor新バージョンの発表がありましたが、Rd...
Mingw でコンパイルできる [[PostgreSQL Native Win32:http...
//
-自分ではちゃんと確認していないのですが、PostgreSQL7.5でW...
#comment
**Rのロゴ表示 [#j50ce055]
>[[谷村]] (2004-05-14 (金) 22:29:19)~
~
ロゴ表示コマンド(addlogo)の紹介がRjpWikiの中に見あたらな...
> library(pixmap)
> x <- read.pnm(system.file("pictures/logo.ppm", package...
> plot(runif(100),runif(100))
> addlogo(x,px=c(0.8,1),py=c(0,0.2))
#ref(addlogo.png,center)~
***応用例 [#dd062563]
方位記号をロゴとして加える
#ref(nagasaki.png,center)~
//
-上の方位記号の画像はフリーのものでしょうか？ -- &new{20...
-ArcViewのものを参考に適当に作成したものです。 -- [[谷村]...
#comment
**日本人の活躍 [#ted6da05]
>[[谷村]] (2004-05-14 (金) 21:47:25)~
~
汎用空間データ関数パッケージ(sp)のヘルプを何気なく見たら~
holepolys Hisaji Ono's hole problem~
を発見！尾野さん、素晴らしいです。RでGISも期待しています。~
//
#comment
**carのrecode(条件反射でコードを書くと) [#t1afe2f2]
>[[なかま]] (2004-05-11 (火) 08:47:51)~
~
ぼやーっと、~
x <- 0:100
y <- recode ( x, '0:19 = "未成年" ; 20:59 = "成人" ; els...
のようにしたものの、18才未満の方がいいかなと、~
recode ( x, 'x[ x < 18] = "未成年" ; x[ x < 60] = "成人"...
とかしてしまい、~
recode ( x, 'x[ x <= 12] = "小児" ; x[ x < 18] = "未成年...
x[ x < 60] = "成人" ; else = "老人" ')
とした時点でパーサーエラー。~
'='でsplitしているのでエラーになったんですね。~
で、条件式を書いてもOKとは書いてない事に最後に気づいてし...
# 良く読め<オレ~
でもそれって、~
x <- 0:10
recode ( x, '0:5 = ";" ; 6:8 = ":" ; else = "x"' )
もエラーになる訳で。。。~
//
-地道にcutとlevelsでやれば，最後の例はエラーにはなりませ...
-ありがとうございます。地道に精進致します。 -- [[なかま]]...
#comment
**dynamicGraphについて（Vineユーザのみなさま) [#j2b6880a]
>[[谷村]] (2004-05-10 (月) 16:43:06)~
~
dynamicGraphはtcltkに依存しています。Vineオリジナル版では...
//
-仮に有効にしてもそのままではPostScriptに吐いてもWinでもL...
-あっなかまさん。VinePlusへR-1.9日本語版をputする件はいか...
-うぐっ、もう少し、心の準備が出来てから。。。(汗; -- [[な...
#comment
**カンサス州知事は R ユーザー？ [#dbd0446b]
>[[間瀬]] (2004-04-29 (木) 14:41:20)~
~
最近の r-help ML にカンサス州知事 Kathleen Sebelius から...
//
#comment
**車輪を二度発明すること [#ve47fe05]
> (2004-04-23 (金) 18:00:47)~
~
splancs に areapl というのがあってポリゴンの面積を計算で...
~
ポリゴンの面積を計算するのは，R だとものすごく短く書ける。~
頂点の座標を時計回りに順にn×2行列にして与える。~
area <- function(xy)
{
n <- length(x <- xy[,1])
sum((c(x[-1], x[1])-c(x[n], x[-n]))*xy[,2])/2
}
# test run. answer=28.5
xy <- matrix(c(9,4,7,5,6,4,7,2,5,1,2,3,3,4,1,5,1,7,3,7,4...
> area(xy)
[1] 28.5
//
-areaplは、Fortranでコンパイルされた関数を呼んでいますが...
-30時間かかる計算が1時間でできると，きっとうれしい。30秒...
#comment
**全世界の R ユーザー数を推定する (r-help 記事より) [#seb...
> (2004-04-20 (火) 08:14:42)~
~
I have been trying to determine the size of the R user ba...
~
1. It appears that there are about 100,000 S-PLUS users.~
~
Rationale: According to Insightful's 2002 Annual Report, ...
~
Conclusion: S-PLUS costs $3500 (Windows) or $4500 (Linux/...
~
Does anyone has any other information that would give som...
~
2. At least one R book has achieved sales of just over 5,...
~
Conclusion: Very few books sell to more than 12% of the ...
~
3. There are a total of 3225 unique subscribers to the th...
//
-r-help はこの件で盛り上がっています。やはり、皆気になっ...
#comment
**R 1.9.0 正式公開 [#v48c431e]
>[[間瀬 ]] (2004-04-12 (月) 22:20:15)~
~
R 1.9.0 が正式に公開されました。各バイナリーの公開はすこ...
//
-lfactorial も同時に設けられたそうです。統計プログラムを...
~
> factorial(10)
[1] 3628800
> lfactorial(10)
[1] 15.10441
> x <-matrix(1:4, nc=2)
> addmargins(x)
[,1] [,2] [,3]
[1,] 1 3 4
[2,] 2 4 6
[3,] 3 7 10
-VineLinux用のrpmを本家にputしました。本家及びミラーに反...
#comment
**emacs-21とessと日本語 [#q0d91554]
>[[なかま@金沢市民]] (2004-04-11 (日) 18:09:50)~
~
ヘルプで日本語が出ない。essでの日本語の扱いで旨く行かない...
ありましたが、解決出来ました。~
~
(set-language-environment "Japanese")
(set-default-coding-systems 'euc-japan-unix)
(set-terminal-coding-system 'euc-japan-unix)
(set-keyboard-coding-system 'euc-japan-unix)
(set-buffer-file-coding-system 'euc-japan-unix)
(require 'ess-site)
(define-key ess-mode-map "?177" 'delete-char)
(setq ess-ask-for-ess-directory nil)
(setq ess-pre-run-hook
'((lambda () (setq S-directory default-directory)
(setq default-process-coding-system '(euc...
)))
//
-setq ess-pre-run-hook に (setq default-process-coding-sy...
#comment
**The R Package for multivariate and spatial analysis [#e...
> (2004-04-07 (水) 07:04:15)~
~
r-help 記事によれば [[The R Package for multivariate and ...
//
-個人的にはGoogleにかかりやすいようにRの名前の方を変えて...
-私は「R 統計」と，"統計"をかませて検索します． -- &new{...
#comment
**円周率はなぜπか (オフトピック) [#ce3cc40f]
>[[間瀬茂]] (2004-04-06 (火) 06:40:47)~
~
「工学のためのデータサイエンス入門」の23頁のコラム「円周...
であると書きましたが、これは不正確でした。確かに π が円周...
オイラーの1748年出版の広く読まれた著書「無限小解析入門」...
既に1647年に π を円周率を表す記号として使った例があるそう...
p とか書いていたそうです。詳しくは「π - 魅力の数」、ジャ...
を見てください。(若干一名から好評でしたのでずにのっておし...
分母が１より大きいすべての有理数 p/q に対し次の不等式
|π-p/q| ＞ 1/q^{20.6}
があるからπはランダムではない、と書いており、びっくりした...
( 10^{-5} ?sum_{n=-?infty}^{?infty} exp(-n^2/10^{10}) )^2
や、C言語でπを2400桁計算できる158語のプログラムなど、面白...
//
-考えてみたら上の無限級数は N(0, 10^{10}/2) の密度関数の ...
-細かいことで申し訳ありませんが、書名、出版社名、発行年が...
-失礼、題名と出版社はおっしゃる通りです。2003年に第二刷が...
#comment
**survey package [#rb516e5c]
>[[あいだま]] (2004-04-03 (土) 16:08:23)~
~
Survey Packageの翻訳が不適切です。↓~
階層化、クラスター化、非一様標本調査に対する要約統計量、...
~
Summary statistics, generalised linear models, and genera...
~
層化抽出、集落抽出、等確率でない抽出の場合の標本のための~
記述統計、一般化線形モデル、pseudolikelihood estimationの...
テイラー近似による分散の推定、replicate weights、~
//
//-自己レスで失礼しますが。replicate weightsの訳は当方知...
-自己レスで失礼しますが。replicate weightsの訳は当方知り...
-どこの訳を指していらっしゃるのでしょうか、情報不足です。...
-replicate weights はただの標本抽出用のウェイトではなくて...
#comment
**ご意見をお聞きします [#i7514b06]
>[[間瀬茂]] (2004-04-03 (土) 10:28:59)~
~
間瀬です。私どもが企画した「工学のためのデータサイエンス...
//
-趣旨に合えば問題ないと思います。実際、個人Wikiサイトに別...
-お騒がわせしましたが、結局複雑な数式を表現するには latex...
-[[工学のためのデータサイエンス入門の正誤表:http://www.is...
#comment
**The R Reference Manual 出版！ [#sdf77e03]
>[[間瀬茂]] (2004-04-03 (土) 00:32:09)~
~
~
R のベースパッケージ中のすべての関数の参照マニュアルを二...
//
#comment
**RのSVMで学習した結果を外部ファイルにエクスポート [#ddf7...
>[[高階知巳]] (2004-03-28 (日) 11:11:54)~
~
ウィーン工科大学の授業で使うために作られたというe1071パッ...
//
-e1071_1.4 が先日CRANでリリースされ、その中に、上述のパッ...
#comment
**rimageのWindows版バイナリが公開されました [#ra12e198]
>[[高階知巳]] (2004-03-28 (日) 11:01:17)~
~
RでJPEG画像ファイルの読み込みと基本的な画像処理機能を提供...
//
- Windows ユーザーです．早速パッケージを頂きました．examp...
-Rでは画像の表示は遅いですが、特にWindowsでは遅い様です。...
-高階さんが飼ってらっしゃる猫でしたか^ ^．すみません，猫...
-高階知巳さん、日本発 R アドオンパッケージ第一号(?)嬉しい...
-間瀬先生、ありがとうございます。私もパッケージ作者が増え...
#comment
**Nepia の対数 （オフトピック） [#j5aaadbf]
>[[間瀬茂]] (2004-03-27 (土) 09:41:08)~
~
また年寄りの懐古趣味記事です。統計学でも必須の対数関数の...
x ＝ 10^7*(1-10^{-7})^y
という関係にある数 x.y を数表化したのだそうです。比例定数...
//
-「大数の法則」のお話共々，興味深く読ませて頂いております...
-おっしゃる通り、私供の本「工学のためのデータサイエンス入...
#comment
**Knoppix/Math/2004 [#a1993610]
>[[間瀬茂]] (2004-03-24 (水) 16:28:29)~
~
R-jpに岡田さんが投稿されていますが、3月の筑波大学の数学会...
//
-Ring Server で公開をはじめました。収録しているRは日本語...
-[[Knoppix/Math/2004:http://unit.aist.go.jp/it/knoppix/is...
#comment
**useR! 2004 プログラム公開 [#g85fe9a8]
> (2004-03-23 (火) 19:32:17)~
~
日本人(?)による講演も三組。 [[プログラム:http://www.ci.tu...
//
#comment
**R 1.9.0 alpha [#a4aba883]
>[[間瀬茂]] (2004-03-21 (日) 07:02:30)~
~
r-help 記事によると R 1.9.0 alpha を試した人から、MCMC 法...
//
-1.9.0ではアンダースコアを変数名に使えるようになるとか，...
#comment
**Mac OS X Panther と Raqua （r-help 記事より） [#j6697ecd]
> (2004-03-12 (金) 17:22:45)~
~
ところで Raqua とは Mac OS X 用の R のことですよね(か？)
I've been nicely surprised, glancing through O'reilly's ...
for Unix Geeks" (by Brian Jepson and Ernest E. Rothman),...
the sixth section of chapter seven is dedicated to Raqua...
a few lines to describe R and Stefano's port to mac OS X...
(graphic windows and R console), and even a small apples...
interact with Raqua.~
BTW, the book is great.~
~
--Mathieu~
//
- => 青木さんの頁で確認して安堵。ところで Raqua というネ...
-aquaはMac OS XのGUIで、これを使ってＲが動くからだと思い...
-老婆心ながら、1.9 以後から単にRと呼ばれるようになりまし...
#comment
**R マニュアル翻訳チーム結成の呼びかけ（ただしフランスの...
> (2004-03-12 (金) 07:15:42)~
~
r-help にフランス人から投稿があり、前から自分が進めている...
//
#comment
~
**大数の法則 (オフトピック) [#z8471938]
>[[間瀬茂]] (2004-03-11 (木) 17:09:54)~
~
確率・統計の本の最初には必ずでてくる「大数の法則」ですが...
//
#comment
**eps ファイルの編集方法 [#tdc076bf]
>[[舟尾]] (2004-03-07 (日) 18:53:38)~
~
本筋から離れるかもしれませんが，皆さんは R で出力した eps...
//
-恥ずかしながら、(e)ps ファイルを編集するなぞ考えもしませ...
-UNIX系のOS (と Mac OS X) には [[Tgif:http://www.tobata.i...
-有料ですが、[[花子2004:http://www.justsystem.co.jp/hanak...
-花子が eps 形式に対応したというのは初耳でした．[[Canvas:...
-凄い技ですね。ぼくは手作業編集したいときはwin.metafile()...
-花子2004 は eps ファイルの編集が出来るんですか！かなり興...
-あと，[[ここの掲示板:http://www.visio.jp/bbs/bbsView.asp...
-ついでですが、Knoppix-jp の Openoffice はこれまで何故か...
-既存の eps も，[[pstoedit:http://www.pstoedit.net/pstoed...
-度々、枠領域を明示・・・ではなくて、asp=??を指定できる関...
#comment
**R on Zaurusにてグラフ出力時のエラー [#wf17f80c]
>[[aki]] (2004-03-07 (日) 06:19:07)~
~
現在、Zaurus上のX/Qtで動作するRを用意しています。~
コンパイルとipkの用意はできたのですが、いざ起動して~
> x<-rnorm(50)
> y<-rnorm(50)
> plot(x,y)
とすると、~
Warning messages:
1: Nonfinite axis limits [GScale(nan,nan,1, .); log=0]
2: relative range of values =hogehoge * EPS, is small (a...
3: Nonfinite axis limits [GScale(-inf,inf,2, .); log=0]
4: relative range of values =hogehoge * EPS, is small (a...
というメッセージと共にへんてこりんなグラフが表示されてし...
検索すると、~
http://lists.debian.org/debian-arm/2003/debian-arm-200311...
http://maths.newcastle.edu.au/~rking/R/devel/03b/1083.html~
が見つかりましたが、結局のところ解決法が分かりません。~
http://sgp.zaurii.net/binaries/index.html#R~
のRは使えるような話が[[LinuxザウルスSL-C700でR]]には載っ...
どなたか解決法が分かる方はいらっしゃいませんか?~
//
-ATLAS http://math-atlas.sourceforge.net/ が必要?かもしれ...
-以前、Ｘ外した１７１postscriptバージョンを作成した時も同...
-度々、枠領域を明示・・・ではなくて、asp=??を指定できる関...
-BLASとLAPACKを用意してみましたが症状変わらず。うーん、あ...
-FORTRANコンパイラをg77からf2cにしてみましたが今度はエラ...
-f2cでコンパイルしましたがエラー症状変わらず。というわけ...
#comment
**高校生向け R [#i2a9515d]
> (2004-03-07 (日) 00:25:48)~
~
なるほど、高校生に R を使わせるとは考えつきませんでした。...
//
-はじめまして。中学・高校現場にいる者です。実はRを利用し...
-よろしくお願いします。 -- &new{2004-03-12 (金) 22:04:11};
-もしかすると、R は統計学以上に高校生用の情報教育のプラッ...
-職場での公開について許諾が得られたので、PDF形式で公開い...
-すごいですね。数学の確率・統計の一部として「R を使えばこ...
-「データを解析する経験をしてみよう。そのときのツールとし...
-だんだん形を整えていって、「記述統計学の基礎（平均・分散...
-全員対象なので「賛否両論」というところだったでしょうか。...
-ちなみに（内容をみてもわかると思うのですが）あらかじめこ...
#comment
**Statistiques avec R [#n573321b]
>[[間瀬茂]] (2004-03-04 (木) 21:08:48)~
~
フランス語ながら「Statistiques avec R」というマニュアルは...
//
-リンク集にも書いてありますが，特に4章のグラフィックスに...
-他の章も一読あれ。やはりグラフィックスがたくさんあります...
-r-help でも好評のようです。一人が反応メイルで arigato go...
-r-help ではまさにニュースになったようです。これの英語訳...
-著者の Zoonekynd さんはどうやら計算機オタクの数学者らし...
-Zoonekynd という人物に興味を持ち少しグーグルしてみました...
--[[TeX用フォント一覧:http://tex.loria.fr/fontes/zoonekyn...
--[[Metapost による各種図形・図式作図例:http://tex.loria....
--[[難易度順漢字一覧表:http://zoonek.free.fr/kanji/list.p...
--[[Gimp マニュアル:http://zoonek2.free.fr/UNIX/23_perl_g...
-Zoonekynd さんはやはり日本語を理解するようです。中級程度...
#comment
**「データサイエンス入門」出版案内 [#ac86cc3d]
>[[間瀬茂]] (2004-03-03 (水) 00:15:21)~
~
既に少しお知らせしましたが、第二の R 本が三月末に刊行予定...
//
-楽しみにしています。目次がわかるwebサイトはまだないので...
-最寄の本屋の在庫が切れないか…，今から心配しております． ...
-追加コメント。原稿作成中、ますます R の魅力にはまりまし...
-「慶応大学のある先生」とはどなたでしょうか？少し興味があ...
-[[慶応大学のある先生:http://web.sfc.keio.ac.jp/~watanabe...
-ありがとうございました -- [[Yano]] &new{2004-03-10 (水) ...
#comment
**表計算ソフトの統計機能 [#lea93164]
> (2004-02-28 (土) 19:20:21)~
~
Excel の統計機能の問題については青木先生のホーム頁に詳し...
//
-OpenOffice.orgのcalcも，相関係数のエラー（[[青木先生のペ...
-やはり駄目ですか。悪いところも含めて Excel にできるだけ...
-例の Excel をこきおろすページ(^_^;)というのは，無理矢理...
#comment
**Vine Linux 版の R のインストール法の記事? [#l3505498]
> (2004-02-25 (水) 17:36:10)~
~
Vine Linux 版の R のインストール法の記事がいつのまにか消...
//
-私が書いたネットワークインストールの記事でしょうか。apt...
# rpm -ivh R-1.8.1-4vl1.i386.rpm
#comment
**S/R で商売 [#q59cc20e]
>[[間瀬茂]] (2004-02-25 (水) 17:11:06)~
~
r-help に流れた広告。既に日本にも例があるようですが、R の...
XSolutions Corp ([1]www.xlsolutions-corp.com) is proud...
a 2-day "Advanced R/Splus programming" taught by R Dev...
Core Team Guru.
*********San Francisco ----- March 22-23
*********Washington DC ---- March 22-23
*********Chicago ----------- March 25-26
*********New York --------- March 22-23
*********Atlanta ----------- March 25-26
*********Boston ------------ TBD
If your location of choice isn't listed, please let us...
Reserve your seat Now at earlybird rates (payment due...
class)
Registration:
Email Sue Turner: [2]sue@xlsolutions-corp.com
Phone: 206-686-1578
Visit us: [3]www.xlsolutions-corp.com/training.htm
Course outline:
- Overview of R/S fundamentals: Syntax and Semantics
- Class and Inheritance in R/S-Plus
- Concepts, Construction and good use of language obje...
- Coercion and efficiency
- Object-oriented programming in R and S-Plus
- Advanced manipulation tools: Parse, Deparse, Substit...
- How to fully take advantage of Vectorization
- Generic and Method Functions; S4 (S-Plus 6)
- Search path, databases and frames Visibility
- Working with large objects
- Handling Properly Recursion and iterative calculations
- Managing loops; For (S-Plus) and for() loops
- Consequences of Lazy Evaluation
- Efficient Code practices for large computations
- Memory management and Resource monitoring
- Writing R/S-Plus functions to call compiled code
- Writing and debugging compiled code for R/S-Plus sys...
- Connecting R/S-Plus to External Data Sources
- Understanding the structure of model fitting functio...
- Designing and Packaging efficiently
Early-bird group research fee: $995!
This course will also deal with lots of S-Plus efficie...
any special topics from participants is welcome.
Please let us know if you and your colleagues are inte...
class to take advantage of group discount. Register no...
your seat!
Cheers,
Elvis Miller, PhD
Manager Training.
XLSolutions Corporation
206 686 1578
[4]www.xlsolutions-corp.com
References
1. http://www.xlsolutions-corp.com/
2. mailto:sue@xlsolutions-corp.com
3. http://www.xlsolutions-corp.com/training.htm
4. http://www.xlsolutions-corp.com/
//
#comment
**Windowsでパッケージを自作 [#qaf45149]
>[[かくれR狂]] (2004-02-25 (水) 04:21:27)~
~
Windowsでパッケージ自作に挑戦していました。~
山口県立大の中澤先生の[[統計処理ソフトウェアRについてのTi...
必要なソフト探しからだったので結構大変でした。~
Tipsに紹介しようかと考えていますが、~
もうすでにどこかで紹介されているのでしょうか？~
//
-日本語での解説は基本マニュアルにかかれている以外は知りま...
-Ｔｉｐｓに載せてみました。どなたか挑戦してみてください。...
-高階さんのrimageがあるので，第二号ですね。ぼくがMsysと角...
#comment
**新パッケージ mvpart (rpart の多変量データ版) 2004.2.22 ...
> (2004-02-23 (月) 07:15:05)~
~
The package mvpart is now available.~
mvpart includes partitioning based on (1) multivariate nu...
(2) dissimilarity matrices.~
~
The package mvpart is a modification of rpart --~
-- authors of original: Terry M Therneau and Beth Atkinso...
Includes some modified routines from vegan -- Jari Oksane...
//
-mvpartと同様の multivariate regression trees のソースコ...
#comment
~
**「Rを味見してみる」 [#h8fcfce3]
>[[舟尾]] (2004-02-21 (土) 15:02:56)~
~
R界で(も)お馴染みの青木先生が，またもや凄まじいコンテンツ...
[[R を味見してみる:http://aoki2.si.gunma-u.ac.jp/R/R0.htm...
これで「Rというソフトがあるらしいがインストールするのは…...
//
-すみませんでした。諸般の事由で，中止致します。おさわがせ...
-最近気づきましたが[[中間さんち:http://r.nakama.ne.jp/Rwe...
-青木先生のコンテンツが無くなってしまってげんなりしていた...
-Javaアプレットが見れる状態にして、しばらく待つと見ること...
-Rweb の html 15行目にエラーがありませんか？ -- &new{200...
-15行目直しました。「R を味見してみる」で覗いてたのは私で...
-周りの人に紹介しています．なかなか好評でありまして，R を...
#comment
**ブッシュ大統領曰く、私は統計学者でも預言者でもない。 [#...
>[[間瀬]] (2004-02-19 (木) 18:44:42)~
~
r-help 記事によれば、ブッシュ大統領のスタッフが、今年末ま...
//
#comment
**舟尾チップス集 [#q1c11030]
>[[間瀬]] (2004-02-13 (金) 22:06:29)~
~
舟尾チップス集にはまだ私も知らない話題がたくさんあります...
//
-私は[ファイルへのデータ出力例]の節での，ファイル出力の際...
-お気づきの「へー」はどうぞ適当な箇所に追加お願いします。...
-ここに逐次書き込みをするのかと勘違いしました。お詫びする...
#comment
**useR! 2004 - The R User Conference [#k4bcc755]
> (2004-02-12 (木) 15:07:36)~
~
2004年5月20日‐22にウィーンのTechnische大学でRユーザカンフ...
開催されます。演題締切は2月15日ですのでお早めに。~
http://www.ci.tuwien.ac.at/Conferences/useR-2004/~
//
#comment
**投稿記事 どうむ について [#s2951c8c]
~
> test(.Machine$double.eps)
TRUE
> test(.Machine$double.eps/2)
TRUE
となりました。
業務でRを使用しているのですが、if(Cond)のCondにそもそもこ...
とりあえずご報告まで。
//
-こうした記事は独立して投稿すると、すぐ行方不明になります...
> .Machine$double.eps/1e10
[1] 2.220446e-26
と .Machine$double.eps 以下の正実数も立派に存在します。次...
> if(1+.Machine$double.eps==1) cat("TRUE?n") else cat("F...
FALSE
> if(1+.Machine$double.eps/2==1) cat("TRUE?n") else cat(...
TRUE
-- &new{2004-02-02 (月) 17:19:39};
-申し訳ございません -- どうむ &new{2004-02-02 (月) 17:39:...
-「if (x) の不審な(?)行動」に対するRESだったのですが当方...
-試してみたら.Machine$double.eps/1e308 が零になりました。...
> if(.Machine$double.eps/1e307) cat("TRUE?n") else cat("...
TRUE
> if(.Machine$double.eps/1e308) cat("TRUE?n") else cat("...
FALSE
-どうむさんのご指摘のように「if (x) の不審な(?)行動」の記...
#comment
**文字列のオブジェクトへの変換方法 => 追加コメントは Q&A ...
>[[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~
つまり、文字列として取り込んだものをオブジェクトとして扱...
いろいろ探してみたのですが、見つかりませんでした。~
どなたか教えてもらえないでしょうか？~
-getを使うとオブジェクトの中身が取り出せる。 money[get(a)...
#comment
**too large data! [#la36eaf9]
>[[間瀬]] (2004-01-29 (木) 17:36:59)~
~
r-help 記事より。R でどれだけ大きなデータが処理できるかは...
R-1.8.1,
compiled as a 64bit application for a Solaris 5.8, 64 bi...
8Gb of RAM を一人で使える
データは 2,000,000 行超で、多くは因子の60変数、一部は連...
ファイルからの読み込みとセーブに5分以上、再読み込みに15...
S-plus では可能だった解析を行なった
一回の反復で45分後に R はすべてのメモリを使い果たし、ス...
メモリ使用は最大を指定している（はず）
実際の解析コードは以下の通り
library(package = "statmod", pos = 2,
lib.loc = "/home/jeg002/R-1.8.1/lib/R/R_LIBS")
qc.B3.tweedie <- glm(formula = pp20B3 ~ ageveh + anpol +
categveh + champion + cie + dossiera +
faq13c + faq5a + kmaff + kmprom + nb...
rabprof + sexeprin + newage,
family = tweedie(var.power = 1.577,
link.power = 0),
etastart = log(rep(mean(qc.b3.sans.o...
'pp20B3']), nrow(qc.b3.sans.occ))),
weights = unsb3t1,
trace = T,
data = qc.b3.sans.occ)
コアメンバーの Dalgaard から次のようなコメント
データはそれ自体で 0.5-1GB になるはず
60変数のうち実際使っているのは15変数だから、最初からそれ...
因子の水準が多過ぎる可能性？
とくに連続変数を因子に変換すれば違った結果が得られるかも...
-再び r-help 記事より。R の現在の基本ライブラリーは最大利...
-上記記事に関連して。６４ビットCPU の Linux 機では１６GB...
#comment
** 不適当記事一件削除 [#abb2ca4d]
**Rnews 2003/3 公開 [#ud22374c]
> (2004-01-22 (木) 18:29:42)~
~
We have published the 2003/3 issue of R News on
http://cran.R-project.org/doc/Rnews
where you can download the newsletter as PDF or Postscri...
will propagate to the CRAN mirrors within a day or two.
Contents of this issue:
Dimensional Reduction for Data Mapping
R as a Simulation Platform in Ecological Modelling
Using R for Estimating Longitudinal Student Achievement ...
lmeSplines
Debugging Without (Too Many) Tears
Th R2HTML Package
R Help Desk: Package Management
Programmer's Niche: Regular Expressions
Recent Events
Upcoming Events
Changes in R 1.8.1
Changes on CRAN
R Foundation News
//
#comment
**直線回帰の信頼区間 [#a8cb8751]
>[[ぶんしょ]] (2004-01-18 (日) 19:21:24)~
~
回帰させた直線の95％信頼区間を知りたいのですが、どのよう...
//
-predict()関数を使うのが便利です。直線回帰でいいなら，exa...
-ありがとうございます。さっそく試してみます。 -- [[ぶんし...
#comment
~
**[質問] 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
-丁寧な回答をありがとうございます。ただ、私の記事は質問コ...
-どういたしまして。依然としてなぜ lm オブジェクトを R プ...
-いえいえ、私のHPは間瀬さんが翻訳されたR入門や他の方々の...
#comment
**ESRI shape ファイル [#ve85d9d7]
>[[なかま＠金沢市民]] (2004-01-02 (金) 18:37:18)~
~
あけまして、おめでとうございます。~
http://www.diva-gis.org/data/DataServer.htm~
に、色んな国のがありました。~
日本のだと、都道府県レベルまででしたが、~
plot(j[j$ADM2 == "Ishikawa"])
とか、最低限は出来そう。~
//
#comment
終了行:
COLOR(blue){SIZE(20){RjpWiki自体へのコメントでも、雑談で...
COLOR(red){注意}：新規記事用の入力欄は以下の目次の直後に...
----
COLOR(red){過去の記事のアーカイブ} 。アーカイブ中の記事へ...
-[[なんでも掲示板アーカイブ（１）]] (元記事が 2003-07-17 ...
-[[なんでも掲示板アーカイブ（２）]] (元記事が 2003-11-27 ...
----
#contents
----
回答を期待する記事は原則 COLOR(red){質問コーナー [[Ｑ＆Ａ...
~
#article
**なるほど！ [#ecceab77]
> (2004-06-30 (水) 23:58:28)~
~
r-help 記事でまた嬉しくなる発見をしました。言われてみれば...
test1 <- function(x) 1+(x-2)*2*(x-2<0)+(x-2)*(-2)*(x-2>=...
plot(test, -1, 5)
こうした論理値の使い方を知らなければ次のように定義したで...
test2 <- function(x) ifelse(x < 2, 1+2*(x-2), 1-2*(x-2))
ついでに奇妙なことに気づきました。test1 は二回の論理判断...
> set.seed(222); system.time(test1(runif(1e6)))
[1] 0.38 0.30 0.69 0.00 0.00
> set.seed(222); system.time(test2(runif(1e6)))
[1] 0.72 0.33 1.06 0.00 0.00
//
-計算速度も速くなるんですねぇ！勉強になりました．例として...
f(x) = 0.6φ(x)+0.4ψ(x)
となる密度 f(x) に従う乱数を生成する関数を定義してみます．
generator <- function(n) {
data1 <- rnorm(n)-1 # φ(x) に従う乱数
data2 <- rnorm(n)+2 # ψ(x) に従う乱数
data3 <- (runif(n) <= 0.6) # 確率0.6でdata1を採択し，
data1*data3+data2*(1-data3) # 確率0.4でdata2を採択する...
}
mydata <- generator(1000)
([[ヒストグラムと密度の推定:http://www.okada.jp.org/RWiki...
-R は lisp で動いているのでしたっけ。知能的な最適化が組み...
上の関数は，以前に拝見して，次のような書き換えをしてみた...
> generator <- function(n) {
+ data1 <- rnorm(n)-1 # φ(x) に従う乱数
+ data2 <- rnorm(n)+2 # ψ(x) に従う乱数
+ data3 <- (runif(n) <= 0.6) # 確率0.6でdata1を採択...
+ data1*data3+data2*(1-data3) # 確率0.4でdata2を採択...
+ }
> system.time(mydata <- generator(100000))
[1] 0.24 0.02 0.28 0.00 0.00
> generator2 <- function(n) {
+ (rnorm(n)-1)+3*(runif(n) > 0.6)
+ }
> system.time(mydata2 <- generator2(100000))
[1] 0.15 0.00 0.15 0.00 0.00
ついでに，先の関数も
Rweb:> test1 <- function(x) 1+(x-2)*2*(x-2<0)+(x-2)*(-2)...
Rweb:> system.time(test1(rnorm(1000000)+2))
[1] 2.79 0.13 5.94 0.00 0.00
Rweb:> test2 <- function(x) 1+2*(x-2)*((x<2)-(x>=2))
Rweb:> system.time(test2(rnorm(1000000)+2))
[1] 1.36 0.09 1.45 0.00 0.00
人手で最適化しても，たいした効果はないようにも見えます。 ...
-R は Scheme を手本にしているようですが、Lisp インタプリ...
-す、素晴らしい！鳥肌が立ちました！関数 generator2() は速...
-上の書き換え (rnorm(n)-1)+3*(runif(n) > 0.6) は元のコー...
-私も最初「間違いなのでは？」と思ったのですが、
--確率0.6でφ(x)に従う乱数を、確率0.4でψ(x)に従う乱数を生...
という手順は、
--まず平均-1，分散1の正規乱数（φ(x)に従う乱数）を生成する
--この乱数を、確率0.4で平均2，分散1の正規乱数（ψ(x)に従う...
ことと同じになると思います。変更手順は単に「φ(x)に従う乱...
-「元のコードの二つの rnorm(n) が同じベクトルという間違い...
-なるほど！ 勘違い失礼。お詫びに（？）平均、分散がともに...
# 2つの正規分布の混合分布乱数
# n: 発生する混合正規乱数の総数
# p1: 最初の正規分布の比率
# N1: 最初の正規分布の平均、標準偏差のベクトル c(mu1, sd1)
# N2: 2番目の正規分布の平均、標準偏差のベクトル c(mu2, s...
rmixnorm <- function(n, p1, N1, N2) {
P <- (runif(n) < p1)
(P*N1[2] + (1-P)*N2[2])*rnorm(n) + (P*N1[1] + (1-P)*N...
}
-[[Rコードの最適化例：混合正規乱数の発生コード]] にありが...
-問題の解析
Rweb:> rmixnorm <- function(n, p1, N1, N2) {
+ P <- (runif(n) < p1)
+ (P*N1[2] + (1-P)*N2[2])*rnorm(n) + (P*N1[1] + (1-P...
+ }
Rweb:> mean(sapply(1:100, function(i) system.time(rmixno...
[1] 0.10932
100回やって平均を求めた。~
単に混合分布のデータを作るだけなら，乱数の発生順までラン...
Rweb:> rmixnorm2 <- function(n, p1, N1, N2)
+ {
+ n1 <- floor(n*p1)
+ c(rnorm(n1)*N1[2]+N1[1], rnorm(n-n1)*N2[2]+N2[1])
+ }
Rweb:> mean(sapply(1:100, function(i) system.time(rmixno...
[1] 0.04744
これじゃ，それぞれの分布の個数が決まってしまうじゃないか...
Rweb:> rmixnorm3 <- function(n, p1, N1, N2)
+ {
+ n1 <- sum(runif(n) < p1)
+ c(rnorm(n1)*N1[2]+N1[1], rnorm(n-n1)*N2[2]+N2[1])
+ }
Rweb:> mean(sapply(1:100, function(i) system.time(rmixno...
[1] 0.06088
みたいな。 -- &new{2004-07-03 (土) 18:28:01};
-(なるほど)^3、感服！ -- &new{2004-07-03 (土) 19:22:16};
-ああ，わたしってほんとにうっかりもの。平均値も標準偏差も...
Rweb:> rmixnorm4 <- function(n, p1, N1, N2)
+ {
+ n1 <- sum(runif(n) < p1)
+ c(rnorm(n1, mean=N1[1], sd=N1[2]), rnorm(n-n1, mean...
+ }
Rweb:> mean(sapply(1:100, function(i) system.time(rmixno...
[1] 0.05924
時間は大して変わらないけど，可読性はずいぶん違う。 -- &n...
-コメント -- &new{2004-07-04 (日)00:01:00};
> mean(sapply(1:100, function(i) system.time(rmixnorm4(1...
[1] 0.05618
ではなく(これは 100x5 行列の平均になってしまう)
> mean( sapply( 1:100, function(i) system.time( rmixnorm...
[1] 0.1298
としないとまずい。
-最適化された関数をいろいろ紹介して頂きましたが，「本当に...
generator <- function(n) {
set.seed(7777); data1 <- rnorm(n)-1
set.seed(7777); data2 <- rnorm(n)+2
set.seed(7777); data3 <- (runif(n) <= 0.6)
data1*data3+data2*(1-data3)
}
mydata <- generator(1000) # ...
plot(density(mydata)) # ...
generator2 <- function(n) {
set.seed(7777); data1 <- rnorm(n)-1
set.seed(7777); data1 <- data1 + 3*(runif(n) > 0.6)
}
mydata <- generator2(1000) # ...
plot(density(mydata)) # ...
rmixnorm <- function(n, p1, N1, N2) {
set.seed(7777); P <- (runif(n) < p1)
set.seed(7777)
(P*N1[2] + (1-P)*N2[2])*rnorm(n) + (P*N1[1] + (1-P)*N2...
}
mydata <- rmixnorm(1000,0.6,c(-1,1),c(2,1)) # ...
plot(density(mydata)) # ...
#ref(density.gif, center)
- 上と同じ動作をするような，rmixnorm2() 〜 rmixnorm4() の...
-こういうコーディング例の検証を考えると，「[[RUnit を使う...
-set.seed(7777); rnorm(20) と set.seed(7777); rnorm(15); ...
set.seed(7777); data1 <- rnorm(n)-1
set.seed(7777); data2 <- rnorm(n)+2
は，
set.seed(7777); data1 <- rnorm(n)-1; data2 <- data1+3
と全く同じということではないでしょうか。実際にやってみる...
> n <- 1000
> set.seed(7777); data1 <- rnorm(n)-1
> set.seed(7777); data2 <- rnorm(n)+2
> set.seed(7777); data11 <- rnorm(n)-1; data12 <- data11+3
> sum(data12==data2)
[1] 905
差はどれくらいかというと
[1] -4.440892e-16 -2.220446e-16 -2.220446e-16 2.220446e...
程度 -- &new{2004-07-04 (日) 16:52:40};
-乱数の不思議~
同種類の乱数だと以下のようになるが
> set.seed(1)
> rnorm(10)
[1] -0.6264538 0.1836433 -0.8356286 1.5952808 0.3295...
[9] 0.5757814 -0.3053884
> rnorm(5, mean=5, sd=10)
[1] 20.117812 8.898432 -1.212406 -17.146999 16.249309
この乱数列の最後の5つは
> set.seed(1)
> rnorm(15, mean=5, sd=10)
[1] -1.264538 6.836433 -3.356286 20.952808 8.295...
[9] 10.757814 1.946116 20.117812 8.898432 -1.212...
の，最後の5つと同じになっている。
しかし，
> set.seed(1); rnorm(5); runif(5); rnorm(5, mean=5, sd=10)
[1] -0.6264538 0.1836433 -0.8356286 1.5952808 0.3295078
[1] 0.2059746 0.1765568 0.6870228 0.3841037 0.7698414
[1] 4.942328 29.046534 12.635935 -2.990092 -6.476570
のように途中で，異なった種類の乱数発生を挟むと，違ったも...
-set.seed(n) は指定した正整数 n から乱数の初期値を適当に...
-「というか何度も同じ系列を使うことになりむしろ有害」と書...
-私の言葉足らずでした．すみません…．私がやりたかったこと...
よく考えると，rmixnorm2() 〜 rmixnorm4() は，生成すべき片...
よく考えると，rmixnorm2() 〜 rmixnorm4() は，生成すべき片...
（追記）generator() の data1 と data2 には同じ乱数列を与...
#comment
**日本人による新パッケージ [#vf5e7758]
> (2004-06-30 (水) 08:21:13)~
~
We would like to announce the release of our software, w...
available through CRAN.
MNP: R Package for Fitting the Multinomial Probit Models
Abstract:
MNP is a publicly available R package that fits the Baye...
probit models via Markov chain Monte Carlo. Along with t...
multinomial probit model, it can also fit models with di...
sets for each observation, and complete or partial order...
available alternatives. The computation is based on the ...
data augmentation algorithm that is developed by Imai an...
``A Bayesian Analysis of the Multinomial Probit Model Us...
Augmentation,'' Journal of Econometrics, Forthcoming.
Kosuke Imai, Department of Politics, Princeton University
Jordan R. Vance, Department of Computer Science, Princet...
David A. van Dyk, Department of Statistics, University o...
//
#comment
**R for Bioinformatics [#s76a2aea]
>[[nami]] (2004-06-28 (月) 07:59:24)~
~
オライリー本家から「R for Bioinformatics」が、11月に出る...
[[amazon.co.jp:http://www.amazon.co.jp/exec/obidos/ASIN/0...
日本語化されないかな〜~
//
-オライリージャパンの方に直接聞いたのですが、現時点では和...
-どうもありがとうございました。なんとかがんばっていただき...
#comment
**R の一様乱数 runif() [#ie70b59f]
>[[間瀬茂]] (2004-06-27 (日) 00:35:43)~
~
r-help にオヤッと思う記事がありました。~
> set.seed(222)
> sum(duplicated(runif(1e5)))
[1] 4
つまり、10万個の一様乱数中にタイが４個(実際は２組のタイ)...
//
-32ビットマシンで，8バイト実数だと，表現できる正の実数は ...
-MT法は内部的にはもっとたくさんのビットを使っています。さ...
-周期と相異なる乱数の個数は無関係でしょう。 -- &new{2004...
-Ripley が注意しているのは
> qbirthday(classes=2^(53))
[1] 111743588
> qbirthday(classes=2^(32))
[1] 77163
つまり３２ビット実数一様乱数ではタイが生じる確率が 0.5 以...
-- &new{2004-06-27 (日) 10:59:34};
#comment
**package [#g9c79caf]
> (2004-06-26 (土) 23:33:46)~
~
パッケージquadprogの関数solve.QP()を利用してQPを解く際、~
そのパッケージ例に当てはまらないような目的関数を解くこと...
//
-これは独り言ですか、質問ですか？質問ならば場所が不適当で...
#comment
**関数 getAnywhere() [#g6b38f3c]
>[[間瀬]] (2004-06-24 (木) 21:09:49)~
~
また不思議な関数を発見しました。getAnywhere()。引数で与え...
//
#comment
**R fortunes [#lebbe91d]
> (2004-06-20 (日) 08:25:24)~
~
パッケージ更新案内で見つけ、R fortunes というパッケージを...
//
-SAS や「お間抜けな」r-help 投稿記事が格好の的になってい...
Some of us feel that type III sum of squares and so-call...
statistical nonsense which should have been left in SAS.
-- Brian D. Ripley
s-news (May 1999)
To paraphrase provocatively, `machine learning is statis...
of models and assumptions'.
-- Brian D. Ripley (about the difference between machi...
statistics)
useR! 2004, Vienna (May 2004)
Can one be a good data analyst without being a half-good...
answer to that is, `No.' The long answer to that is, `No.'
-- Frank Harrell
1999 S-PLUS User Conference, New Orleans (October 1...
#comment
**RCC (R compiler) ! [#s8315127]
>[[間瀬茂]] (2004-06-19 (土) 23:51:47)~
~
r-help 記事に驚くべきものを発見しました。R コンパイラー！...
[[参考１:http://www.cs.rice.edu/Colloquia/garvin-04-21.sh...
[[参考２:http://hipersoft.cs.rice.edu/rcc/]]~
//
-gcc, bcc に続いて rcc ですか！また凄まじい題材で修士論文...
-参考ホームページの内容だけからは良くわからないが、R コー...
#comment
**R は S(-plus) 互換？ [#a1655aa3]
>[[間瀬]] (2004-06-13 (日) 07:16:11)~
~
最近、J. Chambers の「データによるプログラミング」森北出...
//
#comment
**ｚｉｐのパスワード [#ja6562ba]
>[[初心者]] (2004-06-10 (木) 23:21:48)~
~
unix上でzip圧縮した際にパスワードを設定したいのですが、~
どなたかご存知な方いらっしゃいいますか？~
windowsだと基本的にzip圧縮の際にパスワードを設定できるよ...
一応manコマンドで調べると-eという暗号化のオプションが存在...
環境はsolarisです。よろしくお願いします。~
//
-せめて R との関連を申し訳でも書かないと顰蹙をかううこと...
$ zip -e abc abc.tex # abc.zip という名...
Enter password: # パスワード入力
Verify password: # もう一度
updating: abc.tex (deflated 67%)
$ ls
./ abc.tex abc.zip
$ unzip abc.zip
Archive: abc.zip
[abc.zip] abc.tex password: # 設定パスワードを...
replace abc.tex? [y]es, [n]o, [A]ll, [N]one, [r]ename: y
inflating: abc.tex
-少し強引にRの話題に振りますが、圧縮アルゴリズムって公開...
-?zip.file.extract もちろん R で書いてあるわけではなく、...
-なるほど、内部にunzipのコードは持っていて、すべてのプラ...
#comment
**R-jp, RjpWiki の来歴について [#u6fe8b45]
>[[間瀬茂]] (2004-06-10 (木) 21:19:54)~
~
明日は RjpWiki の開設一周年記念日です。管理者の岡田さん、...
//
-Rの日本での普及の歴史をそろそろまとめるといいかもしれま...
#comment
**予想したとおりになるでしょうか [#w264e6ae]
> (2004-06-10 (木) 13:27:27)~
~
以下のコード，出力はどうなると思いますか。~
やってみて，予想したとおりになったでしょうか。~
私は，驚きました。~
for (i in 1:10) {
if (i == 3) i <- i+5
print(i)
}
//
-スコープについて学ぶときにとてもよい教材ですね（print(i)...
-if 文の { } の有無は，私が意外と思ったことには関係ない...
先に挙げたコードは，以下のコードと同じことを表しているの...
for (i in 1:10) {
if (i == 3) {
i <- i+5
}
print(i)
}
ついでに，試してみたらこんなことも。
n <- 10
for (i in 1:n) {
if (i == 5) {
n <- 20
}
print(i)
}
どうしても for を使わなくてはならなかった局面で（笑），出...
~
help を見れば，書いてあることなんですけどね。~
The index 'seq' in a 'for' loop is evaluated at the start...
loop; changing it subsequently does not affect the loop.
-なるほど、forループではforループの最初のところに戻ってこ...
-駄目といわれると、掟破りをしたくなるのが悲しい人間のさが...
-次も不思議ではありませんか。ループ変数が最終値で残ってい...
> rm(i); for (i in 1:10) {i^2}; i
[1] 10
-私のループ変数の理解がぐらついてしまって非常に混乱してお...
-あくまでも推測ですが，for の最初に制御変数の集合を作って...
> x <- c(3,2,3,8,5)
> for (i in x) { if (i==2) x<-c(4:6); print(i)}
[1] 3
[1] 2
[1] 3
[1] 8
[1] 5
> x
[1] 4 5 6
> i
[1] 5
-i^2 は副作用なし（i に付値してませんから）です。 -- &ne...
> rm(i); for (i in 1:10) {i<-i^2}; i
[1] 100
こうしてみる？
> rm(i); sum<-0; for (i in 1:10) {sum<-sum+i^2}; print(i...
[1] 10
[1] 385
-C では確か、ループ変数をループ中で変更できませんでしたっ...
-C でも，BASIC でも，FORTRAN でも，AWK でも可能です。とい...
- 少し R をからかってみました。文句をいうが、へたれないと...
> rm(i); for (i in 1:3) {print(i); rm(i)}
Warning message: )}remove: variable "i" was not found
NULL
[1] 1
[1] 2
[1] 3
//-記事の衝突を調整中舟尾さんの記事をうっかり消してしまっ...
//-キャッシュに残っていたので，復帰させてみました。 -- &...
//-三人が集中してたのね -- &new{2004-06-10 (木) 22:39:11};
-in 以降の式を評価してスタックして、ループの間だ 指定した...
# 呪われた言語と言われるゆえん?
@ary = (1,2,3);
for $i ( (@ary,4) ){
push @ary,5 if ( $i == 2 );
print $i,"?n";
}
print join(",",@ary),"?n";
for $i ( "a","b","c" )
{
$i="c";
print $i,"?n";
}
1
2
3
4
1,2,3,5
Modification of a read-only value attempted at a.pl line...
-確かに以下のようなへんなことが出来る以上 「制御変数」と...
> r <- list(1, "a", matrix(1:4,nc=2))
> for (i in r) print(i)
[1] 1
[1] "a"
[,1] [,2]
[1,] 1 3
[2,] 2 4
#comment
**CRAN更新パッケージ一覧の表示方法 [#e2c27c40]
>[[岡田]] (2004-06-04 (金) 17:37:34)~
~
とりあえず、MenuBarの最新30件に混ざって出るようにしてみま...
//
-別ページ表示に一票 -- &new{2004-06-04 (金) 21:40:46};
#comment
**無相関は独立性を意味しない、例 [#l8e28809]
> (2004-06-01 (火) 22:58:00)~
~
R とは直接関係ない話題で済みません。また学識深い R ユーザ...
//
-正規分布のような裾の広い分布で示すより，一様分布のような...
-一様分布でも原点対称（例えば [-0.5, 0.5] 上の一様分布）...
#comment
**「R-Tips」のPDF版とその原稿 (04/06/07に修正版をup) [#y0...
>[[舟尾]] (2004-02-03 (火) 22:00:00)~
~
竹澤先生の御厚意により、竹澤さんのウェブスペースに[[R-Tip...
//
-素晴らしい。２８０ページの力作！もともとこんなに量があり...
-竹澤（中央農業総合研究センター・農業情報研究部）です。 ...
-改めて全体を見ました。感想は「すごい」。内容の豊富さも素...
-非常に恐縮であります…。「Introduction to R」の日本語訳版...
このPDF(+html)をネットに上げるにあたって、管理等を全てお...
-リンク切れでしょうか？PDF を見ることができないのですが・...
-復帰したようですね。 -- &new{2004-02-08 (日) 20:57:41};
-ですね。リンクの張り間違いではなかったようです。 -- &ne...
-ありがとうございます。これはものすごいPDFが出ましたね。 ...
- 2/6 に竹澤先生にお願いして up させてもらった PDF は，煩...
-内容についてではなく、つまらないことで恐縮です。印刷して...
-こちらのレーザープリンタ・インクジェットプリンタでは印刷...
-私も１２ページだけ印刷できませんでした。軽い版でも駄目で...
-↑この書き込みの者です。原因は分かりませんでしたが，解決...
-大変申し訳ありません．お手数をおかけしております…．上の...
-すみません。最初に12頁が印刷できないと書き込みした者です...
-お役に立てなくて申し訳ありませんでした…．また原因が分か...
#comment
**統計グラフで特許出願！ [#t93e7aa8]
>[[間瀬茂]] (2004-05-18 (火) 10:15:03)~
~
r-help 記事によると、三次元棒グラフの代替グラフである dia...
グラフが特許出願中だそうです。~
~
[[実物はここを参照:http://www.jhsph.edu/Press_Room/Press_...
//
-引用されたグラフ中の数字が読みとれないので，ネットサーフ...
https://stat.ethz.ch/pipermail/r-help/2003-August/036697....
そのほか，http://maths.newcastle.edu.au/~rking/R/help/03b...
-統計グラフで特許申請というアイデアそのものは「目から鱗」...
-ついでに。特許は有用性ではなく、新規性に対して与えられる...
-このグラフの特許問題に関する R-help 記事に Ross Ihaka が...
（R をオープンソースとする決定をしたこと）は私には考えて...
わかりました。私の所属する機関は知的な財産、特にソフトウ...
(diabolical ママ)な方針を持っています。我々がソフトを密...
漏らした時、私は私の作品の自由な公開について承認を求める...
要とし、辞職から給料のカットまでいくつもの脅しを受けまし...
#comment
**Bioconductor 1.4 released [#s38ae704]
> (2004-05-18 (火) 07:08:01)~
~
r-help でBioconductor新バージョンの発表がありましたが、Rd...
Mingw でコンパイルできる [[PostgreSQL Native Win32:http...
//
-自分ではちゃんと確認していないのですが、PostgreSQL7.5でW...
#comment
**Rのロゴ表示 [#j50ce055]
>[[谷村]] (2004-05-14 (金) 22:29:19)~
~
ロゴ表示コマンド(addlogo)の紹介がRjpWikiの中に見あたらな...
> library(pixmap)
> x <- read.pnm(system.file("pictures/logo.ppm", package...
> plot(runif(100),runif(100))
> addlogo(x,px=c(0.8,1),py=c(0,0.2))
#ref(addlogo.png,center)~
***応用例 [#dd062563]
方位記号をロゴとして加える
#ref(nagasaki.png,center)~
//
-上の方位記号の画像はフリーのものでしょうか？ -- &new{20...
-ArcViewのものを参考に適当に作成したものです。 -- [[谷村]...
#comment
**日本人の活躍 [#ted6da05]
>[[谷村]] (2004-05-14 (金) 21:47:25)~
~
汎用空間データ関数パッケージ(sp)のヘルプを何気なく見たら~
holepolys Hisaji Ono's hole problem~
を発見！尾野さん、素晴らしいです。RでGISも期待しています。~
//
#comment
**carのrecode(条件反射でコードを書くと) [#t1afe2f2]
>[[なかま]] (2004-05-11 (火) 08:47:51)~
~
ぼやーっと、~
x <- 0:100
y <- recode ( x, '0:19 = "未成年" ; 20:59 = "成人" ; els...
のようにしたものの、18才未満の方がいいかなと、~
recode ( x, 'x[ x < 18] = "未成年" ; x[ x < 60] = "成人"...
とかしてしまい、~
recode ( x, 'x[ x <= 12] = "小児" ; x[ x < 18] = "未成年...
x[ x < 60] = "成人" ; else = "老人" ')
とした時点でパーサーエラー。~
'='でsplitしているのでエラーになったんですね。~
で、条件式を書いてもOKとは書いてない事に最後に気づいてし...
# 良く読め<オレ~
でもそれって、~
x <- 0:10
recode ( x, '0:5 = ";" ; 6:8 = ":" ; else = "x"' )
もエラーになる訳で。。。~
//
-地道にcutとlevelsでやれば，最後の例はエラーにはなりませ...
-ありがとうございます。地道に精進致します。 -- [[なかま]]...
#comment
**dynamicGraphについて（Vineユーザのみなさま) [#j2b6880a]
>[[谷村]] (2004-05-10 (月) 16:43:06)~
~
dynamicGraphはtcltkに依存しています。Vineオリジナル版では...
//
-仮に有効にしてもそのままではPostScriptに吐いてもWinでもL...
-あっなかまさん。VinePlusへR-1.9日本語版をputする件はいか...
-うぐっ、もう少し、心の準備が出来てから。。。(汗; -- [[な...
#comment
**カンサス州知事は R ユーザー？ [#dbd0446b]
>[[間瀬]] (2004-04-29 (木) 14:41:20)~
~
最近の r-help ML にカンサス州知事 Kathleen Sebelius から...
//
#comment
**車輪を二度発明すること [#ve47fe05]
> (2004-04-23 (金) 18:00:47)~
~
splancs に areapl というのがあってポリゴンの面積を計算で...
~
ポリゴンの面積を計算するのは，R だとものすごく短く書ける。~
頂点の座標を時計回りに順にn×2行列にして与える。~
area <- function(xy)
{
n <- length(x <- xy[,1])
sum((c(x[-1], x[1])-c(x[n], x[-n]))*xy[,2])/2
}
# test run. answer=28.5
xy <- matrix(c(9,4,7,5,6,4,7,2,5,1,2,3,3,4,1,5,1,7,3,7,4...
> area(xy)
[1] 28.5
//
-areaplは、Fortranでコンパイルされた関数を呼んでいますが...
-30時間かかる計算が1時間でできると，きっとうれしい。30秒...
#comment
**全世界の R ユーザー数を推定する (r-help 記事より) [#seb...
> (2004-04-20 (火) 08:14:42)~
~
I have been trying to determine the size of the R user ba...
~
1. It appears that there are about 100,000 S-PLUS users.~
~
Rationale: According to Insightful's 2002 Annual Report, ...
~
Conclusion: S-PLUS costs $3500 (Windows) or $4500 (Linux/...
~
Does anyone has any other information that would give som...
~
2. At least one R book has achieved sales of just over 5,...
~
Conclusion: Very few books sell to more than 12% of the ...
~
3. There are a total of 3225 unique subscribers to the th...
//
-r-help はこの件で盛り上がっています。やはり、皆気になっ...
#comment
**R 1.9.0 正式公開 [#v48c431e]
>[[間瀬 ]] (2004-04-12 (月) 22:20:15)~
~
R 1.9.0 が正式に公開されました。各バイナリーの公開はすこ...
//
-lfactorial も同時に設けられたそうです。統計プログラムを...
~
> factorial(10)
[1] 3628800
> lfactorial(10)
[1] 15.10441
> x <-matrix(1:4, nc=2)
> addmargins(x)
[,1] [,2] [,3]
[1,] 1 3 4
[2,] 2 4 6
[3,] 3 7 10
-VineLinux用のrpmを本家にputしました。本家及びミラーに反...
#comment
**emacs-21とessと日本語 [#q0d91554]
>[[なかま@金沢市民]] (2004-04-11 (日) 18:09:50)~
~
ヘルプで日本語が出ない。essでの日本語の扱いで旨く行かない...
ありましたが、解決出来ました。~
~
(set-language-environment "Japanese")
(set-default-coding-systems 'euc-japan-unix)
(set-terminal-coding-system 'euc-japan-unix)
(set-keyboard-coding-system 'euc-japan-unix)
(set-buffer-file-coding-system 'euc-japan-unix)
(require 'ess-site)
(define-key ess-mode-map "?177" 'delete-char)
(setq ess-ask-for-ess-directory nil)
(setq ess-pre-run-hook
'((lambda () (setq S-directory default-directory)
(setq default-process-coding-system '(euc...
)))
//
-setq ess-pre-run-hook に (setq default-process-coding-sy...
#comment
**The R Package for multivariate and spatial analysis [#e...
> (2004-04-07 (水) 07:04:15)~
~
r-help 記事によれば [[The R Package for multivariate and ...
//
-個人的にはGoogleにかかりやすいようにRの名前の方を変えて...
-私は「R 統計」と，"統計"をかませて検索します． -- &new{...
#comment
**円周率はなぜπか (オフトピック) [#ce3cc40f]
>[[間瀬茂]] (2004-04-06 (火) 06:40:47)~
~
「工学のためのデータサイエンス入門」の23頁のコラム「円周...
であると書きましたが、これは不正確でした。確かに π が円周...
オイラーの1748年出版の広く読まれた著書「無限小解析入門」...
既に1647年に π を円周率を表す記号として使った例があるそう...
p とか書いていたそうです。詳しくは「π - 魅力の数」、ジャ...
を見てください。(若干一名から好評でしたのでずにのっておし...
分母が１より大きいすべての有理数 p/q に対し次の不等式
|π-p/q| ＞ 1/q^{20.6}
があるからπはランダムではない、と書いており、びっくりした...
( 10^{-5} ?sum_{n=-?infty}^{?infty} exp(-n^2/10^{10}) )^2
や、C言語でπを2400桁計算できる158語のプログラムなど、面白...
//
-考えてみたら上の無限級数は N(0, 10^{10}/2) の密度関数の ...
-細かいことで申し訳ありませんが、書名、出版社名、発行年が...
-失礼、題名と出版社はおっしゃる通りです。2003年に第二刷が...
#comment
**survey package [#rb516e5c]
>[[あいだま]] (2004-04-03 (土) 16:08:23)~
~
Survey Packageの翻訳が不適切です。↓~
階層化、クラスター化、非一様標本調査に対する要約統計量、...
~
Summary statistics, generalised linear models, and genera...
~
層化抽出、集落抽出、等確率でない抽出の場合の標本のための~
記述統計、一般化線形モデル、pseudolikelihood estimationの...
テイラー近似による分散の推定、replicate weights、~
//
//-自己レスで失礼しますが。replicate weightsの訳は当方知...
-自己レスで失礼しますが。replicate weightsの訳は当方知り...
-どこの訳を指していらっしゃるのでしょうか、情報不足です。...
-replicate weights はただの標本抽出用のウェイトではなくて...
#comment
**ご意見をお聞きします [#i7514b06]
>[[間瀬茂]] (2004-04-03 (土) 10:28:59)~
~
間瀬です。私どもが企画した「工学のためのデータサイエンス...
//
-趣旨に合えば問題ないと思います。実際、個人Wikiサイトに別...
-お騒がわせしましたが、結局複雑な数式を表現するには latex...
-[[工学のためのデータサイエンス入門の正誤表:http://www.is...
#comment
**The R Reference Manual 出版！ [#sdf77e03]
>[[間瀬茂]] (2004-04-03 (土) 00:32:09)~
~
~
R のベースパッケージ中のすべての関数の参照マニュアルを二...
//
#comment
**RのSVMで学習した結果を外部ファイルにエクスポート [#ddf7...
>[[高階知巳]] (2004-03-28 (日) 11:11:54)~
~
ウィーン工科大学の授業で使うために作られたというe1071パッ...
//
-e1071_1.4 が先日CRANでリリースされ、その中に、上述のパッ...
#comment
**rimageのWindows版バイナリが公開されました [#ra12e198]
>[[高階知巳]] (2004-03-28 (日) 11:01:17)~
~
RでJPEG画像ファイルの読み込みと基本的な画像処理機能を提供...
//
- Windows ユーザーです．早速パッケージを頂きました．examp...
-Rでは画像の表示は遅いですが、特にWindowsでは遅い様です。...
-高階さんが飼ってらっしゃる猫でしたか^ ^．すみません，猫...
-高階知巳さん、日本発 R アドオンパッケージ第一号(?)嬉しい...
-間瀬先生、ありがとうございます。私もパッケージ作者が増え...
#comment
**Nepia の対数 （オフトピック） [#j5aaadbf]
>[[間瀬茂]] (2004-03-27 (土) 09:41:08)~
~
また年寄りの懐古趣味記事です。統計学でも必須の対数関数の...
x ＝ 10^7*(1-10^{-7})^y
という関係にある数 x.y を数表化したのだそうです。比例定数...
//
-「大数の法則」のお話共々，興味深く読ませて頂いております...
-おっしゃる通り、私供の本「工学のためのデータサイエンス入...
#comment
**Knoppix/Math/2004 [#a1993610]
>[[間瀬茂]] (2004-03-24 (水) 16:28:29)~
~
R-jpに岡田さんが投稿されていますが、3月の筑波大学の数学会...
//
-Ring Server で公開をはじめました。収録しているRは日本語...
-[[Knoppix/Math/2004:http://unit.aist.go.jp/it/knoppix/is...
#comment
**useR! 2004 プログラム公開 [#g85fe9a8]
> (2004-03-23 (火) 19:32:17)~
~
日本人(?)による講演も三組。 [[プログラム:http://www.ci.tu...
//
#comment
**R 1.9.0 alpha [#a4aba883]
>[[間瀬茂]] (2004-03-21 (日) 07:02:30)~
~
r-help 記事によると R 1.9.0 alpha を試した人から、MCMC 法...
//
-1.9.0ではアンダースコアを変数名に使えるようになるとか，...
#comment
**Mac OS X Panther と Raqua （r-help 記事より） [#j6697ecd]
> (2004-03-12 (金) 17:22:45)~
~
ところで Raqua とは Mac OS X 用の R のことですよね(か？)
I've been nicely surprised, glancing through O'reilly's ...
for Unix Geeks" (by Brian Jepson and Ernest E. Rothman),...
the sixth section of chapter seven is dedicated to Raqua...
a few lines to describe R and Stefano's port to mac OS X...
(graphic windows and R console), and even a small apples...
interact with Raqua.~
BTW, the book is great.~
~
--Mathieu~
//
- => 青木さんの頁で確認して安堵。ところで Raqua というネ...
-aquaはMac OS XのGUIで、これを使ってＲが動くからだと思い...
-老婆心ながら、1.9 以後から単にRと呼ばれるようになりまし...
#comment
**R マニュアル翻訳チーム結成の呼びかけ（ただしフランスの...
> (2004-03-12 (金) 07:15:42)~
~
r-help にフランス人から投稿があり、前から自分が進めている...
//
#comment
~
**大数の法則 (オフトピック) [#z8471938]
>[[間瀬茂]] (2004-03-11 (木) 17:09:54)~
~
確率・統計の本の最初には必ずでてくる「大数の法則」ですが...
//
#comment
**eps ファイルの編集方法 [#tdc076bf]
>[[舟尾]] (2004-03-07 (日) 18:53:38)~
~
本筋から離れるかもしれませんが，皆さんは R で出力した eps...
//
-恥ずかしながら、(e)ps ファイルを編集するなぞ考えもしませ...
-UNIX系のOS (と Mac OS X) には [[Tgif:http://www.tobata.i...
-有料ですが、[[花子2004:http://www.justsystem.co.jp/hanak...
-花子が eps 形式に対応したというのは初耳でした．[[Canvas:...
-凄い技ですね。ぼくは手作業編集したいときはwin.metafile()...
-花子2004 は eps ファイルの編集が出来るんですか！かなり興...
-あと，[[ここの掲示板:http://www.visio.jp/bbs/bbsView.asp...
-ついでですが、Knoppix-jp の Openoffice はこれまで何故か...
-既存の eps も，[[pstoedit:http://www.pstoedit.net/pstoed...
-度々、枠領域を明示・・・ではなくて、asp=??を指定できる関...
#comment
**R on Zaurusにてグラフ出力時のエラー [#wf17f80c]
>[[aki]] (2004-03-07 (日) 06:19:07)~
~
現在、Zaurus上のX/Qtで動作するRを用意しています。~
コンパイルとipkの用意はできたのですが、いざ起動して~
> x<-rnorm(50)
> y<-rnorm(50)
> plot(x,y)
とすると、~
Warning messages:
1: Nonfinite axis limits [GScale(nan,nan,1, .); log=0]
2: relative range of values =hogehoge * EPS, is small (a...
3: Nonfinite axis limits [GScale(-inf,inf,2, .); log=0]
4: relative range of values =hogehoge * EPS, is small (a...
というメッセージと共にへんてこりんなグラフが表示されてし...
検索すると、~
http://lists.debian.org/debian-arm/2003/debian-arm-200311...
http://maths.newcastle.edu.au/~rking/R/devel/03b/1083.html~
が見つかりましたが、結局のところ解決法が分かりません。~
http://sgp.zaurii.net/binaries/index.html#R~
のRは使えるような話が[[LinuxザウルスSL-C700でR]]には載っ...
どなたか解決法が分かる方はいらっしゃいませんか?~
//
-ATLAS http://math-atlas.sourceforge.net/ が必要?かもしれ...
-以前、Ｘ外した１７１postscriptバージョンを作成した時も同...
-度々、枠領域を明示・・・ではなくて、asp=??を指定できる関...
-BLASとLAPACKを用意してみましたが症状変わらず。うーん、あ...
-FORTRANコンパイラをg77からf2cにしてみましたが今度はエラ...
-f2cでコンパイルしましたがエラー症状変わらず。というわけ...
#comment
**高校生向け R [#i2a9515d]
> (2004-03-07 (日) 00:25:48)~
~
なるほど、高校生に R を使わせるとは考えつきませんでした。...
//
-はじめまして。中学・高校現場にいる者です。実はRを利用し...
-よろしくお願いします。 -- &new{2004-03-12 (金) 22:04:11};
-もしかすると、R は統計学以上に高校生用の情報教育のプラッ...
-職場での公開について許諾が得られたので、PDF形式で公開い...
-すごいですね。数学の確率・統計の一部として「R を使えばこ...
-「データを解析する経験をしてみよう。そのときのツールとし...
-だんだん形を整えていって、「記述統計学の基礎（平均・分散...
-全員対象なので「賛否両論」というところだったでしょうか。...
-ちなみに（内容をみてもわかると思うのですが）あらかじめこ...
#comment
**Statistiques avec R [#n573321b]
>[[間瀬茂]] (2004-03-04 (木) 21:08:48)~
~
フランス語ながら「Statistiques avec R」というマニュアルは...
//
-リンク集にも書いてありますが，特に4章のグラフィックスに...
-他の章も一読あれ。やはりグラフィックスがたくさんあります...
-r-help でも好評のようです。一人が反応メイルで arigato go...
-r-help ではまさにニュースになったようです。これの英語訳...
-著者の Zoonekynd さんはどうやら計算機オタクの数学者らし...
-Zoonekynd という人物に興味を持ち少しグーグルしてみました...
--[[TeX用フォント一覧:http://tex.loria.fr/fontes/zoonekyn...
--[[Metapost による各種図形・図式作図例:http://tex.loria....
--[[難易度順漢字一覧表:http://zoonek.free.fr/kanji/list.p...
--[[Gimp マニュアル:http://zoonek2.free.fr/UNIX/23_perl_g...
-Zoonekynd さんはやはり日本語を理解するようです。中級程度...
#comment
**「データサイエンス入門」出版案内 [#ac86cc3d]
>[[間瀬茂]] (2004-03-03 (水) 00:15:21)~
~
既に少しお知らせしましたが、第二の R 本が三月末に刊行予定...
//
-楽しみにしています。目次がわかるwebサイトはまだないので...
-最寄の本屋の在庫が切れないか…，今から心配しております． ...
-追加コメント。原稿作成中、ますます R の魅力にはまりまし...
-「慶応大学のある先生」とはどなたでしょうか？少し興味があ...
-[[慶応大学のある先生:http://web.sfc.keio.ac.jp/~watanabe...
-ありがとうございました -- [[Yano]] &new{2004-03-10 (水) ...
#comment
**表計算ソフトの統計機能 [#lea93164]
> (2004-02-28 (土) 19:20:21)~
~
Excel の統計機能の問題については青木先生のホーム頁に詳し...
//
-OpenOffice.orgのcalcも，相関係数のエラー（[[青木先生のペ...
-やはり駄目ですか。悪いところも含めて Excel にできるだけ...
-例の Excel をこきおろすページ(^_^;)というのは，無理矢理...
#comment
**Vine Linux 版の R のインストール法の記事? [#l3505498]
> (2004-02-25 (水) 17:36:10)~
~
Vine Linux 版の R のインストール法の記事がいつのまにか消...
//
-私が書いたネットワークインストールの記事でしょうか。apt...
# rpm -ivh R-1.8.1-4vl1.i386.rpm
#comment
**S/R で商売 [#q59cc20e]
>[[間瀬茂]] (2004-02-25 (水) 17:11:06)~
~
r-help に流れた広告。既に日本にも例があるようですが、R の...
XSolutions Corp ([1]www.xlsolutions-corp.com) is proud...
a 2-day "Advanced R/Splus programming" taught by R Dev...
Core Team Guru.
*********San Francisco ----- March 22-23
*********Washington DC ---- March 22-23
*********Chicago ----------- March 25-26
*********New York --------- March 22-23
*********Atlanta ----------- March 25-26
*********Boston ------------ TBD
If your location of choice isn't listed, please let us...
Reserve your seat Now at earlybird rates (payment due...
class)
Registration:
Email Sue Turner: [2]sue@xlsolutions-corp.com
Phone: 206-686-1578
Visit us: [3]www.xlsolutions-corp.com/training.htm
Course outline:
- Overview of R/S fundamentals: Syntax and Semantics
- Class and Inheritance in R/S-Plus
- Concepts, Construction and good use of language obje...
- Coercion and efficiency
- Object-oriented programming in R and S-Plus
- Advanced manipulation tools: Parse, Deparse, Substit...
- How to fully take advantage of Vectorization
- Generic and Method Functions; S4 (S-Plus 6)
- Search path, databases and frames Visibility
- Working with large objects
- Handling Properly Recursion and iterative calculations
- Managing loops; For (S-Plus) and for() loops
- Consequences of Lazy Evaluation
- Efficient Code practices for large computations
- Memory management and Resource monitoring
- Writing R/S-Plus functions to call compiled code
- Writing and debugging compiled code for R/S-Plus sys...
- Connecting R/S-Plus to External Data Sources
- Understanding the structure of model fitting functio...
- Designing and Packaging efficiently
Early-bird group research fee: $995!
This course will also deal with lots of S-Plus efficie...
any special topics from participants is welcome.
Please let us know if you and your colleagues are inte...
class to take advantage of group discount. Register no...
your seat!
Cheers,
Elvis Miller, PhD
Manager Training.
XLSolutions Corporation
206 686 1578
[4]www.xlsolutions-corp.com
References
1. http://www.xlsolutions-corp.com/
2. mailto:sue@xlsolutions-corp.com
3. http://www.xlsolutions-corp.com/training.htm
4. http://www.xlsolutions-corp.com/
//
#comment
**Windowsでパッケージを自作 [#qaf45149]
>[[かくれR狂]] (2004-02-25 (水) 04:21:27)~
~
Windowsでパッケージ自作に挑戦していました。~
山口県立大の中澤先生の[[統計処理ソフトウェアRについてのTi...
必要なソフト探しからだったので結構大変でした。~
Tipsに紹介しようかと考えていますが、~
もうすでにどこかで紹介されているのでしょうか？~
//
-日本語での解説は基本マニュアルにかかれている以外は知りま...
-Ｔｉｐｓに載せてみました。どなたか挑戦してみてください。...
-高階さんのrimageがあるので，第二号ですね。ぼくがMsysと角...
#comment
**新パッケージ mvpart (rpart の多変量データ版) 2004.2.22 ...
> (2004-02-23 (月) 07:15:05)~
~
The package mvpart is now available.~
mvpart includes partitioning based on (1) multivariate nu...
(2) dissimilarity matrices.~
~
The package mvpart is a modification of rpart --~
-- authors of original: Terry M Therneau and Beth Atkinso...
Includes some modified routines from vegan -- Jari Oksane...
//
-mvpartと同様の multivariate regression trees のソースコ...
#comment
~
**「Rを味見してみる」 [#h8fcfce3]
>[[舟尾]] (2004-02-21 (土) 15:02:56)~
~
R界で(も)お馴染みの青木先生が，またもや凄まじいコンテンツ...
[[R を味見してみる:http://aoki2.si.gunma-u.ac.jp/R/R0.htm...
これで「Rというソフトがあるらしいがインストールするのは…...
//
-すみませんでした。諸般の事由で，中止致します。おさわがせ...
-最近気づきましたが[[中間さんち:http://r.nakama.ne.jp/Rwe...
-青木先生のコンテンツが無くなってしまってげんなりしていた...
-Javaアプレットが見れる状態にして、しばらく待つと見ること...
-Rweb の html 15行目にエラーがありませんか？ -- &new{200...
-15行目直しました。「R を味見してみる」で覗いてたのは私で...
-周りの人に紹介しています．なかなか好評でありまして，R を...
#comment
**ブッシュ大統領曰く、私は統計学者でも預言者でもない。 [#...
>[[間瀬]] (2004-02-19 (木) 18:44:42)~
~
r-help 記事によれば、ブッシュ大統領のスタッフが、今年末ま...
//
#comment
**舟尾チップス集 [#q1c11030]
>[[間瀬]] (2004-02-13 (金) 22:06:29)~
~
舟尾チップス集にはまだ私も知らない話題がたくさんあります...
//
-私は[ファイルへのデータ出力例]の節での，ファイル出力の際...
-お気づきの「へー」はどうぞ適当な箇所に追加お願いします。...
-ここに逐次書き込みをするのかと勘違いしました。お詫びする...
#comment
**useR! 2004 - The R User Conference [#k4bcc755]
> (2004-02-12 (木) 15:07:36)~
~
2004年5月20日‐22にウィーンのTechnische大学でRユーザカンフ...
開催されます。演題締切は2月15日ですのでお早めに。~
http://www.ci.tuwien.ac.at/Conferences/useR-2004/~
//
#comment
**投稿記事 どうむ について [#s2951c8c]
~
> test(.Machine$double.eps)
TRUE
> test(.Machine$double.eps/2)
TRUE
となりました。
業務でRを使用しているのですが、if(Cond)のCondにそもそもこ...
とりあえずご報告まで。
//
-こうした記事は独立して投稿すると、すぐ行方不明になります...
> .Machine$double.eps/1e10
[1] 2.220446e-26
と .Machine$double.eps 以下の正実数も立派に存在します。次...
> if(1+.Machine$double.eps==1) cat("TRUE?n") else cat("F...
FALSE
> if(1+.Machine$double.eps/2==1) cat("TRUE?n") else cat(...
TRUE
-- &new{2004-02-02 (月) 17:19:39};
-申し訳ございません -- どうむ &new{2004-02-02 (月) 17:39:...
-「if (x) の不審な(?)行動」に対するRESだったのですが当方...
-試してみたら.Machine$double.eps/1e308 が零になりました。...
> if(.Machine$double.eps/1e307) cat("TRUE?n") else cat("...
TRUE
> if(.Machine$double.eps/1e308) cat("TRUE?n") else cat("...
FALSE
-どうむさんのご指摘のように「if (x) の不審な(?)行動」の記...
#comment
**文字列のオブジェクトへの変換方法 => 追加コメントは Q&A ...
>[[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~
つまり、文字列として取り込んだものをオブジェクトとして扱...
いろいろ探してみたのですが、見つかりませんでした。~
どなたか教えてもらえないでしょうか？~
-getを使うとオブジェクトの中身が取り出せる。 money[get(a)...
#comment
**too large data! [#la36eaf9]
>[[間瀬]] (2004-01-29 (木) 17:36:59)~
~
r-help 記事より。R でどれだけ大きなデータが処理できるかは...
R-1.8.1,
compiled as a 64bit application for a Solaris 5.8, 64 bi...
8Gb of RAM を一人で使える
データは 2,000,000 行超で、多くは因子の60変数、一部は連...
ファイルからの読み込みとセーブに5分以上、再読み込みに15...
S-plus では可能だった解析を行なった
一回の反復で45分後に R はすべてのメモリを使い果たし、ス...
メモリ使用は最大を指定している（はず）
実際の解析コードは以下の通り
library(package = "statmod", pos = 2,
lib.loc = "/home/jeg002/R-1.8.1/lib/R/R_LIBS")
qc.B3.tweedie <- glm(formula = pp20B3 ~ ageveh + anpol +
categveh + champion + cie + dossiera +
faq13c + faq5a + kmaff + kmprom + nb...
rabprof + sexeprin + newage,
family = tweedie(var.power = 1.577,
link.power = 0),
etastart = log(rep(mean(qc.b3.sans.o...
'pp20B3']), nrow(qc.b3.sans.occ))),
weights = unsb3t1,
trace = T,
data = qc.b3.sans.occ)
コアメンバーの Dalgaard から次のようなコメント
データはそれ自体で 0.5-1GB になるはず
60変数のうち実際使っているのは15変数だから、最初からそれ...
因子の水準が多過ぎる可能性？
とくに連続変数を因子に変換すれば違った結果が得られるかも...
-再び r-help 記事より。R の現在の基本ライブラリーは最大利...
-上記記事に関連して。６４ビットCPU の Linux 機では１６GB...
#comment
** 不適当記事一件削除 [#abb2ca4d]
**Rnews 2003/3 公開 [#ud22374c]
> (2004-01-22 (木) 18:29:42)~
~
We have published the 2003/3 issue of R News on
http://cran.R-project.org/doc/Rnews
where you can download the newsletter as PDF or Postscri...
will propagate to the CRAN mirrors within a day or two.
Contents of this issue:
Dimensional Reduction for Data Mapping
R as a Simulation Platform in Ecological Modelling
Using R for Estimating Longitudinal Student Achievement ...
lmeSplines
Debugging Without (Too Many) Tears
Th R2HTML Package
R Help Desk: Package Management
Programmer's Niche: Regular Expressions
Recent Events
Upcoming Events
Changes in R 1.8.1
Changes on CRAN
R Foundation News
//
#comment
**直線回帰の信頼区間 [#a8cb8751]
>[[ぶんしょ]] (2004-01-18 (日) 19:21:24)~
~
回帰させた直線の95％信頼区間を知りたいのですが、どのよう...
//
-predict()関数を使うのが便利です。直線回帰でいいなら，exa...
-ありがとうございます。さっそく試してみます。 -- [[ぶんし...
#comment
~
**[質問] 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
-丁寧な回答をありがとうございます。ただ、私の記事は質問コ...
-どういたしまして。依然としてなぜ lm オブジェクトを R プ...
-いえいえ、私のHPは間瀬さんが翻訳されたR入門や他の方々の...
#comment
**ESRI shape ファイル [#ve85d9d7]
>[[なかま＠金沢市民]] (2004-01-02 (金) 18:37:18)~
~
あけまして、おめでとうございます。~
http://www.diva-gis.org/data/DataServer.htm~
に、色んな国のがありました。~
日本のだと、都道府県レベルまででしたが、~
plot(j[j$ADM2 == "Ishikawa"])
とか、最低限は出来そう。~
//
#comment
ページ名:
![[RjpWiki]](image/../Rlogo.png "[RjpWiki]")
