【R言語】統計解析フリーソフトR 第6章【GNU R】 [無断転載禁止]©2ch.net

1132人目の素数さん2017/08/03(木) 19:23:12.67ID:Hq1blL0O
R は統計計算とグラフィックスのための言語・環境です。
統計計算で重宝するデータ型や、複数要素を処理する演算や関数、
解析結果を表示するグラフィックなど、多彩な機能を提供します。

●関連サイト
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
※前スレ
【R言語】統計解析フリーソフトR 第5章【GNU R】
http://rio2016.2ch.net/test/read.cgi/math/1380168442/

289132人目の素数さん2017/11/25(土) 01:57:31.82ID:sNF3t9bW
>>278
>argsは見れる関数と見れない関数があるのであんまり役に立たない
だったら、例が悪いよ。見られない関数の例をあげないと

> ヘルプで見ようと思ったら書いてなかった
こっちだってわざわざ下請け関数のオプションまでは書かないよ。
ヘルプを見て未解決ならソースを見てくれ。

290◆2VB8wsVUoo 2017/11/25(土) 06:17:05.29ID:kb2pRtxG

291◆2VB8wsVUoo 2017/11/25(土) 06:17:23.25ID:kb2pRtxG

292◆2VB8wsVUoo 2017/11/25(土) 06:17:39.26ID:kb2pRtxG

293◆2VB8wsVUoo 2017/11/25(土) 06:17:55.24ID:kb2pRtxG

294◆2VB8wsVUoo 2017/11/25(土) 06:18:11.86ID:kb2pRtxG

295◆2VB8wsVUoo 2017/11/25(土) 06:18:28.49ID:kb2pRtxG

296◆2VB8wsVUoo 2017/11/25(土) 06:18:44.12ID:kb2pRtxG

297◆2VB8wsVUoo 2017/11/25(土) 06:19:01.27ID:kb2pRtxG

298◆2VB8wsVUoo 2017/11/25(土) 06:19:19.56ID:kb2pRtxG

299◆2VB8wsVUoo 2017/11/25(土) 06:19:36.14ID:kb2pRtxG

300132人目の素数さん2017/11/26(日) 17:41:27.37ID:ahD9OIUm
rm(list=ls())

で R のメモリーをクリア

301132人目の素数さん2017/11/28(火) 13:00:14.06ID:kbsEO2Dt
先週長々と質問させていただいたものです。その後、試行錯誤して行った結果十解決しました。

以前は「summary」としか入れていなかったのですが、「dd$summary」と入力したところ
全体のデータが出力されました。それ以後「dd$summary$dir.area」等で出せるようになりました。
大変初期的な場面で躓いておりましたが、現在何とか進んでおります。
色々とお答えくださった方々、とても参考になりました。本当にありがとうございました。

302132人目の素数さん2017/12/10(日) 13:11:49.22ID:cG4z9akW
N(=100)回コインをなげてn(=5回)以上続けて表がでる確率。

seqn<-function(n=5,N=100,p=0.5){
rn=rbinom(N,1,p)
count=0
for(i in 1:N){
if(rn[i] & count<n){
count=count+1
}
else{
if(count==n) {return(TRUE)}
else{
count=0
}
}
}
return(count==n)
}
mean(replicate(10^5,seqn()))

> mean(replicate(10^5,seqn()))
[1] 0.81085

案外、高い確率になった。

303132人目の素数さん2017/12/10(日) 22:28:22.98ID:WuWuRvn7
ポアソンでやるやつだっけ

304132人目の素数さん2017/12/11(月) 07:47:07.43ID:EURogHTh
pooledVariance <- function(...) {
args = list(...)
n.args=length(args)
ss2=0
df=0
for(i in 1:n.args){
ss2 = ss2 + var(args[[i]])*(length(args[[i]])-1)
df = df + (length(args[[i]])-1)
}
ss2/df
}

effectsize <- function(y1,y2){
diff=mean(y1)-mean(y2)
var=(var(x1)*(length(x1)-1)+ var(x2)*(length(x2)-1))/(length(c(y1,y2))-2)
sd=sqrt(var)
diff/sd
}

library(effsize)
cohen.d()

305132人目の素数さん2017/12/11(月) 07:50:22.91ID:EURogHTh
>>303
単なる二項分布。
コインが5回続けて表がでたら、0.5^5 <0.05なのでイカサマコインといわれちゃいそうなんだが、
100回やってみると案外、5回表が続くので確率を計算しようと思ったが、解析的にできる頭がないので
シミュレーションしてみた。

306132人目の素数さん2017/12/11(月) 07:55:42.42ID:EURogHTh
>302は

> rbinom(100,1,0.5)
[1] 0 1 0 0 1 0 0 0 1 1 1 0 1 0 1 1 0 1 1 0 1 1 0 0 0 1 0 1 1 1 1 0
[33] 0 1 0 0 0 0 0 0 1 0 1 1 0 0 0 1 0 0 1 0 1 1 1 1 1 1 1 1 1 1 1 0
[65] 0 0 1 1 0 0 1 0 1 0 0 0 0 0 0 1 0 0 1 1 1 1 1 0 1 1 0 0 1 1 0 0
[97] 0 0 1 1

