r - 如何在 fitdistrplus 中使用尺度和位置参数拟合 t 分布

Question

如何使用 fitdistrplus 来估计 t 分布的尺度参数的位置？我知道我需要提供初始值(在 MASS 中效果很好)，但在这个包中只允许使用 df 。你有什么解决办法吗？

非常感谢。

Answer 1

基础R中的位置和尺度参数可以根据这个Wikipedia进行扩展文章。

library(fitdistrplus)
x<-rt(100,23)
dt_ls <- function(x, df=1, mu=0, sigma=1) 1/sigma * dt((x - mu)/sigma, df)
pt_ls <- function(q, df=1, mu=0, sigma=1)  pt((q - mu)/sigma, df)
qt_ls <- function(p, df=1, mu=0, sigma=1)  qt(p, df)*sigma + mu
rt_ls <- function(n, df=1, mu=0, sigma=1)  rt(n,df)*sigma + mu
fit.t<-fitdist(x, 't_ls', start =list(df=1,mu=mean(x),sigma=sd(x))) 
summary(fit.t)