【R言語】統計解析フリーソフトR 第5章【GNU R】at MATH
【R言語】統計解析フリーソフトR 第5章【GNU R】 - 暇つぶし2ch918:132人目の素数さん
17/04/26 11:59:53.14 tqLIcCJs.net
>>876
Stanって知らんかった。
便利そうだな使ってみる。
ありがとう!

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

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

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

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

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

924:ま
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がやればいい。わけで、それを有限にしており、できなかった、のだ
だから、元素が悪いのだすべて火、土、水、空に原因がある

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

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

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

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

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

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


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


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


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


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


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


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


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


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


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


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


R Developersより

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

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

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


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


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


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


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


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


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


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


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


952:¥
17/05/26 09:18:09.48 iCP5fMHR.net


953: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が逆になってしまうのはなぜなんでしょう。
整数を代入するだけでは整数だと予測してくれないんでしょうか?

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


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

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

a <- 5L

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

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


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

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


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


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


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


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


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


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


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


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


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


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

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

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

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

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

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

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

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



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


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


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


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


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


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


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


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


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


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


979: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))

よろしくお願いします。

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

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



982: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が最小になっている

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



984: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")

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

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

986: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)

987: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")

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

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


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


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


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


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


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


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


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


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


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


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

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


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

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



1003: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()を使うべき

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

1005:132人目の素数さん
17/07/25 00:41:44.53 BWjSvXL+.net
数学板はなんで荒らし投稿ばかりなの?
ロボットのしわざかな?
常連さんのPC�


1006:ェ乗っ取られているの?



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

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


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


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


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


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


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


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


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


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


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


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

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



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

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



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


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


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


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


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


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


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


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


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


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


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

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


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


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


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


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


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


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


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


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


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


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


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


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


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


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


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

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


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