【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/

343132人目の素数さん2018/02/27(火) 01:06:58.09ID:GkJoHhTZ
Windows環境で使うから

344132人目の素数さん2018/02/27(火) 01:16:41.39ID:hS0OJ3qQ
会社規定なんでしかたない
UNIX環境使えるのが羨ましい

345132人目の素数さん2018/02/27(火) 01:25:59.40ID:34tpClZB
Docker for WindowsでRStudioサーバー動かせば?

346132人目の素数さん2018/02/27(火) 12:19:33.31ID:ZcuB0rsn
なぜPythonスレッドがない?

347132人目の素数さん2018/02/27(火) 14:05:09.17ID:O+8uJ5V+
板名読める?

348132人目の素数さん2018/03/04(日) 20:14:10.91ID:nWbPE6nD
>>342
IME無効になるのは俺だけじゃなかったと知ってほっとした。

349132人目の素数さん2018/03/04(日) 21:28:34.58ID:R7ZPBSuG
ストアアプリも同じ症状でるからRStudio固有の問題でなくWindows環境の不治の病だと思ってる

350132人目の素数さん2018/03/26(月) 20:57:59.83ID:WAyucZm1
Run any R code you like. There are over three thousand R packages preloaded.

https://rdrr.io/snippets/

351132人目の素数さん2018/03/26(月) 21:11:25.16ID:WAyucZm1
>>350
日本語あると動作しないが、
グラフを描いてくれるのはうれしい。


http://imagizer.imageshack.com/img923/2879/m59LVv.png

352132人目の素数さん2018/04/04(水) 07:52:09.98ID:PZp+DZN4
Fisher test検定時に
p<2.2e-16
と表示されるんですが、これより小さい値の指数桁数を正確に表記する方法教えて下さい。
例えば5.8e-35となるようにです。

353132人目の素数さん2018/04/04(水) 08:05:49.72ID:lBVECsOD
返値の中にあるp値を参照しなされ

354132人目の素数さん2018/04/04(水) 10:02:52.31ID:CRvlhZKw
fisher.test関数の返り値はlist型で、その中にp.valueという名前でp値が格納されているから$演算子を使って直接参照するか、broom::tidy関数に返り値を渡してdata.frame形式で出力してやれば見れる

355132人目の素数さん2018/04/04(水) 12:43:50.08ID:LbKgW3kd
>>352
>>354が言う直接的な参照
> fisher.test(matrix(c(1,120,130,2),2))$p.value
[1] 1.691912e-69

356132人目の素数さん2018/04/04(水) 13:06:17.36ID:PZp+DZN4
352です。よく分かりました!
ありがとうございます!

357132人目の素数さん2018/04/04(水) 14:44:33.05ID:LbKgW3kd
>>356
技術的な助言をしたけど、学術的に言えば、
p < 0.01 は全て p < 0.01 として、具体的なp値を考える意味はないと思うよ。
一部の例外的な研究分野を除いて(e.g., 遺伝統計学)。

358132人目の素数さん2018/04/04(水) 15:18:50.72ID:PZp+DZN4
はい、まさにその例外的な分野で使おうとしてます。ありがとうございます。

359◆2VB8wsVUoo 2018/04/07(土) 06:52:38.98ID:yx+HETs3

360◆2VB8wsVUoo 2018/04/07(土) 06:53:00.80ID:yx+HETs3

361◆2VB8wsVUoo 2018/04/07(土) 06:53:20.77ID:yx+HETs3

362◆2VB8wsVUoo 2018/04/07(土) 06:53:40.58ID:yx+HETs3

363◆2VB8wsVUoo 2018/04/07(土) 06:54:01.92ID:yx+HETs3

364◆2VB8wsVUoo 2018/04/07(土) 06:54:24.32ID:yx+HETs3

365◆2VB8wsVUoo 2018/04/07(土) 06:54:47.76ID:yx+HETs3

366◆2VB8wsVUoo 2018/04/07(土) 06:55:08.41ID:yx+HETs3

367◆2VB8wsVUoo 2018/04/07(土) 06:55:30.46ID:yx+HETs3

368◆2VB8wsVUoo 2018/04/07(土) 06:55:54.53ID:yx+HETs3

369132人目の素数さん2018/04/07(土) 10:24:54.55ID:maRES4NF
>>357
多重比較だと意味あるかも

370132人目の素数さん2018/04/19(木) 21:35:35.54ID:GVMUXyX9
Rのガンマ関数はいくつでオーバーフローするかやってみた。

> i=1
> while(gamma(i)!=Inf){
+ i=i+1
+ }
Warning message:
In gamma(i) : value out of range in 'gammafn'
> i
[1] 172
> gamma(172)
[1] Inf
Warning message:
value out of range in 'gammafn'
> gamma(171)
[1] 7.257416e+306

371132人目の素数さん2018/04/29(日) 00:21:51.67ID:5dW+xNwa
matplot()で折れ線グラフ描いたときに、X軸をカテゴリで示したいのですが、
可能でしょうか?

例えばtemp <- c("0時間","8時間","24時間","48時間")として、
matplot()の引数にtempをとるやり方です。
他にもやり方あれば教えてください。

