A=c(79,74,88,80,80,66,65,86,84,80,78,72,71,74,86,96,77,81,76,80, 76,75,78,87,87,74,85,84,76,77,76,74,85,74,76,77,76,74,81,76) B=length(which(A>70)) #大於70的樣本數 n=length(A) #樣本數 p=0.5 mu=n*p #二項分配期望值 s=sqrt(n*p*(1-p)) #二項分配標準差 ZB.test=(B-mu)/s #檢定統計量 (標準化) cat("ZBtest-statistic=",ZB.test,"\n") #在此可不加"/n"(只出現一個資料) #"/n"=>當有多個資料時,會自動換行不會擠在同一行 if(ZB.test>qnorm(1-0.05)){ #如果檢定統計量大於Z(0.05)=1.645 cat("reject H0 \n") #=>顯示reject H0 }else{ #否則 cat("not reject Ho \n") #顯示 => not reject H0 } p.value=1-pnorm(ZB.test) #P值 p(Z>ZBtest)=1-p(Z