Hi,

I am trying to use the Fisher scoring method with a geometric distribution,
with p = .07, 100 observations from the geom distrib,  and 10 iterations.
I cannot quite get the code to work.
Can anyone see the mistake?




n <- 100

 p <- 0.07

 x <- rgeom(n, p)

s  <- sum(x)

 f <- function(x, p) p*(1-p)^x

 L <- function(p)  p^n*(1-p)^s

logL <- function(p)  n*log(p)+s*(log(1-p))

logLprime <- function(p)  (n/p)-(s/(1-p))



I <- n/p^2*(1-p)

iter <- 10

p[1] <- .06

p[2] <- .11





for (i in 1:10)

{

                pnew <- p[i]+logLprime(p[i])/I*(p[i])

}

	[[alternative HTML version deleted]]

