Grstats Forum

estimation in R

Επισκόπηση προηγούμενης Θ.Ενότητας Επισκόπηση επόμενης Θ.Ενότητας Πήγαινε κάτω

estimation in R

Δημοσίευση από sansindy2 Την / Το Δευ 9 Μαρ 2015 - 0:18

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

sansindy2

Posts : 7
Join date : 03/06/2014

Επιστροφή στην κορυφή Πήγαινε κάτω

Απ: estimation in R

Δημοσίευση από barney84 Την / Το Δευ 9 Μαρ 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




avatar
barney84

Posts : 18
Join date : 25/01/2014

Επιστροφή στην κορυφή Πήγαινε κάτω

Επισκόπηση προηγούμενης Θ.Ενότητας Επισκόπηση επόμενης Θ.Ενότητας Επιστροφή στην κορυφή


 
Δικαιώματα σας στην κατηγορία αυτή
Δεν μπορείτε να απαντήσετε στα Θέματα αυτής της Δ.Συζήτησης