【R言語】統計解析フリーソフトR 第5章【GNU R】at MATH
【R言語】統計解析フリーソフトR 第5章【GNU R】 - 暇つぶし2ch901:¥ ◆2VB8wsVUoo
16/10/23 02:16:16.15 3ZIBPF5r.net


902:¥ ◆2VB8wsVUoo
16/10/23 02:16:33.13 3ZIBPF5r.net


903:¥ ◆2VB8wsVUoo
16/10/23 02:16:51.23 3ZIBPF5r.net


904:¥ ◆2VB8wsVUoo
16/10/23 02:17:08.24 3ZIBPF5r.net


905:¥ ◆2VB8wsVUoo
16/10/23 02:17:24.20 3ZIBPF5r.net


906:¥ ◆2VB8wsVUoo
16/10/23 02:17:40.94 3ZIBPF5r.net


907:¥ ◆2VB8wsVUoo
16/10/23 02:17:58.13 3ZIBPF5r.net


908:¥ ◆2VB8wsVUoo
16/10/23 02:18:16.48 3ZIBPF5r.net


909:¥ ◆2VB8wsVUoo
16/10/23 02:18:34.98 3ZIBPF5r.net


910:132人目の素数さん
16/11/13 17:45:07.40 AJbUUUy0V
R-3.3.1   MagaiHin
R-3.2.5   Last

911:132人目の素数さん
16/11/14 05:39:42.16 LliUCPaqU
無料誹謗中傷対策のブッダワークス
URLリンク(www.buddha-works.net)

912:132人目の素数さん
17/03/18 22:10:26.69 9lKmiZM8.net
盛り上げていこうぜ!

913:132人目の素数さん
17/03/19 08:25:18.19 BUbMna52.net
2chはゴミ

914:132人目の素数さん
17/03/30 22:10:28.44 K8rGc65l.net
Rってどれくらい使えたら役に立つんだ?ベイズ階層位はできんとだめ?

915:132人目の素数さん
17/03/31 16:49:02.43 LyGz7dL7.net
目的によるんじゃない? 考えてみ

916:132人目の素数さん
17/04/23 08:56:55.19 qEAq5gtK.net
R 3.4.0
URLリンク(www.r-project.org)

917:132人目の素数さん
17/04/24 23:49:44.87 38+Bilc8.net
>>873
Stanを使えばベイズ階層も苦にならない。
Samplingの待ち時間が苦痛かも。

918:132人目の素数さん
17/04/25 18:50:16.56 M0sS+FAa.net
>>875
おかげでアップデートに気づいたわ

ちなみにyou stupid darknessってどういう意味やろか?

919:132人目の素数さん
17/04/25 19:39:43.28 24RIix2u.net
URLリンク(kurtisscaletta.com)

920:132人目の素数さん
17/04/25 19:46:55.22 24RIix2u.net
1. なにそれ?
2. 暗いと不平を言うよりも、すすんで灯りをつけましょうって聞いたんだ。
3. それはそうだけど、世の常としてそう思わない人もいるだろうね。
4. 暗いぞバカヤロー!

って感じかな?

921:132人目の素数さん
17/04/26 11:52:55.96 tqLIcCJs.net
取引先の、統計をやってるという人と話をした。
名刺には「シニアスタティスティシャン」と書いてある。

話しをしてみると、統計学のダメな入門書に書いてあるような
大ウソを堂々と主張する。

だが相手は取引先であるから、「馬鹿かてめーは!」とか言えない。

言えないからここに書くのだが、
名刺に「シニアスタティスティシャン」とか書くなよ馬鹿たれが!

922:132人目の素数さん
17/04/26 11:59:53.14 tqLIcCJs.net
>>876
Stanって知らんかった。
便利そうだな使ってみる。
ありがとう!

923:869
17/04/26 12:31:23.28 Vq01HD2H.net
>>878
>>879
ありがとう
ピーナッツが読めるぐらいに英語が堪能ならなぁ

924:ま
17/04/28 19:19:21.81 10ctknwL.net
a+b=C

