読者です 読者をやめる 読者になる 読者になる

義務論理について考える

義務論理(Deontic logic)について考えてみます

Rで統計学を勉強

 

統計学入門 (基礎統計学)

統計学入門 (基礎統計学)

 

 Rで統計学の演習をします。教科書は↑など。

 

第3章 練習問題 3.4 「ブートストラップ」

 

X<-c(71,68,66,67,70,71,70,73,72,65,66)
Y<-c(69,64,65,63,65,62,65,64,66,59,62)

11家族内での兄弟(X)と姉妹(Y)の身長(単位:インチ)を定義。

 

RH<-numeric(length=200)

で、200個の相関係数を格納するベクトルを用意。>

 

P<-numeric(length=11)

Q<-numeric(length=11)

ランダムで抽出した重複を許したそれぞれ11人の兄弟・姉妹の身長を格納するベクトルを用意。

 

for(i in 1:200){
+ for(j in 1:11){
+ x<-ceiling(11*runif(1))
+ P[j]<-X[x]
+ Q[j]<-Y[x]
+ }
+ RH[i]<-cor(P,Q)

}

x<ceiling(11*runif(1))で1, 2, ..., 11の間の乱数をxに代入し、x番目の兄弟、姉妹の身長をそれぞれP, Qの成分とする。

 

RH[i]<-cor(P,Q)でP, Qの相関係数をベクトルRHのi番目の成分とし、この操作をi=1から200に対して繰り返す。

 

RHの出力はこんな感じ

[1]  0.55728214  0.57050127  0.54492414  0.12037289  0.62819438 -0.08205745
[7]  0.74958324  0.40779441  0.42724703  0.72669959  0.65698889  0.61878342
[13]  0.60464357  0.40730498  0.72086351  0.32640087  0.40207860  0.13190595
[19]  0.26411104  0.38569492  0.49337943 -0.23736079  0.30412522  0.35966362
[25]  0.60767428  0.74449427  0.79584666  0.84242801  0.66121530  0.71611047
[31]  0.75613186  0.36191582  0.94004157  0.76550517  0.30380274  0.76572215
[37]  0.83173726  0.53547268  0.05194916  0.80062453  0.33953346  0.51372783
[43]  0.51107145  0.28787727  0.25240061  0.68094811  0.61173661  0.68703547
[49]  0.42512473  0.63448400  0.08869934  0.90844215  0.54884769  0.55754439
[55] -0.11557517  0.56502971  0.76126892  0.84011821  0.59438059  0.30364884
[61]  0.49954122  0.50332500  0.39488507  0.26826997 -0.01806696  0.69173252
[67]  0.82930383  0.92035262  0.59145892  0.55456405  0.72038145  0.92113869
[73]  0.57143167  0.61066158  0.31216755  0.98953754  0.57915293  0.63475045
[79]  0.78401783  0.68385851  0.35848477  0.13164007  0.63865792  0.51439051
[85]  0.75689163  0.62079330  0.78901611  0.61226530  0.48528225  0.35077531
[91]  0.49672801  0.55805471  0.50671802  0.31494988  0.63902050  0.81051095
[97]  0.32776074  0.73868273 -0.11409536  0.80490388  0.45082259  0.70130104
[103]  0.44436879  0.57181053  0.53832710  0.91328594  0.50763725  0.80237742
[109]  0.83843312  0.90838152  0.14071698  0.03456969  0.81529516  0.55851781
[115]  0.29718923  0.73648711 -0.02180718  0.66803727 -0.07643805  0.59063650
[121]  0.36623351  0.86795434  0.84570843  0.63632078  0.29771015  0.46119557
[127]  0.34194836  0.79077086 -0.09673017  0.53832710  0.56918698  0.10195922
[133]  0.21449648  0.28286032  0.55546381  0.31360194  0.66407859  0.19146144
[139]  0.55183105  0.85599203  0.59303963  0.41966962  0.72082725  0.58297093
[145]  0.25009899 -0.14907120  0.58583205  0.68703547  0.78438278  0.54631228
[151]  0.79731755  0.25315008  0.53713752  0.87458105  0.43159605  0.19246032
[157]  0.53601911  0.24467110  0.36006791  0.77870084  0.40005754  0.83028467
[163]  0.13252098  0.72700788  0.28832267 -0.12575650  0.51209677  0.28606586
[169]  0.61263651  0.79066556  0.62665158  0.29315098  0.51327194  0.71040512
[175]  0.53206413  0.26842336  0.97010960  0.41973983  0.78167946  0.65687279
[181]  0.73876144  0.63217743  0.47809457  0.65335799  0.25421615  0.24540597
[187]  0.67237908  0.65782904  0.25968909  0.53783625  0.46075611  0.27257372
[193]  0.66655706  0.49531459  0.77023776  0.94028222  0.44632051  0.49883928
[199]  0.60999607  0.71264311
>

最後にhist(RH)でヒストグラムを出力

f:id:umemura_wataru:20170128004421p:plain