【数据分析入门】R语言之广义线性回归与logistics回归

广义线性回归

?glm
data(breslow.dat, package="robust") names(breslow.dat)
summary(breslow.dat[c(6，7，8，10)])
attach(breslow.dat) 

> fit <-glm(sumY ~ Base + Age + Trt,data=breslow.dat, family=poisson(link="log")) 
> summary(fit)

coef(fit)
exp(coef(fit))

logistics回归

data(Affairs, package="AER") summary(Affairs)
table(Affairs$affairs)
prop.table(table(Affairs$affairs)) prop.table(table(Affairs$gender))

Affairs$ynaffair[Affairs$affairs > 0] <-1 Affairs$ynaffair[Affairs$affairs==0] <-0 Affairs$ynaffair <-factor(Affairs$ynaffair,                     levels=c(0,1),labels=c("NO","Yes"))
table(Affairs$ynaffair)

attach(Affairs)
fit<-glm(ynaffair ~ gender + age + yearsmarried + children + 
                   religiousness+ education+occupation +rating, data=Affairs,family=binomial())
summary(fit)

fit1<-glm(ynaffair ~  age + yearsmarried + 
                   religiousness+  +rating, data=Affairs,family=binomial())
summary(fit)

anova(fit,fit1,test="Chisq")

data.frame(rating =c(1,2,3,4，5)
age = mean(Affairs$age),
yearsmarried = mean(Affairs$yearsmarried), rcliaiousness = mean(Affairs$religiousness)
testdata$prob <- predict(fitl, newdata=testdata, type="response")
testdata <-data.frame(rating = mean(Affairs$rating),
age = seq(17,57,10),
yearsmarried = mean(Affairs$yearsmarried), reliaiousness= mean(Affairs$religiousness)