925:ま
17/04/28 20:57:46.81 10ctknwL.net
>>883今の現代数学は間違ってできてるため、すうじでフっ素をあらわしてる。
そのため、人をすいとる。養分の話しか、してない。数学が人を、おしつぶして、時をすいとってしまった
数式をかえよう
X(f)が、あらたな、人を救う。数式。だということになるよね

926:ま
17/04/28 21:05:59.12 10ctknwL.net
4代元素が、すべて、支払いをし、責任をおうべきだと、思う。この世は、人間をつくり、だまし、最後には、すう式で、すいとってしまい、4代元素は次第には、元素になり、死にいたらしめる、などフリーメーソンなどつくった、元素からの安易な暗号の挑戦状なんだと思う

927:ま
17/04/29 02:13:36.09 PXwKZAGU.net
>>883 a+b=abでいい。c+c=ccでいいと思うんだけど…
同じものを+んだから、同じものになればいいわけで。a×90乗=aで、いいいわけじゃない。なんか、変化して元素に、なるなら、回りのものも、全部元素にならないと、いけないのに、自分たちは、固くのこってて、これを、やりたいなら、人が死んでから、やるべき

928:ま
17/04/29 02:39:27.48 PXwKZAGU.net
>>886固くのこされた人の気持ちを考えるべき。だE=MC2にしたって
E=C2で、やるべきなのに、添字にMを、そえて、そのエネルギーの座標に、人間をさしかえてるだけの、話しでね、
時間一次元と空間3次元の4成分だけ、もたされて、それを、吸う、すう、しき、とか、うけとるだけやないか。元素って。最初からE=C2で、この世を回すべき。
E=Mにかきかえるべきだ。それが、正しいと思うし。正当的。アルケーが歩け

元素はまちがった空間をつくっている
E=C2、または、E=M、で、場をつくるべき
相対性理論はまちがっている
死後の世界観を座標に4つ、ひきつける空間をもたせて、座標に吸わせることを、表しており、正当性のかけらもない
死後なら死後で座標をCがやればいい。わけで、それを有限にしており、できなかった、のだ
だから、元素が悪いのだすべて火、土、水、空に原因がある

929:132人目の素数さん
17/04/29 02:58:57.49 OVQ3EEWI.net
馬鹿が現れた

930:132人目の素数さん
17/05/01 17:04:16.76 xPREFfVb.net
-δC2こうで

931:132人目の素数さん
17/05/01 17:09:14.52 xPREFfVb.net
>>883я

932:132人目の素数さん
17/05/01 17:12:13.05 xPREFfVb.net
タダシさ

933:132人目の素数さん
17/05/01 17:19:22.91 xPREFfVb.net
タダシさ=こたえ

934:¥
17/05/08 20:06:42.83 OR+quqWp.net


935:¥
17/05/08 20:14:24.76 OR+quqWp.net


936:¥
17/05/08 20:14:49.35 OR+quqWp.net


937:¥
17/05/08 20:15:09.69 OR+quqWp.net


938:¥
17/05/08 20:15:31.60 OR+quqWp.net


939:¥
17/05/08 20:15:52.66 OR+quqWp.net


940:¥
17/05/08 20:16:14.65 OR+quqWp.net


941:¥
17/05/08 20:16:41.12 OR+quqWp.net


942:¥
17/05/08 20:17:02.39 OR+quqWp.net


943:¥
17/05/08 20:17:24.76 OR+quqWp.net


944:132人目の素数さん
17/05/18 12:12:04.69 b7yGK+hkS
ウィンドゥズにインストールされたRではR検定に合格できません。
R専用のスパコンでなければ、すべての解析ができないからです。
Fortran2008の関数が導入してあるコンパイラーがあれば、検定用Rのパッケージが構築可能だとか
GCCはバージョン4.12.5で導入予定。完成は未定です。


R Developersより

945:132人目の素数さん
17/05/25 19:49:56.38 zLanYWay.net
数値の20170501を2017-05-01ってな感じで日付に変換したいんですが、何か良い方法はないでしょうか。
力技は考えたんですが、100万行くらい変換しないといけないので、力技以外で変換したいんです。

946:132人目の素数さん
17/05/25 20:11:21.15 t+im+FUk.net
as.Date("20170501",format = "%Y%m%d")

947:¥
17/05/26 09:14:55.04 iCP5fMHR.net


