Advertisement
Guest User

Untitled

a guest
Jun 19th, 2019
83
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.68 KB | None | 0 0
  1. require(lmtest)
  2. require(splines)
  3. x=rnorm(20 ,0,1)
  4. y=rep(c(0,1),times=10)
  5.  
  6. fit1=glm(y~x ,family = "binomial")
  7.  
  8. fit2=glm(y~ns(x,df=2) ,family = "binomial")
  9.  
  10. Model 1: y ~ x
  11. Model 2: y ~ ns(x, df = 2)
  12. #Df LogLik Df Chisq Pr(>Chisq)
  13. 1 2 -12.952
  14. 2 3 -12.757 1 0.3912 0.5317
  15.  
  16. summary(fit1)
  17.  
  18. Coefficients:
  19. Estimate Std. Error z value Pr(>|z|)
  20. (Intercept) 0.2367 0.4987 0.475 0.635
  21. x -1.1804 0.9331 -1.265 0.206
  22.  
  23. (Dispersion parameter for binomial family taken to be 1)
  24.  
  25. Null deviance: 27.726 on 19 degrees of freedom
  26. Residual deviance: 25.904 on 18 degrees of freedom
  27. AIC: 29.904
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement