| 所在主题: | |
| 文件名: t.test.rar | |
| 资料下载链接地址: https://bbs.pinggu.org/a-1226124.html | |
本附件包括:
|
|
| 附件大小: | |
|
对汤老师的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、本站为非盈利性质的学术交流网站,鼓励和保护原创作品,拒绝未经版权人许可的上传行为。本站如接到版权人发出的合格侵权通知,将积极的采取必要措施;同时,本站也将在技术手段和能力范围内,履行版权保护的注意义务。 (如有侵权,欢迎举报) |
|
京ICP备16021002号-2 京B2-20170662号
京公网安备 11010802022788号
论坛法律顾问:王进律师
知识产权保护声明
免责及隐私声明