948:¥
17/05/26 09:15:13.28 iCP5fMHR.net


949:¥
17/05/26 09:15:33.73 iCP5fMHR.net


950:¥
17/05/26 09:15:55.04 iCP5fMHR.net


951:¥
17/05/26 09:16:19.41 iCP5fMHR.net


952:¥
17/05/26 09:16:43.19 iCP5fMHR.net


953:¥
17/05/26 09:17:05.17 iCP5fMHR.net


954:¥
17/05/26 09:17:26.87 iCP5fMHR.net


955:¥
17/05/26 09:17:48.95 iCP5fMHR.net


956:¥
17/05/26 09:18


957::09.48 ID:iCP5fMHR.net



958:132人目の素数さん
17/06/16 15:45:46.35 r+tPpOy0.net
R version 3.1.1ユーザです。

> a = 5
> is.double(a)
[1] TRUE
> is.integer(a)
[1] FALSE

のようにTRUEとFALSEが逆になってしまうのはなぜなんでしょう。
整数を代入するだけでは整数だと予測してくれないんでしょうか?

959:¥
17/06/16 16:12:39.43 TDPjrtUP.net


>112 名前:¥ ◇2VB8wsVUoo 2017/06/15(木) 08:39:06.70 ID:7Rm0/VD6
> ★★★数学徒は情熱的な霊感により主観的に暮らし、唯ひたすら自己の世界に沈潜すべき。★★★
>
> 佐藤幹夫を見よ、ラマヌジャンを見よ
>
> 本物の¥
>

960:132人目の素数さん
17/06/16 17:21:40.28 nDRgYZZV.net
>>916
その代入だと実数として扱われる

a <- 5L

とすれば整数として扱われるようになる

961:¥
17/06/16 17:44:03.74 TDPjrtUP.net


>112 名前:¥ ◇2VB8wsVUoo 2017/06/15(木) 08:39:06.70 ID:7Rm0/VD6
> ★★★数学徒は情熱的な霊感により主観的に暮らし、唯ひたすら自己の世界に沈潜すべき。★★★
>
> 佐藤幹夫を見よ、ラマヌジャンを見よ
>
> 本物の¥
>

962:¥
17/06/16 19:31:56.18 TDPjrtUP.net


963:¥
17/06/16 19:32:16.87 TDPjrtUP.net


964:¥
17/06/16 19:32:37.07 TDPjrtUP.net


965:¥
17/06/16 19:32:57.87 TDPjrtUP.net


966:¥
17/06/16 19:33:17.21 TDPjrtUP.net


967:¥
17/06/16 19:33:35.76 TDPjrtUP.net


968:¥
17/06/16 19:33:55.49 TDPjrtUP.net


969:¥
17/06/16 19:34:15.63 TDPjrtUP.net


970:¥
17/06/16 19:34:36.83 TDPjrtUP.net


971:¥
17/06/16 19:34:57.09 TDPjrtUP.net


972:907
17/06/17 10:07:37.47 oboZrWLI.net
>>918
ユーザーが明示的に指定してやらないとすべて実数になっちゃうんですね?
お答えくださってどうも有難うございました。

973:¥
17/06/17 10:26:55.03 x2f6D4gs.net
★★★忖度と処世術に汚染された日本人:権威主義的な支配と損したくない人達★★★
  ~~~芳雄氏が言う『研究者としての基本的態度』とは一体何だろうか~~~

佐藤幹夫:自分自身の素朴な疑問に真剣に耳を傾ける。⇒不滅の金字塔を打ち立てる。
糞父芳雄:人間関係を駆使し他人を操り根回しを行う。⇒ハリボテお教授として君臨。

隠蔽の財務省、嘘吐きの文科省、そして問答無用に屈服させる官邸。コレでも先進国?

