Displaying 2 results from an estimated 2 matches for "freqtab".
2011 Feb 07
2
question mle again
...--------------
lambda.data <- runif(10,0.5,10)
ll <- function(lambda = 1) {
cat("x in ll()",x,"\n")
y.fit <- dpois(x, lambda)
sum( (y - y.fit)^2 )
}
lapply(1:10, FUN = function(x){
raw.data <- rpois(100,lambda.data[x])
freqTab <- count(raw.data)
x <- freqTab$x
y <- freqTab$freq / sum(freqTab$freq)
cat("x in lapply", x,"\n")
fit <- mle(ll)
coef(fit)
})
Can anybody help?
Antje
2011 Feb 01
1
mle question
...# R-code ---------------------------------
lambda.data <- runif(10,0.5,10)
ll <- function(lambda = 1) {
cat("x in ll()",x,"\n")
y.fit <- dpois(x, lambda)
sum( (y - y.fit)^2 )
}
lapply(1:10, FUN = function(x){
raw.data <- rpois(100,lambda.data[x])
freqTab <- count(raw.data)
x <- freqTab$x
y <- freqTab$freq / sum(freqTab$freq)
cat("x in lapply", x,"\n")
fit <- mle(ll)
coef(fit)
})
Can anybody help?
Antje