【R言語】統計解析フリーソフトR 第5章【GNU R】at MATH
【R言語】統計解析フリーソフトR 第5章【GNU R】 - 暇つぶし2ch967: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