(佐藤師がしてたのは本物の研究だ。だが)芳雄氏がしてたのはケケケ、ケンキュウ。
外見を繕って偉そう見せさえすれば何でもヨロシ。ほんで教授になりさえすれば研究の
中身なんて何でもヨロシ。そもそも論文なんてモンは、外国の権威ある雑誌に掲載され
さえすれば、その中身のギロンなんて何でもヨロシ。そやし適当に書いてしまえ~~~
中身がダメだと知ってて、ソレでもSTAP論文を外国に投稿して受理される。発覚したら
適当に言い逃れる醜い態度。オツムのダメな大学院生に「虚偽の良品ラベル」を貼って
世間に出荷するハリボテ大学は詐欺行為そのもの。世間に媚びを売って客商売に徹し、
『売れさえすれば学生の脳の質なんて何でもヨロシ』と居直る大学。そしてブランド名
だけを見て仕入れる世間。●●は一流大学やさかい、きっと優秀なエリートやろwww

中身を何も説明しないで、問答無用に上から押し付ける。ソレをイチャモンで騒いで、
そして邪魔して潰そうとする周囲の下々。大学教員も国会議事堂も、そして馬鹿板人の
遣ってる事も皆同じだ。日本人はバカ民族であり、今は外国にもちゃんとバレてるので
海外からも軽蔑されるだけであり、そのうちにどの国からも信用されなくなるだろう。

近視眼的で打算的な人生観を息子に押し付ける父親と、大脳に栄養が足りてない連中が
跋扈する永田町や霞が関に支配される国に住む不幸、一体どうしてくれるというのか。

☆☆☆数学徒が馬鹿板をしたらダメ。さもないと国家議事堂みたいになります。☆☆☆



974:¥
17/06/17 12:35:00.60 x2f6D4gs.net


975:¥
17/06/17 12:35:20.71 x2f6D4gs.net


976:¥
17/06/17 12:35:41.80 x2f6D4gs.net


977:¥
17/06/17 12:36:00.84 x2f6D4gs.net


978:¥
17/06/17 12:36:21.78 x2f6D4gs.net


979:¥
17/06/17 12:36:42.46 x2f6D4gs.net


980:¥
17/06/17 12:37:03.93 x2f6D4gs.net


981:¥
17/06/17 12:37:23.71 x2f6D4gs.net


982:¥
17/06/17 12:37:59.24 x2f6D4gs.net


983:¥
17/06/17 12:38:23.00 x2f6D4gs.net


984:132人目の素数さん
17/07/23 16:07:33.29 AqYrJa+q.net
下のy式をもっとスマートに書く方法を教えてください
つまり、yの式は実際には+1, +2 …, +100までしなくてはいけないのです。
x <- 1:10
y <- c(sum((x+1)^2), sum((x+2)^2), sum((x+3)^2), …, sum((x+100)^2))

よろしくお願いします。

985:132人目の素数さん
17/07/23 16:33:37.21 FF49HrjI.net
z<-c()
for(i in 1:100){
z<-c(z, sum(x+i)^2))
}
でもできるけどダメ?

986:¥
17/07/23 16:35:16.97 r8UcZ6ry.net
〒〒〒馬鹿板は悪い習慣であり、この行為は脳を悪くする。そやし足を洗いなさい。〒〒〒



987:132人目の素数さん
17/07/23 17:11:12.91 AqYrJa+q.net
>>943
ありがとうございます。一つ解決しました。
しかし、続いて教わった方法で、平均の最小2乗法による最小値を図示しようと思ったのですが
seq(0, 10, .01)を入れると、x軸がなぜか100倍されてしまします。
なぜなのか分かりませんか?おそらく、このforを使った解決策が
その後の応用(最小値を図示する)に適切でないのでしょうが
それがなんでなのか・・・(私はforについて、勉強し直します)

x <- c(3,5,6,8,2,9)
mean(x)

z<-c()
for(i in seq(from=0, to=10, by=.01)){
z<-c(z,sum(x-i)^2)
}
plot(z, type ="l") # x=5.5が最小のはずなのに、550が最小になっている

988:¥
17/07/23 17:45:58.37 r8UcZ6ry.net
〒〒〒馬鹿板は悪い習慣であり、この行為は脳を悪くする。そやし足を洗いなさい。〒〒〒



989:132人目の素数さん
17/07/23 23:11:56.44 k8lcunzY.net
横軸の数値用ベクトルwにサンプリング数値を保存しておいて、
グラフ化時のx軸の値として指定してみてはどうか?

x <- c(3,5,6,8,2,9)
mean(x)

