oneside <- function(a,b,u1,u2,c) {
a1<- qnorm(1-a)
z<- u1-u2
b1<- qnorm(1-b)
((a1-b1)^2*c^2)/z^2
}
a <- 0.05
b <- 0.1
u1 <- 1.4
u2 <- 1.42
c <- 0.04^2
oneside(a,b,u1,u2,c)
我是想输入a,b,u1,u2,c计算((a1-b1)^2*c^2)/z^2 可是结果不对，哪里有问题？我怀疑是a1计算的问题，结果输入了a的值，计算a1报错，说错误找不到a1.求大神帮忙呀，我是初学者