なんでも掲示板アーカイブ（４）
をテンプレートにして作成
[
トップ
] [
新規
|
一覧
|
検索
|
最終更新
|
ヘルプ
]
開始行:
COLOR(blue){SIZE(20){RjpWiki自体へのコメントでも、雑談で...
COLOR(red){注意}：新規記事用の入力欄は以下の目次の直後に...
----
COLOR(red){過去の記事のアーカイブ} 。アーカイブ中の記事へ...
-[[なんでも掲示板アーカイブ（１）]] (元記事が 2003-07-17 ...
-[[なんでも掲示板アーカイブ（２）]] (元記事が 2003-11-27 ...
-[[なんでも掲示板アーカイブ（３）]] (元記事が 2004-01-01 ...
----
#contents
----
回答を期待する記事は原則 COLOR(red){質問コーナー [[Ｑ＆Ａ...
~
#article
**経済学における統計的有意の記号 [#odc45137]
>[[y.a.]] (2004-11-23 (火) 01:01:55)~
~
[[Ｑ＆Ａ (初級者コース)]]にて，経済学の記号のつけ方の慣習...
結論から言うと，人やジャーナルによって異なるものの，特に...
~
ただ，~
cutpoints= c(0, 0.01, 0.05, 0.1), symbols = c("***", "**...
記法は，少なくとも経済学の分野には少なからず存在はすると...
~
以下，経済学のトップジャーナル群の論文，掲載予定論文から...
~
1) t値を示して星はつけない（でもp値はつけない）派
- [[American Economic Reviewの例(p.1396, Tab.6):http://fa...
- [[Quarterly Journal of Economicsの例(p.1298, Tab.IV):ht...
~
2) (1,5,10) -> (***,**,*) 記法前近代派
- [[Review of Economic Studiesの例(p.22, Tab.6):http://ww...
- その他，[[Journal of Development Economics誌:http://www...
//
-統計学の大口ユーザーである心理学等でも類似の慣習があるの...
-うちの学部の経済関係の先生も，ｔ値は書くけどＰ値は書かな...
某誌（医学系）に載ったチュートリアルで
＋ または △： p < 0.1 傾向あり，傾向差
＊： p < 0.05 有意差あり
＊＊： p < 0.01 明らかな差あり
＊＊＊： p < 0.001 著明な差あり
というのがあり，＋とか△はいかがなものか？とおもったら，心...
へ〜。郷に入らば郷に従えか -- &new{2004-11-25 (木) 18:58...
-Rを離れて雑談になってしまうのですが、「統計的検定（は使...
#comment
**R からキーワード検索する関数 (from r-help) [#f64c2f99]
>[[QDU]] (2004-11-20 (土) 20:51:33)~
~
今日の r-help に面白い関数が紹介されていましたので紹介し...
help.search.archive<-function(string){
RURL="http://www.google.com/u/newcastlemaths"
RSearchURL=paste(RURL,"?q=",string,sep='')
browseURL(RSearchURL)
return(invisible(0))
}
たとえば help.search.archive("plot3d") を試してみて下さい。
//
-やはり r-help 投稿の同様の関数： -- [[QDU]] &new{2004-11...
Inspired by the functions that Barry Rawlingson and Dave...
searching Rwiki and R-help archive, I've made up a funct...
search on Prof. Baron's site (Thanks to Prof. Baron's he...
query string!):
RSiteSearch <- function(string, restrict="Rhelp", format...
sortby="score", matchesPerPage=1...
URL <- "http://finzi.psych.upenn.edu/cgi-bin/htsearch"
qstring <- paste(URL, "?config=htdigrun1", sep="")
## replace spaces with "%20" in the query
string <- paste("words=", gsub(" ", "%20", string), ...
mpp <- paste("matchesperpage=", matchesPerPage, sep=...
format <- charmatch(format, c("long", "short"))
if (format == 0) stop("format must be either long or...
format <- paste("format=builtin-", switch(format, "l...
sortby <- charmatch(sortby, c("score", "time", "titl...
if (sortby == 0) stop("wrong sortby specified")
sortby <- paste("sort=",
switch(sortby, "score", "time", "tit...
sep="")
res <- charmatch(restrict, c("Rhelp", "doc", "functi...
if (res == 0) stop("wrong restriction specified")
res <- switch(res, "Rhelp00/archive|Rhelp01/archive|...
"finzi.psych.upenn.edu/R/doc",
"finzi.psych.upenn.edu/R/library",
"finzi.psych.upenn.edu/R/doc|finzi.psy...
res <- paste("restrict=", res, sep="")
qstring <- paste(qstring, res, format, sortby, string...
browseURL(qstring)
invisible(qstring)
}
The options roughly correspond to the options on that se...
some people find this somewhat helpful. Comments/sugges...
improvement are much appreciated.
Best,
Andy
#comment
**The hidden costs of GPL softwares [#f6ce8d13]
>[[間瀬茂]] (2004-11-18 (木) 22:54:56)~
~
今日の R-help は表題の投稿への反応で盛り上がりました。い...
//
- 他のソフトと比較するのは良くないことなのかもしれません...
ヘルプが充実していて使いやすいといえば Mathematica ですが...
ヘルプに関して云えば，R のヘルプは関数の構文と引数の説明...
これは私の勝手な意見ですが，「学習コストが高い」のは「か...
-ある集会で「SASの日本語化をしましょう」という趣旨の講演...
-SASのGUIメニューや基本的な計算のアウトプットは日本語化さ...
-船尾さんがおっしゃるEXCELのクッキング本とつながりますが...
参考までに，上田・ 苅田・本田(2003) EXCEL徹底活用多変量解...
-linux,*BSD,glibc,readline,gzip,bzip,XOrg,gcc,libgoto...R...
#comment
**S-Plus や R の関数のカテゴリー化された検索システム [#bf...
>[[間瀬茂]] (2004-11-18 (木) 17:42:15)~
~
r-help に S-Plus や R の膨大な関数を用途別にカテゴリ化し...
~
[[S-Plus/R Function Finder:http://biostat.mc.vanderbilt.e...
//
#comment
**R 2.0.1 for VineLinux 3.0 [#lc1a3bf1]
>[[谷村]] (2004-11-16 (火) 20:39:12)~
~
本家(CRAN)にR-2.0.1-0vl1.i386.rpmをputしました。これはオ...
//
#comment
**どこかでみたような [#fc6abf97]
>[[なかま]] (2004-11-14 (日) 00:55:59)~
~
[[中文TeX:http://www.ctex.org/forums/index.php?showtopic=...
ユーザー数ってどれくらいいるんでしょ.~
//
-びっくりしましたね。ただ、先月中国の集会に参加した時印象...
-本当にすみませんでした。この文章(Rを薦める100理由)は、私...
-消す必要はないと思いますが(^^;。中国語になってる！のが面...
-どうせこちらも「無断」引用しまくっていますから、お構い無...
-日本語のサイトを中国の人が読めるということは、逆もまた真...
-批評してる(?)のは日本人のような気がしますが...(^_^; -- [...
-よしよしさんの「良くない」っていう発言は、不偏推定値では...
#comment
**Ｒで物理乱数 [#v53f4e78]
>[[垂水共之]] (2004-11-10 (水) 17:21:09)~
~
[[FDK:http://randomstreamer.ntl.co.jp]]という会社がUSB接...
Windows版のRで使えるようにDLLも作ってもらいました。~
Rユーザーの方はサポートCD-ROM特別サンプルも手に入れてくだ...
//
-おお、すばらしいですね! -- &new{2004-11-12 (金) 08:16:3...
-はたして、お小遣いで買える物なのでしょうか...(^^; -- [[...
-小遣いではちと高いですが、3万程度と聞いています -- [[垂...
-set.seed(xxx) は原理的にできないわけですね。 -- &new{20...
-生成速度がHPには公表されていませんがどれくらいなのか一寸...
-ありがとうございます。3万がんばり(なにを?)ます。 -- [[な...
-内蔵してますRPG100の乱数生成速度は250kbits/sのようですよ...
-ちなみに、販売しているのは「日本テクノ・ラボ（株）」見た...
-Random StreamerはSPLUSで使える様です。 -- [[ななし]] &ne...
#comment
**かんたんな入力チェックを入れるのは [#l3f27661]
>[[なかま]] (2004-11-08 (月) 01:21:09)~
~
--- article.inc.php.orig
+++ article.inc.php
@@ -189,7 +189,24 @@
$article_rows = article_ROWS;
$article_cols = article_COLS;
$string = <<<EOD
-<form action="$script" method="post">
+<script>
+ <!--
+ function isNull(elem,jname){
+ if(elem.value==''){
+ elem.focus();
+ alert(jname+'がはいっていると嬉しいかも.');
+ return true;
+ }
+ return false;
+ }
+ function chkArticle(myform){
+ if(isNull(myform.name,'お名前'))return false;
+ if(isNull(myform.subject,'題名'))return false;
+ return true;
+ }
+ -->
+</script>
+<form action="$script" method="post" onSubmit="return c...
<div>
<input type="hidden" name="article_no" value="$articl...
<input type="hidden" name="plugin" value="article" />
~
まあ,気休めですが効果はあるかも.~
//
-ありがとうございます。早速導入してみました:-) -- [[岡田]...
-無題が減るといいですね.:-) -- [[なかま]] &new{2004-11-13...
#comment
**R とアメリカ大統領選挙 [#mac5766b]
> (2004-11-03 (水) 20:00:03)~
~
今回も米大統領選挙は接戦ですね。~
~
今回はオハイオ州が注目の的ですが、R の car パッケージに...
//
#comment
**日本人によるアドオンパッケージ Zelig [#zd08ed8f]
> (2004-10-23 (土) 14:48:44)~
~
今日初めて気がつきましたが MNP というアドオンパッケージの...
//
#comment
**中国のR事情 [#p98b8c1a]
>[[間瀬茂]] (2004-10-19 (火) 23:28:09)~
~
集会で中国にいって来ました。桂林にある広西師範大学が会場...
//
-r-helpにあったとおりGB18030で表記できれば広まると思いま...
-データもグラフ出力も、そのまま中国語(gb2312)は大丈夫です...
-中国：Linuxは普及せず、ウイルスソフトで国産好調＝サーチ...
-うちのサイトに来る怪しげパケット,scanは中国が一番多いの...
-広西師範大学 って GNU なんですね。(^_^; -- [[なかま]] &n...
-そうなんです。Guangxi Normal University。広西チアン自治...
-PS. 学内をおっかなびっくり見物していたら、ひとりだけ茶髪...
-広西なんて中国の中でもかなり田舎なんだ！中国の地方差はす...
-PS.学食の値段も地方によって違うんだよ。上海の学食ならこ...
#comment
**Mac 版 R 2.0.0 [#a50884d8]
>[[青木繁伸]] (2004-10-18 (月) 11:27:03)~
~
まだ筑波大学のmirrorには現れていない~
working directry の設定で，メニューからやると設定できない...
注釈に日本語を使っている関数のリストを表示する（関数名だ...
一部パッケージがまだ対応できていない（MASS などが使えない...
環境設定メニューがない（なくてもいいようにはできていない...
~
というような具合で，R2.0.1の出現を待ちたい。~
//
-MASS などは，使えました。なぜ昨日使えなかったのか不明だ...
#comment
**R 2.0.0 のリリース近し？ => ついに公開されました（祝^祝...
>[[間瀬茂]] (2004-10-02 (土) 00:17:53)~
~
今日の r-help に R の 2.0.0 に関する問い合わせがあり、R ...
Or just wait for the release on Monday October 4th.
//
-素晴らしい！ついに R もバージョン 2 に突入ですねぇ．Wind...
-10/4(日本時間10/5の何時頃?)にリリース予定です。betaは既...
-*nix方面はなんとか構築。うーん、窓用たいへんかも...(T_T)...
-何かメジャーバージョンアップにふさわしい目玉があるのでし...
-表現が適切かはわかりませんが、地味に細かく*沢山*拡張した...
-パッケージ化待たなくともconfigure,make,make installでも...
-[[2.0 SERIES NEWS :https://svn.r-project.org/R/trunk/NEW...
-あのー、日本語バージョンはいつ頃に。(だから金沢に向かっ...
-どうせすぐ 2.0.1 もしくは 2.1 に代わりますから慌てる必要...
-出ました。http://cran.r-project.org/src/base/R-2 -- [[ak...
-リリースニュースによれば、この版は R に対する革新的な変...
As was the case with R 1.0.0, this new version represent...
more than a radical change to R. We do plan to celebrate...
with press releases and such.
-r-help でコアメンバーの Dalgaard 氏が次のようにコメント...
No, I don't think we have arrived (would we want to?). W...
have unresolved issues in byte compiling and event loops...
formal methods techniques are only just there now. There...
Bayesianfront, too - important even for us hardcore freq...
-大変遅くなりましたが、R-2.0.0-0vl1.i386.rpm for VineLinu...
-いつもデバッグ用に(ふと、やっちまった?と思う場合があるの...
#comment
**pukiwiki-mode [#xa4a4f21]
>[[谷村]] (2004-09-29 (水) 21:57:08)~
~
今日、[[pukiwiki-mode:http://www.bookshelf.jp/pukiwiki/pu...
~
ブラウザから書き込むのと比べて効率がよさそうなのですが、...
//
#comment
**Rwebに入力された驚くべき物 [#f5935bf2]
>[[なかま]] (2004-09-27 (月) 23:40:17)~
~
たまーにビックリするような物が実行されたりしますんですが...
~
[18時 20分 40秒] [緑茶] 1 個を手に入れました。
[18時 20分 41秒] [緑茶] 1 個を手に入れました。
[18時 20分 42秒] [銀杏] 1 個を手に入れました。
[18時 20分 45秒] [緑茶] 1 個を手に入れました。
[18時 20分 52秒] 経験値が 64 上がりました。
[18時 21分 13秒] 経験値が 104 上がりました。
[18時 21分 24秒] コンボで５度目の攻撃に成功すると、コン...
[18時 21分 25秒] 経験値が 64 上がりました。
[18時 21分 27秒] 経験値が 104 上がりました。
~
これをいったいどうしようと思ったのか非常に気になります(笑...
//
-ひょっとして Rweb という名前のオンラインゲームがあるので...
#comment
**sink 関数のネスト [#i28f3260]
>[[QDU]] (2004-09-27 (月) 21:15:06)~
~
たまたま 2ch の R 関連のやりとりを見ていて次のことに気づ...
~
'sink()' or 'sink(file=NULL)' ends the last diversio...
specified type). As from R version 1.3.0 there is a...
diversions for normal output, so output reverts to t...
diversion (if there was one). The stack is of up to...
connections (20 diversions).
~
# 参考例
> example(lm)
> sink("test.summary")
> sink("test2.summary") # sink 関数のネスト
> summary(resid(lm.D9) - resid(lm.D90)) # test2.summary ...
> sink() #
> summary(resid(lm.D9) - resid(lm.D90)) # test.summary ...
> sink()
> summary(resid(lm.D9) - resid(lm.D90)) # 標準出力に戻った
Min. 1st Qu. Median Mean 3rd Qu. ...
-7.772e-16 -6.245e-17 2.776e-17 1.405e-17 1.145e-16 ...
//
-出力を条件ごとに切替えて使うのによいかと。ネストと言うよ...
#comment
**Excel からコピーしたセル範囲を R に読み込む [#x99a1f1f]
>[[青木繁伸]] (2004-09-19 (日) 14:56:45)~
~
エクセルのセル範囲を R に読み込むために，今まではわざわざ...
from.excel <- function(nc)
{
matrix(scan(""), byrow=TRUE, nc=nc)
}
エクセルで必要な部分をコピーします。
x <- from.excel(2)
などとして，プロンプトが出たらペーストします。~
入力の終わりはMacintoshだとcommand+.~
汎用関数にする価値もないのでこのまま。
//
- ?scan によれば Linux だと CTRL+D, MS Windows では CTRL+...
-これはいいですね。 -- [[林 啓一]] &new{2004-09-27 (月) ...
- 参考になりました -- [[田中]] &new{2007-08-05 (日) 09:55...
#comment
//
**CRAN のすべてのパッケージを含むドキュメントからの全文検...
>[[間瀬茂]] (2004-09-12 (日) 22:35:51)~
~
r-help に「これこれするには何を使えばよいでしょう」タイプ...
[[R site search:http://finzi.psych.upenn.edu/search.html]...
//
#comment
**Rの宣伝講演をします [#cad3e233]
>[[谷村]] (2004-09-08 (水) 23:57:14)~
~
2004年9月10日東京大学本郷キャンパス山上会館で「空間疫学分...
http://www.csis.u-tokyo.ac.jp/japanese/research_activitie...
//
-素晴らしい！どのような反響があったのか，是非お聞かせ下さ...
-慶應三田で開かれたSVG Open と重なって行けなかったのが残...
-遅くなりました。「Rは名前は知っているけど空間情報科学分...
-空間分析ではかつては Splus と組んでいたESRI社の Spatial...
-ところで、リンク集にある谷村先生のご著書『保健医療のため...
-発行日は10月5日になっていますが、実質は9月27日頃の予定。...
-谷村先生、ご回答ありがとうございます。発売を期待しており...
-谷村さんに講演を依頼した者です。Rを知らない人にも「何が...
-今ごろ言っても遅いのですが、関西オープンソース2004で15分...
#comment
**大男総身に知恵が周りかね？ [#y0f55429]
>[[間瀬茂]] (2004-09-08 (水) 19:37:45)~
~
最新の r-help 記事に 64 bit の Sun Sparc マシンで R を動...
It uses more memory at all times and so gc() takes longe...
slowdown, for example 90 vs 80 secs for a run of R-devel...
either Sun's Forte 7 or gcc 3.4.1 compilers). But `25-3...
unexpected and needs investigation.
The only difference between 32-bit and 64-bit versions o...
be the size of the pointers and (probably) less efficien...
no reason to expect a performance boost with 64-bit appl...
to do more work and are only worthwhile if you need the ...
memory or also on disc as 64-bit applications use large ...
//
-このコメントの通り、単純な６４ビット化は性能向上にはつな...
-DNAマイクロアレイの解析(bioconductor)では64bitの威力を実...
#comment
//不埒な書き込みがあったので消去しておきます
//タイトルが「**お小遣い稼ぎ！」だから，内容は推して知る...
**データセットの引用 [#hfc3cf59]
>[[間瀬茂]] (2004-08-26 (木) 17:38:14)~
~
「工学のためのデータサイエンス入門」執筆の際、執筆者の間...
~
Dear Arne,
I'm not a lawyer, but I believe that data per se are not
copyrightable.
Regards,
John
~
>> -----Original Message-----
>> From: r-help-bounces@stat.math.ethz.ch
>> [mailto:r-help-bounces@stat.math.ethz.ch] On Behalf O...
>> Sent: Wednesday, August 25, 2004 12:06 PM
>> To: r-help@stat.math.ethz.ch
>> Cc: jfox@mcmaster.ca
>> Subject: [R] License for including datasets in packages
>>
>> Dear All,
>>
>> I would like to publish a function for 'heckit' estim...
>> together with two examples from Greene's and Wooldrid...
>> econometric textbooks.
>> These examples use the dataset of Mroz (1987) that is...
>> available in John Fox' "car" package. However, not all
>> variables that are used in my examples are available ...
>> "car" package. Therefore, I want to put the full Mroz...
>> in my package. This full dataset can be downloaded fr...
>> internet sites (e.g.
>> pages.stern.nyu.edu/~wgreene/Text/econometricanalysis...
>> www.stata.com/texts/eacsap/), but I did not find any ...
>> information.
>>
>> Whom do I have to ask wether I am allowed to put thes...
>> in a GPLed R package? (Mr. Mroz, Mr. Greene, STATA,
>> University of Michigan (PSID), . . . )
>>
>> @John: Do you have any license information on these d...
>>
>> Thanks,
>> Arne
//
-GPLはプログラムについてのみが適用範囲かと。データについ...
-CreativeCommonsってデータには適用できないんでしょうか? -...
-CreativeCommons とは初耳ですが何のことでしょうか？ -- [[...
-http://www.creativecommons.jp/ にも説明があるようです（...
-米国のレッシグ教授が始めたプロジェクトです。自由利用の範...
#comment
**女子マラソンの金メダリスト [#l1406a16]
> (2004-08-26 (木) 00:52:23)~
~
アテネオリンピック女子マラソンの金メダリストの所属会社が...
//
-[[いろいろあるのかも:http://www.isdnet.co.jp/~Ombudsman/...
#comment
**IBM developerWorks の記事 [#d000819c]
>[[nanashi]] (2004-08-14 (土) 11:47:22)~
~
ちょっと古いですが、Rの紹介記事です。興味があったら読んで...
~
http://www-6.ibm.com/jp/developerworks/linux/040116/j_l-s...
//
-[[コメントアーカイブ]]（2004-01-29 (木) 20:54:36）にて紹...
#comment
**sas2R [#y57bb720]
> (2004-08-13 (金) 12:13:34)~
~
(いささか古い r-help 記事によると) SAS から R プログラム...
//
-http://www.biostatistics.dk/sas2r/ -- [[学生]] &new{2004...
-四月馬鹿ですか、だまされた。でも本気にした、本気にしたい...
-あったらいいな，とおもうpackageですよね．一生懸命調べち...
-「一生懸命調べちゃいました」<- ごめんなさい。 -- [[元記...
#comment
**車輪で特許をとった男の話（オフトピック） [#w33dfbac]
>[[間瀬茂]] (2004-08-09 (月) 14:56:10)~
~
暑いですね、を口実に仕事もせずに無駄知識を紹介します。「...
~
(1) 相続税の低下の見込みが延命効果を持つことの研究~
(2) 落下するバタートーストの力学的研究~
(3) 米国の各州毎の地獄へ行く人数の推計~
(4) ミハイル・ゴルバチョフが反キリストでない確率の算出~
(5) ビスケットを紅茶に浸す理想的な方法の考察
L^2 = (γ x D x t)/(4 x η) ここで
L は紅茶にビスケットを浸す幅、γ は紅茶の表面力、D はビス...
t は浸す時間、η は紅茶の粘性率
また英国規格協会（BSI）には紅茶の正式ないれ方に関する６頁...
//
#comment
**現在の CRAN 登録パッケージ数 [#ad4bfc6e]
> (2004-08-05 (木) 01:40:16)~
~
数えてみたら 2004.08.05 現在 CRAN 登録パッケージ数は 381 ...
//
#comment
**BLAS for Pentium4 512kB(Windows, Linux) [#vc93e740]
>[[後藤]] (2004-08-03 (火) 22:41:42)~
~
現在改良版の BLAS のリリースを準備しているのですが、とり...
~
http://www.cs.utexas.edu/users/kgoto/libraries/libgoto_p4...
http://www.cs.utexas.edu/users/kgoto/win2k/libgoto_p4_512...
~
にありますので、wget で取ってください。なお、R のインスト...
//
-complexも動作しました。CPUが違うにもかかわらず(^^;,ATLAS...
#comment
**数学者は如何にしてパン屋の不正を見抜いたか [#d1162628]
>[[間瀬茂]] (2004-08-03 (火) 20:55:57)~
~
(最近読んだ本から紹介)~
~
数学者ポアンカレは、毎日買っているパン屋の1kgのパンがしば...
~
R でシミュレーションすると次のようになるでしょうか。y は...
~
x <- rnorm(365, mean=950, sd=30)
y <- sapply( 1:365, function(i) max(rnorm(500, mean=950,...
par(mfrow=c(1,2))
hist(x); hist(y)
#ref(poincare2.png)
//
-あ，この話，ガウスではなかったのですか。すっかり，誤解し...
-今手元に本が無いので出典をかけませんが、ポアンカレと書い...
-[[この話ですね:http://www.anecdotage.com/index.php?aid=3...
-ご指摘ありがとうございます。「この話」の解説は、正規分布...
-出典は「確率統計で世界を読む」B.K. ホランド著、林大訳、...
#comment
**お手軽 WinR で ATLASクッキング [#sba4b38f]
>[[なかま]] (2004-07-28 (水) 10:02:42)~
~
セレロンの人でも早くなると思います.１世代前のATLASがいい...
あまり知らない指令が来てもCPUは遅くなるでしょうし.
http://toucan.stats.ox.ac.uk/R/RWin/ATLAS/
から適当なRblas.dllをダウンロードしてR_HOME/bin にコピー...
//
-Windowsでもできたんですね。早速入れました。最初にATLAS入...
-下記のtest2を実行したところ [1] 1.21 0.01 1.41 NA NA...
-辛口な教授の労力には脱帽の一言です.今夜から頭を向けて寝...
-標準状態に戻して計り直してみました。結果は[1] 6.48 0.06 ...
-ちなみにマシンはAthlonXP2500+(Barton)@2700+(2GHz)+nForce...
-ATLASって何ですか？ぐぐってみたけどよく分かりませんでし...
- Debian パッケージの解説より引用。automatically の意味が...
Automatically Tuned Linear Algebra Software,generic stat...
an approach for the automatic generation and optimizatio...
software. Currently ATLAS supplies optimized versions fo...
set of linear algebra kernels known as the Basic Linear ...
Subroutines (BLAS), and a subset of the linear algebra r...
LAPACK library. . This package includes the static libra...
symbolic links needed for program development. .
Homepage: http://math-atlas.sourceforge.net/
-ありがとうございます。Rを加速するためのライブラリで、Rの...
-古代は人々がCPUのキャッシュに合わせて人手でBLASの最適化...
-ATLASはSMP対応版も作成できるようですのでそれとマルチプロ...
-LAPACKまでとなると、考えただけでも面倒そうな気がしますが...
-すかさず宣伝。Debian では ATLAS 化された BLAS, LAPACK が...
-FreeBSDのportsでもコンパイルは自動化されてますよ。math/a...
-というわけで最適化済BLASを参照するようにLAPACKを作成して...
-LAPACKの高速化はRにはあまり効果が無いと言うのは本当だっ...
-まだ奥の手があるんですか！ -- [[間瀬]] &new{2004-07-29 (...
-http://www.aspect-sdm.org/Parallel-R/ で背後はLAM/MPI(rs...
-中間さんの test2 関数で再度調べてみました。libgoto なし...
libgoto 無し mean = 1.5024 sd = 0.1505
fivenum(xnogoto)
[1] 1.200 1.400 1.480 1.645 2.100
libgoto あり mean = 1.3056 sd = 0.2866
fivenum(xgoto)
[1] 0.000 1.230 1.325 1.475 1.670
以下のボックスプロットの左は libgoto 無し、右はありです。
libgoto ありで時間 0 というはずれ値は素性が良くわかりま...
#ref(goto.png, left)
-並列処理で一番簡単なのはsnowをSOCKでrshパスワード無しに...
-マルチプロセッサやクラスタリングに逝ってしまうとマシンの...
-もしかしたらPentium4最適化(SSE2を使う)方が効果大きいかも...
-Intel Compiler v8 だと *FLAGS="-mp -pc64" を指定。更にif...
#comment
**お手軽Rで後藤クッキング [#a1d83731]
>[[なかま]] (2004-07-27 (火) 22:33:31)~
~
http://www.cs.utexas.edu/users/kgoto/ から自分のCPUにあった
libgoto と xerbla.c をダウンロード.
gzip -d libgoto_p4_128-r0.94.so.gz
gcc -c xerbla.c
gcc -shared -o libgoto.so ~/libgoto_p4_128-r0.94.so xerb...
Rgoto と言うシェルスクリプト作成
#!/bin/bash
export LD_PRELOAD=~/libgoto.so:/usr/lib/libm.so
${RPROG:-R} -q $*
起動
$ ./Rgoto
で, R を実行すれば標準BLASやATLASの変わりに後藤(日本語訳...
元もとATLAS使っているなら劇的には早くはならないと思います...
//
-中間さんの処方箋通りにしてみました。1000 by 1000 の行列 ...
-FreeBSD派はいかにすれば・・・。というかATLASより速いんで...
-どうも goto ライブラリは CPU 毎に最適化チューンナップを...
-CPU最適化するのはATLASも同じでは。 -- [[aki]] &new{2004-...
-そうなんですか。そういえば Debian の atlas3 パッケージは...
-ATLASを最大性能(キャッシュのSIZEやCPU毎のオプテイマイズ)...
-FreeBSDのportsでは、ATLASもコンパイルからやれます。ただ...
-/proc/cpuinfo からL2ひろって atlasのsrpmを作れば嬉しい? ...
-cpuid から動的に最適なルーチンを呼び出すということは可能...
-もしかして噂の後藤さんご本人でしょうか。わざわざご投稿光...
-です。ライブラリに対する希望がありましたらお知らせくださ...
-よもや御本人から投稿していただけるとは考えず(豪雨汗)、失...
-なんかとんでもないトンチンカンなことを書いていたのでは？...
-今のところ、再配布は遠慮していただいています。ライセンス...
-LD_PRELOAD はすべてのプロセスの起動に影響を与えるため(プ...
-そうですね。Rgoto等のシェルスクリプトにくるむように書い...
-ぐわ、ちょっと見てない間に・・・。えー、FreeBSD用がある...
-ソースといっても、中身はほとんどアセンブラなので、FreeBS...
-げ、そうなんですか。それはまいった。うーん、Linuxに戻る...
-ほとんど違いはないので、あとは趣味の問題かと。FreeBSD 用...
#comment
**FreeBSDのportsから日本語化済Rをインストールする [#k371b...
>[[aki]] (2004-07-27 (火) 14:17:01)~
~
http://www.fifthdimension.jp/index.php?R%2FTips~
に半自動化スクリプトも含めて書いておきました。~
余計なこともやってますがとりあえずパッチの適用までは確認...
動作確認が取れましたらこちらのTipsにでもミラーしていただ...
//
-devPS.cのコンパイルでコケますね・・・(滝汗)。 -- [[aki]]...
-スクリプトの方がタコでパッチがちゃんと当たってませんでし...
-MAIN_CFLAGSがdevPS.cのコンパイル時に反映されてない? よう...
-R-1.9.1.l10n.psxfigtex.20040625.patchも微妙に改変しない...
-あ、すみません。devPS.cは直してあるんですが、公開し忘れ...
-では、うちのやつは消してしまっても問題ないのでしょうか。...
-今にして思えば日本語か掲示板の方が適当だったか・・・。 -...
-お手数おかけします。パッチ置いときました。ATLAS良いです...
-なかまさんのここのところの宣伝で気になってきていますが、...
-ldd usr/lib/R/bin/R.bin or usr/lib/R/bin/libR.so(共有で...
-そのままでは libgoto.so は見つかりませんでしたから、Debi...
-gotoはソースが無いからなんでは? -- [[なかま]] &new{2004-...
-FreeBSDマシンでも下に出ているtest2を走らせてみました。結...
-遅すぎるなーと思ってLinuxエミュを切ったら一つ目が1.15に...
-この結果はWinXPマシン同様2GHzにクロックアップした結果で...
#comment
**ELT [#wa7dca9b]
> (2004-07-27 (火) 12:37:19)~
~
率直に言って、ＲはＥＬＴ（※）には向いていないと思うのです...
?*nixな人だとAWK,PERL等になるのでしょうか？~
さすがに、SASをELTに使ってRで分析する人はいないですよね。~
~
※ELT~
Extraction（抽出）、Load（ロード）、Transformation（変換）~
//
-ELTってなんですか？ -- &new{2004-07-27 (火) 19:32:03};
-Every Little Thing でしょう(^_^;) -- &new{2004-07-27 (火...
#comment
**予約キーワード [#ube48076]
>[[青木繁伸]] (2004-07-27 (火) 01:23:05)~
~
名前付きの引数の指定の時，名前は前方から数文字で指定できる~
しかし，以下のようなこともある~
> x <- c(3.1, 2.4, 1.3, 4.2, 5.3)
> hist(x, breaks=1:6) # フルネーム（当然ながら，ちゃんと...
> hist(x, b=1:6) # b で始まる引数が2つあるので，識別不能...
Error in hist.default(x, b = 1:6) : argument 2 matches m...
> hist(x, br=1:6) # ちゃんと動く
> hist(x, bre=1:6) # ちゃんと動く
> hist(x, brea=1:6) # ちゃんと動く
> hist(x, break=1:6) # break は予約キーワードなので，エ...
Error: syntax error
ちょっと釈然としない。~
自作関数でも当然ながら，同じ問題は生じる。
> a <- function(fortran=0) print(fortran)
> a(f=3)
[1] 3
> a(fo=4)
[1] 4
> a(for=5)
Error: syntax error
> a(fort=6)
[1] 6
//
-全く解決策ではないですが、a(`for`=5) とバックチックでく...
#comment
**Eclipse Plug-in [#oedfdd0d]
> (2004-07-26 (月) 20:06:05)~
~
とうとう出ました。r-help に情報がありました。~
~
まだ試していません。~
~
Eclipse 3.0用で以下の機能があるようです。~
~
シンタックスのハイライト、括弧のマッチングチェック、IDE...
~
以下より入手可能です。~
~
~
http://www.walware.de/rplugin.zip~
//
-どうもXPではまったく動かない。どなたか実行に成功された方...
-[[ここ:http://www.walware.de/go?to=statet]]にプラグイン...
#comment
**TV や映画に出ていたR言語 [#f4cd2e51]
> (2004-07-24 (土) 01:51:01)~
~
てなことないでしょうね、~
//
-昔、S-PLUSの無料紹介ビデオ（英語・字幕なし）があったなぁ...
// -何でもとはいいながら、あまりに情報量のない話題はいか...
#comment
**浮動小数演算の怪 [#ud3aaae9]
>[[QDU]] (2004-07-21 (水) 11:01:10)~
~
これまた r-help 記事ネタ。次の一見無邪気な引き算の結果は...
> 0.6-0.3-0.3
[1] 0
> 0.6-0.4-0.2
[1] -5.551115e-17
> 0.6-0.5-0.1
[1] -2.775558e-17
> 0.6-0.2-0.4
[1] -5.551115e-17
> (0.6-0.2)-0.4
[1] -5.551115e-17
> 0.6-(0.2+0.4)
[1] -1.110223e-16
//
-2進表現で循環小数にならない10進数以外は計算機中の実数は...
Excel は，=3.6-1.2-1.4 でも 1 を表示しますので，単に「0に...
-R では次のようになります。この辺の機微もよう分からん。--...
> 3.6-1.2-1.4
[1] 1
> 1-(3.6-1.2-1.4)
[1] -4.440892e-16
-3.6-1.2-1.4 を 1-4.440892e-16 の値として保持するには 1 ...
#comment
**バージョン４のS言語はコンパイル機能を持つ！ [#m616a9c3]
> (2004-07-18 (日) 23:40:54)~
~
たまたま[[ベル研:http://cm.bell-labs.com/cm/ms/department...
//
#comment
**THE R BOOK の続編【R入門＋セミコマンドリファレンス】 [#...
>[[舟尾]] (2004-07-17 (土) 00:00:53)~
~
来年あたりに出版して頂けるであろう THE R BOOK の続編を私...
書籍として刊行するということで，ある程度は差別化を図る予...
[[http://cse.naro.affrc.go.jp/takezawa/r-tips/r.html]]~
です．そこで，もし良ければ皆様に R-Tips のダメ出しをお願...
//
-出版社が了解するかどうか怪しいですが、できた分から公開し...
-ご意見ありがとうございます．それはさすがに出版社の方に怒...
-日本語パッチで拡張された機能（日本語の変数等）も含まれる...
-日本語パッチで拡張されたのは日本語の為の設定ファイル(X11...
-貴重なご意見ありがとうございます．par 関係のオプションは...
-R-Tipsをベースにしたリファレンス本ということで、楽しみに...
-RjpWikiはご自由にお使いください。~
ありがとうございます！リファレンスを書く際には RjpWiki は...
-9月いっぱいでご意見の受付を終了いたしました．青木さん，...
#comment
**source() 関数のネスト [#afea9436]
>[[QDU]] (2004-07-12 (月) 22:08:24)~
~
r-help におもしろい話題がありました。普通 source() 関数を...
# ファイル "a.R"
a <- function(x) x^2
source("b.R")
~
# ファイル "b.R"
cat(a(10),"?n") # 単に a(10) では何故か何も表示されない？
~
# 実行結果
> source("a.R") # 中で source("b.R") を実行し、a(10) を...
100
#comment
**goto [#k4833353]
>[[なかま]] (2004-07-08 (木) 19:26:32)~
~
「ゴーツー」では無く「ごとう」なんですね。(^_^;
test2 <- function ( n=500 )
{
A<-array(rnorm(n^2), dim=c(n,n))
B<-array(rnorm(n^2), dim=c(n,n))
C<-array(rnorm(n^2), dim=c(n,n))
D<-array(rnorm(n^2), dim=c(n,n))
BA <- B%*%A
return (system.time(A%*%solve(t(BA)%*%BA+C)%*%BA%*%D ))
}
貧乏な私のマシン(ここでRのビルドとかしてるんですが)
> test2()
[1] 23.69 0.07 23.99 0.00 0.00
貧乏なRweb-jpマシン
> test2()
[1] 10.91 0.11 11.08 0.00 0.00
某所R-1.8.1
> test2()
[1] 9.822266 0.015625 9.839844 0.000000 0.000000
某所R-1.9.1 --with-blas=-lgoto
> test2()
[1] 1.01757812 0.02050781 1.04199219 0.00000000 0.00000000
ディスク欲しさにCPUをCeleronに落としてはいけない訳ですね...
//
-嗚呼、南北問題！ R.Ihaka 氏がどこかに「私はワークステー...
-[[libgoto:http://www.cs.utexas.edu/users/flame/goto/]] ...
-むーん、でもよくよく観察してみると警告が(環境?)。。。ATL...
-当家のDNS,IPSEC-SGはp5-133ですが、test2()は130程でした。...
-Rweb-jp の R をBLASにATLAS利用に変えてみました。[[Rweb-j...
#comment
**RjpWiki トップページについて [#w9418581]
>[[間瀬茂]] (2004-07-08 (木) 17:42:18)~
~
ここ数日私等がトップページをいじくっております。どうして...
//
#comment
**R で データマイニング [#d6c8e120]
> (2004-07-06 (火) 18:21:43)~
~
r-help に以下のサイトの情報が流れていました。~
~
http://www.liacc.up.pt/~ltorgo/DataMiningWithR/~
//
#comment
**漸近回帰モデル [#bb4ccc28]
>[[間瀬茂]] (2004-07-05 (月) 14:42:49)~
~
最近、R の stats パッケージに漸近回帰モデル (asymptotic r...
//
#comment
**誕生日のパラドックス [#e6150096]
> (2004-07-02 (金) 23:43:29)~
~
pbirthday は正しい値（？）を返さない~
これぞ，まさしく，パラドックス(^_^;)~
//
-実例は？ ？birthday によれば近似計算を使っているとの注意...
-以下を実行のこと。~
pbirthday2 <- function(n)
{
1-exp(lfactorial(365)-n*log(365)-lfactorial(365-n))
}
n <- 2:50
result <- cbind(n, sapply(n, pbirthday), sapply(n, pbirt...
colnames(result) <- c("n", "pbirthday", "pbirthday2")
result
実行結果（n が大きくても，小さくても，中くらいでも，不正...
n pbirthday pbirthday2
[1,] 2 0.005479452 0.002739726
[2,] 3 0.016438356 0.008204166
[3,] 4 0.026301370 0.016355912
[4,] 5 0.034246575 0.027135574
:
[21,] 22 0.499315068 0.475695308
[22,] 23 0.530136986 0.507297234
[23,] 24 0.531506849 0.538344258
[24,] 25 0.583542799 0.568699704
:
[47,] 48 0.957357026 0.960597973
[48,] 49 0.961562289 0.965779609
[49,] 50 0.967534247 0.970373580
関数本体はもちろん，1-prod((365:(365-n+1))/rep(365,n)) ...
//
-?birthday で説明されているように、ここで使われている近似...
-所詮，最終的な答え（確率）は，0〜1の範囲の値しか取らない...
-組み合わせ問題では(少なくとも途中では)すぐに超天文学的な...
-具体例があれば，きぼんぬ。 -- &new{2004-07-06 (火) 22:5...
#comment
**64 bit 版 R ！ [#e0d79fe1]
> (2004-07-02 (金) 23:01:19)~
~
r-help 記事より。お金持ちには６４ビット版 R で ４ Gigabyt...
Did you compile R as 64-bit executable on the Irix?
If not, R will be subjected to the 4GB limit of 32-bit s...
//
-RBookのbioconductorのところでも言及しましたが、同じ4GBで...
#comment
**anti-R vitriol [#c4c7cf18]
>[[間瀬茂]] (2004-07-01 (木) 17:10:54)~
~
昨日の r-help に表題のようなタイトル(vitriol というのは硫...
A colleague is receiving some data from another person. ...
SAS and it takes 30s and uses 64k RAM. That person then ...
and it takes 10 minutes and uses a gigabyte of RAM. Per...
It's not that I think SAS is such great software,
it's not. But I really hate badly designed
software. R is designed by committee. Worse,
it's designed by a committee of statisticians.
They tend to confuse numerical analysis with
computer science and don't have any idea about
software development at all. The result is R.
I do hope [your colleague] won't have to waste time do...
[this analysis] in an outdated and poorly designed piece
of software like R.
Would any of the "committee" like to respond to this? Or...
collective forehead and wonder how someone could get su...
>system.time(spam <- read.table("data_dmc2003_train.tx...
header=T,
colClasses=c(rep("numeric", 83...
"character...
[1] 15.92 0.09 16.80 NA NA
> system.time(spam <- read.table("data_dmc2003_train.tx...
[1] 187.29 0.60 200.19 NA NA
//
#comment
**なんで runif はランダムなの？ [#ka5cdc11]
> (2004-07-01 (木) 13:20:06)~
~
意味不明なタイトルと思われるでしょうが、つぎの結果をどう...
(つまり百万個一様乱数を発生する手間は同じはず、念には念を...
test <- function (i) {set.seed(222); system.time(runif(1...
x <- sapply(1:1000, test)
fivenum(xx)
[1] 0.24 0.28 0.30 0.30 0.47
//
-ほかに何のプログラムも走っていなければ同じ時間になるので...
Rweb:> system.time(for (i in 1:10000000) x <- 1)
[1] 7.01 0.10 7.13 0.00 0.00
Rweb:> system.time(for (i in 1:10000000) x <- 1)
[1] 5.85 0.07 5.92 0.00 0.00
これは，Ｒの問題ではないでしょう。UNIX の time コマンドも...
-なるほど、すると system.time でコード速度をはかる時は、...
#comment
終了行:
COLOR(blue){SIZE(20){RjpWiki自体へのコメントでも、雑談で...
COLOR(red){注意}：新規記事用の入力欄は以下の目次の直後に...
----
COLOR(red){過去の記事のアーカイブ} 。アーカイブ中の記事へ...
-[[なんでも掲示板アーカイブ（１）]] (元記事が 2003-07-17 ...
-[[なんでも掲示板アーカイブ（２）]] (元記事が 2003-11-27 ...
-[[なんでも掲示板アーカイブ（３）]] (元記事が 2004-01-01 ...
----
#contents
----
回答を期待する記事は原則 COLOR(red){質問コーナー [[Ｑ＆Ａ...
~
#article
**経済学における統計的有意の記号 [#odc45137]
>[[y.a.]] (2004-11-23 (火) 01:01:55)~
~
[[Ｑ＆Ａ (初級者コース)]]にて，経済学の記号のつけ方の慣習...
結論から言うと，人やジャーナルによって異なるものの，特に...
~
ただ，~
cutpoints= c(0, 0.01, 0.05, 0.1), symbols = c("***", "**...
記法は，少なくとも経済学の分野には少なからず存在はすると...
~
以下，経済学のトップジャーナル群の論文，掲載予定論文から...
~
1) t値を示して星はつけない（でもp値はつけない）派
- [[American Economic Reviewの例(p.1396, Tab.6):http://fa...
- [[Quarterly Journal of Economicsの例(p.1298, Tab.IV):ht...
~
2) (1,5,10) -> (***,**,*) 記法前近代派
- [[Review of Economic Studiesの例(p.22, Tab.6):http://ww...
- その他，[[Journal of Development Economics誌:http://www...
//
-統計学の大口ユーザーである心理学等でも類似の慣習があるの...
-うちの学部の経済関係の先生も，ｔ値は書くけどＰ値は書かな...
某誌（医学系）に載ったチュートリアルで
＋ または △： p < 0.1 傾向あり，傾向差
＊： p < 0.05 有意差あり
＊＊： p < 0.01 明らかな差あり
＊＊＊： p < 0.001 著明な差あり
というのがあり，＋とか△はいかがなものか？とおもったら，心...
へ〜。郷に入らば郷に従えか -- &new{2004-11-25 (木) 18:58...
-Rを離れて雑談になってしまうのですが、「統計的検定（は使...
#comment
**R からキーワード検索する関数 (from r-help) [#f64c2f99]
>[[QDU]] (2004-11-20 (土) 20:51:33)~
~
今日の r-help に面白い関数が紹介されていましたので紹介し...
help.search.archive<-function(string){
RURL="http://www.google.com/u/newcastlemaths"
RSearchURL=paste(RURL,"?q=",string,sep='')
browseURL(RSearchURL)
return(invisible(0))
}
たとえば help.search.archive("plot3d") を試してみて下さい。
//
-やはり r-help 投稿の同様の関数： -- [[QDU]] &new{2004-11...
Inspired by the functions that Barry Rawlingson and Dave...
searching Rwiki and R-help archive, I've made up a funct...
search on Prof. Baron's site (Thanks to Prof. Baron's he...
query string!):
RSiteSearch <- function(string, restrict="Rhelp", format...
sortby="score", matchesPerPage=1...
URL <- "http://finzi.psych.upenn.edu/cgi-bin/htsearch"
qstring <- paste(URL, "?config=htdigrun1", sep="")
## replace spaces with "%20" in the query
string <- paste("words=", gsub(" ", "%20", string), ...
mpp <- paste("matchesperpage=", matchesPerPage, sep=...
format <- charmatch(format, c("long", "short"))
if (format == 0) stop("format must be either long or...
format <- paste("format=builtin-", switch(format, "l...
sortby <- charmatch(sortby, c("score", "time", "titl...
if (sortby == 0) stop("wrong sortby specified")
sortby <- paste("sort=",
switch(sortby, "score", "time", "tit...
sep="")
res <- charmatch(restrict, c("Rhelp", "doc", "functi...
if (res == 0) stop("wrong restriction specified")
res <- switch(res, "Rhelp00/archive|Rhelp01/archive|...
"finzi.psych.upenn.edu/R/doc",
"finzi.psych.upenn.edu/R/library",
"finzi.psych.upenn.edu/R/doc|finzi.psy...
res <- paste("restrict=", res, sep="")
qstring <- paste(qstring, res, format, sortby, string...
browseURL(qstring)
invisible(qstring)
}
The options roughly correspond to the options on that se...
some people find this somewhat helpful. Comments/sugges...
improvement are much appreciated.
Best,
Andy
#comment
**The hidden costs of GPL softwares [#f6ce8d13]
>[[間瀬茂]] (2004-11-18 (木) 22:54:56)~
~
今日の R-help は表題の投稿への反応で盛り上がりました。い...
//
- 他のソフトと比較するのは良くないことなのかもしれません...
ヘルプが充実していて使いやすいといえば Mathematica ですが...
ヘルプに関して云えば，R のヘルプは関数の構文と引数の説明...
これは私の勝手な意見ですが，「学習コストが高い」のは「か...
-ある集会で「SASの日本語化をしましょう」という趣旨の講演...
-SASのGUIメニューや基本的な計算のアウトプットは日本語化さ...
-船尾さんがおっしゃるEXCELのクッキング本とつながりますが...
参考までに，上田・ 苅田・本田(2003) EXCEL徹底活用多変量解...
-linux,*BSD,glibc,readline,gzip,bzip,XOrg,gcc,libgoto...R...
#comment
**S-Plus や R の関数のカテゴリー化された検索システム [#bf...
>[[間瀬茂]] (2004-11-18 (木) 17:42:15)~
~
r-help に S-Plus や R の膨大な関数を用途別にカテゴリ化し...
~
[[S-Plus/R Function Finder:http://biostat.mc.vanderbilt.e...
//
#comment
**R 2.0.1 for VineLinux 3.0 [#lc1a3bf1]
>[[谷村]] (2004-11-16 (火) 20:39:12)~
~
本家(CRAN)にR-2.0.1-0vl1.i386.rpmをputしました。これはオ...
//
#comment
**どこかでみたような [#fc6abf97]
>[[なかま]] (2004-11-14 (日) 00:55:59)~
~
[[中文TeX:http://www.ctex.org/forums/index.php?showtopic=...
ユーザー数ってどれくらいいるんでしょ.~
//
-びっくりしましたね。ただ、先月中国の集会に参加した時印象...
-本当にすみませんでした。この文章(Rを薦める100理由)は、私...
-消す必要はないと思いますが(^^;。中国語になってる！のが面...
-どうせこちらも「無断」引用しまくっていますから、お構い無...
-日本語のサイトを中国の人が読めるということは、逆もまた真...
-批評してる(?)のは日本人のような気がしますが...(^_^; -- [...
-よしよしさんの「良くない」っていう発言は、不偏推定値では...
#comment
**Ｒで物理乱数 [#v53f4e78]
>[[垂水共之]] (2004-11-10 (水) 17:21:09)~
~
[[FDK:http://randomstreamer.ntl.co.jp]]という会社がUSB接...
Windows版のRで使えるようにDLLも作ってもらいました。~
Rユーザーの方はサポートCD-ROM特別サンプルも手に入れてくだ...
//
-おお、すばらしいですね! -- &new{2004-11-12 (金) 08:16:3...
-はたして、お小遣いで買える物なのでしょうか...(^^; -- [[...
-小遣いではちと高いですが、3万程度と聞いています -- [[垂...
-set.seed(xxx) は原理的にできないわけですね。 -- &new{20...
-生成速度がHPには公表されていませんがどれくらいなのか一寸...
-ありがとうございます。3万がんばり(なにを?)ます。 -- [[な...
-内蔵してますRPG100の乱数生成速度は250kbits/sのようですよ...
-ちなみに、販売しているのは「日本テクノ・ラボ（株）」見た...
-Random StreamerはSPLUSで使える様です。 -- [[ななし]] &ne...
#comment
**かんたんな入力チェックを入れるのは [#l3f27661]
>[[なかま]] (2004-11-08 (月) 01:21:09)~
~
--- article.inc.php.orig
+++ article.inc.php
@@ -189,7 +189,24 @@
$article_rows = article_ROWS;
$article_cols = article_COLS;
$string = <<<EOD
-<form action="$script" method="post">
+<script>
+ <!--
+ function isNull(elem,jname){
+ if(elem.value==''){
+ elem.focus();
+ alert(jname+'がはいっていると嬉しいかも.');
+ return true;
+ }
+ return false;
+ }
+ function chkArticle(myform){
+ if(isNull(myform.name,'お名前'))return false;
+ if(isNull(myform.subject,'題名'))return false;
+ return true;
+ }
+ -->
+</script>
+<form action="$script" method="post" onSubmit="return c...
<div>
<input type="hidden" name="article_no" value="$articl...
<input type="hidden" name="plugin" value="article" />
~
まあ,気休めですが効果はあるかも.~
//
-ありがとうございます。早速導入してみました:-) -- [[岡田]...
-無題が減るといいですね.:-) -- [[なかま]] &new{2004-11-13...
#comment
**R とアメリカ大統領選挙 [#mac5766b]
> (2004-11-03 (水) 20:00:03)~
~
今回も米大統領選挙は接戦ですね。~
~
今回はオハイオ州が注目の的ですが、R の car パッケージに...
//
#comment
**日本人によるアドオンパッケージ Zelig [#zd08ed8f]
> (2004-10-23 (土) 14:48:44)~
~
今日初めて気がつきましたが MNP というアドオンパッケージの...
//
#comment
**中国のR事情 [#p98b8c1a]
>[[間瀬茂]] (2004-10-19 (火) 23:28:09)~
~
集会で中国にいって来ました。桂林にある広西師範大学が会場...
//
-r-helpにあったとおりGB18030で表記できれば広まると思いま...
-データもグラフ出力も、そのまま中国語(gb2312)は大丈夫です...
-中国：Linuxは普及せず、ウイルスソフトで国産好調＝サーチ...
-うちのサイトに来る怪しげパケット,scanは中国が一番多いの...
-広西師範大学 って GNU なんですね。(^_^; -- [[なかま]] &n...
-そうなんです。Guangxi Normal University。広西チアン自治...
-PS. 学内をおっかなびっくり見物していたら、ひとりだけ茶髪...
-広西なんて中国の中でもかなり田舎なんだ！中国の地方差はす...
-PS.学食の値段も地方によって違うんだよ。上海の学食ならこ...
#comment
**Mac 版 R 2.0.0 [#a50884d8]
>[[青木繁伸]] (2004-10-18 (月) 11:27:03)~
~
まだ筑波大学のmirrorには現れていない~
working directry の設定で，メニューからやると設定できない...
注釈に日本語を使っている関数のリストを表示する（関数名だ...
一部パッケージがまだ対応できていない（MASS などが使えない...
環境設定メニューがない（なくてもいいようにはできていない...
~
というような具合で，R2.0.1の出現を待ちたい。~
//
-MASS などは，使えました。なぜ昨日使えなかったのか不明だ...
#comment
**R 2.0.0 のリリース近し？ => ついに公開されました（祝^祝...
>[[間瀬茂]] (2004-10-02 (土) 00:17:53)~
~
今日の r-help に R の 2.0.0 に関する問い合わせがあり、R ...
Or just wait for the release on Monday October 4th.
//
-素晴らしい！ついに R もバージョン 2 に突入ですねぇ．Wind...
-10/4(日本時間10/5の何時頃?)にリリース予定です。betaは既...
-*nix方面はなんとか構築。うーん、窓用たいへんかも...(T_T)...
-何かメジャーバージョンアップにふさわしい目玉があるのでし...
-表現が適切かはわかりませんが、地味に細かく*沢山*拡張した...
-パッケージ化待たなくともconfigure,make,make installでも...
-[[2.0 SERIES NEWS :https://svn.r-project.org/R/trunk/NEW...
-あのー、日本語バージョンはいつ頃に。(だから金沢に向かっ...
-どうせすぐ 2.0.1 もしくは 2.1 に代わりますから慌てる必要...
-出ました。http://cran.r-project.org/src/base/R-2 -- [[ak...
-リリースニュースによれば、この版は R に対する革新的な変...
As was the case with R 1.0.0, this new version represent...
more than a radical change to R. We do plan to celebrate...
with press releases and such.
-r-help でコアメンバーの Dalgaard 氏が次のようにコメント...
No, I don't think we have arrived (would we want to?). W...
have unresolved issues in byte compiling and event loops...
formal methods techniques are only just there now. There...
Bayesianfront, too - important even for us hardcore freq...
-大変遅くなりましたが、R-2.0.0-0vl1.i386.rpm for VineLinu...
-いつもデバッグ用に(ふと、やっちまった?と思う場合があるの...
#comment
**pukiwiki-mode [#xa4a4f21]
>[[谷村]] (2004-09-29 (水) 21:57:08)~
~
今日、[[pukiwiki-mode:http://www.bookshelf.jp/pukiwiki/pu...
~
ブラウザから書き込むのと比べて効率がよさそうなのですが、...
//
#comment
**Rwebに入力された驚くべき物 [#f5935bf2]
>[[なかま]] (2004-09-27 (月) 23:40:17)~
~
たまーにビックリするような物が実行されたりしますんですが...
~
[18時 20分 40秒] [緑茶] 1 個を手に入れました。
[18時 20分 41秒] [緑茶] 1 個を手に入れました。
[18時 20分 42秒] [銀杏] 1 個を手に入れました。
[18時 20分 45秒] [緑茶] 1 個を手に入れました。
[18時 20分 52秒] 経験値が 64 上がりました。
[18時 21分 13秒] 経験値が 104 上がりました。
[18時 21分 24秒] コンボで５度目の攻撃に成功すると、コン...
[18時 21分 25秒] 経験値が 64 上がりました。
[18時 21分 27秒] 経験値が 104 上がりました。
~
これをいったいどうしようと思ったのか非常に気になります(笑...
//
-ひょっとして Rweb という名前のオンラインゲームがあるので...
#comment
**sink 関数のネスト [#i28f3260]
>[[QDU]] (2004-09-27 (月) 21:15:06)~
~
たまたま 2ch の R 関連のやりとりを見ていて次のことに気づ...
~
'sink()' or 'sink(file=NULL)' ends the last diversio...
specified type). As from R version 1.3.0 there is a...
diversions for normal output, so output reverts to t...
diversion (if there was one). The stack is of up to...
connections (20 diversions).
~
# 参考例
> example(lm)
> sink("test.summary")
> sink("test2.summary") # sink 関数のネスト
> summary(resid(lm.D9) - resid(lm.D90)) # test2.summary ...
> sink() #
> summary(resid(lm.D9) - resid(lm.D90)) # test.summary ...
> sink()
> summary(resid(lm.D9) - resid(lm.D90)) # 標準出力に戻った
Min. 1st Qu. Median Mean 3rd Qu. ...
-7.772e-16 -6.245e-17 2.776e-17 1.405e-17 1.145e-16 ...
//
-出力を条件ごとに切替えて使うのによいかと。ネストと言うよ...
#comment
**Excel からコピーしたセル範囲を R に読み込む [#x99a1f1f]
>[[青木繁伸]] (2004-09-19 (日) 14:56:45)~
~
エクセルのセル範囲を R に読み込むために，今まではわざわざ...
from.excel <- function(nc)
{
matrix(scan(""), byrow=TRUE, nc=nc)
}
エクセルで必要な部分をコピーします。
x <- from.excel(2)
などとして，プロンプトが出たらペーストします。~
入力の終わりはMacintoshだとcommand+.~
汎用関数にする価値もないのでこのまま。
//
- ?scan によれば Linux だと CTRL+D, MS Windows では CTRL+...
-これはいいですね。 -- [[林 啓一]] &new{2004-09-27 (月) ...
- 参考になりました -- [[田中]] &new{2007-08-05 (日) 09:55...
#comment
//
**CRAN のすべてのパッケージを含むドキュメントからの全文検...
>[[間瀬茂]] (2004-09-12 (日) 22:35:51)~
~
r-help に「これこれするには何を使えばよいでしょう」タイプ...
[[R site search:http://finzi.psych.upenn.edu/search.html]...
//
#comment
**Rの宣伝講演をします [#cad3e233]
>[[谷村]] (2004-09-08 (水) 23:57:14)~
~
2004年9月10日東京大学本郷キャンパス山上会館で「空間疫学分...
http://www.csis.u-tokyo.ac.jp/japanese/research_activitie...
//
-素晴らしい！どのような反響があったのか，是非お聞かせ下さ...
-慶應三田で開かれたSVG Open と重なって行けなかったのが残...
-遅くなりました。「Rは名前は知っているけど空間情報科学分...
-空間分析ではかつては Splus と組んでいたESRI社の Spatial...
-ところで、リンク集にある谷村先生のご著書『保健医療のため...
-発行日は10月5日になっていますが、実質は9月27日頃の予定。...
-谷村先生、ご回答ありがとうございます。発売を期待しており...
-谷村さんに講演を依頼した者です。Rを知らない人にも「何が...
-今ごろ言っても遅いのですが、関西オープンソース2004で15分...
#comment
**大男総身に知恵が周りかね？ [#y0f55429]
>[[間瀬茂]] (2004-09-08 (水) 19:37:45)~
~
最新の r-help 記事に 64 bit の Sun Sparc マシンで R を動...
It uses more memory at all times and so gc() takes longe...
slowdown, for example 90 vs 80 secs for a run of R-devel...
either Sun's Forte 7 or gcc 3.4.1 compilers). But `25-3...
unexpected and needs investigation.
The only difference between 32-bit and 64-bit versions o...
be the size of the pointers and (probably) less efficien...
no reason to expect a performance boost with 64-bit appl...
to do more work and are only worthwhile if you need the ...
memory or also on disc as 64-bit applications use large ...
//
-このコメントの通り、単純な６４ビット化は性能向上にはつな...
-DNAマイクロアレイの解析(bioconductor)では64bitの威力を実...
#comment
//不埒な書き込みがあったので消去しておきます
//タイトルが「**お小遣い稼ぎ！」だから，内容は推して知る...
**データセットの引用 [#hfc3cf59]
>[[間瀬茂]] (2004-08-26 (木) 17:38:14)~
~
「工学のためのデータサイエンス入門」執筆の際、執筆者の間...
~
Dear Arne,
I'm not a lawyer, but I believe that data per se are not
copyrightable.
Regards,
John
~
>> -----Original Message-----
>> From: r-help-bounces@stat.math.ethz.ch
>> [mailto:r-help-bounces@stat.math.ethz.ch] On Behalf O...
>> Sent: Wednesday, August 25, 2004 12:06 PM
>> To: r-help@stat.math.ethz.ch
>> Cc: jfox@mcmaster.ca
>> Subject: [R] License for including datasets in packages
>>
>> Dear All,
>>
>> I would like to publish a function for 'heckit' estim...
>> together with two examples from Greene's and Wooldrid...
>> econometric textbooks.
>> These examples use the dataset of Mroz (1987) that is...
>> available in John Fox' "car" package. However, not all
>> variables that are used in my examples are available ...
>> "car" package. Therefore, I want to put the full Mroz...
>> in my package. This full dataset can be downloaded fr...
>> internet sites (e.g.
>> pages.stern.nyu.edu/~wgreene/Text/econometricanalysis...
>> www.stata.com/texts/eacsap/), but I did not find any ...
>> information.
>>
>> Whom do I have to ask wether I am allowed to put thes...
>> in a GPLed R package? (Mr. Mroz, Mr. Greene, STATA,
>> University of Michigan (PSID), . . . )
>>
>> @John: Do you have any license information on these d...
>>
>> Thanks,
>> Arne
//
-GPLはプログラムについてのみが適用範囲かと。データについ...
-CreativeCommonsってデータには適用できないんでしょうか? -...
-CreativeCommons とは初耳ですが何のことでしょうか？ -- [[...
-http://www.creativecommons.jp/ にも説明があるようです（...
-米国のレッシグ教授が始めたプロジェクトです。自由利用の範...
#comment
**女子マラソンの金メダリスト [#l1406a16]
> (2004-08-26 (木) 00:52:23)~
~
アテネオリンピック女子マラソンの金メダリストの所属会社が...
//
-[[いろいろあるのかも:http://www.isdnet.co.jp/~Ombudsman/...
#comment
**IBM developerWorks の記事 [#d000819c]
>[[nanashi]] (2004-08-14 (土) 11:47:22)~
~
ちょっと古いですが、Rの紹介記事です。興味があったら読んで...
~
http://www-6.ibm.com/jp/developerworks/linux/040116/j_l-s...
//
-[[コメントアーカイブ]]（2004-01-29 (木) 20:54:36）にて紹...
#comment
**sas2R [#y57bb720]
> (2004-08-13 (金) 12:13:34)~
~
(いささか古い r-help 記事によると) SAS から R プログラム...
//
-http://www.biostatistics.dk/sas2r/ -- [[学生]] &new{2004...
-四月馬鹿ですか、だまされた。でも本気にした、本気にしたい...
-あったらいいな，とおもうpackageですよね．一生懸命調べち...
-「一生懸命調べちゃいました」<- ごめんなさい。 -- [[元記...
#comment
**車輪で特許をとった男の話（オフトピック） [#w33dfbac]
>[[間瀬茂]] (2004-08-09 (月) 14:56:10)~
~
暑いですね、を口実に仕事もせずに無駄知識を紹介します。「...
~
(1) 相続税の低下の見込みが延命効果を持つことの研究~
(2) 落下するバタートーストの力学的研究~
(3) 米国の各州毎の地獄へ行く人数の推計~
(4) ミハイル・ゴルバチョフが反キリストでない確率の算出~
(5) ビスケットを紅茶に浸す理想的な方法の考察
L^2 = (γ x D x t)/(4 x η) ここで
L は紅茶にビスケットを浸す幅、γ は紅茶の表面力、D はビス...
t は浸す時間、η は紅茶の粘性率
また英国規格協会（BSI）には紅茶の正式ないれ方に関する６頁...
//
#comment
**現在の CRAN 登録パッケージ数 [#ad4bfc6e]
> (2004-08-05 (木) 01:40:16)~
~
数えてみたら 2004.08.05 現在 CRAN 登録パッケージ数は 381 ...
//
#comment
**BLAS for Pentium4 512kB(Windows, Linux) [#vc93e740]
>[[後藤]] (2004-08-03 (火) 22:41:42)~
~
現在改良版の BLAS のリリースを準備しているのですが、とり...
~
http://www.cs.utexas.edu/users/kgoto/libraries/libgoto_p4...
http://www.cs.utexas.edu/users/kgoto/win2k/libgoto_p4_512...
~
にありますので、wget で取ってください。なお、R のインスト...
//
-complexも動作しました。CPUが違うにもかかわらず(^^;,ATLAS...
#comment
**数学者は如何にしてパン屋の不正を見抜いたか [#d1162628]
>[[間瀬茂]] (2004-08-03 (火) 20:55:57)~
~
(最近読んだ本から紹介)~
~
数学者ポアンカレは、毎日買っているパン屋の1kgのパンがしば...
~
R でシミュレーションすると次のようになるでしょうか。y は...
~
x <- rnorm(365, mean=950, sd=30)
y <- sapply( 1:365, function(i) max(rnorm(500, mean=950,...
par(mfrow=c(1,2))
hist(x); hist(y)
#ref(poincare2.png)
//
-あ，この話，ガウスではなかったのですか。すっかり，誤解し...
-今手元に本が無いので出典をかけませんが、ポアンカレと書い...
-[[この話ですね:http://www.anecdotage.com/index.php?aid=3...
-ご指摘ありがとうございます。「この話」の解説は、正規分布...
-出典は「確率統計で世界を読む」B.K. ホランド著、林大訳、...
#comment
**お手軽 WinR で ATLASクッキング [#sba4b38f]
>[[なかま]] (2004-07-28 (水) 10:02:42)~
~
セレロンの人でも早くなると思います.１世代前のATLASがいい...
あまり知らない指令が来てもCPUは遅くなるでしょうし.
http://toucan.stats.ox.ac.uk/R/RWin/ATLAS/
から適当なRblas.dllをダウンロードしてR_HOME/bin にコピー...
//
-Windowsでもできたんですね。早速入れました。最初にATLAS入...
-下記のtest2を実行したところ [1] 1.21 0.01 1.41 NA NA...
-辛口な教授の労力には脱帽の一言です.今夜から頭を向けて寝...
-標準状態に戻して計り直してみました。結果は[1] 6.48 0.06 ...
-ちなみにマシンはAthlonXP2500+(Barton)@2700+(2GHz)+nForce...
-ATLASって何ですか？ぐぐってみたけどよく分かりませんでし...
- Debian パッケージの解説より引用。automatically の意味が...
Automatically Tuned Linear Algebra Software,generic stat...
an approach for the automatic generation and optimizatio...
software. Currently ATLAS supplies optimized versions fo...
set of linear algebra kernels known as the Basic Linear ...
Subroutines (BLAS), and a subset of the linear algebra r...
LAPACK library. . This package includes the static libra...
symbolic links needed for program development. .
Homepage: http://math-atlas.sourceforge.net/
-ありがとうございます。Rを加速するためのライブラリで、Rの...
-古代は人々がCPUのキャッシュに合わせて人手でBLASの最適化...
-ATLASはSMP対応版も作成できるようですのでそれとマルチプロ...
-LAPACKまでとなると、考えただけでも面倒そうな気がしますが...
-すかさず宣伝。Debian では ATLAS 化された BLAS, LAPACK が...
-FreeBSDのportsでもコンパイルは自動化されてますよ。math/a...
-というわけで最適化済BLASを参照するようにLAPACKを作成して...
-LAPACKの高速化はRにはあまり効果が無いと言うのは本当だっ...
-まだ奥の手があるんですか！ -- [[間瀬]] &new{2004-07-29 (...
-http://www.aspect-sdm.org/Parallel-R/ で背後はLAM/MPI(rs...
-中間さんの test2 関数で再度調べてみました。libgoto なし...
libgoto 無し mean = 1.5024 sd = 0.1505
fivenum(xnogoto)
[1] 1.200 1.400 1.480 1.645 2.100
libgoto あり mean = 1.3056 sd = 0.2866
fivenum(xgoto)
[1] 0.000 1.230 1.325 1.475 1.670
以下のボックスプロットの左は libgoto 無し、右はありです。
libgoto ありで時間 0 というはずれ値は素性が良くわかりま...
#ref(goto.png, left)
-並列処理で一番簡単なのはsnowをSOCKでrshパスワード無しに...
-マルチプロセッサやクラスタリングに逝ってしまうとマシンの...
-もしかしたらPentium4最適化(SSE2を使う)方が効果大きいかも...
-Intel Compiler v8 だと *FLAGS="-mp -pc64" を指定。更にif...
#comment
**お手軽Rで後藤クッキング [#a1d83731]
>[[なかま]] (2004-07-27 (火) 22:33:31)~
~
http://www.cs.utexas.edu/users/kgoto/ から自分のCPUにあった
libgoto と xerbla.c をダウンロード.
gzip -d libgoto_p4_128-r0.94.so.gz
gcc -c xerbla.c
gcc -shared -o libgoto.so ~/libgoto_p4_128-r0.94.so xerb...
Rgoto と言うシェルスクリプト作成
#!/bin/bash
export LD_PRELOAD=~/libgoto.so:/usr/lib/libm.so
${RPROG:-R} -q $*
起動
$ ./Rgoto
で, R を実行すれば標準BLASやATLASの変わりに後藤(日本語訳...
元もとATLAS使っているなら劇的には早くはならないと思います...
//
-中間さんの処方箋通りにしてみました。1000 by 1000 の行列 ...
-FreeBSD派はいかにすれば・・・。というかATLASより速いんで...
-どうも goto ライブラリは CPU 毎に最適化チューンナップを...
-CPU最適化するのはATLASも同じでは。 -- [[aki]] &new{2004-...
-そうなんですか。そういえば Debian の atlas3 パッケージは...
-ATLASを最大性能(キャッシュのSIZEやCPU毎のオプテイマイズ)...
-FreeBSDのportsでは、ATLASもコンパイルからやれます。ただ...
-/proc/cpuinfo からL2ひろって atlasのsrpmを作れば嬉しい? ...
-cpuid から動的に最適なルーチンを呼び出すということは可能...
-もしかして噂の後藤さんご本人でしょうか。わざわざご投稿光...
-です。ライブラリに対する希望がありましたらお知らせくださ...
-よもや御本人から投稿していただけるとは考えず(豪雨汗)、失...
-なんかとんでもないトンチンカンなことを書いていたのでは？...
-今のところ、再配布は遠慮していただいています。ライセンス...
-LD_PRELOAD はすべてのプロセスの起動に影響を与えるため(プ...
-そうですね。Rgoto等のシェルスクリプトにくるむように書い...
-ぐわ、ちょっと見てない間に・・・。えー、FreeBSD用がある...
-ソースといっても、中身はほとんどアセンブラなので、FreeBS...
-げ、そうなんですか。それはまいった。うーん、Linuxに戻る...
-ほとんど違いはないので、あとは趣味の問題かと。FreeBSD 用...
#comment
**FreeBSDのportsから日本語化済Rをインストールする [#k371b...
>[[aki]] (2004-07-27 (火) 14:17:01)~
~
http://www.fifthdimension.jp/index.php?R%2FTips~
に半自動化スクリプトも含めて書いておきました。~
余計なこともやってますがとりあえずパッチの適用までは確認...
動作確認が取れましたらこちらのTipsにでもミラーしていただ...
//
-devPS.cのコンパイルでコケますね・・・(滝汗)。 -- [[aki]]...
-スクリプトの方がタコでパッチがちゃんと当たってませんでし...
-MAIN_CFLAGSがdevPS.cのコンパイル時に反映されてない? よう...
-R-1.9.1.l10n.psxfigtex.20040625.patchも微妙に改変しない...
-あ、すみません。devPS.cは直してあるんですが、公開し忘れ...
-では、うちのやつは消してしまっても問題ないのでしょうか。...
-今にして思えば日本語か掲示板の方が適当だったか・・・。 -...
-お手数おかけします。パッチ置いときました。ATLAS良いです...
-なかまさんのここのところの宣伝で気になってきていますが、...
-ldd usr/lib/R/bin/R.bin or usr/lib/R/bin/libR.so(共有で...
-そのままでは libgoto.so は見つかりませんでしたから、Debi...
-gotoはソースが無いからなんでは? -- [[なかま]] &new{2004-...
-FreeBSDマシンでも下に出ているtest2を走らせてみました。結...
-遅すぎるなーと思ってLinuxエミュを切ったら一つ目が1.15に...
-この結果はWinXPマシン同様2GHzにクロックアップした結果で...
#comment
**ELT [#wa7dca9b]
> (2004-07-27 (火) 12:37:19)~
~
率直に言って、ＲはＥＬＴ（※）には向いていないと思うのです...
?*nixな人だとAWK,PERL等になるのでしょうか？~
さすがに、SASをELTに使ってRで分析する人はいないですよね。~
~
※ELT~
Extraction（抽出）、Load（ロード）、Transformation（変換）~
//
-ELTってなんですか？ -- &new{2004-07-27 (火) 19:32:03};
-Every Little Thing でしょう(^_^;) -- &new{2004-07-27 (火...
#comment
**予約キーワード [#ube48076]
>[[青木繁伸]] (2004-07-27 (火) 01:23:05)~
~
名前付きの引数の指定の時，名前は前方から数文字で指定できる~
しかし，以下のようなこともある~
> x <- c(3.1, 2.4, 1.3, 4.2, 5.3)
> hist(x, breaks=1:6) # フルネーム（当然ながら，ちゃんと...
> hist(x, b=1:6) # b で始まる引数が2つあるので，識別不能...
Error in hist.default(x, b = 1:6) : argument 2 matches m...
> hist(x, br=1:6) # ちゃんと動く
> hist(x, bre=1:6) # ちゃんと動く
> hist(x, brea=1:6) # ちゃんと動く
> hist(x, break=1:6) # break は予約キーワードなので，エ...
Error: syntax error
ちょっと釈然としない。~
自作関数でも当然ながら，同じ問題は生じる。
> a <- function(fortran=0) print(fortran)
> a(f=3)
[1] 3
> a(fo=4)
[1] 4
> a(for=5)
Error: syntax error
> a(fort=6)
[1] 6
//
-全く解決策ではないですが、a(`for`=5) とバックチックでく...
#comment
**Eclipse Plug-in [#oedfdd0d]
> (2004-07-26 (月) 20:06:05)~
~
とうとう出ました。r-help に情報がありました。~
~
まだ試していません。~
~
Eclipse 3.0用で以下の機能があるようです。~
~
シンタックスのハイライト、括弧のマッチングチェック、IDE...
~
以下より入手可能です。~
~
~
http://www.walware.de/rplugin.zip~
//
-どうもXPではまったく動かない。どなたか実行に成功された方...
-[[ここ:http://www.walware.de/go?to=statet]]にプラグイン...
#comment
**TV や映画に出ていたR言語 [#f4cd2e51]
> (2004-07-24 (土) 01:51:01)~
~
てなことないでしょうね、~
//
-昔、S-PLUSの無料紹介ビデオ（英語・字幕なし）があったなぁ...
// -何でもとはいいながら、あまりに情報量のない話題はいか...
#comment
**浮動小数演算の怪 [#ud3aaae9]
>[[QDU]] (2004-07-21 (水) 11:01:10)~
~
これまた r-help 記事ネタ。次の一見無邪気な引き算の結果は...
> 0.6-0.3-0.3
[1] 0
> 0.6-0.4-0.2
[1] -5.551115e-17
> 0.6-0.5-0.1
[1] -2.775558e-17
> 0.6-0.2-0.4
[1] -5.551115e-17
> (0.6-0.2)-0.4
[1] -5.551115e-17
> 0.6-(0.2+0.4)
[1] -1.110223e-16
//
-2進表現で循環小数にならない10進数以外は計算機中の実数は...
Excel は，=3.6-1.2-1.4 でも 1 を表示しますので，単に「0に...
-R では次のようになります。この辺の機微もよう分からん。--...
> 3.6-1.2-1.4
[1] 1
> 1-(3.6-1.2-1.4)
[1] -4.440892e-16
-3.6-1.2-1.4 を 1-4.440892e-16 の値として保持するには 1 ...
#comment
**バージョン４のS言語はコンパイル機能を持つ！ [#m616a9c3]
> (2004-07-18 (日) 23:40:54)~
~
たまたま[[ベル研:http://cm.bell-labs.com/cm/ms/department...
//
#comment
**THE R BOOK の続編【R入門＋セミコマンドリファレンス】 [#...
>[[舟尾]] (2004-07-17 (土) 00:00:53)~
~
来年あたりに出版して頂けるであろう THE R BOOK の続編を私...
書籍として刊行するということで，ある程度は差別化を図る予...
[[http://cse.naro.affrc.go.jp/takezawa/r-tips/r.html]]~
です．そこで，もし良ければ皆様に R-Tips のダメ出しをお願...
//
-出版社が了解するかどうか怪しいですが、できた分から公開し...
-ご意見ありがとうございます．それはさすがに出版社の方に怒...
-日本語パッチで拡張された機能（日本語の変数等）も含まれる...
-日本語パッチで拡張されたのは日本語の為の設定ファイル(X11...
-貴重なご意見ありがとうございます．par 関係のオプションは...
-R-Tipsをベースにしたリファレンス本ということで、楽しみに...
-RjpWikiはご自由にお使いください。~
ありがとうございます！リファレンスを書く際には RjpWiki は...
-9月いっぱいでご意見の受付を終了いたしました．青木さん，...
#comment
**source() 関数のネスト [#afea9436]
>[[QDU]] (2004-07-12 (月) 22:08:24)~
~
r-help におもしろい話題がありました。普通 source() 関数を...
# ファイル "a.R"
a <- function(x) x^2
source("b.R")
~
# ファイル "b.R"
cat(a(10),"?n") # 単に a(10) では何故か何も表示されない？
~
# 実行結果
> source("a.R") # 中で source("b.R") を実行し、a(10) を...
100
#comment
**goto [#k4833353]
>[[なかま]] (2004-07-08 (木) 19:26:32)~
~
「ゴーツー」では無く「ごとう」なんですね。(^_^;
test2 <- function ( n=500 )
{
A<-array(rnorm(n^2), dim=c(n,n))
B<-array(rnorm(n^2), dim=c(n,n))
C<-array(rnorm(n^2), dim=c(n,n))
D<-array(rnorm(n^2), dim=c(n,n))
BA <- B%*%A
return (system.time(A%*%solve(t(BA)%*%BA+C)%*%BA%*%D ))
}
貧乏な私のマシン(ここでRのビルドとかしてるんですが)
> test2()
[1] 23.69 0.07 23.99 0.00 0.00
貧乏なRweb-jpマシン
> test2()
[1] 10.91 0.11 11.08 0.00 0.00
某所R-1.8.1
> test2()
[1] 9.822266 0.015625 9.839844 0.000000 0.000000
某所R-1.9.1 --with-blas=-lgoto
> test2()
[1] 1.01757812 0.02050781 1.04199219 0.00000000 0.00000000
ディスク欲しさにCPUをCeleronに落としてはいけない訳ですね...
//
-嗚呼、南北問題！ R.Ihaka 氏がどこかに「私はワークステー...
-[[libgoto:http://www.cs.utexas.edu/users/flame/goto/]] ...
-むーん、でもよくよく観察してみると警告が(環境?)。。。ATL...
-当家のDNS,IPSEC-SGはp5-133ですが、test2()は130程でした。...
-Rweb-jp の R をBLASにATLAS利用に変えてみました。[[Rweb-j...
#comment
**RjpWiki トップページについて [#w9418581]
>[[間瀬茂]] (2004-07-08 (木) 17:42:18)~
~
ここ数日私等がトップページをいじくっております。どうして...
//
#comment
**R で データマイニング [#d6c8e120]
> (2004-07-06 (火) 18:21:43)~
~
r-help に以下のサイトの情報が流れていました。~
~
http://www.liacc.up.pt/~ltorgo/DataMiningWithR/~
//
#comment
**漸近回帰モデル [#bb4ccc28]
>[[間瀬茂]] (2004-07-05 (月) 14:42:49)~
~
最近、R の stats パッケージに漸近回帰モデル (asymptotic r...
//
#comment
**誕生日のパラドックス [#e6150096]
> (2004-07-02 (金) 23:43:29)~
~
pbirthday は正しい値（？）を返さない~
これぞ，まさしく，パラドックス(^_^;)~
//
-実例は？ ？birthday によれば近似計算を使っているとの注意...
-以下を実行のこと。~
pbirthday2 <- function(n)
{
1-exp(lfactorial(365)-n*log(365)-lfactorial(365-n))
}
n <- 2:50
result <- cbind(n, sapply(n, pbirthday), sapply(n, pbirt...
colnames(result) <- c("n", "pbirthday", "pbirthday2")
result
実行結果（n が大きくても，小さくても，中くらいでも，不正...
n pbirthday pbirthday2
[1,] 2 0.005479452 0.002739726
[2,] 3 0.016438356 0.008204166
[3,] 4 0.026301370 0.016355912
[4,] 5 0.034246575 0.027135574
:
[21,] 22 0.499315068 0.475695308
[22,] 23 0.530136986 0.507297234
[23,] 24 0.531506849 0.538344258
[24,] 25 0.583542799 0.568699704
:
[47,] 48 0.957357026 0.960597973
[48,] 49 0.961562289 0.965779609
[49,] 50 0.967534247 0.970373580
関数本体はもちろん，1-prod((365:(365-n+1))/rep(365,n)) ...
//
-?birthday で説明されているように、ここで使われている近似...
-所詮，最終的な答え（確率）は，0〜1の範囲の値しか取らない...
-組み合わせ問題では(少なくとも途中では)すぐに超天文学的な...
-具体例があれば，きぼんぬ。 -- &new{2004-07-06 (火) 22:5...
#comment
**64 bit 版 R ！ [#e0d79fe1]
> (2004-07-02 (金) 23:01:19)~
~
r-help 記事より。お金持ちには６４ビット版 R で ４ Gigabyt...
Did you compile R as 64-bit executable on the Irix?
If not, R will be subjected to the 4GB limit of 32-bit s...
//
-RBookのbioconductorのところでも言及しましたが、同じ4GBで...
#comment
**anti-R vitriol [#c4c7cf18]
>[[間瀬茂]] (2004-07-01 (木) 17:10:54)~
~
昨日の r-help に表題のようなタイトル(vitriol というのは硫...
A colleague is receiving some data from another person. ...
SAS and it takes 30s and uses 64k RAM. That person then ...
and it takes 10 minutes and uses a gigabyte of RAM. Per...
It's not that I think SAS is such great software,
it's not. But I really hate badly designed
software. R is designed by committee. Worse,
it's designed by a committee of statisticians.
They tend to confuse numerical analysis with
computer science and don't have any idea about
software development at all. The result is R.
I do hope [your colleague] won't have to waste time do...
[this analysis] in an outdated and poorly designed piece
of software like R.
Would any of the "committee" like to respond to this? Or...
collective forehead and wonder how someone could get su...
>system.time(spam <- read.table("data_dmc2003_train.tx...
header=T,
colClasses=c(rep("numeric", 83...
"character...
[1] 15.92 0.09 16.80 NA NA
> system.time(spam <- read.table("data_dmc2003_train.tx...
[1] 187.29 0.60 200.19 NA NA
//
#comment
**なんで runif はランダムなの？ [#ka5cdc11]
> (2004-07-01 (木) 13:20:06)~
~
意味不明なタイトルと思われるでしょうが、つぎの結果をどう...
(つまり百万個一様乱数を発生する手間は同じはず、念には念を...
test <- function (i) {set.seed(222); system.time(runif(1...
x <- sapply(1:1000, test)
fivenum(xx)
[1] 0.24 0.28 0.30 0.30 0.47
//
-ほかに何のプログラムも走っていなければ同じ時間になるので...
Rweb:> system.time(for (i in 1:10000000) x <- 1)
[1] 7.01 0.10 7.13 0.00 0.00
Rweb:> system.time(for (i in 1:10000000) x <- 1)
[1] 5.85 0.07 5.92 0.00 0.00
これは，Ｒの問題ではないでしょう。UNIX の time コマンドも...
-なるほど、すると system.time でコード速度をはかる時は、...
#comment
ページ名:
![[RjpWiki]](image/../Rlogo.png "[RjpWiki]")
