ISL书的抽奖昨天顺利结束，由于很多评论的点赞数没有达到参与抽奖要求的30个，所以最后只有两位读者进入抽奖备选名单：
我们可以用简单的伯努利试次来进行抽奖，但是由于这是一个统计学的公众号，所以我们要用一点高大上的抽奖方式（毕竟简单事物复杂化是学术研究者最擅长的事情）。我们使用单样本z检验来进行抽奖，大家可以顺便根据以下内容来测验一下自己的假设检验知识！
首先，我们进行零假设为总体均数等于0的单样本z检验，总体标准差为5，我们需要将z检验的power控制在50%，由此，如果获得显著结果，则网名为“腹有诗书气自华”的读者中奖；如果获得不显著结果，则网名为“Sapphire^^”的读者中奖，这样二人就均有50%的概率获奖了！
为了将power控制在50%，我们需要对真实均数进行一些设定，具体检验参数见如下R代码（mu为真实均数）：
n=100sd=5SE=sd/sqrt(n)mu=qnorm(0.95,mean=0)*SE
通过如上R代码进行10000次取样和z检验，我们发现显著结果的比例（power）确实为50%左右，说明我们的设定没问题：
sum(sig_res)/10000 # power应为50%左右，所以结果应该靠近0.5[1] 0.4982
接下来，我们将随机种子设为今天大年初一的日期20250129，再跑一遍全程的代码。最后，我们这10000次检验中第2025次检验的结果作为本次抽奖中奖人的结果，那么第2025次检验的结果是：
> sig_res[2025][1] FALSE
结果为不显著，恭喜网名为“Sapphire^^”的读者获奖！该名读者可以在后台联系我，商量后续的奖励书籍发放！本次抽奖完全公开透明，符合开放科学的可重复原则，欢迎大家复制粘贴以下完整代码来复现我们的抽奖结果！
set.seed(20250129)# 显著-腹有诗书气自华；不显著-Sapphire^^n=100sd=5SE=sd/sqrt(n)mu=qnorm(0.95,mean=0)*SEsig_res=rep(NA,10000)for (i in 1:10000){data=rnorm(n,mu,sd)z_score=mean(data)/SEsig_res[i]=z_score>qnorm(0.95)}sum(sig_res)/10000 # power 应为50%左右，所以结果应该靠近0.5sig_res[2025]
最后，祝贺大家新年快乐！也预祝大家在新的一年沉迷统计，统计水平更上一层楼！
