Simulated Sample from Poisson Distribution
The data set is used to understand the sampling variation of the score function. The simulated data is available in Pawitan (2001).
data(ps)
A data frame with 10 observations on the following 20 variables.
Sample.1a numeric vector
Sample.2a numeric vector
Sample.3a numeric vector
Sample.4a numeric vector
Sample.5a numeric vector
Sample.6a numeric vector
Sample.7a numeric vector
Sample.8a numeric vector
Sample.9a numeric vector
Sample.10a numeric vector
Sample.11a numeric vector
Sample.12a numeric vector
Sample.13a numeric vector
Sample.14a numeric vector
Sample.15a numeric vector
Sample.16a numeric vector
Sample.17a numeric vector
Sample.18a numeric vector
Sample.19a numeric vector
Sample.20a numeric vector
Pawitan, Y. (2001). In All Likelihood. Oxford Science Publications.
Pawitan, Y. (2001). In All Likelihood. Oxford Science Publications.
data(ps) n <- 10 sample_means <- colMeans(ps) poisson_score_fn <- function(theta,xbar) n*(xbar-theta)/theta theta <- seq(from=2,to=8,by=0.2) plot(theta,sapply(theta,poisson_score_fn,xbar=sample_means[1]),"l",xlab= expression(lambda),ylab=expression(S(lambda)),ylim=c(-5,15)) title(main="B: Score Function Plot of the Poisson Model") for(i in 2:20) lines(theta,sapply(theta,poisson_score_fn,xbar=sample_means[i]),"l") abline(v=4) abline(h=0)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.