Hello!
I've fitted a time series to the exponential power density (Subbotin) with
the following function:
library(MASS)
dep<- function(x,location,scale,tail) # Exponential power (=Subbotin)
{
const<- 2*scale*tail^(1/tail) *gamma(1+1/tail)
z<- (x-location)/scale
exp(-1/tail*abs(z)^tail)/const
}
x <- rt(100,df=3)
a <- fitdistr(x, dep, start=list(location=0, scale=1, tail=1.5))
print(a$estimate)
print(a$sd)
Now i want to test the goodness of fit with the Kolmogorow-Smirnov Test! How can
i do this?
thank you very much!
Sara
--