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.求大神帮忙呀,我是初学者


雷达卡


京公网安备 11010802022788号







