Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- set.seed(1)
- #-------------------------
- X = rpois(n=10000,lambda=100)
- B1 = 0.74723
- B2 = 0.09646
- r.error = rpois(n=10000,lambda=0.05)
- Y = exp(B1+B2*X)+r.error
- F = data.frame(Y,X)
- ##GLM fitting-------------------------
- F.test = data.frame(X=seq(90,120,1))
- F.fit.p = glm(Y~X,data=F,family=poisson)
- F.pred.p = predict(F.fit.p,newdata = F.test,se.fit=T)
- F.pred.p$se.fit
- #lm fitting---------------------------
- F.fit.g = glm(log(Y)~X,data=F,family=gaussian)
- F.pred.g = predict(F.fit.g,newdata = F.test,se.fit=T)
- F.pred.g$se.fit
Add Comment
Please, Sign In to add comment