TOP カテ一覧 スレ一覧 100〜終まで 2ch元 削除依頼
Visual Studio 2019
Java Web Application Framework総合 ver2
【Lua】組み込み系言語総合 その7【Squirrel】
OpenMPプログラミング
文字コード総合スレ Part12
【PHP】下らねぇ質問はここに書き込みやがれ 2
Visual Studio Code / VSCode Part7
【DI】Java Spring Frameworkを語るスレ 5.0
Visual Studio 2012 Part8
Ruby>>>>>Java
統計解析R たぶんpart3くらい
- 1 :2012/06/22 〜 最終レス :2020/06/01
- あっても良さそうなんですが,需要ないの?
過去スレのうちのひとつ
統計解析R
http://toro.2ch.sc/test/read.cgi/tech/1285597198/
- 2 :
- このスレッドは天才チンパンジー「アイちゃん」が
言語訓練のために立てたものです。
アイと研究員とのやり取りに利用するスレッドなので、
関係者以外は書きこまないで下さい。
京都大学霊長類研究所
- 3 :
- Rのスレは別の板にあったと思うけど
- 4 :
- 数学板にあったけど、あっちはサーバーが飛んじゃったから、こっちに立ててもいいんじゃないかな。
- 5 :
- R ってデータベース?
- 6 :
- R言語
- 7 :
- SQLみたいなもの
- 8 :
- Rはゴミカス向けだと断言する
俺はゴミカスだがエリートゴミカスだ
お前らのような下級ゴミカスとは格が違う
- 9 :
- 数学版から誘導されてきたけど、
Rjpwikiと裏Rjpwikiのヲチはここでよいの?
- 10 :
- いいんじゃないの
- 11 :
- これのクラスとかのオブジェクト指向言語的な機能って
なんか目指すところが普通の言語と違ってる感じだよな
- 12 :
- >>11
具体的にどのように?
- 13 :
- programmerというよりはstatisticianが使いやすいことを念頭にしているから?
- 14 :
- >>13
R言語って、統計解析ソフトのマクロ言語みたいなものだから、
一般のプログラミング言語と比較する方が間違っていると思う。
Rの基本関数はCやfortranで書かれたものを呼び出しているに過ぎない。
- 15 :
- なんでこの板にスレ立てたの?
- 16 :
- Rを使ってる奴等はなにに使っているのか
- 17 :
- 俺はゴミカスだがエリートゴミカスでR
お前らのような下級ゴミカスとは格が違うのでR
- 18 :
- ゴミ
- 19 :
- そういえば S のスレはなんでないの?
- 20 :
- Mのスレならあるぞ
ちなみに俺はMだ
- 21 :
- >>16
10個くらいのパラメータあるときに、
pairsで、総当りの相関見るのは便利
- 22 :
- これってマスマティカみたいな使い方できるの?
- 23 :
- >>22
つmaxima
- 24 :
- ありがとうございます。
maximaも荒らしが凄いのねプログラム板に移動したいですね。
- 25 :
- 2ch の見過ぎで ステマ ティカに見えた。
- 26 :
- ●関連サイト
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
- 27 :
- 【R言語】統計解析フリーソフトR 第4章【GNU R】
http://uni.2ch.sc/test/read.cgi/math/1294561909/
- 28 :
- データを分割してスペクトルを測定しようとして次のようなプログラムを書いてみたんですが、
「以下にエラー spec.pgram(x, ...) : 添え字が許される範囲外です 」
と出て動かないですが何がおかしいんでしょうか?
行ごとに実行して動いているんですけど。
function(x)
{
n <- length(x)
l <- 40000
m <- floor(n/l)
for (i in 0:(m -1))
{
spectrum(ts(x[i*l+1 : (i+1)*l],deltat=.0025))
}
}
- 29 :
- >>28
カッコ追加してみた。
spectrum(ts(x[(i*l+1) : ((i+1)*l)],deltat=.0025))
- 30 :
- 馬鹿には無理
- 31 :
- RなんかよりもIncanterを使うべき
- 32 :
- >>31
>Incanter
知らなかったので調べてみたら、
本当にExcel統計に負けるほどの基本関数しかなさそう。
経済、医学、生態とか分野別になってもいないし、
ベイズ推定とかデジタル地図の空間演算とか何もなさそう。
> RなんかよりもIncanterを使うべき
お薦めの点はどういうところ?
- 33 :
- 俺が作った(きりっ
- 34 :
- lispで書ける。それだけだが、それこそが何ものにもかえがたい。
- 35 :
- groovyからR呼べないの?Javaって機械学習周りのOSSツール沢山あるのに
最近、Apache族がGroovyをデファクトで採用しそうな感あるよね
- 36 :
- >>31
と言っている本人も使っていないくせに。まだpytnonでnumpyとscikits statsmodels & pandasを使う方がまだ実用性がある。
- 37 :
- Pythonにはthink.statsて物もあるからね。
普及はしてナイっぽいが。
- 38 :
- Excelとか信用できない
- 39 :
- タイトルに次からはR言語って入れろや
次スレは5年は後だろうけど
- 40 :
- >>39
RのことをR言語というは2重の意味で反対って、
御大が仰っていたが。
- 41 :
- 人いないの?
- 42 :
- >>41
いるけど、何故か質問は数学版の方にしか来ない
- 43 :
- なる
- 44 :
- >>42
>>41
数学版
http://uni.2ch.sc/test/read.cgi/math/1380168442/
- 45 :
- ↓こんなんあるんだね。なかなか良さそう。
http://www.rstudio.com/
- 46 :
- >>45
ESSを使えよ
- 47 :
- え?
- 48 :
- ESSって何?
- 49 :
- >>48
Emacs Speaks Statistics
世界の主なR開発者とユーザが使っているEmacs用のインターフェイス
昔はRStudioとかなくて、pmgとかRKWardとかはあったけど、
あまり使い物にならなかったので、事実上、ESSしか選択肢がなかった。
だから、Rの古参ヘビーユーザはほとんどESSユーザ。
# Emacsを団塊世代の化石と揶揄するのはなしの方向で
- 50 :
- ESSで、なにげにJuliaも使える
便利
- 51 :
- 何故、EmacsはVimより重いのか
- 52 :
- こんなとこでEmacsとVimの論争はやめてよ!
- 53 :
- なにこれ言語?
- 54 :
- Rって汎用言語よりも遥かに短く式書けていいよね
- 55 :
- >>53
行列演算をサポートするスクリプト言語。
おもに統計の手法開発向けに考えられたS言語というのがあって
それの兄弟言語。
無償なんで結構つかっている人は多い。
- 56 :
- gsにならなくて良かったと思う
- 57 :
- 確かに
- 58 :
- 今、本を読んで、
多くの常識はずれなことに当惑している
ここはプログラム板だから、
プログラムの観点から言わせてもらうと、
まず行ではなく、列方向に、
データがまとめられていることに、
ショックを受けた
次に、C言語の配列がなく、
リストを配列と呼んでいるらしい
もう、ついていけん!これは文系脳?
- 59 :
- 論理的に述べよ
- 60 :
- Lisp脳からすると配列とか言われると頭がモヤモヤする。
- 61 :
- >>58
Rの目的は、自分が考えたオレ様統計量やモデルを簡単に実装すること。
既存の言語に過度に振り回されず、簡潔明瞭に統計モデルが記述できればそれでよい。
> 次に、C言語の配列がなく、
> リストを配列と呼んでいるらしい
arrayとlistは別物だけど、そこは理解している?"["関数は理解している?
> a <- array(1:9, dim = c(3, 3))
> a[2, 3]
[1] 8
> b <- list(x = matrix(1:6, 2), y = letters[1:5])
> b
$x
[,1] [,2] [,3]
[1,] 1 3 5
[2,] 2 4 6
$y
[1] "a" "b" "c" "d" "e"
> b[[1]][2, 2]
[1] 4
- 62 :
- あっ、逆か。
Cのリスト構造がRの配列に相当すると主張しているのか?
でも、そうすると余計に意味不明。
データフレームが変数ごとにまとめられているのは、
データベース的観点から言えばあたりまえ。
- 63 :
- もともと、Sの世界でデータフレームというものはなかった。
あくまで手法開発ツールなので、それでよしということだったが
データ処理という観点で必要性が高いため追加された。
ただし、データ管理という点では今だ無いに等しいので、ビジネスで
使うには他ソフトとのコラボレーションがどうしても必要。
- 64 :
- データベース的には行指向のほうが普通じゃない?
最近は列指向DBもあるけど
構造体とかオブジェクトとかとは逆だけどそれはそういうもんだと思うしかないのでは
- 65 :
- >データベース的には行指向のほうが普通じゃない?
実世界のデータ自体がそういう形だからね。
数理系の目的以外で使用するのはちと厳しい。
- 66 :
- >>58
行列データの物理的な格納方式のことを言っているのだったら、
C言語は行列の格納は行指向(row-major)だけど、
Fortranだと列指向(column-major)なので、Rが列指向なのはむしろ自然。
- 67 :
- RDBに慣れた人がKVSに文句言ってるくらい可笑しい話
- 68 :
- >次に、C言語の配列がなく、
>リストを配列と呼んでいるらしい
漏れの勘違い
Rのベクトルは、C言語の静的配列で、
同じデータ型の要素が並ぶもの
Rのリストは、C言語の構造体で、
異なるデータ型の要素が並ぶもの
Rのデータフレームは、Rのリストの集積
「アート・オブ・R・プログラミング」を読んでいるが、
オライリーの本は外人なので、
わかりにくいし、疲れる。時間のムダ
アンクの絵本みたいに、
図で表示すればすぐわかるのに
- 69 :
- >>68
文句ばっかりw
しかも間違っているし。
- 70 :
- >>58
数値計算専用のFortranは列指向なわけだが
- 71 :
- 68みたいに何の素養もないボケがRを使うような機会ってあるの?
- 72 :
- ひょっとして手当たりしだいに履歴書だせば、Rでデータ解析するような仕事に付けるの?
- 73 :
- >>72
試してみれば?
- 74 :
- >>74
http://detail.chiebukuro.yahoo.co.jp/qa/question_detail/q14130946449
- 75 :
- やはりそのようなパッケージはないということでしょうか?
- 76 :
- >>75
Tucker法とかなら、ThreeWayパッケージ
- 77 :
- ご返答ありがとうございます
three Wayパッケージについて調べてみたいと思います。
- 78 :
- ・グラフを重ねて描画
・グラフを並べて描画
この2つの両立は可能でしょうか?
以下の様な流れのコードで、可能かと思ったのですが、
ひとつのグラフエリアに全グラフを出力してしまい、実現できませんでした。
よろしくお願いします。
//描画領域を1行4列に
jpg('test.jpg')
par(mfcol=c(1,4))
//4回繰り返す
while(4回繰り返す){
//1つのグラフに、3本の線を描画
while(3回繰り返す){
par(new=T)
hogehoge = 描画したいデータ
plot(hogehoge,type="l")
}
}
//画像ファイル出力
dev.off()
- 79 :
- >>78
> ・グラフを重ねて描画
> ・グラフを並べて描画
> この2つの両立は可能でしょうか?
可能と言えば可能。
ただし、1つのグラフに、3本の線≠グラフを重ねること
複数の折れ線グラフを描画するときには、
plot()で軸などを描画した後に、lines()で線を描くのが基本。
matplot()を使うのもあり。
plot(runif(10), runif(10), type = "l")
lines(runif(10), runif(10), lty = 2, col = 2)
lines(runif(10), runif(10), lty = 3, col = 3)
これを1行4列に
op <- par(mfcol=c(1,4))
sapply(1:4, function(i){
hogehoge <- as.data.frame(x = runif(30), y = runif(30))
plot(hogehoge[1:10, ], type = "l")
lines(hogehoge[11:20, ], lty = 2, col = 2)
lines(hogehoge[21:30, ], lty = 3, col = 3)
})
par(op)
みないな感じにする。
なお、// のコメント禁止。代入に=を使うのは禁止。
初心者はpar(new=T)の使用禁止。
while()も真にやむを得ない場合以外は避ける(回数が決まっているループに使うのはもってのほか)。
- 80 :
- 画像ファイル形式に、表(行列)を出力する方法ってありますか?
無ければデータフレームをwriteでcsvに書き出して、スクショとってペイント切り貼りなので、
できれば避けたいです・・・
- 81 :
- gridExtraパッケージのgrid.table()でいける
- 82 :
- >>81
質問者じゃないけど、参考になった。
複雑な図にたくさんの数字ラベルをつけて、
横に数字とラベル名の対照表を簡単につけられないかと長年考えていたが、
これなら簡単そうだな。
op <- par(mfrow = c(1,2))
x <- data.frame(x = runif(10), y = runif(10))
plot(x, cex=2)
i <- 1:10
text(x, as.character(i))
grid.table(data.frame(Month = month.name[i]))
par(op)
だが、古典グラフィックとGridの共存がうまくいかないw
- 83 :
- 突っ込まれる前に追記。
もちろん、下記のようにすれば、ある程度実現することは分かっているけど、
ちゃんとレイアウトされた表にすることが難しい
x <- data.frame(x = runif(10), y = runif(10))
plot(x, xlim = c(0, 1.5), cex=2)
i <- 1:10
text(x, as.character(i))
legend("topright", paste(formatC(1:10, width=2),
month.name[i]), title = "Month", bty = "n")
- 84 :
- >>81
ありがとうございます。
無事表示できました
贅沢を言えば、
>>82の方が言うように
既存のplotと共存させてなくなってきました。
ライブラリを読みましたが、既存のplotとの共存はできないのでしょうか
@A
@プロットしたグラフ
Aその値の表
- 85 :
- >>84
古典グラフィックと混載するのではなく、
全てをGrid化して制御すれば、自由自在に描画できると思うよ。
とりあえず、viewport()を勉強してくれ
ちょっと長くなるけど、下記のコードは期待通りになった。
library(gridExtra)
grid.newpage()
pushViewport(viewport(x = 0, width = 0.8, just="left"))
pushViewport(plotViewport())
grid.rect(); grid.xaxis(); grid.yaxis()
x <- runif(10); y <- runif(10)
grid.points(x, y, size = unit(1.5, "char"))
grid.text(as.character(1:10), x = x, y = y)
grid.text("Grid test plot", y = unit(1, "npc") + unit(2, "lines"))
popViewport()
pushViewport(viewport(x = 1.2, width = .2, just = "right"))
grid.table(data.frame(Month = month.name[1:10]))
- 86 :
- 【動物】外へ出たネコはどこへ行くのか? 1000匹のデータから解明するプロジェクト [8/8]
http://daily.2ch.sc/test/read.cgi/newsplus/1407513849/
- 87 :
- RのクラスってC言語上でクラスを作るみたいなもんだな
- 88 :
- pygmentsを使ってみたが、
$ pygmentize -N tmp.R
splus
$ pygmentize -N tmp.r
rebol
拡張子が小文字だとREBOLと認識される模様。
- 89 :
- 代入にほとんど<-が使われていますが、
=を<-と同じとして使ってよいのでしょうか?
- 90 :
- >>89
だめ。意味が違うし、結果が異なる場合もある。
ベテランでも<-にすべきところで=を使っている人もいるので、
混乱すると思うけど、S言語の仕様に従った方が間違いがない。
<-と=で結果が変わってしまう例:
> (x <- 9:11)
[1] 9 10 11
> mean(x = 1:3); mean(x)
[1] 2
[1] 10
> mean(x <- 1:3); mean(x)
[1] 2
[1] 2
> x
[1] 1 2 3
引数に<-を使うと、mean(x)の値が予期せず変わってしまう。
9,10,11の平均値を2として論文に書いてしまう人はいないだろうが、
複雑な計算だったら、気がつかずにそのまま解析結果にしてしまうかも知れない。
- 91 :
- >>90
ありがとうございました
最近勉強始めたのですが、教えていただけなければ全部=で覚えるところでした
また機会ありましたらよろしくおねがいします
- 92 :
- ちょっとお聞きしたいのですが、
x y z f
1 1 1 f1
1 1 2 f2
1 1 3 ・・
1 2 1
1 2 2
1 2 3
1 3 1
1 3 2
・・
のようなデータがあって、
例えばzが1だけのx,y,fだけで3次元サーフェスプロットを書きたいのですが、
Rだとそういったことはもしかしてできたりしますでしょうか?
(可能なら他のz=2,z=3でのサーフェスプロットも同時に表示するなど・・)
突然ですいませんがもしよろしければご教示いただけると幸いです
- 93 :
- subsetで抜き出せばいいんじゃないの?
- 94 :
- aggregate
- 95 :
- >>92
> dat <- data.frame(expand.grid(1:3, 1:3, 1:2)[,3:1], as.factor(paste0("f",1:18)))
> names(dat) <- c("x", "y", "z", "f")
> dat
x y z f
1 1 1 1 f1
2 1 1 2 f2
3 1 1 3 f3
4 1 2 1 f4
5 1 2 2 f5
6 1 2 3 f6
[中略]
17 2 3 2 f17
18 2 3 3 f18
というデータがあったとして、
> (dat1 <- dat[dat$z == 1, ])
x y z f
1 1 1 1 f1
4 1 2 1 f4
7 1 3 1 f7
10 2 1 1 f10
13 2 2 1 f13
16 2 3 1 f16
とすればzが1だけのデータになる。これを行をx、列をyとしたときのfの行列を作成してpersp()で表示すればよい。
> m <- matrix(NA, nrow = max(dat1$x), ncol = max(dat1$y))
> for(i in 1:nrow(dat1)){m[dat1[i, 1], dat1[i, 2]] <- dat1[i, 4]}
> persp(m)
factor型のfを強制的に整数型にしているけど、適宜、自分の思うとおりにアレンジしてくれ。
- 96 :
- >>93-96
ありがとうございました
特に具体的なやり方まで教えていただき、感謝です
試させていただこうと思います
- 97 :
- マスゴミ・売国奴・医療業界が隠そうとする真実---------------------安楽死---------------------奴隷に勝手に死なれては困る
安楽死旅行企画が大人気|竹田恒泰チャンネル
https://www.youtube.com/watch?v=XmP1TRsAe88
武田邦彦:安楽死と大麻、そして売春・・・オランダに学ぶ
https://www.youtube.com/watch?v=nWV8YOY39tw
安楽死党
https://www.youtube.com/watch?v=8nU2UaSlGx0
自殺は後遺症が怖い!だから-----------------------------------安楽死制度-------------------------------------安心して生きるために
- 98 :
- まだあったのかこのスレ
- 99 :
- bash スクリプト内でRのスクリプトを起動して統計処理させてその結果を
1か0かで戻り値($?)として返すにはどうすれば良いでしょうか?
- 100 :
- シェルスクリプト test.sh の内容
#!/bin/bash
r --vanilla --slave < test.R
# コマンドの直後に、終了コードを見る
echo $?
100〜のスレッドの続きを読む
任天堂「今後C++は捨てJavaScriptで開発していく」
JAVAってこんなことも出来ないの?
BASICの宿題はお前にまかせた
Android Studio Part3
1行ずつC言語を書いてくスレ(目標なし)
【計測】LabVIEW相談室【制御】
日下部陽一著 作ってわかるCプログラミング(第6版)
人工知能はゼルダの伝説をクリアできないだろうな
Deep learning
【論理】Prolog【初心者】
--------------------
BTSの雑談スレ1326
レーザー方式メガネスレ opt.1
わくまゆこと和久田麻由子アナウンサー vol.19
【QMA】ニュースクイズで使われそうな時事ネタ書き込むスレ3
ウイイレ2013の不具合・不満・要望をぶちまけるスレ
HuGっと!プリキュア 第31話「時よ、すすめ!メモリアルキュアクロック誕生!」 [757453285]
💪ダンベルで筋トレNo.6
エレガンス系ブログ・インスタヲチ32
【米国】アメリカのオフィスでは、意外にも社内恋愛が生まれにくい
【MHW】ランススレ 94突き目
ゲーム開発者「ソニーにPCとXbox Oneのクロスプレイを入れないように指示された」
【植物】日中韓の「桜の起源」論争に終止符?=中国専門家「奪うな、桜は日本のものだ」[03/26]
【難病】膠原病全般 その4
クッカー総合スレ Part33
GMスレッドPart70
東大理三女子のハロウィンの仮装がこちら [309927646]
ガチャガチャは昔は一回20円〜だった。
断酒ネット@アルコール依存症『断酒継続・97』
御朱印ブロガーウォッチスレ25
【9984】ソフトバンクグループ【627】
TOP カテ一覧 スレ一覧 100〜終まで 2ch元 削除依頼