で5回以上1が連続するときTRUEを返す関数なのだが
もっと簡単にやれないかなぁとは思っている。

rep(1,5) %in% rbinom(100,1,0.5)は 1個ずつ評価されてTRUEが5個返ってくるだけ。

文字列にしてgrepを使うとなんとかなりそうな気がしないでもないのだけど....

307132人目の素数さん2017/12/22(金) 17:21:16.51ID:U1Bq8bdb
すまん、教えてほしいだけど
分析するために初めてRをインストールしようと思って、このスレのあるように公式サイト行ったら、esetが「JS/Redirector.NAV トロイの木馬」を検知したんだが…;
どうしたらいいだ…

308132人目の素数さん2017/12/22(金) 20:09:03.17ID:ol6XbcBE
そっちじゃなくてCRAN行け

309132人目の素数さん2017/12/22(金) 20:42:31.02ID:U1Bq8bdb
ホームページダメとかわけわかんねぇ・・・とんでもねぇな
ありがとう

310132人目の素数さん2017/12/22(金) 20:50:20.52ID:uzZ0hejh
>>307
エロサイトにアクセスしてないw?

311132人目の素数さん2017/12/22(金) 21:27:01.31ID:RNoKRtjI
>>309
CRANはここ
https://cran.r-project.org


プロジェクトの方はESETの誤検知っぽいんだよな

312132人目の素数さん2017/12/23(土) 09:28:44.66ID:n0SBd+bp
>>310
してねぇよww
totalvirusで調べると一件引っかかるし、Redirector検知だから、しょうがないね

313132人目の素数さん2017/12/24(日) 20:42:51.55ID:CT/NKMd7
col=rgb(runif(1),runif(1),runif(1),runif(1))で色指定すると

走らせるたびに色がちがっておもしろい。

hist(rnorm(100),col=rgb(runif(1),runif(1),runif(1),runif(1)))

314132人目の素数さん2017/12/31(日) 17:23:19.35ID:14tdpK/Y
stanやJAGSのコードでgamma関数を使おうとして
y = gamma(x)
と、書いたらエラーになった。

stanだと y=tgamma(x)、JAGSだとy=exp(loggam(x))で動作した。

315132人目の素数さん2018/01/02(火) 08:46:50.95ID:qdmBZ37O
ある大学の入学者男女の比率は1であるという帰無仮説を検定する課題が花子と太郎に課された。

花子は50人を調査できたら終了として入学者を50人をみつけて18人が女子であるという結果を得た。
帰無仮説のもとで
50人中18人が女子である確率は 0.01603475
これ以下になるのは50人中0〜18人と32〜50人が女子の場合なので
両側検定して
> sum(dbinom(c(0:18,32:50),50,0.5))
[1] 0.06490865
> binom.test(18,50,0.5)$p.value
[1] 0.06490865
で帰無仮説は棄却できないと結論した。
http://i.imgur.com/XDIp9rM.png

一方、本番と十八番が好きな太郎は一人ずつ調べて18人めの女子がみつかったところで調査を終えることにした。
18人めがみつかったのは花子と同じく50人めであった。
帰無仮説のもとで
18人がみつかるのが50人めである確率は0.005772512
これ以下になるのは23人以下50人以上番めで女子18人めがみつかった場合なので
両側検定して
pnb=dnbinom(0:999,18,0.5)
> 1 - sum(pnb[-which(pnb<=dnbinom(50-18,18,0.5))]) # < 0.05
[1] 0.02750309
http://i.imgur.com/K3T7utr.png
で帰無仮説は棄却される。

どちらの検定が正しいか、どちらも正しくないか?
検定する意図によってp値が変わるのは頻度主義統計の欠陥といえるか?

花子の横軸は裏口入学者数、太郎の横軸はサンプル数なので
サンプルでの裏口入学率を横軸にして95%信頼区間を示す。
花子の検定での信頼区間は0.36〜0.72で18/50を含む、p=0.06491
http://i.imgur.com/SeTLk8K.jpg
太郎の検定での信頼区間は0.375〜0.72で18/50を含まない、p= 0.0275
http://i.imgur.com/tNzlfxe.jpg
主観である、検定の中止の基準の差でp値や信頼区間が変化するのは変だという批判である。

316132人目の素数さん2018/01/03(水) 10:35:48.54ID:YJfyxrv+
(訂正)
ある大学の入学者男女の比率は1であるという帰無仮説を検定する課題が花子と太郎に課された。

花子は50人を調査できたら終了として入学者を50人をみつけて18人が女子であるという結果を得た。
帰無仮説のもとで
50人中18人が女子である確率は 0.01603475
これ以下になるのは50人中0〜18人と32〜50人が女子の場合なので
両側検定して
> sum(dbinom(c(0:18,32:50),50,0.5))
[1] 0.06490865
> binom.test(18,50,0.5)$p.value
[1] 0.06490865
で帰無仮説は棄却できないと結論した。
http://i.imgur.com/XDIp9rM.png

