如何在R中编写学生的t分布？

Question

mu <- 0
sigma <- 1
x <- seq(from = -3, to  =3, by = .1)
nu<-10
y.1 <- gamma*((nu+1)/2)/(sqrt(pi*nu)*gamma*(nu/2))*(1+(t^2)/nu)^-((nu+1)/2)

这是结果：

伽玛错误*（（nu + 1）/ 2）：二进制的非数字参数 运算符

我想写学生的t分布。如何解决这个问题？

Answer 1

gamma是一种功能：如果您想将其应用于仅写的内容：

gamma((nu + 1)/2)

没有*。 另外，您还必须将t替换为x： 所以这应该工作：

y.1 <- gamma((nu+1)/2) / (sqrt(pi*nu) * gamma(nu/2)) * (1+(x^2)/nu)^-((nu+1)/2)

无论如何，t分布的密度函数已经在R中实现。 参见?dt