Displaying 2 results from an estimated 2 matches for "loglik_sum".
2002 Apr 22
3
glm() function not finding the maximum
...9.98, 554.13, 104.88, 170.15, 945.1, 143.52)
#Fits data to a gamma distribution using glm()
gamma1_function(data){
n_length(data)
m_summary(glm(data~1, family=Gamma(link=identity)))
shape_1/as.numeric(m$disp)
scale_as.numeric(m$coeff[1]*m$disp)
dev.res_-2*log(dgamma(data,shape=shape,scale=scale))
loglik_sum(dev.res) #actually -2 * log like
list(loglik=loglik,par=c(shape,scale))
}
#Fits data to a gamma distribution "by hand" using optim()
gamma2_function(data){
n_length(data)
m_summary(glm(data~1, family=Gamma))
shape_1/as.numeric(m$disp)
#L = -Log likelihood
L_function(x){-(-n*log(gamma(...
2013 Nov 03
1
Comparison of two weibull distributions
Hello,
How can I do a test of two weibull distributions?
I have two weibull distribution sets from two wind datasets in order to
check whether they are same.
I thought 2 sample t-test would be applicable but I couldn't find any ways
to do that on the Internet.
Does anyone know what type of test is applicable to my purpose? and what R
function can you recommend?
Plus, if it turned out that