一方、十八という数字が好きな太郎は一人ずつ調べて18人めの女子がみつかったところで調査を終えることにした。
18人めがみつかったのは花子と同じく50人めであった。
帰無仮説のもとで
18人がみつかるのが50人めである確率は0.005772512
これ以下になるのは23人以下50人以上番めで女子18人めがみつかった場合なので
両側検定して
pnb=dnbinom(0:999,18,0.5)
> 1 - sum(pnb[-which(pnb<=dnbinom(50-18,18,0.5))]) # < 0.05
[1] 0.02750309
http://i.imgur.com/K3T7utr.png
で帰無仮説は棄却される。

どちらの検定が正しいか、どちらも正しくないか?
検定する意図によってp値が変わるのは頻度主義統計の欠陥といえるか?

花子の横軸は女子数、太郎の横軸はサンプル数なので
サンプルでの女子の割合を横軸にして95%信頼区間を示す。
花子の検定での信頼区間は0.36〜0.72で18/50を含む、p=0.06491
http://i.imgur.com/SeTLk8K.jpg
太郎の検定での信頼区間は0.375〜0.72で18/50を含まない、p= 0.0275
http://i.imgur.com/tNzlfxe.jpg
主観である、検定の中止の基準の差でp値や信頼区間が変化するのは変だという批判である。

317132人目の素数さん2018/01/15(月) 16:30:41.76ID:wJofbCL/
kainokousiki<-function(a,b,c){return (-b+sqrt(b^2-4*a*c))/(2*a)} #解の公式
kainokousiki(1,-5,6)

でrunすると3じゃなくて6を返すんだけど、どこが間違ってる?

318132人目の素数さん2018/01/15(月) 16:35:00.20ID:wJofbCL/
自己解決
かっこが足りなかった

319132人目の素数さん2018/01/19(金) 12:12:23.00ID:k3h1PrrK
Pythonのスレはないのか

320◆2VB8wsVUoo 2018/01/21(日) 20:11:34.35ID:oUqQkvBY

321◆2VB8wsVUoo 2018/01/21(日) 20:11:56.34ID:oUqQkvBY

322◆2VB8wsVUoo 2018/01/21(日) 20:12:22.56ID:oUqQkvBY

323◆2VB8wsVUoo 2018/01/21(日) 20:12:47.97ID:oUqQkvBY

324◆2VB8wsVUoo 2018/01/21(日) 20:13:05.14ID:oUqQkvBY

325◆2VB8wsVUoo 2018/01/21(日) 20:13:34.57ID:oUqQkvBY

326◆2VB8wsVUoo 2018/01/21(日) 20:14:02.66ID:oUqQkvBY

327◆2VB8wsVUoo 2018/01/21(日) 20:14:16.99ID:oUqQkvBY

328◆2VB8wsVUoo 2018/01/21(日) 20:14:35.47ID:oUqQkvBY

329◆2VB8wsVUoo 2018/01/21(日) 20:15:07.01ID:oUqQkvBY

330132人目の素数さん2018/02/05(月) 19:08:34.64ID:AIJT+apj
機械学習をきっかけにPythonに逆転された感じだね

331132人目の素数さん2018/02/06(火) 18:35:39.90ID:tAZA/Fp/
『Rを使った〜』だとPythonじゃないのかよって思うよね

332132人目の素数さん2018/02/06(火) 21:37:32.10ID:tUqX17n9
ずっとRだけでPython触ったこと無いけど、覚え直す価値ある?
環境構築からもう面倒なイメージ

333132人目の素数さん2018/02/06(火) 22:43:23.82ID:qQMNyZjW
Python自体は:と直後のインデントさえ気を付ければかなり簡単
3系は数が全て小数扱いなので楽

Anacondaというパッケージでインストールすれば、今流行りのJupyter Notebookという開発環境で対話的にコーディングできる(Rも使える)

334132人目の素数さん2018/02/07(水) 20:40:49.48ID:mgaw9oVw
アナコンダてのがRぽくできるのね、ありがとう
dplyrやggplot2みたいに素人でも簡単便利だといいんだけど
Pythonはオブジェクト志向ぽいしすぐ諦めそう

335132人目の素数さん2018/02/12(月) 16:19:11.21ID:NSJ4iUa4
ブラウザ環境なくなっちゃったの?

336132人目の素数さん2018/02/12(月) 17:06:56.33ID:hMO/kWPg
>>335
誤爆?
何をブラウズする環境が無くなったの?

337132人目の素数さん2018/02/13(火) 18:21:54.45ID:1hic99Cx
ブラウザでプログラミングする環境

338132人目の素数さん2018/02/13(火) 21:09:53.07ID:EFzybxrC
>>337
RStudio Serverの話?

339132人目の素数さん2018/02/13(火) 22:27:26.57ID:Qyl0hNg7
RStudio Severあるで。Dockerで使うのがいいんじゃないかな?!

新着レスの表示
レスを投稿する