I have a data set $Y$, $Y$ is assumed to follow a poisson distribution with mean = $\lambda$ .
However, each element $y_i$ in $Y$ is accompanied with a covariate $x_i$. So here I want to use poisson regression to model this dataset. The link function is assumed to be log($\lambda(i)$)=$\beta_0$+$\beta_1x_i$.
I then want to use glm(Y~X,family='poisson') to fit the model and get $\beta_0$ and $\beta_1$.
After get $\beta_0$ and $\beta_1$, can I calculate the probability mass function given new $y$ and $x$.
I want to know if all of my above thought make sense?