朱莉娅相当于R的qnorm()？

Question

我试图将这些行从R翻译成Julia:

n <- 100
mean <- 0
sd <- 1
x <- qnorm(seq(1 / n, 1 - 1 / n, length.out = n), mean, sd)

但是,我在使用qnorm函数时遇到了麻烦.我搜索了"分位数函数"并找到了quantile()函数.但是,R的版本返回长度为100的向量,而Julia的版本返回长度为5的向量.

这是我的尝试:

import Distributions
n = 100
x = Distributions.quantile(collect(range(1/n, stop=1-1/n, length=n))) 

Answer 1

在Julia 1.0下你应该广播这样的调用quantile:

quantile.(Normal(0, 1), range(1/n, 1-1/n, length = n))

或这个

quantile.(Normal(0, 1), range(1/n, 1-1/n, length = n))

关键部分是quantile没有定义广播行为(还),所以你必须把它包装成支持的东西quantile.