372132人目の素数さん2018/04/29(日) 01:55:59.75ID:PngyHOQZ
>>371
matplot(..., xaxt="n")
axis(1, at=seq(along=temp), lab=temp)

373132人目の素数さん2018/04/30(月) 16:50:51.65ID:t3vhzyao
>>372
遅くなりましたがありがとうございました。
できました!

374◆2VB8wsVUoo 2018/04/30(月) 23:55:02.74ID:y1TqbRSE

375◆2VB8wsVUoo 2018/04/30(月) 23:55:22.50ID:y1TqbRSE

376◆2VB8wsVUoo 2018/04/30(月) 23:55:43.19ID:y1TqbRSE

377◆2VB8wsVUoo 2018/04/30(月) 23:56:02.89ID:y1TqbRSE

378◆2VB8wsVUoo 2018/04/30(月) 23:56:23.17ID:y1TqbRSE

379◆2VB8wsVUoo 2018/04/30(月) 23:56:40.16ID:y1TqbRSE

380◆2VB8wsVUoo 2018/04/30(月) 23:56:53.49ID:y1TqbRSE

381◆2VB8wsVUoo 2018/04/30(月) 23:57:13.43ID:y1TqbRSE

382◆2VB8wsVUoo 2018/04/30(月) 23:57:32.05ID:y1TqbRSE

383◆2VB8wsVUoo 2018/04/30(月) 23:57:53.02ID:y1TqbRSE

384132人目の素数さん2018/05/01(火) 18:57:32.34ID:iUBwAKWd
特定の長方形の中に複数の長方形を最小面積で敷き詰める平面充填に関するパッケージってありませんかね

385132人目の素数さん2018/05/01(火) 19:35:44.26ID:iUBwAKWd
やっぱ自分でアルゴリズムを書いてみます

386132人目の素数さん2018/05/02(水) 00:09:39.04ID:xpXz1bJY
Package`deldir'
これどうなの?

387132人目の素数さん2018/05/06(日) 22:15:43.42ID:BK1CxH7U
# jonckheereテストを書いてみた

jonckheere <- function(L,
alternative = c("two.sided", "increasing", "decreasing"),
cat=TRUE){
# L : list of vectors A1,A2,...,Ak, with assumed tendency
How.Many.Greater.Pairs <- function(A,B){ # How many pairs of A[i] > B[j], count as 0.5 when equal,
n.a = length(A)
n.b = length(B)
how.many.greater.pairs = 0
for(i in 1:n.a){
for(j in 1:n.b){
how.many.greater.pairs = how.many.greater.pairs+ifelse(A[i]==B[j],0.5,A[i]>B[j])
}
}
return(how.many.greater.pairs)
}
Sum.of.Greater.Pairs <- function(L){ #L=list(A1,,,,Ak),A1 < A2 < A3,..,< Ak : vector
k = length(L)
comb = combn(1:k,2) # possible combinaition of pairs to compare
n.comb = ncol(comb) # how many combinations
J = 0 # sum of greater pairs
for(i in 1:n.comb){
J = J + How.Many.Greater.Pairs(L[[comb[1,i]]],L[[comb[2,i]]])
}
return(J)
}
J = Sum.of.Greater.Pairs(L)
n = sapply(L,length)
N = sum(n)
EJ = (N^2-sum(n^2))/4
VJ = (N^2*(2*N+3)-sum(n^2*(2*n+3)))/72
Z = (J-EJ)/sqrt(VJ)
alternative = match.arg(alternative)
p.value = switch(alternative, 'two.sided' = 2 * min(pnorm(Z), pnorm(-Z), 0.5),
'increasing' = pnorm(Z),
'decreasing' = pnorm(-Z))
if(cat){
cat( 'p.value = ', p.value,'\n')
cat('alternative hypothesis: ' ,alternative,'\n')
}
invisible(p.value)
}

388132人目の素数さん2018/05/22(火) 20:34:14.53ID:8iM7y57s
代入は「<-」派?「=」派?

389132人目の素数さん2018/05/22(火) 21:16:34.39ID:iB1pjrmI
>>388
実態調査か何か?
<-と=は挙動が違う場合があるので、使い分けていますが、
代入はどっちかと問われたら、無論 <- または ->

なお、
> 1 -> x
これはエラーにならないけど、
> 1 = x
1 = x でエラー: 代入の左辺が不正 (do_set) です
これはエラー

390132人目の素数さん2018/05/23(水) 11:17:56.08ID:OSJ/4EBd
>>389
俺は基本=派。

関数の定義は
z.test <- function(x,n=16,sigma=1){
z=sqrt(n)*mean(x)/sigma
2*pnorm(abs(z),lower=FALSE)
}
と書いている。

391132人目の素数さん2018/05/23(水) 11:24:57.67ID:MGQGuwX9
>1 = x でエラー

当たり前
そんな使い方なんてするかよ

他言語と同じく=一文字の方がすっきりしてイイ

392132人目の素数さん2018/05/24(木) 14:28:38.60ID:ExPgBsbL
こういうのが紛らわしいから、俺は = 推奨。

x <- 1
if(x <- 1) print('YES')
if(x < -1) print('YES')

393132人目の素数さん2018/05/24(木) 14:36:00.11ID:EQ5K0CF7
だよねぇ
<-良くない

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