Λέσχη Φίλων Στατιστικής - GrStats forum
estimation in R Forumgrstats

Join the forum, it's quick and easy

Λέσχη Φίλων Στατιστικής - GrStats forum
estimation in R Forumgrstats
Λέσχη Φίλων Στατιστικής - GrStats forum
Would you like to react to this message? Create an account in a few clicks or log in to continue.
Για προβλήματα εγγραφής και άλλες πληροφορίες επικοινωνήστε με : grstats.forum@gmail.com ή grstats@stat-athens.aueb.gr

Go down
avatar
sansindy2
Posts : 7
Join date : 2014-06-03

estimation in R Empty estimation in R

Mon 9 Mar 2015 - 0:18
Καλησπέρα πως θα μπορούσα να εκτιμήσω τις παραμέτρους α,β,γ της παρακάτω σχέσης με την μέθοδο μέγιστης πιθανοφάνειας (mle) στην R?
[1-exp(-α-β*χ)]*exp(-γ*χ)
barney84
barney84
Posts : 19
Join date : 2014-01-25

estimation in R Empty Re: estimation in R

Mon 9 Mar 2015 - 18:10


Υποθέτω εννοείς ότι έχεις κάποια δεδομένα Χ που ακολουθούν την κατανομή με σ.π.π. [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




Back to top
Permissions in this forum:
You cannot reply to topics in this forum