なんでも掲示板アーカイブ（１）
をテンプレートにして作成
[
トップ
] [
新規
|
一覧
|
検索
|
最終更新
|
ヘルプ
]
開始行:
COLOR(blue){SIZE(20){なんでも掲示板アーカイブ（１）}}~
SIZE(20){COLOR(red){注意}}：ここに保存した記事へのコメン...
元記事日付 2003-07-17 から 2003-11-26 分 ~
#contents
~
**行列データの抜き出し方(時系列) [#ef3afdea]
>[[さかもと]] (2003-11-26 (水) 10:08:12)~
~
ある時系列データセット(1980年から2002年)を~
>ts(data="data set name",start=1980,end=2002,frequency=1)...
として読み込んだ後，ある変数の1980年のデータを取り出した...
>D[1980,] # ← 1980年のデータすべてを表示としたいのですが...
確認できていることは，~
>D[1,] # ← だったら1980年のデータを参照することができ...
年の指定ができるとそうさもしやすいのですが･･･．なにとぞご...
//
-window(D,start=1980,end=1980)ではどうでしょう。 -- [[な...
-アドホックですが、頻繁に年度で参照したければ次のように専...
> library(ts)
> data(sunspot)
> sunspot
Time Series: Start = 1936 End = 1972 Frequency = 1
[1] 40 115 100 80 60 40 23 10 10 25 75 145 130 ...
[20] 10 60 190 180 175 120 50 35 20 10 15 30 60...
> ofs <- function(n) {n+1-attributes(sunspot)$tsp[[1]]} ...
> x[ofs(1936)]
[1] 40
-使いやすいかどうか分かりませんがラベルをつけるという手も...
> names(x)=1936:1972
> x["1936"]
1936
40
-早速のリプライ多謝でございます．早速試させていただきます...
#comment
**tests/p-r-random-tests.R [#x6d4cee4]
>[[なかま]] (2003-11-25 (火) 23:48:51)~
~
どうも乱数部分に問題?があるような気がするんですが目がショ...
追えません。(爆)~
for ((i=0;i<=100;i++))
do
R -q --vanilla < R-1.8.1/tests/p-r-random-tests....
done
とか
for ((i=0;i<=100;i++))
do
Rterm.exe -q --vanilla < R-1.8.1/tests/p-r-rando...
done
~
等とすると、テストがFAILする場合がたまにあります。~
Rの中のMTはsrc/main/RNG.cの中にあって、~
A C-program for MT19937: Real number version([0,1)-inte...
(1999/10/28)
なバージョンです。~
どなたか、ハックしませんか。(^^;~
# 明日から客先なんで。。。~
//
-行く前に。一晩寝てみて、これってFAILEDする事もあるって論...
#comment
**CODAとBOA どちらがおすすめ？ [#we982c22]
> (2003-11-25 (火) 18:39:33)~
~
BUGSのサンプルデータで作成した出力ファイル(*.out,*.ind)を...
(coda_0.5-14.tar.gz)とBOA(boa-1_0_1R.tgz)の両方を試してみ...
どちらも同じような感じだったのですが、それぞれの特徴的な...
うか？どちらがおすすめでしょうか？~
どちらも最新版の日付は2003年1月で開発が止まっているように...
//
#comment
**SOMパッケージの違い [#h23c02cb]
>[[amaki]] (2003-11-25 (火) 06:17:59)~
~
Rには"class"と"som"の2種類のsomパッケージが入っていますが...
どのように違うのか教えてください。~
"class"がコホネンのsom_pakを基にしたもので、~
"som"は別物と考えてよろしいのでしょうか？~
~
今、100000ベクトルほどあるデータをsomパッケージでマップ化...
とても時間がかかってしまい(学習回数、1000,10000で半日)、~
困っています。~
学習効率、プロットの精度の観点から両者の違い等ありました...
どなたか教えて下さい。~
~
デザインの領域のモノで初心者です。~
ご助言、よろしくお願いします。~
//
-少し表現が分かりにくかと思うので訂正します。現在使ってい...
#comment
**cygwinでRはできる？ [#wac2a965]
>[[taku]] (2003-11-24 (月) 20:17:51)~
~
ｃｙｇｗｉｎでＲがやりたいとおもい、~
http://datamining.tama.ac.jp/~yama/R/install.html~
を参考にインストールしようとしていましたが、~
ｍａｋｅで失敗してしまいます。~
make[3]:*** [devices.o] Error1~
make[2]:*** [R] Error2~
make[1]:*** [R] Error1~
make:*** [R] Error1~
cygwinではRはできないのでしょうか？~
~
よろしくお願いします。~
//
-確か、Cygwin は現在サポートしていないのではないですか？W...
-そうなのですか。ありがとうございます。cygwinではできない...
-以下の参考ページを見つけております。~
http://www1.accsnet.ne.jp/~gen/R-rev/cygwin-R.html~
http://www1.accsnet.ne.jp/~gen/R-rev/cygwin-R.pdf~
私自身まだうまくいってませんので、うまくインストールされ...
-CygwinとMingwではいくつかのヘッダーの位置が異なるのと、d...
-いろいろありがとうございます。参考ページのとおりにまずや...
-とりあえず[[これ:http://www1.accsnet.ne.jp/~gen/R-rev/cy...
#comment
**RGUI内でWindows の環境変数の設定できますか？ [#te479ecd]
>[[SJava]] (2003-11-23 (日) 20:16:18)~
~
現在SJava をWindows内で使用しています。~
~
SJava を動作させるには、環境変数"JAVA_HOME"にJREのパス...
~
現在マシーンにはJavaSDKを導入していて、このSDKのパスをJ...
~
SJAVAを使用するときだけ、JAVA_HOMEのパスを変更したいの...
//
-1.バッチでくるむ。2.Sys.putenv("JAVA_HOME"="hoge")。3."....
-なかまさん、ありがとうございます。 -- &new{2003-11-27 (...
#comment
**Rでベイズ統計入門 [#r10cd017]
> (2003-11-21 (金) 03:18:38)~
~
Rを使ってベイズ統計を初歩から学べるようなサイトってある...
//
-ベイズ統計の意味によりますが、NIST にベイズ統計の実際的...
#comment
**多変量解析の演習門題 [#v681cae3]
>[[宮島 のぼる]] (2003-11-20 (木) 14:38:48)~
~
(X,Y){X=xi,Y=yj}でX,Yが独立でない時ρxy＝０となる例は？~
（ただし、i,j=1,2,3,‥‥）~
//
-これはもしかして学校のレポートの課題ではないでしょうね。...
-少々きついコメントを書いてしまいましたが、もし真面目な質...
#comment
**Rgraphviz on R for Win [#l0dcb804]
> (2003-11-19 (水) 15:22:42)~
~
BioConductor のgraphvizを RのWin版で動かされている方、...
~
MinGWでgraphvizをビルドできなくて、現在指をくわえたまま...
//
-cygwinでビルドできませんか -- &new{2003-11-19 (水) 21:5...
-開発版(>=R 2.1.0)が http://www.bioconductor.org/packages...
-[[正規版:http://www.bioconductor.org/packages/bioc/1.7/s...
#comment
**統計数理研究所ワークショップの発表内容のWebでの公開 [#w...
> (2003-11-17 (月) 22:37:40)~
~
統計数理研究所のワークショップですが、平日なので拝聴す...
~
このワークショップでの発表はPDFやPPTで、Webで公表されな...
//
-代表者に話はしておきますが、そうした形で発表者が原稿を用...
#comment
**誕生月でスポーツ能力が決まる？！ [#ua44e33f]
>[[間瀬]] (2003-11-17 (月) 22:25:52)~
~
今日の日経新聞に面白い記事がでていました。今期の日本の J...
> x=c(23,20,21,51,57,43,46,42,46,38,26,19)
> chisq.test(x)
Chi-squared test for given probabilities
data: x
X-squared = 55.3889, df = 11, p-value = 6.578e-08
~
となりやはり強く有意なようです。これで我が家の男三人がス...
//
-二乗では？出生数の分布も一様でないそうですが、それを調整...
-自乗も二乗も特に区別されずに使われているのでは。私はむし...
-おお！わたしも長年の謎(自乗か二乗か)が解けてすっきりしま...
-ワルノリして更に。昭和４５年以降月別出生数は殆んど同一と...
-カイ二乗が約2,830件,カイ自乗が約722件でした。-- google &...
-これほどの量をどのように集計したのですか。まさか一つ一つ...
-私の愛用のウェッブ英和・和英辞典のスペースアルクによると...
-常用漢字ではニをジ(漢語)と読まないので国語的には自乗が正...
-ちなみに、二乗派はどう発音しているのでしょうか？「にじょ...
#comment
**CRAN での自作パッケージの登録方法 [#f9382483]
> (2003-11-16 (日) 16:09:14)~
~
CRAN で自作のRパッケージを登録したい場合、どのような手...
//
-日本人によるパッケージがあればいいなと思っていました。も...
#comment
**[[Rから他言語利用]] の gmrf.c のフルソースコード希望 [#...
> (2003-11-14 (金) 14:44:19)~
~
現在、オープンソースの C++ で作成されたパッケージをRよ...
~
ところで、[[Rから他言語利用]]にある、 gmrf.c のフルソー...
//
-構いませんが、RjpWiki に載せるには少し長過ぎるかと。御希...
#comment
**回帰分析に関する日本語文献 [#a7caada0]
>[[間瀬茂]] (2003-11-11 (火) 20:05:53)~
~
教えて下さい。重回帰分析 and/or 非線形回帰分析に関する良...
//
-間瀬先生，こっちでもお世話になります．え??，非線形回帰に...
-もし Amazon でのキーワード検索結果が完璧だとすると、現在...
-英語で理論の本と意味の本とでは、どんな本が良いのでしょう...
-すみません，上の書き込みで"え??"となったところは，文字化...
-それはさておき，確かに回帰分析を勉強する本はないですね．...
#comment
**x <- c(....) vs x <- read.csv(...,...) [#b25991fe]
>[[まさ]] (2003-11-11 (火) 14:13:28)~
~
prec <- c(3.2, 3.1, 5.1, 4.8, 8.3, 9.8, 8.3, 6.6, 5.1, 3....
prec <- read.csv("tempdata.csv",header=T)として以下のこと...
op <- par(mfrow=c(2,2))~
#ヒストグラム~
hist(prec,main="Histgram")~
を実行すると~
Error in hist.default(prec, main = "Histgram") : ~
`x' must be numeric
と返ってきます。~
listで確認したところ~
prec <- c(...)では~
[[1]]~
[1] 3.2 3.1 5.1 4.8 8.3 9.8 8.3 6.6 5.1 3.8 5...
[16] 7.4 6.2 7.6 7.5 7.0 6.2 6.0 8.8 10.7 8.5~
となり、グラフが作成できました。~
prec <- read.csv("tempdata.csv",header=T)では~
[[1]]~
temp
1 3.2~
2 3.1~
3 5.1~
4 4.8~
5 8.3~
6 9.8~
7 8.3~
8 6.6~
9 5.1~
10 3.8~
11 5.2~
12 5.6~
13 6.5~
14 5.7~
15 5.7~
16 7.4~
17 6.2~
18 7.0~
19 6.7~
20 5.7~
21 6.2~
22 6.0~
23 8.8~
24 10.7~
25 8.5~
となり、グラフが作成できません。~
prec <- c(3.2, 3.1, 5.1, 4.8, 8.3, 9.8, 8.3, 6.6, 5.1, 3....
temp~
3.2~
3.1~
5.1~
4.8~
8.3~
9.8~
8.3~
6.6~
5.1~
3.8~
5.2~
5.6~
6.5~
5.7~
5.7~
7.4~
6.2~
7~
6.7~
5.7~
6.2~
6~
8.8~
10.7~
8.5~
のデータで実現するにはどうしたらよいでしょうか？~
//
-read.csvで付値したらデータフレームになるので，変数名も指...
-長くなって他のトピが見づらくなるので、編集お願いします。...
#comment
**エクセルの四捨五入 [#ie6f07e7]
>[[かずっち]] (2003-11-11 (火) 11:05:43)~
~
エクセルで202.85-201.4=1.45で小数１桁でroundをかけると~
1.5にならなければいけないのに1.4になります何度計算をして...
何か設定があるのでしょうか？~
//
-Rでもそうですね。round(1.35,1)は1.4になるし，round(1.55,...
-R の help(round) によれば, round は 'IEEE standard' (IEE...
-かずっちさんの事例は Excel VBA 関数での現象と思いますが,...
-Excel 固有の疑問はこの Wiki にはふさわしくないかと愚考し...
-なるほど。1.45から1.4までと1.5までは差が同じなので，末尾...
#comment
**RMySQLをwindowsで使いたい [#q2037c66]
>[[sakamoto]] (2003-11-06 (木) 12:59:50)~
~
RMySQLがwindowsで使えるようにした方はいらっしゃいませんか...
//
-RMySQL/inst/{INSTALL.win,README.windows}に従えば良いよう...
-なかまさんありがとうございます。早速試してみたいと思いま...
#comment
**Rで0-1整数計画問題を解く [#x18848a5]
> (2003-11-01 (土) 23:07:53)~
~
[[Rで数理計画]]を参考にさせていただいておりますが、とこ...
//
-パッケージlpSolveのlpは整数計画問題も対応できるはずです...
-ken さん、ご教授ありがとうございます -- &new{2003-11-04...
#comment
**SAS ー＞ R [#ief8272c]
> (2003-10-29 (水) 12:05:04)~
~
R-jp の記事を転載します。R＋SAS両刀使いの方いませんか。~
~
RをWindowsで勉強しています。さてよくSASを使うとTypeII SS...
//
-carパッケージに含まれています -- &new{2003-10-29 (水) 2...
//
-中澤さんの [[Rtips:http://phi.ypu.jp/swtips/R.html]] に...
#comment
**Rgui（非日本語パッチR for Win）のフォント指定をコマンド...
> (2003-10-29 (水) 10:48:04)~
~
日本語パッチを当てていない、Windows用のR（1.8.0）で、"R...
//
-R_HOME/etc/Rconsole -- [[なかま]] &new{2003-10-29 (水) 1...
-なかまさん、早速のレスありがとうございます。ところで、コ...
-もしかして、Rのコマンドラインから変更してダイレクトに反...
#comment
**[質問]名前に利用できる文字 [#c806f207]
>[[_（アンダーバー）が使えると良いな]] (2003-10-29 (水) 1...
~
質問というよりは、疑問なのですが...~
なぜ、名前につかえる文字に _（アンダーバー）がなく~
.（ドット）が含まれているのですか？~
歴史的な理由でもあるのでしょうか？~
//
-付値(代入)の演算子だったからです。 -- [[なかま]] &new{20...
-R 1.8 からは _ は使えなくなったそうですが、ついでに変数...
-変数名に"_"は使えません。essを使えば、"_"は"<-"に入力時...
-1.8.0から`が加わり、QuotedSymbolValueの中はStringValue(c...
-最近の r-help 記事によると数カ月後には名前の一部としてア...
#comment
**Rの列数の限界値について [#k8e82294]
>[[sallyh]] (2003-10-28 (火) 16:09:52)~
~
こんにちは。~
Rの列数の限界値について調べています。~
read.tableでデータを読み込ませようとした時、列数が30000列...
いくらメモリを増やしても読み込むことが出来ませんでした。~
~
そこで、Rの列数の上限を調べているのですが、調べ方が悪いら...
良く分かりません。~
（エクセルは確か256列でしたよね？）~
~
もしかして、ものすごく基礎のことかもしれませんが、~
ご存知の方いらっしゃいましたら、もし良ければ教えて下さい。~
できれば、Windows版とUnix版で違いがあるのか等も知りたいで...
~
よろしくお願いします。~
//
-答えではないのですが，help(read.table)すると，read.table...
-x=matrix(1.0,c(57000000,57000000)) はOK(すごいけど、おそ...
-列数に制限は無いです。read.tableは列名を得るために何度も...
-ありがとうございます。助かりました。参考にします！！ -- ...
#comment
**Rでマルチ・エージェント [#y282da6f]
>[[M.A.S.]] (2003-10-23 (木) 00:46:19)~
~
マルチ・エージェントというと、SWARMが有名ですが、Rから...
~
CRANでは見当たりませんが？~
//
#comment
**RExcel【(D)COM】の使い方に関して [#u98e2623]
>[[激・初心者です]] (2003-10-21 (火) 17:13:18)~
~
R Start をクリックすると、~
『Fatal error:unable to open the base package』⇒ ~
『Could not start RServer』⇒~
『実行時エラー'462'リモートサーバがないか使用できる状態で...
というエラーが出てしまい、何もできません。~
~
原因がお判りになるかた、レスを頂ければ幸いです。~
//
-ちなみに、RExcel.xlaはきちんとアドインされているようです...
#comment
**[質問] Rの分散分析について [#g15514e3]
>[[R初心者]] (2003-10-20 (月) 21:57:10)~
~
Rの初心者です．つまずいていることがあります．~
~
分散分析をする際に，~
「母数モデル・変量モデル・混合モデル」の~
いずれかを指定したいのですが，~
変量モデルを指定する関数（SASでいうVARCOMP）は，あるので...
//
-私自身は使う必要がないので詳細はわかりませんが、こうした...
and Ripley (邦訳あり)に関連する話題がありそうですから見て...
On Wed, 29 Aug 2001, Troels Ring wrote:
> Dear friends.
>
> I know this has been discussed before but as lively as...
> things might have changed. In Venables and Ripley(199)...
> of variance components including robust methods is dis...
> anything like that presently but might have overlooked...
> do it ? Is it obsolete ?
It is not there in R and not planned. (You can check pla...
developer.r-project.org.) It is still in S-PLUS, so not ...
-- &new{2003-10-20 (月) 22:10:45};
-分散分析 randam and mixed については、以下の URL にある...
-おお．レスがこんなについてたとは・・・nlmeに当たってみま...
#comment
**text関数を使って‥ [#pda2309a]
>[[R初心者]] (2003-10-18 (土) 15:29:31)~
~
text関数を使って描いたグラフのある位置に、求めた重要なパ...
具体例をあげますと‥~
text(100,2,parameter)とすると当然画面上には変数parameter...
//
- 例えば次の例が参考になりますか。paste 関数を使って単一...
> x=runif(10)
> plot(1:length(x), x)
> legend(2, 0.5, paste("no of point", 10, " ")) # 最...
> legend(2, 0.7, paste("no of point", length(x), " "))
#ref(test.png, left)
-- &new{2003-10-18 (土) 18:25:19};
#comment
**[解決]table()で総和を出す方法 [#tba12ab6]
>[[初心者2]] (2003-10-18 (土) 05:34:19)~
~
を考えています。~
table(c(1,1,2,2,3))
とかで、一番最後にTotal 5と表示させたいのですが、いい方法...
//
- テーブル化された結果から元のデータ数を計算したいという...
> x<-table(c(1,1,2,2,3))
> sum(x)
[1] 5
-table関数の仕様を変えたいという意味です。 -- [[初心者2]]...
-ttab <- function(x) { c(table(x),sum(table(x))) } として...
-ttab(floor(runif(10000,1,7)))でもちゃんと一番右に10000と...
-ありがとうございます。ttabを使わせてもらいます。 -- &ne...
-NAを数えたい方は、-- &new{2003-10-21 (火) 11:39:01};
ttab <- function(x) { c(x<-table(x,exclude=NULL),sum=sum...
#comment
**lsfit関数について [#x1128534]
>[[R初心者]] (2003-10-17 (金) 16:04:08)~
~
lsfit(x,y)$coefで回帰係数であるaとbがX，Interceptとして得...
//
-すみません。過去の「最初の一歩：必要な出力だけ欲しい時」...
#comment
**R と Ruby [#m1a1835d]
> (2003-10-15 (水) 10:08:18)~
~
RPy (R from Python) http://rpy.sourceforge.net/about.html~
このようなのがあるのですね。~
Ruby にもあるのでしょうか？ 知っていたら教えてください。~
//
-[[Rとインターフェースのあるアプリ]]にリンク先があります...
-R Math library を Ruby からつかうものは実装があります。R...
#comment
**コホネンのSOM_PAKをRから利用できるパッケージ [#p843ef29]
> (2003-10-10 (金) 18:56:10)~
~
コホネンのSOM_PAKですが、これをRで呼び出して利用できる...
~
//
-自己レスです。class パッケージにSOMはありますね。 -- &n...
-som(旧名 GeneSOM) というパッケージもありますね。 -- &n...
#comment
**[解決] 配列に収めたデータを正規化したいと思っているので...
>[[激・初心者です]] (2003-10-09 (木) 11:40:20)~
~
全くどうしたら良いかわかりません。~
列ごとに平均・標準偏差を求めて正規化を行いたいのですが、...
階層的クラスタリングに用いるデータで、その後の処理(Euclid...
VisualBasicで自作したツールでは不便でRを利用しようと考え...
どなたか、気が向いた方、返答頂ければ幸いです。~
//
-sweep()（http://finzi.psych.upenn.edu/R/library/base/htm...
-関数を簡単に定義できるので、関数を自作し、applyすればい...
-早速、ありがとうございます。 -- &new{2003-10-09 (木) 16...
-上記返答を参考に模索してみます(すぐ利用出来るほど賢くな...
-例えば
> x <- matrix(runif(20), c(5,4))
> xx <- scale(x, apply(x,2,mean), apply(x,2,sd))
行について操作したければ、y <- t(x) と転置しておいて、最...
-素晴らしいです、出来ました！ご丁寧にありがとうございます...
#comment
**中澤先生の「Rによる統計解析の基礎」の出版社 [#obb2faac]
> (2003-10-06 (月) 13:23:35)~
~
件名は１０月に発行されるようですが、Amazonではヒットし...
~
どこの出版社から出るのでしょうか？~
//
-株式会社ピアソン・エデュケーションのようです。 -- &new{...
-中澤 港：「Rによる統計解析の基礎」ピアソン・エデュケー...
#comment
**カーネル密度推定のわかりやすい文献 [#ib1a5b12]
> (2003-10-01 (水) 11:22:20)~
~
R言語はいくつかのパッケージで多くの種類のカーネル密度推...
//
//-この御質問に対する答えとしては、 A. Bowman and A. Azza...
//-書き込みの衝突が起きたそうなので、もう一度書きます。
-この御質問に対する答えとしては、 A. Bowman and A. Azzali...
#comment
**[解決]二つのグラフにまたがるタイトル [#b4ecaf4e]
>林 (2003-09-30 (火) 06:04:46)~
~
お世話になります。
#ref(graphs.png)
のようなふたつのグラフにまたがるようにタイトルをつけたい...
//
-layoutは使った事が無いんですがだいたい以下のようにしてま...
par(mfrow=c(1,2),oma=c(0,0,2,0))
deg<-seq(-180,180)
plot(deg,sin(deg/180*pi),type="l",main="sin")
plot(deg,cos(deg/180*pi),type="l",main="cos")
mtext(side=3,line=0,cex=1.4,font=2,outer=TRUE,"sin cos")
-なるほどうまい。ついでですが、グラフィックス実例集で lay...
-おっ、有難うございました。いつもながら、やはり関数がある...
-上手く行かないのはpngだけであればpngを疑うしかないかとwi...
-なぜ X11 デバイスには出力されるのに png デバイスだと複数...
test() # なかまさんの例を X11 デバイスに出力
pp <- recordPlot() # X11画面への作図手順を記録
png("test.png") # png デバイスを開く
replayPlot(pp) # png デバイスに出力
dev.off() # デバイスを閉じる
-- [[間瀬]] &new{2003-10-01 (水) 21:33:13};
-間瀬先生=debianと解釈。libpngのバージョンはいくつなんで...
-ローカルなトラブルですか。r-help を探しても報告は見つか...
-うちのだと、system(paste("ldd ",R.home(),"/modules/R_X11...
#comment
**[解決]リストの名前タグ [#qf3828bd]
>[[dimnamesのように]] (2003-09-21 (日) 11:01:04)~
~
リストにあとから名前タグをつけるには、どうすればいいので...
//
- names 関数を使います(変更にも使えます。一部だけ使うには...
> x=list(1:3,matrix(1:4,c(2,2)))
> x
[[1]]
[1] 1 2 3
[[2]]
[,1] [,2]
[1,] 1 3
[2,] 2 4
> names(x) <- c("vector","matrix")
> x$vector
[1] 1 2 3
> x$matrix
[,1] [,2]
[1,] 1 3
[2,] 2 4
> x[["vector"]]
[1] 1 2 3
> names(x)[1] <- "vector" # 一部指定にはこのようなやり方...
> names(x)[1] <- "" # 名前タグを(見かけ上)取り去る
-x$matrixには、colnames関数やrownames関数も使えました。あ...
#comment
**[解決]関数lagを使いたいのですが [#hc5921db]
>[[lager]] (2003-09-19 (金) 16:02:15)~
~
時系列分析をやっているのですが、~
RはSPlusのうち時系列分析関係の関数がかなり省かれているよ...
（たとえばtsplotなども）~
lagを使いたいのですが、どこかにないものでしょうか。~
//
-すみません、パッケージでtsを入れたらできました。 -- [[la...
#comment
**[解決]メモリの問題 [#eba4994a]
>林 啓一 (2003-09-18 (木) 03:20:51)~
~
メモリが512Mのマシンを再起動後、Rのみを立ち上げて~
7万＊60変数ぐらいのデータフレームをつくり(下のやつです)~
他のソフトで分析するためにwrite.tableで書き出そうとすると~
メモリが足りないといわれてしまします。~
~
メモリを変えるor小さく分けて書き出す以外になにか解決方法~
ありませんでしょうか？~
//
-メッセージにもよりますが既に下の方の、Reached total allo...
-write.tableがメモリを使うのは x <- as.matrix(x) でmatrix...
-memory.limit(size=1024)で17Mの.txtファイルの書き出しに成...
http://www.geocities.co.jp/SiliconValley-SanJose/3220/XP-...
を読むといらないサービスばっかりですね＞M$
#comment
**CRANのakimaパッケージ内のinterp関数を入れたんですが…。 ...
>[[wash-and-go]] (2003-09-17 (水) 18:21:17)~
~
うまくいきません。途中でエラー出てしまいます。~
~
>Error in .Fortran(as.single(x), as.single(y), as.single(...
> as.integer(nx), : ~
> function name must be a string (of length 1)~
~
といった感じです。他にもパッケージが入っておらずエラーが...
でしょうか？どれを入れていいかもわかりません。~
//
-お使いの環境はなにか分かりませんが~
~
R1.7.1のWin版では、~
library(akima)~
test<-interp(rnorm(10),rnorm(10),rnorm(10))~
image(test$x,test$y,test$z)~
~
で、なんら問題は出ませんが。 -- &new{2003-09-17 (水) 1...
-このパッケージは Fortran コードを使っています。もし Linu...
-R から Fortran を呼び出す場合、.Fortran("関数名文字列",...
-R1.7.0のWin版です。Fortranのコンパイラは入っていないと思...
-たしか、ここから落としてきました。http://www.bioconducto...
-R_HOME/library/akima/DESCRIPTION の Buildには何と書いて...
-http://www.bioconductor.org/CRAN/の1.7のディレクトリの中...
-ftp://ftp.u-aizu.ac.jp/pub/lang/R/CRAN/bin/windows/contr...
-インストール時に古いakimaは消されたのでしょうか、akimaの...
-消しましたがうまくいきません。一度、バージョンを新しくし...
-バージョンを1.8.0にしてakimaを入れてみたらうまくいきまし...
#comment
**NP問題 [#a24cd16e]
>[[なかま]] (2003-09-17 (水) 14:05:38)~
~
どこかにNP問題の一覧みたいなのって無いんでしょうか。~
Rで紐解いてるのがあるともっと嬉しいです。~
//
-一覧があるかどうかは知りません。R の base パッケージには...
-素晴らし過ぎます。でも、キーワードが思い付かないと本当に...
-ga は http://www.dst.unive.it/~claudio/R/index.html#ga ...
-CRAN パッケージ gafit は曲線当てはめという特殊目的のよう...
#comment
**[解決] 効率のいいデータの変形(stack) [#l5cbf661]
>林 啓一 (2003-09-16 (火) 13:13:43)~
一人6回観察されている1万2千件のレコードを1万2千*6回という...
> memory.size(30000000)
[1] 829953112
> for(id in 1: nrow(data)){ for(t in c(15,22,33,46,50,51...
newdata<-rbind(newdata,cbind(t(data[id,]),outcome=...
Error in rbind(newdata, cbind(t(data[id, ]), outc..: Un...
memory (this request is for 2725680 bytes, 830035032 byt...
//
-いまひとつイメージがつかめませんが、t(data[id,])はt(data...
-もちろん perl 等を使ってデータファイルを preprocessing ...
-loopは出来るだけ使わない、cbindやrbindの繰り返しでデータ...
S-PLUS for UNIX/Linux Version 6 Programmer's Guide ~
http://www.msi.co.jp/splus/support/download/ ~
のp.509-14に書いてありました。必要なデータセットは
y<-c(15,22,33,46,50,51)
system.time(stackdata<-matrix(rep(as.vector(t(as.matrix(...
outcome<-as.vector(as.matrix(data[,y]))
time<-sort(rep(1:6,12686))
cbind(stackdata,outcome,time,unique=1:12686*6)
で作成できました。お騒がせしました。2行目のn*mの行列をた...
-単純にデータフレームを縦に6個並べると言う意味でしょうか?...
-なかまさんありがとうございます。
stackdata<-rbind(data,data,data,data,data,data)
が一番早そうですね。[[Rでいいコードを書く:http://plaza.um...
//なかまさんレスありがとうございます。あの件はどうなりま...
#comment
**SASの番組 [#mf75b10b]
> (2003-09-15 (月) 20:19:54)~
~
Rと直接関係がありませんが、いま、NTVで有名な福利厚生につ...
//
-番組を見ていないでの意味不明です。SAS(あの統計ソフトです...
-「世界まる見え！テレビ特捜部」（20:00-20:54）の中の「仰...
-R関係と言う事で、Rの日本語化をしている人の福利厚生を紹介...
#comment
**データフレームからの要素の取り出し [#n6f8c8c8]
>[[竹澤]] (2003-09-10 (水) 11:32:56)~
~
竹澤（中央農業総合研究センター・農業情報研究部）です。~
~
以下のＲオブジェクトを実行しました。
function ()
{
xx <- data.frame(x1=c("a1","b1","b1","d1"), x2=c("a2","b...
print("xx=")
print(xx)
yy <- xx[,1]
print("yy=")
print(yy)
zz <- xx[,2]
qq <- c(yy,zz)
print("qq=")
print(qq)
rr <- c(as.vector(yy), as.vector(zz))
print("rr=")
print(rr)
}
結果は以下のものです。
[1] "xx="
x1 x2
1 a1 a2
2 b1 b2
3 b1 b2
4 d1 b2
[1] "yy="
[1] a1 b1 b1 d1
Levels: a1 b1 d1
[1] "qq="
[1] 1 2 2 3 1 2 2 2
[1] "rr="
[1] "a1" "b1" "b1" "d1" "a2" "b2" "b2" "b2"
データフレームからの要素の取り出しには注意が必要なのです...
ってなことは自明だ、という方が多いのでしょうか。~
//
-おっしゃりたいことを分かりにくいので、もう少し説明しても...
-http://buran.u-gakugei.ac.jp/~mori/LEARN/R/R-intro-jp.17...
-なるほど。 -- [[竹澤]] &new{2003-10-03 (金) 08:55:22};
-「6.3 データフレーム」の「数値ベクトルと論理値、因子はそ...
#comment
**2chぽいですが [#j85f7958]
> (2003-09-08 (月) 19:53:26)~
~
http://epio.jpinfo.ne.jp/1057/index.html~
グラフがExcelのようですが、Excelで求めたんですかね。~
//
#comment
**[質問]--vsize --nsize オプションの使い方について… [#n88...
>[[INA]] (2003-09-04 (木) 11:32:47)~
~
初心者ですが、経済分野の「産業連関分析」で、Rを利用させ...
ただ私の環境では、少し大きなデータを読み込ませるとすぐ...
~
どなたか似たような経験をなさった方はおられないでしょう...
~
「産業連関分析」では作業の手始めとして「工業統計調査」...
つまり~
~
品目1 品目2 品目3… 計
産業分類1~
産業分類2~
産業分類3~
：~
計~
~
という具合です。~
~
データの形式は~
~
'産業分類','品目分類','出荷額'←実際は英語です~
'1843','184312',33000~
'3494','349419',1125~
'1931','193111',10979~
：
~
の形で続いており、約59万レコードあります。（サイズは７MB...
~
これを~
~
>data <- read.csv("c:?Hogehoge.csv",header=TRUE,sep=",")~
~
で読み込もうとしているのですが、どうしても~
~
$"Reached total allocation of 126Mb: see help(memory.size...
~
のメッセージが出て、最後まで読み込むことができません。~
~
念のためvsizeオプションを使い、コマンドラインから~
~
>rgui --vsize 36M ~
~
として走らせて見ましたが、やはり読みきれません。~
~
~
ちなみに同じようなデータ、約18万レコードは問題なく読み込...
（その後、>tapply(data$出荷額,list(data$産業分類,data$品...
ご存知の方、どうぞご指導のほどよろしくお願い致します。~
//
-GlobalMemoryStatus から、R_max_memory = min(1024 * Mega,...
-R上からでも、memory.limit()で認識している物理メモリサイ...
-ご指導有り難うございました。雷の影響かチップがおかしくな...
-どういうデータなのかにもよると思いますが、ちょうど126M以...
-PCをパワーアップして1Gまでメモリを積んだところ、問題なく...
#comment
**最初の一歩：必要な出力だけ欲しい時 [#r26634af]
>[[角煮]] (2003-09-02 (火) 10:38:14)~
~
初心者の私が、必要な出力だけ欲しい時、悩んでしまいました...
~> x<-c(1,3,5) ~
~> y<-c(3,5,4) ~
~> wout<-lsfit(x,y)$coef~
~> wout[1]~
Intercept
3.25
~> wout［［1］］~
[1] 3.25
~> cat(wout[1])~
3.25
これで、やっと数値3.25だけ得られました。（ここまで来るの...
//
-lsfit(x,y)$coef["Intercept"] -- &new{2003-09-02 (火) 11...
-私もてこずりました。リスト型式は慣れると便利です。 ~
Betaの書き出し http://plaza.umin.ac.jp/~epi/cgi-bin/fswi...
にも書きましたが、SASやSTATAでは簡単に書き出せない統計量...
(違ったらすみません、)書き出しやすいのがR(S)のいいところ...
１．str() でリストの中身 (階層構造も) が参照できる~
２．リストは$名前で参照できるとともに、
[[順番]][[その階層での順番]]
でも参照できるということがポイントだと思います。私なら
str(lsfit(x,y))
lsfit(x,y)[[1]][[1]]
でその後の解析に使います。cat()までして数字したいという状...
#comment
**樹形モデル [#sfa9a1e8]
> (2003-08-30 (土) 11:56:24)~
~
S-plusのtree関数に相当するものは、Rのパッケージに存在する...
アドバイスよろしくお願いします。~
//
-http://cran.r-project.org/src/contrib/PACKAGES.html#tree...
#comment
**software patent in Europe [#i6b09b58]
> (2003-08-27 (水) 19:27:59)~
~
9月１日、"software patent "の導入について欧州議会で投票が...
//
-http://swpat.ffii.org/index.ja.html の事でしょうか? -- [...
#comment
**Splusでのinterp()の関数はRにはありますか？ [#k079aa19]
>[[wash-and-go]] (2003-08-25 (月) 16:15:24)~
~
影響関数を描くのに必要な関数なんですが、~
Rにはinterp()という関数はあるでしょうか？~
//
-CRANのakimaパッケージ内のinterp関数（不規則分布のポイン...
#comment
**gstatの日本語マニュアル [#m9966bc7]
>[[丸山]] (2003-08-22 (金) 13:47:42)~
~
CRANに多変量地球統計学のパッケージとして提供されているgst...
レフェレンスマニュアルを和訳しました．まだ未完成ですが，...
ある方は御覧下さい．~
~
http://home.csis.u-tokyo.ac.jp/~maruyama/gstat/gstatmanjp...
~
作り方のエッセンス(vine linuxの場合)~
・/usr/lib/R/share/texmf/Rd.styを?usepackage{Rd.sty}~
で指定する．~
・/usr/lib/R/library/gstat/latex/*.texを和訳・編集して~
つなぎ合わせる．~
・コンパイルしてdviファイルを作り，dvipdfm~
//
-地理学ご出身の丸山先生に、なんですが、３１ページの｢メウ...
-すみません．改訂時に修正します．この川を知らなかったので...
#comment
#comment
*R Commanderのメニューについて [#nd8861b6]
> (2003-08-03 (日) 07:50:20)
~
>R Commanderのメニューって、デフォルトの関数以外でも指定...
-0.90ではetcに*.Rを入れればできると書かれています。もっと...
-日本語化R.dllを入れると、utf-8が読めないので日本語menuの...
-あっ、元もと国際化版のTcl/Tkですから、なんか弄れば出そう...
-0.90の問題ではなかったんですね。R.dllを元に戻してから林...
-Vineのtcltkを引き摺り出して8.0.5のWinのtcltk(SJIS)をこさ...
-今ひとつ分かりませんが、0.90は日本語化R.dllでは動作しな...
-メニュ‐ファイルをRで読み込むんですが、SJISをきちんと読み...
-ここだと判りにくいと思うので詳細？は[[Rと日本語]]にまと...
#comment
*Rのライセンスについて [#m2f627bd]
>(2003-08-03 (日) 03:32:38)
~
>ところで、RはよりゆるやかなLGPLなどのラインセンスには移...
-将来のことは R Foundation メンバーが決めるのでしょうが、...
-他のGPLのソースとか(LGPLのソースも)もRには含まれています...
#comment
*夏時間 [#w1e8068e]
>[[なかま]] (2003-08-01 (金) 11:08:15)
~
>検証しようとする方は、~
export TZ=est5edt # 夏時間
export TZ=est5 # 夏時間無し
>日本にもそういえば戦前、二つの時間があった時期がありまし...
# 北米におけるtimezone は以下参照~
http://cvs.parisc-linux.org/glibc/timezone/northamerica?r...
日本も明治6年より前は太陰暦だし(閏3月とか表現は難しい)、...
商用ソフトでも、西暦9999年まで大丈夫と言いながら、西暦0年...
# そのうち実装する人が現れるのだろうか?~
#comment
*VineLinuxユーザ向けESSのrpmこさえました [#j23a64a8]
>[[なかま]] (2003-07-31 (木) 18:36:50)
~
>http://www.nakama.ne.jp/memo/cran-R/ess/ からどうぞ。~
ところで、VineLinux使ってる人ってどの位いるんでしょうか。~
-素晴らしい!!!早速ダウンロードさせてもらいました。インス...
-おまけに、私のW2K上のMeadowにessを入れた例?も追加してみ...
#comment
*[質問]if(0){}を色分けしてくれるエディタ [#q268515e]
>[[コメント多用するので、、、]] (2003-07-30 (水) 11:07:09)
~
>Rでの複数行コメント~
if(0){~
コメント~
}~
を色分けしてくれるエディタというとどのようなエディタでし...
-これは本来のコメントではなく、正規の構文の特性を生かした...
-裏技ということは承知しています。正規表現とかでうまくいく...
-ヒントになるかどうかわかりませんが、Xemacs では対応括弧...
-Windows 上の editor でよければ。 -- &new{2003-09-11 (木...
-済ません、上の書き込みの続きです。真魚[まな] www.vector...
#comment
*Eclipse プラグイン [#m70202cd]
> (2003-07-30 (水) 07:08:12)
~
>R-helpで流れているように、R用のEclipseプラグインの開発を...
-7/30から今までに何のレスも付いていないところを見ると、お...
#comment
*四捨五入 [#q79b293f]
>[[なかま]] (2003-07-29 (火) 09:42:52)
~
>私は銀行屋でもなく、科学技術計算をする訳でもないので、ro...
いささか都合が悪いです。~
oround <- function(f,n=0)
{
five <- 0.5
signed <- sign(f)
signed[signed == 0] <- 1
r <- f
p <- 10**(n)
r <- abs(r) * p
r <- floor(r + five)
r / 10**(n) * signed
}
>~
と、とりあえず濁して書きましたが、どなたか車輪を既に発明...
-だいぶ前に、Excel の四捨五入は IEEE 方式でないと聞きまし...
-無知で申し訳ないのですが、IEEE 方式って何ですか？ -- SI...
-google で四捨五入をキーワードに検索するとわかります。普...
-[[JIS,ISO式四捨五入]] を参照 -- SIZE(10){2003-07-29 (火...
-よく、消費税の小数点以下を四捨五入と言うのを見ます。Exce...
-[[JIS,ISO式四捨五入]] にシミュレーション結果
-Banker's rounding(五捨五入)とround off(四捨五入)の説明す...
-書き込むべき所が分かりませんでしたので、こちらに書きます...
- oroundにもエラーが出るようです。> oround(0.135,2) [1] 0...
#comment
*不思議な関数 [#q057dd73]
r-help の記事から (2003-07-28 (月) 13:36:01)
> "+<-" <- function(lhs, rhs, value) value - rhs
> x + 1 <- 3
> x
[1] 2
つまり x + 1 == 3 が成り立つように x の値が決まる。面白い...
-1.8.0からは出来なくなりました。(バグでしょやっぱ) -- &n...
#comment
*R の使い方 [#kbca8f39]
> (2003-07-26 (土) 21:05:31)
~
>「うがっ！ここは R の出番だ！」~
さぁ，貴方はどうする？~
(1) RGui 起動！~
(2) コマンドラインから R 起動~
(3) お気に入りのエディタ起動（含・MS Word）で，関数を書き...
(4) お気に入りのエディタから R 起動 (ex. ESS)~
~
普段どのように使っているかも知りたい所です．~
私は，Emacs から ESS ですが，小数派でしょうか？~
#comment
*curve関数 [#q966bc62]
> (2003-07-26 (土) 19:10:28)
~
>curve(sin,-100,100)とxの範囲を広くすると見るも無残な曲線...
これは仕様ですか？~
-デフォールトがこうなるという意味では仕様です。そもそもほ...
-curve のhelp を見たら，n という引数が指定できますね。
> curve(sin,-100,100,n=10000)
とか書けば（出力デバイスの解像度が許せばだけど）ちったあ...
#ref(sin1.png)
-あまり横に長いものをおかれると他の記事が見にくくなります...
-n という引数を学習しました。ありがとうございます。画像小...
-質問者です。Helpにもしっかり載ってますね・・・お蔭様で理...
#comment
*ラグ変数の作り方？ [#h90e6f2b]
> (2003-07-25 (金) 12:11:44)
~
>教えてください。~
Rにおいて、回帰分析をする際に1期前のラグ変数をどうやって...
また、時系列データをRで扱う方法として、CSV形式ファイルをr...
何か間違っているのでしょうか？よろしくご指導ください。~
//移動しました、まとまってから新しいページを作るというこ...
-一期前の変数ですか。
> x <- 1:10
> lag.x <- c(NA, x)[-length(x)-1]
かな？NA を取るなどはその後で -- SIZE(10){2003-07-25 (金...
- ts ライブラリーにはそのものずばり lag という関数があり...
Description:
Computed a lagged version of a time series, shifting...
back by a given number of observations.
Usage:
lag(x, k = 1, ...)
Arguments:
x: A vector or matrix or univariate or multivariate ...
k: The number of lags (in units of observations).
#comment
* matrix(1:15, nc=3, byrow=TRUE) に対して，各行に 11:13 ...
-5x3 行列を 5x6 行列にするという意味ですか。rbind、cbind ...
-cbind(matrix(1:15, nc=3, byrow=TRUE), matrix(rep(11:13,5...
-エレガントではありませんが、、、
cbind(matrix(1:15,nc=3,b=T),11,12,13)
はどうでしょうか？足りない要素は繰り返されるというのがポ...
-x <- matrix(1:15, nc=3, byrow=TRUE) とすると，t(t(x)+11:...
-Tips なのかな -- SIZE(10){2003-07-25 (金) 11:49:25}
-そっかt()か・・・ -- SIZE(10){2003-07-25 (金) 12:10:16}
#comment
*数値ベクトル x(x(i), i = 1,2,...,n）があり，xx(1) = 0,j=...
-たとえば，x <- 1:10; n <- length(x); c(0, cumsum(x))[-(n...
-同じことですが c(0,cumsum(x)[-length(x)])。いずれにして...
#comment
*viバインドのコマンド編集 [#v929bf49]
>[[~emacs]] (2003-07-24 (木) 12:40:23)
~
>Rはreadlineライブラリを使ってるようなのでemacsバインドで...
ご存知の方がいらっしゃったら教えてください。~
~
よろしくお願いします。~
-readlineだから、(M-C-j)でvi、(C-e)でemacs。 -- [[なかま]...
-レスありがとうございます。M-C-jとは（Meta-Control-j）で...
-Metaは何なのでしょうか。。。。EscでもAltでもできませんで...
-[Alt]+[Ctrl]+[j] です。man readline に書いてあります。in...
-.inputrcで実現できました。ありがとうございました。 -- [[...
#comment
* RGui の help() で検索はできる？ [#n0c6287d]
> (2003-07-23 (水) 15:18:08)
~
>MS Windows 阪の RGui で ''help()'' を使うと別ウィンドウ...
CUI 阪のコマンドラインですとページャ (ex. ''less'') の機...
-help("hoge",htmlhelp=T)ではダメでしょうか。どのみち別窓...
-help("hoge",htmlhelp=T) これっていいですね。私はいままで...
-help.start() しておけばよいのですね．ってか何故検索でき...
-emac/essやexyzzyなどhelpの結果を検索できるのではないでし...
#comment
*RSSからHTMLの生成 [#za3c11b6]
>[[橘]] (2003-07-19 (土) 14:57:47)
~
>うちのサイトでこちらのWikiのRSSを利用したページ~
http://www.ikushimo.com/sal/diaryandwiki.html ~
を作成しました.~
サーバに負荷がかかるとやめるかもしれませんがとりあえずご...
このページのソースに書いてあるタグを打ち込んだページを用...
ただし, うちのサーバかこちらのサーバのどちらかが止まって...
//-リンクURLの末尾にチルダが付いていて跳べませんね…. 再掲...
-便利に使っています。--チルダの位置を変えました。
#comment
* [質問] R の因子分析法について -> 希望の手法は現在の R ...
> (2003-07-18 (金) 06:32:48)
~
>SPSSなどではメインである主因子法を、なぜRはサポートしな...
-ある分野では中心的である手法も、統計学全体から見れば、む...
-レスありがとうございます。principal component analysis ...
-とんちんかんなコメント失礼しました。
-主因子法より最尤法が明らかに優れているので，あえて主因子...
-質問者です、懇切丁寧なレスありがとうございます。因子分析...
-「色々な手法が試せるほうが、過去の研究の手法と比較しなが...
-オープンソースなのですから、どうしても必要としている方が...
-biased sampling かも知れません。旧来の方法に不満を持つ人...
#comment
*[質問] データ読み込み時に列や行を選ぶことは出来ますか？ ...
>SASから移行です (2003-07-18 (金) 03:5:42)
~
>大きなデータを扱っているもので、上記のことができるといい...
-データがファイルに平文で書かれていると仮定します。また、...
-SAS ではそうした作業を Proc Data で行いやすいので，そう...
-補足ありがとうございます。そのようなことです、上記のよう...
-UNIXで列を選ぶだけなら、pipeコマンドを使って awk をR内部...
-R に読み込めるデータの限界はどの程度なのでしょうか？ テ...
-Rに読み込めるデータ量は、R起動時に --max-vsize とか --ma...
-help("Memory") and help("memory.profile") を参考に調整し...
-scan() を使うのかなぁ？ -- SIZE(10){2003-07-20 (日) 00:...
-ありがとうございます。?scanに書いてあるようにskipやnline...
scan("ex.data", skip = 1, nlines=1)# only 1 line after t...
のように読み込みたい(途中の)複数行を選択できますね、--SAS...
-?read.table読み返したらskipはread.tableでも使えます。大...
#comment
*[質問] perlからRの関数を呼ぶ方法は？ -> とりあえず解決？...
>[[福重貴雄]] (2003-07-17 (木) 18:48:03)
~
>(j)perlを使っていますが、perlからRの関数を呼べるようにす...
RSperlというものがあるらしいのですが、ネットで検索しても...
具体的には qbeta()を呼び出したいのですが。~
よろしくお願いします。~
-こういう時は CRAN の検索エンジンで R の ML の過去記事か...
-Tips 集 [[PerlスクリプトからRを起動する]] が参考になりま...
-[[Rとインターフェースのあるアプリ]]に"RSperl"のリンク張...
-ありがとうございました。実行時間は急がないのですが、コー...
#comment
*SAS JMP [#jc479d3e]
> (2003-07-17 (木) 17:27:21)
~
> http://www.jmp.com/japan/news_and_event/jmp_press_relea...
[[Rを動かしているプラットフォーム]]が影響ってことは流石...
-そもそもRを使う動機は、Linux上で統計パッケージを使いたか...
#comment
*[質問] それぞれのステップにかかった時間を知る方法は？ ->...
>[[SASから移行です。]] (2003-07-17 (木) 03:08:27)
~
>SASのようにプラグラムを走らせるのにかかった時間を出す関...
-[[Rの関数定義の基本]]を参照して下さい -- SIZE(10){2003-...
#comment
----
終了行:
COLOR(blue){SIZE(20){なんでも掲示板アーカイブ（１）}}~
SIZE(20){COLOR(red){注意}}：ここに保存した記事へのコメン...
元記事日付 2003-07-17 から 2003-11-26 分 ~
#contents
~
**行列データの抜き出し方(時系列) [#ef3afdea]
>[[さかもと]] (2003-11-26 (水) 10:08:12)~
~
ある時系列データセット(1980年から2002年)を~
>ts(data="data set name",start=1980,end=2002,frequency=1)...
として読み込んだ後，ある変数の1980年のデータを取り出した...
>D[1980,] # ← 1980年のデータすべてを表示としたいのですが...
確認できていることは，~
>D[1,] # ← だったら1980年のデータを参照することができ...
年の指定ができるとそうさもしやすいのですが･･･．なにとぞご...
//
-window(D,start=1980,end=1980)ではどうでしょう。 -- [[な...
-アドホックですが、頻繁に年度で参照したければ次のように専...
> library(ts)
> data(sunspot)
> sunspot
Time Series: Start = 1936 End = 1972 Frequency = 1
[1] 40 115 100 80 60 40 23 10 10 25 75 145 130 ...
[20] 10 60 190 180 175 120 50 35 20 10 15 30 60...
> ofs <- function(n) {n+1-attributes(sunspot)$tsp[[1]]} ...
> x[ofs(1936)]
[1] 40
-使いやすいかどうか分かりませんがラベルをつけるという手も...
> names(x)=1936:1972
> x["1936"]
1936
40
-早速のリプライ多謝でございます．早速試させていただきます...
#comment
**tests/p-r-random-tests.R [#x6d4cee4]
>[[なかま]] (2003-11-25 (火) 23:48:51)~
~
どうも乱数部分に問題?があるような気がするんですが目がショ...
追えません。(爆)~
for ((i=0;i<=100;i++))
do
R -q --vanilla < R-1.8.1/tests/p-r-random-tests....
done
とか
for ((i=0;i<=100;i++))
do
Rterm.exe -q --vanilla < R-1.8.1/tests/p-r-rando...
done
~
等とすると、テストがFAILする場合がたまにあります。~
Rの中のMTはsrc/main/RNG.cの中にあって、~
A C-program for MT19937: Real number version([0,1)-inte...
(1999/10/28)
なバージョンです。~
どなたか、ハックしませんか。(^^;~
# 明日から客先なんで。。。~
//
-行く前に。一晩寝てみて、これってFAILEDする事もあるって論...
#comment
**CODAとBOA どちらがおすすめ？ [#we982c22]
> (2003-11-25 (火) 18:39:33)~
~
BUGSのサンプルデータで作成した出力ファイル(*.out,*.ind)を...
(coda_0.5-14.tar.gz)とBOA(boa-1_0_1R.tgz)の両方を試してみ...
どちらも同じような感じだったのですが、それぞれの特徴的な...
うか？どちらがおすすめでしょうか？~
どちらも最新版の日付は2003年1月で開発が止まっているように...
//
#comment
**SOMパッケージの違い [#h23c02cb]
>[[amaki]] (2003-11-25 (火) 06:17:59)~
~
Rには"class"と"som"の2種類のsomパッケージが入っていますが...
どのように違うのか教えてください。~
"class"がコホネンのsom_pakを基にしたもので、~
"som"は別物と考えてよろしいのでしょうか？~
~
今、100000ベクトルほどあるデータをsomパッケージでマップ化...
とても時間がかかってしまい(学習回数、1000,10000で半日)、~
困っています。~
学習効率、プロットの精度の観点から両者の違い等ありました...
どなたか教えて下さい。~
~
デザインの領域のモノで初心者です。~
ご助言、よろしくお願いします。~
//
-少し表現が分かりにくかと思うので訂正します。現在使ってい...
#comment
**cygwinでRはできる？ [#wac2a965]
>[[taku]] (2003-11-24 (月) 20:17:51)~
~
ｃｙｇｗｉｎでＲがやりたいとおもい、~
http://datamining.tama.ac.jp/~yama/R/install.html~
を参考にインストールしようとしていましたが、~
ｍａｋｅで失敗してしまいます。~
make[3]:*** [devices.o] Error1~
make[2]:*** [R] Error2~
make[1]:*** [R] Error1~
make:*** [R] Error1~
cygwinではRはできないのでしょうか？~
~
よろしくお願いします。~
//
-確か、Cygwin は現在サポートしていないのではないですか？W...
-そうなのですか。ありがとうございます。cygwinではできない...
-以下の参考ページを見つけております。~
http://www1.accsnet.ne.jp/~gen/R-rev/cygwin-R.html~
http://www1.accsnet.ne.jp/~gen/R-rev/cygwin-R.pdf~
私自身まだうまくいってませんので、うまくインストールされ...
-CygwinとMingwではいくつかのヘッダーの位置が異なるのと、d...
-いろいろありがとうございます。参考ページのとおりにまずや...
-とりあえず[[これ:http://www1.accsnet.ne.jp/~gen/R-rev/cy...
#comment
**RGUI内でWindows の環境変数の設定できますか？ [#te479ecd]
>[[SJava]] (2003-11-23 (日) 20:16:18)~
~
現在SJava をWindows内で使用しています。~
~
SJava を動作させるには、環境変数"JAVA_HOME"にJREのパス...
~
現在マシーンにはJavaSDKを導入していて、このSDKのパスをJ...
~
SJAVAを使用するときだけ、JAVA_HOMEのパスを変更したいの...
//
-1.バッチでくるむ。2.Sys.putenv("JAVA_HOME"="hoge")。3."....
-なかまさん、ありがとうございます。 -- &new{2003-11-27 (...
#comment
**Rでベイズ統計入門 [#r10cd017]
> (2003-11-21 (金) 03:18:38)~
~
Rを使ってベイズ統計を初歩から学べるようなサイトってある...
//
-ベイズ統計の意味によりますが、NIST にベイズ統計の実際的...
#comment
**多変量解析の演習門題 [#v681cae3]
>[[宮島 のぼる]] (2003-11-20 (木) 14:38:48)~
~
(X,Y){X=xi,Y=yj}でX,Yが独立でない時ρxy＝０となる例は？~
（ただし、i,j=1,2,3,‥‥）~
//
-これはもしかして学校のレポートの課題ではないでしょうね。...
-少々きついコメントを書いてしまいましたが、もし真面目な質...
#comment
**Rgraphviz on R for Win [#l0dcb804]
> (2003-11-19 (水) 15:22:42)~
~
BioConductor のgraphvizを RのWin版で動かされている方、...
~
MinGWでgraphvizをビルドできなくて、現在指をくわえたまま...
//
-cygwinでビルドできませんか -- &new{2003-11-19 (水) 21:5...
-開発版(>=R 2.1.0)が http://www.bioconductor.org/packages...
-[[正規版:http://www.bioconductor.org/packages/bioc/1.7/s...
#comment
**統計数理研究所ワークショップの発表内容のWebでの公開 [#w...
> (2003-11-17 (月) 22:37:40)~
~
統計数理研究所のワークショップですが、平日なので拝聴す...
~
このワークショップでの発表はPDFやPPTで、Webで公表されな...
//
-代表者に話はしておきますが、そうした形で発表者が原稿を用...
#comment
**誕生月でスポーツ能力が決まる？！ [#ua44e33f]
>[[間瀬]] (2003-11-17 (月) 22:25:52)~
~
今日の日経新聞に面白い記事がでていました。今期の日本の J...
> x=c(23,20,21,51,57,43,46,42,46,38,26,19)
> chisq.test(x)
Chi-squared test for given probabilities
data: x
X-squared = 55.3889, df = 11, p-value = 6.578e-08
~
となりやはり強く有意なようです。これで我が家の男三人がス...
//
-二乗では？出生数の分布も一様でないそうですが、それを調整...
-自乗も二乗も特に区別されずに使われているのでは。私はむし...
-おお！わたしも長年の謎(自乗か二乗か)が解けてすっきりしま...
-ワルノリして更に。昭和４５年以降月別出生数は殆んど同一と...
-カイ二乗が約2,830件,カイ自乗が約722件でした。-- google &...
-これほどの量をどのように集計したのですか。まさか一つ一つ...
-私の愛用のウェッブ英和・和英辞典のスペースアルクによると...
-常用漢字ではニをジ(漢語)と読まないので国語的には自乗が正...
-ちなみに、二乗派はどう発音しているのでしょうか？「にじょ...
#comment
**CRAN での自作パッケージの登録方法 [#f9382483]
> (2003-11-16 (日) 16:09:14)~
~
CRAN で自作のRパッケージを登録したい場合、どのような手...
//
-日本人によるパッケージがあればいいなと思っていました。も...
#comment
**[[Rから他言語利用]] の gmrf.c のフルソースコード希望 [#...
> (2003-11-14 (金) 14:44:19)~
~
現在、オープンソースの C++ で作成されたパッケージをRよ...
~
ところで、[[Rから他言語利用]]にある、 gmrf.c のフルソー...
//
-構いませんが、RjpWiki に載せるには少し長過ぎるかと。御希...
#comment
**回帰分析に関する日本語文献 [#a7caada0]
>[[間瀬茂]] (2003-11-11 (火) 20:05:53)~
~
教えて下さい。重回帰分析 and/or 非線形回帰分析に関する良...
//
-間瀬先生，こっちでもお世話になります．え??，非線形回帰に...
-もし Amazon でのキーワード検索結果が完璧だとすると、現在...
-英語で理論の本と意味の本とでは、どんな本が良いのでしょう...
-すみません，上の書き込みで"え??"となったところは，文字化...
-それはさておき，確かに回帰分析を勉強する本はないですね．...
#comment
**x <- c(....) vs x <- read.csv(...,...) [#b25991fe]
>[[まさ]] (2003-11-11 (火) 14:13:28)~
~
prec <- c(3.2, 3.1, 5.1, 4.8, 8.3, 9.8, 8.3, 6.6, 5.1, 3....
prec <- read.csv("tempdata.csv",header=T)として以下のこと...
op <- par(mfrow=c(2,2))~
#ヒストグラム~
hist(prec,main="Histgram")~
を実行すると~
Error in hist.default(prec, main = "Histgram") : ~
`x' must be numeric
と返ってきます。~
listで確認したところ~
prec <- c(...)では~
[[1]]~
[1] 3.2 3.1 5.1 4.8 8.3 9.8 8.3 6.6 5.1 3.8 5...
[16] 7.4 6.2 7.6 7.5 7.0 6.2 6.0 8.8 10.7 8.5~
となり、グラフが作成できました。~
prec <- read.csv("tempdata.csv",header=T)では~
[[1]]~
temp
1 3.2~
2 3.1~
3 5.1~
4 4.8~
5 8.3~
6 9.8~
7 8.3~
8 6.6~
9 5.1~
10 3.8~
11 5.2~
12 5.6~
13 6.5~
14 5.7~
15 5.7~
16 7.4~
17 6.2~
18 7.0~
19 6.7~
20 5.7~
21 6.2~
22 6.0~
23 8.8~
24 10.7~
25 8.5~
となり、グラフが作成できません。~
prec <- c(3.2, 3.1, 5.1, 4.8, 8.3, 9.8, 8.3, 6.6, 5.1, 3....
temp~
3.2~
3.1~
5.1~
4.8~
8.3~
9.8~
8.3~
6.6~
5.1~
3.8~
5.2~
5.6~
6.5~
5.7~
5.7~
7.4~
6.2~
7~
6.7~
5.7~
6.2~
6~
8.8~
10.7~
8.5~
のデータで実現するにはどうしたらよいでしょうか？~
//
-read.csvで付値したらデータフレームになるので，変数名も指...
-長くなって他のトピが見づらくなるので、編集お願いします。...
#comment
**エクセルの四捨五入 [#ie6f07e7]
>[[かずっち]] (2003-11-11 (火) 11:05:43)~
~
エクセルで202.85-201.4=1.45で小数１桁でroundをかけると~
1.5にならなければいけないのに1.4になります何度計算をして...
何か設定があるのでしょうか？~
//
-Rでもそうですね。round(1.35,1)は1.4になるし，round(1.55,...
-R の help(round) によれば, round は 'IEEE standard' (IEE...
-かずっちさんの事例は Excel VBA 関数での現象と思いますが,...
-Excel 固有の疑問はこの Wiki にはふさわしくないかと愚考し...
-なるほど。1.45から1.4までと1.5までは差が同じなので，末尾...
#comment
**RMySQLをwindowsで使いたい [#q2037c66]
>[[sakamoto]] (2003-11-06 (木) 12:59:50)~
~
RMySQLがwindowsで使えるようにした方はいらっしゃいませんか...
//
-RMySQL/inst/{INSTALL.win,README.windows}に従えば良いよう...
-なかまさんありがとうございます。早速試してみたいと思いま...
#comment
**Rで0-1整数計画問題を解く [#x18848a5]
> (2003-11-01 (土) 23:07:53)~
~
[[Rで数理計画]]を参考にさせていただいておりますが、とこ...
//
-パッケージlpSolveのlpは整数計画問題も対応できるはずです...
-ken さん、ご教授ありがとうございます -- &new{2003-11-04...
#comment
**SAS ー＞ R [#ief8272c]
> (2003-10-29 (水) 12:05:04)~
~
R-jp の記事を転載します。R＋SAS両刀使いの方いませんか。~
~
RをWindowsで勉強しています。さてよくSASを使うとTypeII SS...
//
-carパッケージに含まれています -- &new{2003-10-29 (水) 2...
//
-中澤さんの [[Rtips:http://phi.ypu.jp/swtips/R.html]] に...
#comment
**Rgui（非日本語パッチR for Win）のフォント指定をコマンド...
> (2003-10-29 (水) 10:48:04)~
~
日本語パッチを当てていない、Windows用のR（1.8.0）で、"R...
//
-R_HOME/etc/Rconsole -- [[なかま]] &new{2003-10-29 (水) 1...
-なかまさん、早速のレスありがとうございます。ところで、コ...
-もしかして、Rのコマンドラインから変更してダイレクトに反...
#comment
**[質問]名前に利用できる文字 [#c806f207]
>[[_（アンダーバー）が使えると良いな]] (2003-10-29 (水) 1...
~
質問というよりは、疑問なのですが...~
なぜ、名前につかえる文字に _（アンダーバー）がなく~
.（ドット）が含まれているのですか？~
歴史的な理由でもあるのでしょうか？~
//
-付値(代入)の演算子だったからです。 -- [[なかま]] &new{20...
-R 1.8 からは _ は使えなくなったそうですが、ついでに変数...
-変数名に"_"は使えません。essを使えば、"_"は"<-"に入力時...
-1.8.0から`が加わり、QuotedSymbolValueの中はStringValue(c...
-最近の r-help 記事によると数カ月後には名前の一部としてア...
#comment
**Rの列数の限界値について [#k8e82294]
>[[sallyh]] (2003-10-28 (火) 16:09:52)~
~
こんにちは。~
Rの列数の限界値について調べています。~
read.tableでデータを読み込ませようとした時、列数が30000列...
いくらメモリを増やしても読み込むことが出来ませんでした。~
~
そこで、Rの列数の上限を調べているのですが、調べ方が悪いら...
良く分かりません。~
（エクセルは確か256列でしたよね？）~
~
もしかして、ものすごく基礎のことかもしれませんが、~
ご存知の方いらっしゃいましたら、もし良ければ教えて下さい。~
できれば、Windows版とUnix版で違いがあるのか等も知りたいで...
~
よろしくお願いします。~
//
-答えではないのですが，help(read.table)すると，read.table...
-x=matrix(1.0,c(57000000,57000000)) はOK(すごいけど、おそ...
-列数に制限は無いです。read.tableは列名を得るために何度も...
-ありがとうございます。助かりました。参考にします！！ -- ...
#comment
**Rでマルチ・エージェント [#y282da6f]
>[[M.A.S.]] (2003-10-23 (木) 00:46:19)~
~
マルチ・エージェントというと、SWARMが有名ですが、Rから...
~
CRANでは見当たりませんが？~
//
#comment
**RExcel【(D)COM】の使い方に関して [#u98e2623]
>[[激・初心者です]] (2003-10-21 (火) 17:13:18)~
~
R Start をクリックすると、~
『Fatal error:unable to open the base package』⇒ ~
『Could not start RServer』⇒~
『実行時エラー'462'リモートサーバがないか使用できる状態で...
というエラーが出てしまい、何もできません。~
~
原因がお判りになるかた、レスを頂ければ幸いです。~
//
-ちなみに、RExcel.xlaはきちんとアドインされているようです...
#comment
**[質問] Rの分散分析について [#g15514e3]
>[[R初心者]] (2003-10-20 (月) 21:57:10)~
~
Rの初心者です．つまずいていることがあります．~
~
分散分析をする際に，~
「母数モデル・変量モデル・混合モデル」の~
いずれかを指定したいのですが，~
変量モデルを指定する関数（SASでいうVARCOMP）は，あるので...
//
-私自身は使う必要がないので詳細はわかりませんが、こうした...
and Ripley (邦訳あり)に関連する話題がありそうですから見て...
On Wed, 29 Aug 2001, Troels Ring wrote:
> Dear friends.
>
> I know this has been discussed before but as lively as...
> things might have changed. In Venables and Ripley(199)...
> of variance components including robust methods is dis...
> anything like that presently but might have overlooked...
> do it ? Is it obsolete ?
It is not there in R and not planned. (You can check pla...
developer.r-project.org.) It is still in S-PLUS, so not ...
-- &new{2003-10-20 (月) 22:10:45};
-分散分析 randam and mixed については、以下の URL にある...
-おお．レスがこんなについてたとは・・・nlmeに当たってみま...
#comment
**text関数を使って‥ [#pda2309a]
>[[R初心者]] (2003-10-18 (土) 15:29:31)~
~
text関数を使って描いたグラフのある位置に、求めた重要なパ...
具体例をあげますと‥~
text(100,2,parameter)とすると当然画面上には変数parameter...
//
- 例えば次の例が参考になりますか。paste 関数を使って単一...
> x=runif(10)
> plot(1:length(x), x)
> legend(2, 0.5, paste("no of point", 10, " ")) # 最...
> legend(2, 0.7, paste("no of point", length(x), " "))
#ref(test.png, left)
-- &new{2003-10-18 (土) 18:25:19};
#comment
**[解決]table()で総和を出す方法 [#tba12ab6]
>[[初心者2]] (2003-10-18 (土) 05:34:19)~
~
を考えています。~
table(c(1,1,2,2,3))
とかで、一番最後にTotal 5と表示させたいのですが、いい方法...
//
- テーブル化された結果から元のデータ数を計算したいという...
> x<-table(c(1,1,2,2,3))
> sum(x)
[1] 5
-table関数の仕様を変えたいという意味です。 -- [[初心者2]]...
-ttab <- function(x) { c(table(x),sum(table(x))) } として...
-ttab(floor(runif(10000,1,7)))でもちゃんと一番右に10000と...
-ありがとうございます。ttabを使わせてもらいます。 -- &ne...
-NAを数えたい方は、-- &new{2003-10-21 (火) 11:39:01};
ttab <- function(x) { c(x<-table(x,exclude=NULL),sum=sum...
#comment
**lsfit関数について [#x1128534]
>[[R初心者]] (2003-10-17 (金) 16:04:08)~
~
lsfit(x,y)$coefで回帰係数であるaとbがX，Interceptとして得...
//
-すみません。過去の「最初の一歩：必要な出力だけ欲しい時」...
#comment
**R と Ruby [#m1a1835d]
> (2003-10-15 (水) 10:08:18)~
~
RPy (R from Python) http://rpy.sourceforge.net/about.html~
このようなのがあるのですね。~
Ruby にもあるのでしょうか？ 知っていたら教えてください。~
//
-[[Rとインターフェースのあるアプリ]]にリンク先があります...
-R Math library を Ruby からつかうものは実装があります。R...
#comment
**コホネンのSOM_PAKをRから利用できるパッケージ [#p843ef29]
> (2003-10-10 (金) 18:56:10)~
~
コホネンのSOM_PAKですが、これをRで呼び出して利用できる...
~
//
-自己レスです。class パッケージにSOMはありますね。 -- &n...
-som(旧名 GeneSOM) というパッケージもありますね。 -- &n...
#comment
**[解決] 配列に収めたデータを正規化したいと思っているので...
>[[激・初心者です]] (2003-10-09 (木) 11:40:20)~
~
全くどうしたら良いかわかりません。~
列ごとに平均・標準偏差を求めて正規化を行いたいのですが、...
階層的クラスタリングに用いるデータで、その後の処理(Euclid...
VisualBasicで自作したツールでは不便でRを利用しようと考え...
どなたか、気が向いた方、返答頂ければ幸いです。~
//
-sweep()（http://finzi.psych.upenn.edu/R/library/base/htm...
-関数を簡単に定義できるので、関数を自作し、applyすればい...
-早速、ありがとうございます。 -- &new{2003-10-09 (木) 16...
-上記返答を参考に模索してみます(すぐ利用出来るほど賢くな...
-例えば
> x <- matrix(runif(20), c(5,4))
> xx <- scale(x, apply(x,2,mean), apply(x,2,sd))
行について操作したければ、y <- t(x) と転置しておいて、最...
-素晴らしいです、出来ました！ご丁寧にありがとうございます...
#comment
**中澤先生の「Rによる統計解析の基礎」の出版社 [#obb2faac]
> (2003-10-06 (月) 13:23:35)~
~
件名は１０月に発行されるようですが、Amazonではヒットし...
~
どこの出版社から出るのでしょうか？~
//
-株式会社ピアソン・エデュケーションのようです。 -- &new{...
-中澤 港：「Rによる統計解析の基礎」ピアソン・エデュケー...
#comment
**カーネル密度推定のわかりやすい文献 [#ib1a5b12]
> (2003-10-01 (水) 11:22:20)~
~
R言語はいくつかのパッケージで多くの種類のカーネル密度推...
//
//-この御質問に対する答えとしては、 A. Bowman and A. Azza...
//-書き込みの衝突が起きたそうなので、もう一度書きます。
-この御質問に対する答えとしては、 A. Bowman and A. Azzali...
#comment
**[解決]二つのグラフにまたがるタイトル [#b4ecaf4e]
>林 (2003-09-30 (火) 06:04:46)~
~
お世話になります。
#ref(graphs.png)
のようなふたつのグラフにまたがるようにタイトルをつけたい...
//
-layoutは使った事が無いんですがだいたい以下のようにしてま...
par(mfrow=c(1,2),oma=c(0,0,2,0))
deg<-seq(-180,180)
plot(deg,sin(deg/180*pi),type="l",main="sin")
plot(deg,cos(deg/180*pi),type="l",main="cos")
mtext(side=3,line=0,cex=1.4,font=2,outer=TRUE,"sin cos")
-なるほどうまい。ついでですが、グラフィックス実例集で lay...
-おっ、有難うございました。いつもながら、やはり関数がある...
-上手く行かないのはpngだけであればpngを疑うしかないかとwi...
-なぜ X11 デバイスには出力されるのに png デバイスだと複数...
test() # なかまさんの例を X11 デバイスに出力
pp <- recordPlot() # X11画面への作図手順を記録
png("test.png") # png デバイスを開く
replayPlot(pp) # png デバイスに出力
dev.off() # デバイスを閉じる
-- [[間瀬]] &new{2003-10-01 (水) 21:33:13};
-間瀬先生=debianと解釈。libpngのバージョンはいくつなんで...
-ローカルなトラブルですか。r-help を探しても報告は見つか...
-うちのだと、system(paste("ldd ",R.home(),"/modules/R_X11...
#comment
**[解決]リストの名前タグ [#qf3828bd]
>[[dimnamesのように]] (2003-09-21 (日) 11:01:04)~
~
リストにあとから名前タグをつけるには、どうすればいいので...
//
- names 関数を使います(変更にも使えます。一部だけ使うには...
> x=list(1:3,matrix(1:4,c(2,2)))
> x
[[1]]
[1] 1 2 3
[[2]]
[,1] [,2]
[1,] 1 3
[2,] 2 4
> names(x) <- c("vector","matrix")
> x$vector
[1] 1 2 3
> x$matrix
[,1] [,2]
[1,] 1 3
[2,] 2 4
> x[["vector"]]
[1] 1 2 3
> names(x)[1] <- "vector" # 一部指定にはこのようなやり方...
> names(x)[1] <- "" # 名前タグを(見かけ上)取り去る
-x$matrixには、colnames関数やrownames関数も使えました。あ...
#comment
**[解決]関数lagを使いたいのですが [#hc5921db]
>[[lager]] (2003-09-19 (金) 16:02:15)~
~
時系列分析をやっているのですが、~
RはSPlusのうち時系列分析関係の関数がかなり省かれているよ...
（たとえばtsplotなども）~
lagを使いたいのですが、どこかにないものでしょうか。~
//
-すみません、パッケージでtsを入れたらできました。 -- [[la...
#comment
**[解決]メモリの問題 [#eba4994a]
>林 啓一 (2003-09-18 (木) 03:20:51)~
~
メモリが512Mのマシンを再起動後、Rのみを立ち上げて~
7万＊60変数ぐらいのデータフレームをつくり(下のやつです)~
他のソフトで分析するためにwrite.tableで書き出そうとすると~
メモリが足りないといわれてしまします。~
~
メモリを変えるor小さく分けて書き出す以外になにか解決方法~
ありませんでしょうか？~
//
-メッセージにもよりますが既に下の方の、Reached total allo...
-write.tableがメモリを使うのは x <- as.matrix(x) でmatrix...
-memory.limit(size=1024)で17Mの.txtファイルの書き出しに成...
http://www.geocities.co.jp/SiliconValley-SanJose/3220/XP-...
を読むといらないサービスばっかりですね＞M$
#comment
**CRANのakimaパッケージ内のinterp関数を入れたんですが…。 ...
>[[wash-and-go]] (2003-09-17 (水) 18:21:17)~
~
うまくいきません。途中でエラー出てしまいます。~
~
>Error in .Fortran(as.single(x), as.single(y), as.single(...
> as.integer(nx), : ~
> function name must be a string (of length 1)~
~
といった感じです。他にもパッケージが入っておらずエラーが...
でしょうか？どれを入れていいかもわかりません。~
//
-お使いの環境はなにか分かりませんが~
~
R1.7.1のWin版では、~
library(akima)~
test<-interp(rnorm(10),rnorm(10),rnorm(10))~
image(test$x,test$y,test$z)~
~
で、なんら問題は出ませんが。 -- &new{2003-09-17 (水) 1...
-このパッケージは Fortran コードを使っています。もし Linu...
-R から Fortran を呼び出す場合、.Fortran("関数名文字列",...
-R1.7.0のWin版です。Fortranのコンパイラは入っていないと思...
-たしか、ここから落としてきました。http://www.bioconducto...
-R_HOME/library/akima/DESCRIPTION の Buildには何と書いて...
-http://www.bioconductor.org/CRAN/の1.7のディレクトリの中...
-ftp://ftp.u-aizu.ac.jp/pub/lang/R/CRAN/bin/windows/contr...
-インストール時に古いakimaは消されたのでしょうか、akimaの...
-消しましたがうまくいきません。一度、バージョンを新しくし...
-バージョンを1.8.0にしてakimaを入れてみたらうまくいきまし...
#comment
**NP問題 [#a24cd16e]
>[[なかま]] (2003-09-17 (水) 14:05:38)~
~
どこかにNP問題の一覧みたいなのって無いんでしょうか。~
Rで紐解いてるのがあるともっと嬉しいです。~
//
-一覧があるかどうかは知りません。R の base パッケージには...
-素晴らし過ぎます。でも、キーワードが思い付かないと本当に...
-ga は http://www.dst.unive.it/~claudio/R/index.html#ga ...
-CRAN パッケージ gafit は曲線当てはめという特殊目的のよう...
#comment
**[解決] 効率のいいデータの変形(stack) [#l5cbf661]
>林 啓一 (2003-09-16 (火) 13:13:43)~
一人6回観察されている1万2千件のレコードを1万2千*6回という...
> memory.size(30000000)
[1] 829953112
> for(id in 1: nrow(data)){ for(t in c(15,22,33,46,50,51...
newdata<-rbind(newdata,cbind(t(data[id,]),outcome=...
Error in rbind(newdata, cbind(t(data[id, ]), outc..: Un...
memory (this request is for 2725680 bytes, 830035032 byt...
//
-いまひとつイメージがつかめませんが、t(data[id,])はt(data...
-もちろん perl 等を使ってデータファイルを preprocessing ...
-loopは出来るだけ使わない、cbindやrbindの繰り返しでデータ...
S-PLUS for UNIX/Linux Version 6 Programmer's Guide ~
http://www.msi.co.jp/splus/support/download/ ~
のp.509-14に書いてありました。必要なデータセットは
y<-c(15,22,33,46,50,51)
system.time(stackdata<-matrix(rep(as.vector(t(as.matrix(...
outcome<-as.vector(as.matrix(data[,y]))
time<-sort(rep(1:6,12686))
cbind(stackdata,outcome,time,unique=1:12686*6)
で作成できました。お騒がせしました。2行目のn*mの行列をた...
-単純にデータフレームを縦に6個並べると言う意味でしょうか?...
-なかまさんありがとうございます。
stackdata<-rbind(data,data,data,data,data,data)
が一番早そうですね。[[Rでいいコードを書く:http://plaza.um...
//なかまさんレスありがとうございます。あの件はどうなりま...
#comment
**SASの番組 [#mf75b10b]
> (2003-09-15 (月) 20:19:54)~
~
Rと直接関係がありませんが、いま、NTVで有名な福利厚生につ...
//
-番組を見ていないでの意味不明です。SAS(あの統計ソフトです...
-「世界まる見え！テレビ特捜部」（20:00-20:54）の中の「仰...
-R関係と言う事で、Rの日本語化をしている人の福利厚生を紹介...
#comment
**データフレームからの要素の取り出し [#n6f8c8c8]
>[[竹澤]] (2003-09-10 (水) 11:32:56)~
~
竹澤（中央農業総合研究センター・農業情報研究部）です。~
~
以下のＲオブジェクトを実行しました。
function ()
{
xx <- data.frame(x1=c("a1","b1","b1","d1"), x2=c("a2","b...
print("xx=")
print(xx)
yy <- xx[,1]
print("yy=")
print(yy)
zz <- xx[,2]
qq <- c(yy,zz)
print("qq=")
print(qq)
rr <- c(as.vector(yy), as.vector(zz))
print("rr=")
print(rr)
}
結果は以下のものです。
[1] "xx="
x1 x2
1 a1 a2
2 b1 b2
3 b1 b2
4 d1 b2
[1] "yy="
[1] a1 b1 b1 d1
Levels: a1 b1 d1
[1] "qq="
[1] 1 2 2 3 1 2 2 2
[1] "rr="
[1] "a1" "b1" "b1" "d1" "a2" "b2" "b2" "b2"
データフレームからの要素の取り出しには注意が必要なのです...
ってなことは自明だ、という方が多いのでしょうか。~
//
-おっしゃりたいことを分かりにくいので、もう少し説明しても...
-http://buran.u-gakugei.ac.jp/~mori/LEARN/R/R-intro-jp.17...
-なるほど。 -- [[竹澤]] &new{2003-10-03 (金) 08:55:22};
-「6.3 データフレーム」の「数値ベクトルと論理値、因子はそ...
#comment
**2chぽいですが [#j85f7958]
> (2003-09-08 (月) 19:53:26)~
~
http://epio.jpinfo.ne.jp/1057/index.html~
グラフがExcelのようですが、Excelで求めたんですかね。~
//
#comment
**[質問]--vsize --nsize オプションの使い方について… [#n88...
>[[INA]] (2003-09-04 (木) 11:32:47)~
~
初心者ですが、経済分野の「産業連関分析」で、Rを利用させ...
ただ私の環境では、少し大きなデータを読み込ませるとすぐ...
~
どなたか似たような経験をなさった方はおられないでしょう...
~
「産業連関分析」では作業の手始めとして「工業統計調査」...
つまり~
~
品目1 品目2 品目3… 計
産業分類1~
産業分類2~
産業分類3~
：~
計~
~
という具合です。~
~
データの形式は~
~
'産業分類','品目分類','出荷額'←実際は英語です~
'1843','184312',33000~
'3494','349419',1125~
'1931','193111',10979~
：
~
の形で続いており、約59万レコードあります。（サイズは７MB...
~
これを~
~
>data <- read.csv("c:?Hogehoge.csv",header=TRUE,sep=",")~
~
で読み込もうとしているのですが、どうしても~
~
$"Reached total allocation of 126Mb: see help(memory.size...
~
のメッセージが出て、最後まで読み込むことができません。~
~
念のためvsizeオプションを使い、コマンドラインから~
~
>rgui --vsize 36M ~
~
として走らせて見ましたが、やはり読みきれません。~
~
~
ちなみに同じようなデータ、約18万レコードは問題なく読み込...
（その後、>tapply(data$出荷額,list(data$産業分類,data$品...
ご存知の方、どうぞご指導のほどよろしくお願い致します。~
//
-GlobalMemoryStatus から、R_max_memory = min(1024 * Mega,...
-R上からでも、memory.limit()で認識している物理メモリサイ...
-ご指導有り難うございました。雷の影響かチップがおかしくな...
-どういうデータなのかにもよると思いますが、ちょうど126M以...
-PCをパワーアップして1Gまでメモリを積んだところ、問題なく...
#comment
**最初の一歩：必要な出力だけ欲しい時 [#r26634af]
>[[角煮]] (2003-09-02 (火) 10:38:14)~
~
初心者の私が、必要な出力だけ欲しい時、悩んでしまいました...
~> x<-c(1,3,5) ~
~> y<-c(3,5,4) ~
~> wout<-lsfit(x,y)$coef~
~> wout[1]~
Intercept
3.25
~> wout［［1］］~
[1] 3.25
~> cat(wout[1])~
3.25
これで、やっと数値3.25だけ得られました。（ここまで来るの...
//
-lsfit(x,y)$coef["Intercept"] -- &new{2003-09-02 (火) 11...
-私もてこずりました。リスト型式は慣れると便利です。 ~
Betaの書き出し http://plaza.umin.ac.jp/~epi/cgi-bin/fswi...
にも書きましたが、SASやSTATAでは簡単に書き出せない統計量...
(違ったらすみません、)書き出しやすいのがR(S)のいいところ...
１．str() でリストの中身 (階層構造も) が参照できる~
２．リストは$名前で参照できるとともに、
[[順番]][[その階層での順番]]
でも参照できるということがポイントだと思います。私なら
str(lsfit(x,y))
lsfit(x,y)[[1]][[1]]
でその後の解析に使います。cat()までして数字したいという状...
#comment
**樹形モデル [#sfa9a1e8]
> (2003-08-30 (土) 11:56:24)~
~
S-plusのtree関数に相当するものは、Rのパッケージに存在する...
アドバイスよろしくお願いします。~
//
-http://cran.r-project.org/src/contrib/PACKAGES.html#tree...
#comment
**software patent in Europe [#i6b09b58]
> (2003-08-27 (水) 19:27:59)~
~
9月１日、"software patent "の導入について欧州議会で投票が...
//
-http://swpat.ffii.org/index.ja.html の事でしょうか? -- [...
#comment
**Splusでのinterp()の関数はRにはありますか？ [#k079aa19]
>[[wash-and-go]] (2003-08-25 (月) 16:15:24)~
~
影響関数を描くのに必要な関数なんですが、~
Rにはinterp()という関数はあるでしょうか？~
//
-CRANのakimaパッケージ内のinterp関数（不規則分布のポイン...
#comment
**gstatの日本語マニュアル [#m9966bc7]
>[[丸山]] (2003-08-22 (金) 13:47:42)~
~
CRANに多変量地球統計学のパッケージとして提供されているgst...
レフェレンスマニュアルを和訳しました．まだ未完成ですが，...
ある方は御覧下さい．~
~
http://home.csis.u-tokyo.ac.jp/~maruyama/gstat/gstatmanjp...
~
作り方のエッセンス(vine linuxの場合)~
・/usr/lib/R/share/texmf/Rd.styを?usepackage{Rd.sty}~
で指定する．~
・/usr/lib/R/library/gstat/latex/*.texを和訳・編集して~
つなぎ合わせる．~
・コンパイルしてdviファイルを作り，dvipdfm~
//
-地理学ご出身の丸山先生に、なんですが、３１ページの｢メウ...
-すみません．改訂時に修正します．この川を知らなかったので...
#comment
#comment
*R Commanderのメニューについて [#nd8861b6]
> (2003-08-03 (日) 07:50:20)
~
>R Commanderのメニューって、デフォルトの関数以外でも指定...
-0.90ではetcに*.Rを入れればできると書かれています。もっと...
-日本語化R.dllを入れると、utf-8が読めないので日本語menuの...
-あっ、元もと国際化版のTcl/Tkですから、なんか弄れば出そう...
-0.90の問題ではなかったんですね。R.dllを元に戻してから林...
-Vineのtcltkを引き摺り出して8.0.5のWinのtcltk(SJIS)をこさ...
-今ひとつ分かりませんが、0.90は日本語化R.dllでは動作しな...
-メニュ‐ファイルをRで読み込むんですが、SJISをきちんと読み...
-ここだと判りにくいと思うので詳細？は[[Rと日本語]]にまと...
#comment
*Rのライセンスについて [#m2f627bd]
>(2003-08-03 (日) 03:32:38)
~
>ところで、RはよりゆるやかなLGPLなどのラインセンスには移...
-将来のことは R Foundation メンバーが決めるのでしょうが、...
-他のGPLのソースとか(LGPLのソースも)もRには含まれています...
#comment
*夏時間 [#w1e8068e]
>[[なかま]] (2003-08-01 (金) 11:08:15)
~
>検証しようとする方は、~
export TZ=est5edt # 夏時間
export TZ=est5 # 夏時間無し
>日本にもそういえば戦前、二つの時間があった時期がありまし...
# 北米におけるtimezone は以下参照~
http://cvs.parisc-linux.org/glibc/timezone/northamerica?r...
日本も明治6年より前は太陰暦だし(閏3月とか表現は難しい)、...
商用ソフトでも、西暦9999年まで大丈夫と言いながら、西暦0年...
# そのうち実装する人が現れるのだろうか?~
#comment
*VineLinuxユーザ向けESSのrpmこさえました [#j23a64a8]
>[[なかま]] (2003-07-31 (木) 18:36:50)
~
>http://www.nakama.ne.jp/memo/cran-R/ess/ からどうぞ。~
ところで、VineLinux使ってる人ってどの位いるんでしょうか。~
-素晴らしい!!!早速ダウンロードさせてもらいました。インス...
-おまけに、私のW2K上のMeadowにessを入れた例?も追加してみ...
#comment
*[質問]if(0){}を色分けしてくれるエディタ [#q268515e]
>[[コメント多用するので、、、]] (2003-07-30 (水) 11:07:09)
~
>Rでの複数行コメント~
if(0){~
コメント~
}~
を色分けしてくれるエディタというとどのようなエディタでし...
-これは本来のコメントではなく、正規の構文の特性を生かした...
-裏技ということは承知しています。正規表現とかでうまくいく...
-ヒントになるかどうかわかりませんが、Xemacs では対応括弧...
-Windows 上の editor でよければ。 -- &new{2003-09-11 (木...
-済ません、上の書き込みの続きです。真魚[まな] www.vector...
#comment
*Eclipse プラグイン [#m70202cd]
> (2003-07-30 (水) 07:08:12)
~
>R-helpで流れているように、R用のEclipseプラグインの開発を...
-7/30から今までに何のレスも付いていないところを見ると、お...
#comment
*四捨五入 [#q79b293f]
>[[なかま]] (2003-07-29 (火) 09:42:52)
~
>私は銀行屋でもなく、科学技術計算をする訳でもないので、ro...
いささか都合が悪いです。~
oround <- function(f,n=0)
{
five <- 0.5
signed <- sign(f)
signed[signed == 0] <- 1
r <- f
p <- 10**(n)
r <- abs(r) * p
r <- floor(r + five)
r / 10**(n) * signed
}
>~
と、とりあえず濁して書きましたが、どなたか車輪を既に発明...
-だいぶ前に、Excel の四捨五入は IEEE 方式でないと聞きまし...
-無知で申し訳ないのですが、IEEE 方式って何ですか？ -- SI...
-google で四捨五入をキーワードに検索するとわかります。普...
-[[JIS,ISO式四捨五入]] を参照 -- SIZE(10){2003-07-29 (火...
-よく、消費税の小数点以下を四捨五入と言うのを見ます。Exce...
-[[JIS,ISO式四捨五入]] にシミュレーション結果
-Banker's rounding(五捨五入)とround off(四捨五入)の説明す...
-書き込むべき所が分かりませんでしたので、こちらに書きます...
- oroundにもエラーが出るようです。> oround(0.135,2) [1] 0...
#comment
*不思議な関数 [#q057dd73]
r-help の記事から (2003-07-28 (月) 13:36:01)
> "+<-" <- function(lhs, rhs, value) value - rhs
> x + 1 <- 3
> x
[1] 2
つまり x + 1 == 3 が成り立つように x の値が決まる。面白い...
-1.8.0からは出来なくなりました。(バグでしょやっぱ) -- &n...
#comment
*R の使い方 [#kbca8f39]
> (2003-07-26 (土) 21:05:31)
~
>「うがっ！ここは R の出番だ！」~
さぁ，貴方はどうする？~
(1) RGui 起動！~
(2) コマンドラインから R 起動~
(3) お気に入りのエディタ起動（含・MS Word）で，関数を書き...
(4) お気に入りのエディタから R 起動 (ex. ESS)~
~
普段どのように使っているかも知りたい所です．~
私は，Emacs から ESS ですが，小数派でしょうか？~
#comment
*curve関数 [#q966bc62]
> (2003-07-26 (土) 19:10:28)
~
>curve(sin,-100,100)とxの範囲を広くすると見るも無残な曲線...
これは仕様ですか？~
-デフォールトがこうなるという意味では仕様です。そもそもほ...
-curve のhelp を見たら，n という引数が指定できますね。
> curve(sin,-100,100,n=10000)
とか書けば（出力デバイスの解像度が許せばだけど）ちったあ...
#ref(sin1.png)
-あまり横に長いものをおかれると他の記事が見にくくなります...
-n という引数を学習しました。ありがとうございます。画像小...
-質問者です。Helpにもしっかり載ってますね・・・お蔭様で理...
#comment
*ラグ変数の作り方？ [#h90e6f2b]
> (2003-07-25 (金) 12:11:44)
~
>教えてください。~
Rにおいて、回帰分析をする際に1期前のラグ変数をどうやって...
また、時系列データをRで扱う方法として、CSV形式ファイルをr...
何か間違っているのでしょうか？よろしくご指導ください。~
//移動しました、まとまってから新しいページを作るというこ...
-一期前の変数ですか。
> x <- 1:10
> lag.x <- c(NA, x)[-length(x)-1]
かな？NA を取るなどはその後で -- SIZE(10){2003-07-25 (金...
- ts ライブラリーにはそのものずばり lag という関数があり...
Description:
Computed a lagged version of a time series, shifting...
back by a given number of observations.
Usage:
lag(x, k = 1, ...)
Arguments:
x: A vector or matrix or univariate or multivariate ...
k: The number of lags (in units of observations).
#comment
* matrix(1:15, nc=3, byrow=TRUE) に対して，各行に 11:13 ...
-5x3 行列を 5x6 行列にするという意味ですか。rbind、cbind ...
-cbind(matrix(1:15, nc=3, byrow=TRUE), matrix(rep(11:13,5...
-エレガントではありませんが、、、
cbind(matrix(1:15,nc=3,b=T),11,12,13)
はどうでしょうか？足りない要素は繰り返されるというのがポ...
-x <- matrix(1:15, nc=3, byrow=TRUE) とすると，t(t(x)+11:...
-Tips なのかな -- SIZE(10){2003-07-25 (金) 11:49:25}
-そっかt()か・・・ -- SIZE(10){2003-07-25 (金) 12:10:16}
#comment
*数値ベクトル x(x(i), i = 1,2,...,n）があり，xx(1) = 0,j=...
-たとえば，x <- 1:10; n <- length(x); c(0, cumsum(x))[-(n...
-同じことですが c(0,cumsum(x)[-length(x)])。いずれにして...
#comment
*viバインドのコマンド編集 [#v929bf49]
>[[~emacs]] (2003-07-24 (木) 12:40:23)
~
>Rはreadlineライブラリを使ってるようなのでemacsバインドで...
ご存知の方がいらっしゃったら教えてください。~
~
よろしくお願いします。~
-readlineだから、(M-C-j)でvi、(C-e)でemacs。 -- [[なかま]...
-レスありがとうございます。M-C-jとは（Meta-Control-j）で...
-Metaは何なのでしょうか。。。。EscでもAltでもできませんで...
-[Alt]+[Ctrl]+[j] です。man readline に書いてあります。in...
-.inputrcで実現できました。ありがとうございました。 -- [[...
#comment
* RGui の help() で検索はできる？ [#n0c6287d]
> (2003-07-23 (水) 15:18:08)
~
>MS Windows 阪の RGui で ''help()'' を使うと別ウィンドウ...
CUI 阪のコマンドラインですとページャ (ex. ''less'') の機...
-help("hoge",htmlhelp=T)ではダメでしょうか。どのみち別窓...
-help("hoge",htmlhelp=T) これっていいですね。私はいままで...
-help.start() しておけばよいのですね．ってか何故検索でき...
-emac/essやexyzzyなどhelpの結果を検索できるのではないでし...
#comment
*RSSからHTMLの生成 [#za3c11b6]
>[[橘]] (2003-07-19 (土) 14:57:47)
~
>うちのサイトでこちらのWikiのRSSを利用したページ~
http://www.ikushimo.com/sal/diaryandwiki.html ~
を作成しました.~
サーバに負荷がかかるとやめるかもしれませんがとりあえずご...
このページのソースに書いてあるタグを打ち込んだページを用...
ただし, うちのサーバかこちらのサーバのどちらかが止まって...
//-リンクURLの末尾にチルダが付いていて跳べませんね…. 再掲...
-便利に使っています。--チルダの位置を変えました。
#comment
* [質問] R の因子分析法について -> 希望の手法は現在の R ...
> (2003-07-18 (金) 06:32:48)
~
>SPSSなどではメインである主因子法を、なぜRはサポートしな...
-ある分野では中心的である手法も、統計学全体から見れば、む...
-レスありがとうございます。principal component analysis ...
-とんちんかんなコメント失礼しました。
-主因子法より最尤法が明らかに優れているので，あえて主因子...
-質問者です、懇切丁寧なレスありがとうございます。因子分析...
-「色々な手法が試せるほうが、過去の研究の手法と比較しなが...
-オープンソースなのですから、どうしても必要としている方が...
-biased sampling かも知れません。旧来の方法に不満を持つ人...
#comment
*[質問] データ読み込み時に列や行を選ぶことは出来ますか？ ...
>SASから移行です (2003-07-18 (金) 03:5:42)
~
>大きなデータを扱っているもので、上記のことができるといい...
-データがファイルに平文で書かれていると仮定します。また、...
-SAS ではそうした作業を Proc Data で行いやすいので，そう...
-補足ありがとうございます。そのようなことです、上記のよう...
-UNIXで列を選ぶだけなら、pipeコマンドを使って awk をR内部...
-R に読み込めるデータの限界はどの程度なのでしょうか？ テ...
-Rに読み込めるデータ量は、R起動時に --max-vsize とか --ma...
-help("Memory") and help("memory.profile") を参考に調整し...
-scan() を使うのかなぁ？ -- SIZE(10){2003-07-20 (日) 00:...
-ありがとうございます。?scanに書いてあるようにskipやnline...
scan("ex.data", skip = 1, nlines=1)# only 1 line after t...
のように読み込みたい(途中の)複数行を選択できますね、--SAS...
-?read.table読み返したらskipはread.tableでも使えます。大...
#comment
*[質問] perlからRの関数を呼ぶ方法は？ -> とりあえず解決？...
>[[福重貴雄]] (2003-07-17 (木) 18:48:03)
~
>(j)perlを使っていますが、perlからRの関数を呼べるようにす...
RSperlというものがあるらしいのですが、ネットで検索しても...
具体的には qbeta()を呼び出したいのですが。~
よろしくお願いします。~
-こういう時は CRAN の検索エンジンで R の ML の過去記事か...
-Tips 集 [[PerlスクリプトからRを起動する]] が参考になりま...
-[[Rとインターフェースのあるアプリ]]に"RSperl"のリンク張...
-ありがとうございました。実行時間は急がないのですが、コー...
#comment
*SAS JMP [#jc479d3e]
> (2003-07-17 (木) 17:27:21)
~
> http://www.jmp.com/japan/news_and_event/jmp_press_relea...
[[Rを動かしているプラットフォーム]]が影響ってことは流石...
-そもそもRを使う動機は、Linux上で統計パッケージを使いたか...
#comment
*[質問] それぞれのステップにかかった時間を知る方法は？ ->...
>[[SASから移行です。]] (2003-07-17 (木) 03:08:27)
~
>SASのようにプラグラムを走らせるのにかかった時間を出す関...
-[[Rの関数定義の基本]]を参照して下さい -- SIZE(10){2003-...
#comment
----
ページ名:
![[RjpWiki]](image/../Rlogo.png "[RjpWiki]")
