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分布。如何解决这个问题?
答案 0 :(得分:2)
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