z<-c()
w<-c()
for(i in seq(from=0, to=10, by=.01)){
z<-c(z,sum(x-i)^2)
w<-append(w, i)
}
plot(w, z, type ="l")

990:132人目の素数さん
17/07/23 23:14:46.35 k8lcunzY.net
w<-append(w, i)

w<-c(w, i)
の方が簡潔で上との整合性も保てるか。

991:132人目の素数さん
17/07/24 00:56:57.68 7iMASKV5.net
>>942
z <- outer(x, 1:100, function(x, y) (x + y)^2)
apply(z, 2, sum)

992:132人目の素数さん
17/07/24 01:07:34.75 7iMASKV5.net
>>945
y軸の値しか指定しないから、x軸の値は1から始まる整数になる。
xv <- seq(0, 10, by=0.01)
z <- outer(x, xv, function(x, y) (x - y)^2)
yv <- apply(z, 2, sum)
plot(xv, yv, type="l")

993:132人目の素数さん
17/07/24 05:19:21.68 0E8vffzB.net
>>947-950
たくさんの返信ありがとうございます。
仕事から帰ったら、試してみます。

994:¥
17/07/24 09:36:16.88 /2fQM0DG.net


995:¥
17/07/24 09:36:35.46 /2fQM0DG.net


996:¥
17/07/24 09:36:55.11 /2fQM0DG.net


997:¥
17/07/24 09:37:12.58 /2fQM0DG.net


998:¥
17/07/24 09:37:31.25 /2fQM0DG.net


999:¥
17/07/24 09:37:48.24 /2fQM0DG.net


1000:¥
17/07/24 09:38:05.80 /2fQM0DG.net


1001:¥
17/07/24 09:38:23.56 /2fQM0DG.net


1002:¥
17/07/24 09:38:48.61 /2fQM0DG.net


1003:¥
17/07/24 09:39:05.54 /2fQM0DG.net


1004:描
17/07/24 09:53:25.77 3peLAFj6.net
$

1005:¥
17/07/24 11:30:57.34 /2fQM0DG.net


1006:132人目の素数さん
17/07/24 17:15:25.78 RJfWHfZU.net
この熱いあぼーん達はなんなんだろ… <


1007:¥
17/07/24 17:33:19.17 /2fQM0DG.net
〒〒〒馬鹿板は悪い習慣であり、この行為は脳を悪くする。そやし足を洗いなさい。〒〒〒



1008:132人目の素数さん
17/07/24 21:41:07.75 J3GDpKdM.net
>>950
添削。
> system.time(for(i in 1:10){yv <- apply(z, 2, sum)})
ユーザ システム 経過
0.032 0.005 0.037
> system.time(for(i in 1:10){yv <- colSums(z)})
ユーザ システム 経過
0.002 0.000 0.001
ということで、apply()ではなくcolSums()を使うべき

1009:132人目の素数さん
17/07/25 00:40:01.16 w3drJ+gw.net
>>966
初心者には、速度を追及するよりは、応用の効くやりかたを学んでもらったほうが良いのでは?

1010:132人目の素数さん
17/07/25 00:41:44.53 BWjSvXL+.net
数学板はなんで荒らし投稿ばかりなの?
ロボットのしわざかな?
常連さんのPCが乗っ取られているの?

1011:132人目の素数さん
17/07/25 06:50:18.42 7svw9fzv.net
>>947, 939
ありがとうございます。x軸の値の指定が間違っていたのですね。
y軸の値をx軸に指定していても図の形が同じなため気がつきませんでした。
しかし、y軸の値をx軸にあてはめても図の形が同じなのはなぜだろうか?
というか、y軸を指定しないでx軸だけ指定したのに、エラーが出ず
それなりの図(plot)が書けているのはなぜだ?
世の中分からんことばかりです。

>>949, 941, 957
outer()とapply()も実は使ったことなくて、いろいろいじってみて勉強します。
この分野は計算速度とか測って極めていくのか。
colSums()とか初めて見ました。

>>967
働き出してから、統計を勉強し直している初心者です。
ご迷惑おかけします。
基本的には論文(心理学・医療系)を理解するための勉強です。
ふと、このスレで、とんちんかんな質問をしていたら私かもしれません。

1012:¥
17/07/25 08:34:59.00 1OMr9h78.net


