質問編集履歴
3
タグを追加しました。
test
CHANGED
File without changes
|
test
CHANGED
File without changes
|
2
追記、編集しました。
test
CHANGED
File without changes
|
test
CHANGED
@@ -1,5 +1,5 @@
|
|
1
|
-
|
1
|
+
遺伝統計学の基礎という書籍に書かれているにR言語のコードで、理解の難しい箇所がございます。
|
2
|
-
コード
|
2
|
+
コードは無料公開されております。
|
3
3
|
また、Rに関する質問の場合、タグは何をつけるのが適切でしょうか?
|
4
4
|
R,R言語と入れても、該当するものが出てこないようなのですが・・・。
|
5
5
|
|
@@ -28,4 +28,6 @@
|
|
28
28
|
sibIdValueで、0:2kとしているということは、0,1,2,...,2kの2k+1個のデータ?が作成されていると思うのですが、
|
29
29
|
さらにforループを23回回す中で、identityにmumchの1から23までを入力してますが、これってつまり、
|
30
30
|
forが1の時は0:2、2の時は0:4、3の時は0:6みたいに増えていってしまい、それを毎回plotしていると
|
31
|
-
かなり重複する事になると思うのですが
|
31
|
+
かなり重複する事になると思うのですが・・・
|
32
|
+
|
33
|
+
これは実際に重複してプロットしているのでしょうか??
|
1
コードを追記いたしました、コードは公開されており、だれでもDLできるものです。
test
CHANGED
File without changes
|
test
CHANGED
@@ -2,3 +2,30 @@
|
|
2
2
|
コードの理解の難しい箇所がございます、こちらで質問させていただいても良いでしょうか?
|
3
3
|
また、Rに関する質問の場合、タグは何をつけるのが適切でしょうか?
|
4
4
|
R,R言語と入れても、該当するものが出てこないようなのですが・・・。
|
5
|
+
|
6
|
+
```R
|
7
|
+
sibIdValue<-function(k=1){
|
8
|
+
(0:(2*k))/(2*k) # i/(2k)
|
9
|
+
}
|
10
|
+
sibIdProb<-function(k=1){
|
11
|
+
dbinom(0:(2*k),2*k,0.5)
|
12
|
+
}
|
13
|
+
numch<-1:23
|
14
|
+
means<-vars<-rep(0,length(numch))
|
15
|
+
for(i in 1:length(numch)){
|
16
|
+
identity<-sibIdValue(numch[i])
|
17
|
+
prob<-sibIdProb(numch[i])
|
18
|
+
if(i>1){par(new=TRUE)}
|
19
|
+
plot(identity,prob,type="b",ylim=c(0,1))
|
20
|
+
means[i]<-momentProb(identity,prob,order=1,center=FALSE)
|
21
|
+
vars[i]<-momentProb(identity,prob,order=2,center=TRUE)
|
22
|
+
}
|
23
|
+
plot(means)
|
24
|
+
plot(vars)
|
25
|
+
```
|
26
|
+
|
27
|
+
具体的なコードを1つ貼らせて頂きます、たとえばこれなのですが・・・
|
28
|
+
sibIdValueで、0:2kとしているということは、0,1,2,...,2kの2k+1個のデータ?が作成されていると思うのですが、
|
29
|
+
さらにforループを23回回す中で、identityにmumchの1から23までを入力してますが、これってつまり、
|
30
|
+
forが1の時は0:2、2の時は0:4、3の時は0:6みたいに増えていってしまい、それを毎回plotしていると
|
31
|
+
かなり重複する事になると思うのですが、どのような構造になっているのでしょうか?
|