Guest User

Untitled

a guest
Jan 13th, 2018
79
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.51 KB | None | 0 0
  1. set.seed(1)
  2. #-------------------------
  3. X = rpois(n=10000,lambda=100)
  4. B1 = 0.74723
  5. B2 = 0.09646
  6. r.error = rpois(n=10000,lambda=0.05)
  7. Y = exp(B1+B2*X)+r.error
  8. F = data.frame(Y,X)
  9.  
  10. ##GLM fitting-------------------------
  11. F.test = data.frame(X=seq(90,120,1))
  12. F.fit.p = glm(Y~X,data=F,family=poisson)
  13. F.pred.p = predict(F.fit.p,newdata = F.test,se.fit=T)
  14. F.pred.p$se.fit
  15.  
  16. #lm fitting---------------------------
  17. F.fit.g = glm(log(Y)~X,data=F,family=gaussian)
  18. F.pred.g = predict(F.fit.g,newdata = F.test,se.fit=T)
  19. F.pred.g$se.fit
Add Comment
Please, Sign In to add comment