1013:¥
17/07/25 08:35:16.52 1OMr9h78.net


1014:¥
17/07/25 08:35:33.91 1OMr9h78.net


1015:¥
17/07/25 08:35:50.65 1OMr9h78.net


1016:¥
17/07/25 08:36:07.97 1OMr9h78.net


1017:¥
17/07/25 08:36:25.20 1OMr9h78.net


1018:¥
17/07/25 08:36:41.70 1OMr9h78.net


1019:¥
17/07/25 08:36:58.72 1OMr9h78.net


1020:¥
17/07/25 08:37:19.20 1OMr9h78.net


1021:¥
17/07/25 08:37:35.98 1OMr9h78.net


1022:132人目の素数さん
17/07/25 21:48:57.22 PpGdU5wz.net
>>969
x軸の値を指定せずにplotすると、横軸はデータのインデックス(順番)が表示される。
インデックスに対応する、本当の値をx軸の値としてplotしても、当然同じ図形になる。

1023:¥
17/07/25 22:08:54.95 1OMr9h78.net
〒〒〒馬鹿板は悪い習慣であり、この行為は脳を悪くする。そやし足を洗いなさい。〒〒〒



1024:132人目の素数さん
17/07/25 22:22:42.60 ylT3F2qf.net
次スレは>>990あたり?

1025:¥
17/07/25 22:24:57.98 1OMr9h78.net
〒〒〒馬鹿板は悪い習慣であり、この行為は脳を悪くする。そやし足を洗いなさい。〒〒〒



1026:¥
17/07/26 05:47:15.01 e+L8tFA6.net


1027:¥
17/07/26 05:47:35.53 e+L8tFA6.net


1028:¥
17/07/26 05:47:52.90 e+L8tFA6.net


1029:¥
17/07/26 05:48:08.88 e+L8tFA6.net


1030:¥
17/07/26 05:48:25.31 e+L8tFA6.net


1031:¥
17/07/26 05:48:42.11 e+L8tFA6.net


1032:¥
17/07/26 05:48:59.88 e+L8tFA6.net


1033:¥
17/07/26 05:49:26.67 e+L8tFA6.net


1034:¥
17/07/26 05:49:42.86 e+L8tFA6.net


1035:¥
17/07/26 05:50:00.74 e+L8tFA6.net


1036:132人目の素数さん
17/08/02 19:58:43.15 CTkpktyp.net
自由度10のt分布に従う乱数を50個作りたい時は
rt(50,10)
でokでしょうか

1037:¥
17/08/02 20:01:22.80 6XC/wDpJ.net


1038:¥
17/08/02 22:21:30.33 6XC/wDpJ.net


1039:¥
17/08/02 22:21:48.18 6XC/wDpJ.net


1040:¥
17/08/02 22:22:05.33 6XC/wDpJ.net


1041:¥
17/08/02 22:22:21.58 6XC/wDpJ.net


1042:¥
17/08/02 22:22:39.25 6XC/wDpJ.net


1043:¥
17/08/02 22:22:57.70 6XC/wDpJ.net


1044:¥
17/08/02 22:23:16.14 6XC/wDpJ.net


1045:¥
17/08/02 22:23:36.18 6XC/wDpJ.net


1046:¥
17/08/02 22:23:53.86 6XC/wDpJ.net


1047:¥
17/08/02 22:24:11.53 6XC/wDpJ.net


1048:¥
17/08/02 22:24:30.81 6XC/wDpJ.net


1049:¥
17/08/02 22:24:47.55 6XC/wDpJ.net


1050:¥
17/08/02 22:25:05.15 6XC/wDpJ.net


1051:猫
17/08/02 22:25:35.27 6XC/wDpJ.net


1052:1001
Over 1000 Thread.net
このスレッドは1000を超えました。
もう書けないので、新しいスレッドを立ててくださいです。。。
life time: 1406日 9時間 18分 13秒

1053:過去ログ ★
[過去ログ]
■ このスレッドは過去ログ倉庫に格納されています


最新レス表示
レスジャンプ
類似スレ一覧
スレッドの検索
話題のニュース
おまかせリスト
オプション
しおりを挟む
スレッドに書込
スレッドの一覧
暇つぶし2ch