Problem 2, Øving 10

7 views
Skip to first unread message

TS

unread,
Jun 8, 2011, 3:44:38 AM6/8/11
to ST2304 Statistical Modelling for Biologists/Biotechnologists

Hei!

Jeg lurer på hvorfor og når man skal skrive par når man lager
funksjoner? Feks, i problem 2 i øving 10 brukes dette:

> lnLP <- function(par,x,n,time) {
+ beta0 <- par[1]
+ beta1 <- par[2]
+ p <- pnorm(beta0+beta1*time)
+ -sum(dbinom(x,size=n,prob=p,log=T)
+ }

Hvorfor brukes par her feks.?

Jarle Tufto

unread,
Jun 8, 2011, 4:16:32 AM6/8/11
to ST2304 Statistical Modelling for Biologists/Biotechnologists
Hvis du ser på hjelpesiden til optim ser du at funksjonen som optim
skal minimalisere må være skrevet på en bestemt måte:

fn: A function to be minimized (or maximized), with first
argument the vector of parameters over which minimization is
to take place. It should return a scalar result.

Første argument må være en vektor som inneholder parameterne vi ønsker
å maksimalisere funskjonen m.h.t. Dette er bare en konvensjon som er
valgt slik at optim kan snakke med vår likelihoodfunksjon... Om vi
kaller dette argument par eller noe annet spiller ingen rolle.
Reply all
Reply to author
Forward
0 new messages