Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- require(lmtest)
- require(splines)
- x=rnorm(20 ,0,1)
- y=rep(c(0,1),times=10)
- fit1=glm(y~x ,family = "binomial")
- fit2=glm(y~ns(x,df=2) ,family = "binomial")
- Model 1: y ~ x
- Model 2: y ~ ns(x, df = 2)
- #Df LogLik Df Chisq Pr(>Chisq)
- 1 2 -12.952
- 2 3 -12.757 1 0.3912 0.5317
- summary(fit1)
- Coefficients:
- Estimate Std. Error z value Pr(>|z|)
- (Intercept) 0.2367 0.4987 0.475 0.635
- x -1.1804 0.9331 -1.265 0.206
- (Dispersion parameter for binomial family taken to be 1)
- Null deviance: 27.726 on 19 degrees of freedom
- Residual deviance: 25.904 on 18 degrees of freedom
- AIC: 29.904
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement