Displaying 1 result from an estimated 1 matches for "tp16084593p16093575".
2008 Mar 16
2
(no subject)
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)