estimation in R

Καλησπέρα πως θα μπορούσα να εκτιμήσω τις παραμέτρους α,β,γ της παρακάτω σχέσης με την μέθοδο μέγιστης πιθανοφάνειας (mle) στην R?
[1-exp(-α-β*χ)]*exp(-γ*χ)

Υποθέτω εννοείς ότι έχεις κάποια δεδομένα Χ που ακολουθούν την κατανομή με σ.π.π. [1-exp(-α-β*χ)]*exp(-γ*χ)
άρα η loglikehood θα ναι η Σlog([1-exp(-α-β*χ)]*exp(-γ*χ))

Αρχικά φτιάχνεις την -loglikehood
loglike <- function(a,b,c){-sum(log(1-exp(-a-b*x))-c*x)}


Κοίτα γενικά το
?mle (θέλει το πακέτο stats4)
όπως και το ?mle2 (θέλει το πακέτο bbmle)
με το mle θα πρέπει μάλλον να είναι σχετικά κοντά οι start = list()
στο mle2 δεν είναι αναγκαίο, αλλά δεν ξέρω και τι ακριβώς κάνει



μετά τρέχεις κάτι σαν αυτό

fit0 <- mle2(loglike, start = list(a=1,b=1,c=2), method = "BFGS")
summary(fit0)

πρέπει να έχει κάνει install και να χεις φορτώσει το πακέτο bbmle