搜索
人大经济论坛 附件下载

附件下载

所在主题:
文件名:  t.test.rar
资料下载链接地址: https://bbs.pinggu.org/a-1226124.html
本附件包括:
  • t.test.R
附件大小:
对汤老师的p-value有点小小的疑惑?R语言与统计分析-by汤银才-p171
此题是方差sigma已知的情况下, 对假设进行检验,并对均值区间进行估计
z.test<-function(x,n,sigma,alpha,u0=0,alternative="two.sided"){
options(digits=4)
result<-list( )
mean<-mean(x)
z<-(mean-u0)/(sigma/sqrt(n))
p<-pnorm(z,lower.tail=FALSE)
result$mean<-mean
result$z<-z
result$p.value<-p
if(alternative=="two.sided"){
p<-2*p
result$p.value<-p
}
else if (alternative == "greater"|alternative =="less" ){
result$p.value<-p
}
else return("your input is wrong")
result$conf.int<- c(mean-sigma*qnorm(1-alpha/2,mean=0, sd=1,lower.tail = TRUE)/sqrt(n),
mean+sigma*qnorm(1-alpha/2,mean=0, sd=1,lower.tail = TRUE)/sqrt(n))
result
}
红色部分,是疑惑的地方,p.value是取得 lower.tail=F,是上侧的面积,那么 less时即假设H0:u>= u0时,p.value =p是不是错误的呢?另外,双侧检验的时候,z如果是负数,那么2*p(上侧面积)会不会出现 大于 1的情况!!!!

也尝试修改下,大大们帮忙看看有没有统计知识的问题撒?
z.test <- function(x,n,sigma,alpha,u0=0,alternative="two.sided"){
options(digits=4)
result <- list()
mean <- mean(x)
z <-(mean-u0)/(sigma/sqrt(n))
result$mean <- mean
result$z <- z
ifelse(alternative=="two.sided",result$p.value <- 2*(1-pnorm(abs(z))), #双侧检测
ifelse(alternative=="less",result$p.value <- pnorm(z,lower.tail=T), #下侧检验
ifelse(alternative=="greater",result$p.value <-pnorm(z,lower.tail=F), #上侧检验
return("your input is wrong"))))
result$conf.int <- c(mean-qnorm(1-alpha/2)*sigma/sqrt(n),
mean+qnorm(1-alpha/2)*sigma/sqrt(n)) #有误
result
}

刚学到 假设检验, p-value 又是太重要了
自学统计, 没老师。。。
请大大们帮忙指点统计知识的问题撒 谢谢

谢谢 epoh老师的指导 之前把alternative 给搞反了,p.value 的面积真是有点绕


    熟悉论坛请点击新手指南
下载说明
1、论坛支持迅雷和网际快车等p2p多线程软件下载,请在上面选择下载通道单击右健下载即可。
2、论坛会定期自动批量更新下载地址,所以请不要浪费时间盗链论坛资源,盗链地址会很快失效。
3、本站为非盈利性质的学术交流网站,鼓励和保护原创作品,拒绝未经版权人许可的上传行为。本站如接到版权人发出的合格侵权通知,将积极的采取必要措施;同时,本站也将在技术手段和能力范围内,履行版权保护的注意义务。
(如有侵权,欢迎举报)
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

GMT+8, 2025-12-31 20:29