El ajuste de una distribución t de Student en R usando fitdistr () produce el error "valor de diferencia finita no finita"
Un ejemplo reproducible que proporcionará el código de error mencionado cada vez es: (Tenga en cuenta que incluso sin set.seed, el error aparece cada vez)
library(MASS)
set.seed(seed = 1)
data<-rnorm(n = 10000,mean = 0.0002,sd = 0.001)
fitdistr(x = data,densfun = "t")
El mensaje de error es:
Error in stats::optim(x = c(-0.000426453810742332, 0.000383643324222082, :
non-finite finite-difference value [2]
In addition: Warning message:
In log(s) : NaNs produced
El problema es el "valor de diferencia finita no finito". Fitdistr no me da un resultado.
Mi conocimiento: investigué y aparentemente esto podría significar que un parámetro es negativo durante la iteración. Y que la solución podría ser proporcionar un valor inicial mejor o al menos diferente. Pero no pude averiguar cómo hacer esto y no estoy seguro de si este es el problema.
MI PREGUNTA: a) ¿Por qué recibo este mensaje de error yb) cómo puedo solucionarlo en R, de modo que pueda ajustar la distribución de t de estudiante a mis datos distribuidos normalmente?