论坛好冷清啊。
循环次数有点多,下面的程序可能慢一点点。懒得写函数了,但是基本就是这么个思路了。
- #set.seed(100)
- #Q1:
- a1=rnorm(30,10,5)
- a2=rnorm(30,14,5)
- t_test=t.test(a1,a2)
- t_test$p.value
- #Q2:
- b1=rnorm(30,10,5)
- b2=rnorm(30,10,5)
- t_test=t.test(b1,b2)
- t_test$p.value
- #Q3
- market_sample=NULL;
- for(i in 1:100000){
- market_type=sample(c("bull","bear","stagnant"),1,prob=c(0.15,0.8,0.05))
- market_sample=c(market_sample,market_type)
- }
- (result_vector=table(market_sample))
复制代码
|