我很乐意定期发行,比如模拟一个公平的模具.
但是,如果我想模拟以下内容怎么办:
当n为奇数时,X_n = 2(概率为0.6)或-1(概率为0.4).
当n为偶数时,X_n = -2(概率为0.6)或1(概率为0.4).
有什么建议?
试试这个:
sample(c(2,-1), size=100, replace=TRUE, prob=c(0.6,0.4))
sample(c(-2,1), size=100, replace=TRUE, prob=c(0.6,0.4))> n<-20
Run Code Online (Sandbox Code Playgroud)
前20个数字的一个例子:
n<-20
odd_n<-seq(1,n,by=2)
even_n<-seq(2,n,by=2)
odd<-sample(c(2,-1), size=n/2, replace=TRUE, prob=c(0.6,0.4))
even<-sample(c(-2,1), size=n/2, replace=TRUE, prob=c(0.6,0.4))
db<-data.frame(value=c(odd,even),number=c(odd_n,even_n))
row.names(db) <- NULL
db
value number
1 -1 1
2 -2 2
3 2 3
4 -2 4
5 2 5
6 -2 6
7 -1 7
8 1 8
9 2 9
10 -2 10
11 2 11
12 1 12
13 2 13
14 -2 14
15 -1 15
16 -2 16
17 -1 17
18 1 18
19 2 19
20 1 20
Run Code Online (Sandbox Code Playgroud)
| 归档时间: |
|
| 查看次数: |
224 次 |
| 最近记录: |