# 如何合并两个区间
# 比如合并GBD两个年龄组的区间,两个年龄组数据如下
# 0-4岁率val=0.2，lower=0.1,upper=0.4
# 5-9岁率val=0.4，lower=0.25,upper=0.56

set.seed(1234)

# 第一步，计算两组的se,se=(upper-lower)/(1.96*2)
se1=(0.4-0.1)/(1.96*2)
se2=(0.56-0.25)/(1.96*2)

val1=0.2
val2=0.4

# 第二步，生成均值为val，方差为se服从正态分布的随机数随便多少个
a <- val1+se1*rnorm(100000) #rnorm(10000)表示生成10000个均值为0的标准正态分布
b <- val2+se2*rnorm(100000)

c <- a+b

# 第三步有放回抽样：
temp <- sample(c,1000,replace = T)

# 第四部，选择位于2.5,97.5分位数的值，为新的区间
quantile(temp,c(0.025,0.5,0.975))