【R言語】統計解析フリーソフトＲ 第４章【GNU R】
■ このスレッドは過去ログ倉庫に格納されています
0001１３２人目の素数さん
2011/01/09(日) 17:31:49統計計算で重宝するデータ型や、複数要素を処理する演算や関数、
解析結果を表示するグラフィックなど、多彩な機能を提供します。
●関連サイト
The R Project
http://www.r-project.org/
RjpWiki
http://www.okada.jp.org/RWiki/
リンク集
http://www.okada.jp.org/RWiki/?%A5%EA%A5%F3%A5%AF%BD%B8
0093１３２人目の素数さん
2011/06/24(金) 12:32:27.69なんでEが平成なんだろう
0094１３２人目の素数さん
2011/06/24(金) 13:01:48.18マニュアルには、
%E Modifier: use alternative format
と書いてあるよ。つまり日本語ローケルだと、日本語化した形式にするという意味。
>>92を書き込んだ後、MacOSXとWindowsで試してみたが、
どちらも和暦にならなかった。GNUだけかも。ヘルプには環境依存って書いてあるしね。
0095あんでぃはストーカー ◆AdkZFxa49I
2011/06/24(金) 16:34:26.88あんでぃ
0096R入門
2011/06/26(日) 00:02:26.05Rについても勉強したいと思っています。
Rのサイトをみると、結構、「使い方」的なサイトが多く、
汎用言語の視点から書かれた、サイト、本が少ない(？)ような気がします。
たとえば、ほかの言語と比較しながら書かれたような
本、サイトなどあれば、教えていただきたいのですが。
ありますでしょうか。
0097１３２人目の素数さん
2011/06/26(日) 00:19:39.460098１３２人目の素数さん
2011/06/26(日) 01:03:33.08Rの基礎とプログラミング技法
http://www.amazon.co.jp/dp/4431712186
http://d.hatena.ne.jp/sesejun/20090331/p2
Rと他の言語を比較した本は知りませんが、
汎用言語としてのR言語を勉強するならこの本がお勧めです。
0099１３２人目の素数さん
2011/06/27(月) 09:37:34.48そういえば、Rユーザ会で、韓国から来た演者が汎用言語の視点から、Rの問題点などを
発表していたな。興味がなかったので内容は覚えていない。
汎用言語の視点で書かれた書籍はないと思う。
しいてあげるなら、>>98が薦めるリゲスさんの本しかない。
010085
2011/07/04(月) 10:41:00.37xtsをインストールしたんですが
データフレームをxtsにコンバートできません。
> x.xts<-as.xts(x)
以下にエラー as.POSIXlt.character(x, tz, ...) :
character string is not in a standard unambiguous format
Zooオブジェクトに経由させることもできません
> x.xts<-as.xts(read.zoo(x))
以下にエラー xts(coredata(x), order.by = order.by, frequency = frequency, :
order.by requires an appropriate time-based object
データフレームの中身は
日付 株価
20040101 1500
20040102 1450
といった具合に２列のデータなのですが？
0101１３２人目の素数さん
2011/07/04(月) 11:05:02.37ここに書くのではなくて、
それを書いた人に凸してみたら？
> library(xts)
要求されたパッケージ zoo をロード中です
> dat <- data.frame(日付=c(20040101,20040102),株価=c(1500,1450))
> dat
日付 株価
1 20040101 1500
2 20040102 1450
> as.xts(dat)
以下にエラー as.POSIXlt.character(x, tz, ...) :
character string is not in a standard unambiguous format
> rownames(dat) <- c("2004-01-01","2004-01-02")
> as.xts(dat)
日付 株価
2004-01-01 20040101 1500
2004-01-02 20040102 1450
まぁ、参考になれば。
0102１３２人目の素数さん
2011/07/14(木) 14:34:30.11相補誤差関数 erfc を qnorm 関数から定義すると、 qnorm が実数に対してしか定義されていないためエラーになってしまいます。
0103１３２人目の素数さん
2011/07/15(金) 02:10:37.99＞ 汎用言語の視点から書かれた、サイト、本が少ない(？)ような気がします。
Rの本を書いているのは、大学の先生が多いからね。
そして、大学の先生ってのはコンピュータのことや、
プログラミングのことをわかってないから。
Rの使い方だけの説明になってしまうのは
しかたがないと思いますよ。
0104１３２人目の素数さん
2011/07/15(金) 11:23:45.87○ 情報科学が専門ではなく疫学やら統計屋が仕事の教員
0105１３２人目の素数さん
2011/07/16(土) 00:45:00.68○ 大学の先生
Ｘ 情報科学が専門ではなく疫学やら統計屋が仕事の教員
特に情報科学の先生や院生などは、知ったかの馬鹿ばっかりです。
国立情報学研究所を見ればよくわかるでしょ？
税金を使ってＰＣで遊んでるだけ。
全員がコネ採用。
0106１３２人目の素数さん
2011/07/16(土) 01:04:55.28唯のお飾りか？
0107１３２人目の素数さん
2011/07/16(土) 16:16:24.06税金を無駄遣いしながら、馬鹿同士楽しく
PCで遊んでいる組織ですね。
どいつもこいつも馬鹿ばっかり。
コネがないと絶対に採用されません。
公募する前から採用者は決まってます。
0108１３２人目の素数さん
2011/07/16(土) 16:49:55.560109１３２人目の素数さん
2011/07/16(土) 20:02:28.47統計数理研究所
この二つは廃止で結構です。税金の無駄です。
0110１３２人目の素数さん
2011/07/16(土) 20:09:53.71バカばっかりやなw あれで数学やっているつもりのもいるからw
頓珍漢なことばかりを言って来るw
0111猫は嘘吐き ◆MuKUnGPXAY
2011/07/16(土) 20:20:25.680112１３２人目の素数さん
2011/07/16(土) 20:58:33.45しかし、単なる就職先のひとつとなっており、
学歴とコネのみで判断されるから。
実力で研究者になれるところって、少ないと思う。
0113１３２人目の素数さん
2011/07/16(土) 21:06:58.590114１３２人目の素数さん
2011/07/16(土) 21:19:33.030115猫は無職 ◆MuKUnGPXAY
2011/07/17(日) 03:35:47.40猫
0116１３２人目の素数さん
2011/07/17(日) 13:55:02.20最近の調査では勤労者の平均年収が
３７０万円ぐらいらしいんですよ。
ただし最頻値は２５０万円ぐらいですから、
一部の高給取りが平均を上げていると思います。
また、公務員の平均は５００万円ぐらいらしいですが、
公務員だけの平均年収は公表されていません。
研究と称してPCで遊んでいるだけですから、
年収１５０万円でも十分だと思います。
0117猫の夢と痴漢の埋葬 ◆MuKUnGPXAY
2011/07/17(日) 14:08:35.01『研究と称してＰＣで遊んでる人』だけをそういう待遇にしたら良いと
思います。但しソレをどうやって判定スルのかは現実には難しく、なの
で『虚偽院生と優秀院生の分離』と同様の困難が伴い兼ねませんね。
では質問ですが、虚偽院生に対してはどういう扱いが宜しいのでしょう
かね？ やはり追放ですかね？
猫
0118１３２人目の素数さん
2011/07/17(日) 14:12:09.83統計数理研究所では女の子の飲み物調査という論文もあったぞw
0119１３２人目の素数さん
2011/07/17(日) 14:24:35.05虚偽院生とは、論文を書かない院生のことですか？
全く論文を書く気のない院生は追放でしかたないと思います。
少なくとも論文を何本か書いていて、研究費も助成金も１円も
受け取っておらず、アルバイトをして授業料を納めている奴は、
他の人の研究の邪魔をしないかぎりは、追放でなくて放置でいいと思います。
つまり、まあおれのことですけどね。
追放すべきは高い給料、研究費や助成金をたんまりもらっているのに
結果を出さない馬鹿な奴です。
これは院生だけではなく、助手も助教、準教授、教授も含まれます。
それらは能力と業績だけで判定されるべきだと思います。
その業績というのをどう判定するのかが難しいわけですけどね。
0120猫の夢と痴漢の埋葬 ◆MuKUnGPXAY
2011/07/17(日) 14:29:59.14全くその通りだと思います。私が言う虚偽院生とは：
★★★『基礎体力も無ければ動機も無く、唯大学院に入学してしまったゾンビ』★★★
を指し示す言葉です。ですから恐らくは貴方の事ではないでしょう。放置
されていても自分から何かが出来る人は決してゾンビではアリマセンから。
猫
0121１３２人目の素数さん
2011/07/17(日) 14:41:00.64つまり>>119の発言にまったく同意しているという解釈でよいですか？
0123１３２人目の素数さん
2011/07/17(日) 14:51:47.93眼醒めたらゾンビみたいだったが、観てたら生気がもどってきた
この辺のクオリティの高さは邦楽にはない、流石だとおもう。
0124１３２人目の素数さん
2011/07/17(日) 23:18:34.72周囲と意見交換が出来るということは、非常に大きな能力であり財産だと思うよ
そして研究室への貢献度という点で言えば
論文は書かないけどPCで遊べる程度の実力はあって、周囲と意見交換が出来る奴 ＞ 一人黙々と論文を書く奴
だと自分は思う
0125１３２人目の素数さん
2011/07/17(日) 23:41:22.12月１０万でボーナス（一ヶ月)も出て国家的身分は安泰
やりがいなんかまったくないと思うかもしれないが、言うことなしの職場じゃないか
0126１３２人目の素数さん
2011/07/17(日) 23:45:41.25それこそ、そこのなんとか研究所に問い合わせをして詳細資料（統計データ)請求すればいいんじゃないの？
ただ気になったんだけど、rlang使ってるなら「加重平均」ってなんのことかはしってるよね。
0127１３２人目の素数さん
2011/07/18(月) 02:46:14.960128猫は爺惨状の蚊軍 ◆MuKUnGPXAY
2011/07/18(月) 03:33:15.34その考え方は決して認められない。大学は馬鹿がママゴトして遊ぶ場所
では断じてない。当然に：
★★★『黙々と論文を書く奴＞論文は書かない奴』★★★
であり、PCで遊ぶかどうかや、周囲と意見交換が出来るかどうかなんて、
そんな事はどうでもヨロシ。
論文を書くかどうかだけが問題である。勘違いも甚だしい。ケシカラン。
猫
0129１３２人目の素数さん
2011/07/18(月) 03:34:40.12まずお前が実績を出してから批判しような
0130１３２人目の素数さん
2011/07/18(月) 04:44:10.88「左右非対称分布の平均値は分布を的確に示せない」
のは統計の基本常識です。
集計方法と指標の定義: さまざまな平均値/中央値/最頻値
http://web-tan.forum.impressrd.jp/e/2008/07/11/3428
0131１３２人目の素数さん
2011/07/18(月) 05:55:50.190132１３２人目の素数さん
2011/07/18(月) 06:00:29.170133１３２人目の素数さん
2011/07/18(月) 07:56:35.430134１３２人目の素数さん
2011/07/18(月) 07:59:15.180135１３２人目の素数さん
2011/07/18(月) 10:43:44.13どこの大学でも統計の先生の論文は下らんものばかり
0136１３２人目の素数さん
2011/07/18(月) 10:54:16.080137１３２人目の素数さん
2011/07/18(月) 11:20:01.810138１３２人目の素数さん
2011/07/18(月) 11:22:33.75女子高校生の飲み物調査
女子中学生の飲み物調査
女子小学生の飲み物調査
女子幼稚園生の飲み物調査
２０代OLの飲み物調査
３０代のOLの飲み物調査
以下、女と男にかえても論文が書けます 統計数理研究所の教授がこうした論文を書いています
0139１３２人目の素数さん
2011/07/18(月) 11:24:51.54同上 ２ ２部上場企業の場合
同上 ３ ヘラクレス上場企業の場合
同上 ４ 町役場の場合
以下 延々と書けます。
0140１３２人目の素数さん
2011/07/18(月) 11:30:15.870141１３２人目の素数さん
2011/07/18(月) 15:21:16.300142１３２人目の素数さん
2011/07/18(月) 18:25:28.32等高線を書きたいのですが、情報量が少ないため、contour関数では書けません。現在与えられているのは、20個の点に対しての座標(x,y,z)です。2次元上の布置まではできているのですが、等高線を用いてz座標を表現することができません。
どの関数を扱えばいいかなど、わかる方がいらっしゃいましたら教えてください。
0143１３２人目の素数さん
2011/07/18(月) 18:59:47.110144１３２人目の素数さん
2011/07/18(月) 22:51:46.35つか、授業聞いてなかったの？
全部説明したはずだけどね。
geoRで作成しなさいって言ったでしょ？
単位あげないよ！
0145１３２人目の素数さん
2011/07/31(日) 19:04:21.64Rコマンダーでcsv経由。
RGUIはMSGhosicに設定。
コマンドでフォント指定以外に
対策あれば教えてください。
0146１３２人目の素数さん
2011/08/01(月) 00:30:07.52Rというよりそのcsvファイルがくせもの
解決策は2つ。
1. csvの文字コードを文字化けしないものに変更してからRに読み込む
2. Rの中でiconv()で文字コードを変換する
好きな方を選べ
0147１３２人目の素数さん
2011/08/01(月) 05:59:12.09了解です。
クリップボードを経由してみたら、ますます
酷くなりました。
特に列の変数名がボロボロです。
これはなにゆえでしょうか？
0148１３２人目の素数さん
2011/08/01(月) 11:35:11.11全く別のアプローチを提示してみます。
> dat <- read.csv(file.choose())
> dat
bake bakeda
1 32 re
2 242 ga
ここで列名が文字化けしていたとします。
本当に化けている状態を示せないので仮に"bake"とか"bakeda"とかにしています。
> names(dat) <- c("日本語","大丈夫")
> dat
日本語 大丈夫
1 32 re
2 242 ga
というようにnames()で、列名を自由に再設定できます。
列名が、文字化けしていて、なおかつ文字コードを変更する能力がないなら、
上記のような方法もあります。
0149１３２人目の素数さん
2011/08/08(月) 14:36:57.39東大,
弁護士,
Re,
0150１３２人目の素数さん
2011/08/08(月) 14:38:15.29TS10,SBR,VFK10,TKK,VF1,LCCR,SINX,
VF1M4,VF1L2,VF1H2,EMPC,MPE,
4231,4213,3331,3313,145,53A3,6236,
EMPCB,EMJ,LP,CJ,F4,LC,DNA,RNA,
SINT,JEL,23458,2348,DBT,GMO,
AB,APLWJKSJ,PES,WE,CA,RR,ASL,
EPH,ITU,261036,CBS6,1358,G1,AS3,M5,
0151１３２人目の素数さん
2011/08/11(木) 22:25:16.5220個の変数があり、各変数は0 から 0.2 まで 0.05 ずつの値をとることが可能
また、この変数の総和は 1 となる。
この時の組み合わせをすべて求めたいのですが、
どうすれば良いでしょうか?
よろしくお願いいたします
0152１３２人目の素数さん
2011/08/12(金) 01:53:48.61> seq(from=0,to=0.2,by=0.05)
[1] 0.00 0.05 0.10 0.15 0.20
> factorial(5+20-1)/factorial(20)/factorial(5-1)
[1] 10626
この程度の組み合わせなら全て計算して合計が1のものを選んだ方が早いかも。
0153152
2011/08/12(金) 02:20:56.85> n <- seq(from=0,to=0.2,by=0.05)
> m <- expand.grid(n,n,n,n,n,n)
> m.s <- apply(m,1,sum)
> a <- m[m.s==1,]
> head(a)
Var1 Var2 Var3 Var4 Var5 Var6
3125 0.20 0.20 0.20 0.20 0.20 0.00
5625 0.20 0.20 0.20 0.20 0.15 0.05
6125 0.20 0.20 0.20 0.15 0.20 0.05
6225 0.20 0.20 0.15 0.20 0.20 0.05
6245 0.20 0.15 0.20 0.20 0.20 0.05
6249 0.15 0.20 0.20 0.20 0.20 0.05
0154１３２人目の素数さん
2011/08/12(金) 03:46:54.73パターンは次の108通り。合計が20未満の場合は、足りない数だけ1を加える。変数が20個未満の場合は、足りない分だけ0を加える。
例：432→合計が9なので、1を11個加え、43211111111111、さらに数字の数が14個なので、0を6つ加え、「43211111111111000000」が真のパターン
44444 , 44443 , 444422 , 44442 , 4444 , 444332 , 44433 , 444322 , 44432 ,
4443 , 4442222 , 444222 , 44422 , 4442 , 444 , 443333 , 443332 , 44333 ,
4433222 , 443322 , 44332 , 4433 , 4432222 , 443222 , 44322 , 4432 , 443 ,
44222222 , 4422222 , 442222 , 44222 , 4422 , 442 , 44 , 433333 , 4333322 ,
433332 , 43333 , 4333222 , 433322 , 43332 , 4333 , 43322222 , 4332222 ,
433222 , 43322 , 4332 , 433 , 43222222 , 4322222 , 432222 , 43222 , 4322 ,
432 , 43 , 422222222 , 42222222 , 4222222 , 422222 , 42222 , 4222 ,
422 , 42 , 4 , 3333332 , 333333 , 3333322 , 333332 , 33333 , 33332222 ,
3333222 , 333322 , 33332 , 3333 , 33322222 , 3332222 , 333222 , 33322 ,
3332 , 333 , 332222222 , 33222222 , 3322222 , 332222 , 33222 , 3322 , 332 , 33 ,
322222222 , 32222222 , 3222222 , 322222 , 32222 , 3222 , 322 , 32 , 3 ,
2222222222 , 222222222 , 22222222 , 2222222 , 222222 , 22222 , 2222 , 222 , 22 , 2 ,1
各パターン毎に、変数の名前の入れ替えを考えると、合計35,561,166,195通りあると思われる。
0155１３２人目の素数さん
2011/08/16(火) 17:26:31.500156１３２人目の素数さん
2011/08/16(火) 20:01:09.46http://sites.google.com/site/geogebrajp/
0157１３２人目の素数さん
2011/08/17(水) 15:07:10.21式は最大で4つほど、xとyに関しては整数を想定しています。
z=ax+by … (i)
z=cx+dy … (ii)
単純に総当りで(x,y)に対して、zを最小化する式を選ぶだけなら適当に関数を作ればいけそうですが（スマートではないですね…）、その境界をどうやったらグラフに引けるかがわかりません。
どなたか教えてくだされば幸いです。
よろしくお願いします。
0158１３２人目の素数さん
2011/08/17(水) 18:16:26.16>z=ax+by … (i)
>z=cx+dy … (ii)
この式でzが等しくなるax+by=cx+dyを引けばよいです。
>式は最大で4つほど
組み合わせが6組あるから6本引いてください。
>xとyに関しては整数を想定しています。
境界線を引くときは整数という条件を無視し、
グラフを見る時に整数が線のどちら側にあるか判断してください。
0159１３２人目の素数さん
2011/09/01(木) 15:41:46.210160１３２人目の素数さん
2011/09/05(月) 16:54:51.220161１３２人目の素数さん
2011/09/05(月) 17:18:37.28Ｒ入門書ってないでしょうか
0162１３２人目の素数さん
2011/09/05(月) 20:58:24.38(Use R) Springer-Verlag
Jan Graffelman (2010) Multivariate Analysis With Matlab And R CRC Press Online
0163１３２人目の素数さん
2011/10/02(日) 01:52:40.21ワークスペースを保存しようとすると1kbのファイルしか保存されず、保存に失敗します。
原因がわかればよろしくお願いします。
0164163
2011/10/02(日) 02:08:38.40スレ汚しすいません。
0165１３２人目の素数さん
2011/10/04(火) 02:54:00.77plot(g)
こんな感じで半径1の円を描きたいのです。
教えてください。
0166１３２人目の素数さん
2011/10/04(火) 09:53:36.95symbols(0,0,circles=1,inches=FALSE)
お望みはこんな感じ？
0167１３２人目の素数さん
2011/10/04(火) 11:45:17.26symbols(0,0,circles=1,inches=F, xlim=c(0,2), ylim=c(0,2))
図示されたイメージの縦横比(大きさ)を変えると
円の大きさも変化してしまうのですが…
できれば、sqrt(x^2+y^2)の式を用いて描きたいです
0168152
2011/10/04(火) 12:39:57.11>sqrt(x^2+y^2)の式を用いて描きたい
じゃあそうすればいいじゃん。なぜためらう？
> x <- seq(-1,1,length.out=100)
> y <- sqrt(1^2-x^2)
> plot(x,y,type="l")
これで半円だろ。
> plot(c(x,x),c(y,-y),type="l")
これで円になるだろ。
0169１３２人目の素数さん
2011/10/04(火) 12:40:50.470170１３２人目の素数さん
2011/10/04(火) 12:49:40.400171１３２人目の素数さん
2011/10/04(火) 12:54:05.14念のために、ちょっと改良
> plot(c(x,rev(x)),c(y,-y),type="l")
用途によるけど、円を描画したければ、symbolsを使った方が、
アスペクト比に影響されずに正円になるのでよいと思うが。
xy座標で円を描くと、asp=1にでもしない限り、普通は楕円になってしまう。
0172１３２人目の素数さん
2011/10/12(水) 10:45:20.34のようにAICやBICを使うと、
selectlagsをFixedしてラグを変えて一通り試しても
AICやBICを用いた計算に一致する結果がない、
自分で計算したBICとでは選択されたラグが一致しないなど
どこか間違っている気がします。
ソースを見ても
critRes <-AIC(result, k = switch(selectlags, AIC = 2, BIC = log(length(z.diff))))
lags <- which.min(critRes)
となっていて、当方の利用の仕方に間違いがあるのかどうかも分かりません。
どなたか、この点についてご教示いただければと思います。
0173１３２人目の素数さん
2011/10/23(日) 19:16:11.40ロジスティック回帰のような、説明変数に対して単調増加する線形じゃなくて
多分ガウス曲線みたいな感じになると思うのだけれども
なんか適当なパッケージないですかね？
0174１３２人目の素数さん
2011/10/24(月) 21:24:37.50都道府県 学年 テスト点数
京都 1 85
大阪 2 90
京都 3 70
奈良 2 60
東京 _3 90
・・・
というようなデータフレームDFが存在します。
これを各都道府県別に学年毎のテストの平均点を算出したいのですが、
どうすれば良いのでしょうか?
最初に都道府県別にデータフレームを分割すれば
tapply関数などで集約できるのですが、いかんせんfor loopで回すので
時間がかかります。ループを利用しないような簡単な計算方法はないでしょうか?
府県名 <- sort(unique(DF[,1]))
for ( i in 府県名 ){
}
0175174
2011/10/24(月) 21:27:18.06途中で送信してしまいました。
現在は
府県名 <- sort(unique(DF[,1]))
for ( i in 府県名 ){
temp <- sub(DF, DF$府県名==i)
result <- tapply(temp, temp$テスト点数, mean)
}
のようにしてますがいかんせんループで時間がかかります。
よろしくお願いいたします。
0176１３２人目の素数さん
2011/10/24(月) 22:11:15.94整数型配列を二つ用意します。大きさは両方とも6*47=282です。
名前はsum、counter。最初は全て0を入れておきます。
配列の１番目は、北海道の１年に対応させ、２番目は青森の１年．．．４７番は沖縄の一年。
４８番は北海道の２年．．．２８２番目は沖縄の６年。
あるいは、１番は北海道の１年、２番は北海道の２年、．．．、６番は北海道の６年、．．．
２７７番が沖縄の１年、．．．２８２番は沖縄の６年。
どちらでもかまいません。
一人目のデータを読み込み、京都１年に対応するデータのsumを点数の分だけ増やし、counterは1増やす。
これを全データで行います。一度スキャンするだけです。
この目的には、ソートは必要ありません。
0177174
2011/10/24(月) 22:25:22.55ありがとうございます。
Cならそれで十分なのですが、Rの場合ベクトルで扱うのでループ処理は遅くなるらしいのです。
もう一点質問です。
CODE CD VAL
0001 101 3
0001 201 4
0001 333 5
0002 201 2
0002 333 7
・・・
というデータフレームがある場合に
CODE VAL[101] VAL[102] ........ VAL[201] .............
0001 3 0 ........ 4 ...............
というようなデータフレームに変換するにはどうすればよいでしょうか。
これも現在はループを利用して処理をしております
よろしくお願いいたします。
0178１３２人目の素数さん
2011/10/24(月) 22:51:31.71pref grade score
kyoto 1
kyoto 1 85
oosaka 2 90
kyoto 3 70
nara 2 60
tokyo _3 90
・・・
として、
tapply(DF$score, list(grade=DF$grade, pref=DF$pref), mean, na.rm=TRUE)
みたいのではいかんのですか？
いや、私もよくわからんので、Rコマンダーにお伺いをたててみたんですがね。
0179１３２人目の素数さん
2011/10/25(火) 02:06:48.72ヘルプを見るとインデックスはlist of factorS とあるから、複数のベクトルをいれれるんだね。勉強になりました。
tapply(X, INDEX, FUN = NULL, ..., simplify = TRUE)
Arguments
X
an atomic object, typically a vector.
INDEX
list of factors, each of same length as X. The elements are coerced to factors by as.factor.
0180１３２人目の素数さん
2011/10/25(火) 02:12:24.16http://cran.r-project.org/web/packages/HSAUR2/vignettes/Ch_analysing_longitudinal_dataI.pdf
の三ページのreshape()かな
0181１３２人目の素数さん
2011/10/25(火) 02:20:33.81＞CODE CD VAL
＞0001 101 3
＞0001 201 4
＞0001 333 5
＞0002 201 2
＞0002 333 7
＞・・・
＞
＞というデータフレームがある場合に
＞CODE VAL[101] VAL[102] ........ VAL[201] .............
＞0001 3 0 ........ 4 ...............
reshape(data, idvar="CODE", v.name="VAL", time="CD")
でどうでしょう。
0182１３２人目の素数さん
2011/10/25(火) 02:28:43.43reshape(data, idvar="CODE", v.name="VAL", timevar="CD", direction="wide")
でした。
0183１３２人目の素数さん
2011/10/25(火) 09:18:31.89精度を上げるための特徴選択を行おうと思い、
そのためのパッケージとして、"penalizedSVM"が
あると知りました（元々BioInformatics用みたいですが）。
http://r.789695.n4.nabble.com/Stepwise-SVM-Variable-selection-td3178718.html
http://cran.r-project.org/web/packages/penalizedSVM/index.html
しかし日本語の解説やリファレンスが皆無で、
英語でも分かりやすい解説がなく、
パッケージの公式リファレンスでは自分にはイミフでした。
http://cran.r-project.org/web/packages/penalizedSVM/penalizedSVM.pdf
どなたか実際にお使いの方がいれば、具体的に
どのように使用するか教えていただけないでしょうか。
018417
2011/10/25(火) 22:44:03.742つとも無事うまく解決できました。
やっぱりよくヘルプを読む必要がありますね。
精進します。
0185１３２人目の素数さん
2011/10/26(水) 00:50:54.901対1の全組み合わせのt.testを一気に行う方法があるでしょうか?
(loopで回す関数を自作したりしない出来合の方法で)
0186１３２人目の素数さん
2011/10/26(水) 07:28:10.01それは、R以前の問題として、不適切だと思うのですが。
あるいは駄目な例としてやるってことかな？
二元配置の分散分析か、何らかの多重比較ではだめなんでしょうか？
0187１３２人目の素数さん
2011/10/31(月) 21:41:20.58＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿
AAAAA 122 NA 2.5 3.8 NA
AAAAA 154 2.4 1 4 5
AAAAA 166 3.4 4.0 4 5
AAAAA 178 NA 33 NA NA
・
・
BBBBB 122 4 3.4 4.3 1
BBBBB 123 NA 2 4 4
BBBBB 154 2 5 5 6
・
・
というようなデータをもつデータフレーム
これをCODE_2毎に相関を算出したいと考えております。
現在は
allcode <- sort(unique(DF$CODE_2))
for ( c_code in 1:length(allcode)) {
TEMPDF <- subset(DF, DF$CODE_2==allcode[c_code])
cor(TEMPDF[, -1][, -1], use = "pairwise.complete.obs")
・・・
}
でおこなっているのですが、いかんせんCODE_2の数が多く、
ループ数が非常におおくなります。
解決策などありますでしょうか？
よろしくお願いします
0188１３２人目の素数さん
2011/11/01(火) 08:20:00.30code_2を因子化、tapply cor
0189１３２人目の素数さん
2011/11/01(火) 10:39:24.41でいいんじゃないの？
0190１３２人目の素数さん
2011/11/01(火) 11:23:41.66>>187じゃないけど、by()ってこんなに便利だったのか。
> V=c(NA,round(runif(10)*10,1))
> dat <- data.frame(C1=sample(LETTERS,50,replace=TRUE),
+ C2=sample(100:105,50,replace=TRUE),
+ V1=sample(V,50,replace=TRUE),
+ V2=sample(V,50,replace=TRUE),
+ V3=sample(V,50,replace=TRUE),
+ V4=sample(V,50,replace=TRUE))
> by(dat[,3:6],dat$C2,cor)
dat$C2: 100
V1 V2 V3 V4
V1 1 NA NA NA
V2 NA 1 NA NA
V3 NA NA 1 NA
V4 NA NA NA 1
------------------------------------------------------------
dat$C2: 101
V1 V2 V3 V4
V1 1 NA NA NA
V2 NA 1 NA NA
V3 NA NA 1.0000000 -0.4458286
V4 NA NA -0.4458286 1.0000000
------------------------------------------------------------
dat$C2: 102
V1 V2 V3 V4
V1 1.0000000 NA 0.1731169 NA
V2 NA 1 NA NA
V3 0.1731169 NA 1.0000000 NA
V4 NA NA NA 1
[以下略]
0191１３２人目の素数さん
2011/11/01(火) 11:47:20.47すまんこれ教えてくれ
0192１３２人目の素数さん
2011/11/01(火) 12:23:35.93見れねーよ
■ このスレッドは過去ログ倉庫に格納されています