出差回来了，回个文给你。
## 建立一组虚假资料并进行 2 x 4 设计的 logistic regression
y <- c(rep(0, 1), rep(1, 19), rep(0, 16), rep(1, 4))
fa <- gl(2, 20)
fb <- gl(4, 1, 40)
fit <- glm(y ~ fa + fb, family = binomial)
summary(fit)
## 小数点后想看很多位数
print(summary(fit)$coefficients, digits=10)
## 针对 fa 和 fb 的 Wald test
library(lmtest)
waldtest(fit, . ~ . - fa, test = "Chisq")
waldtest(fit, . ~ . - fb, test = "Chisq")
请研究上面的例子。有问题欢迎再讨论。
※ 引述《lepin2001 (排年代的票真闷)》之铭言：
: [问题类型]:
: 程式咨询(我想用R 做某件事情，但是我不知道要怎么用R 写出来)
: [软件熟悉度]:
: 入门(写过其他程式，只是对语法不熟悉)
: [问题叙述]:
: 我用r 想做到spss的logistic regression
: 如图
: http://ppt.cc/SgUa
: 但缺黄色的和红色框线部份
: 不知该如何做
: [程式范例]:
: R的CODE 与ＳＰＳＳ
: data_test<-read.csv("test.csv",header=T)
: attach(data_test)
: agelevel<-factor(AgeGroup,levels=c(1,2,3,4))
: ublevel<-factor(Urbanization,levels=c(1,2,3,4))
: test_lm<-glm(Result~ExpGroup+relevel(agelevel,1)+relevel(ublevel,1),
: + data=data_test,family=binomial)
: summary(test_lm)
: R code 与spss 如图
: http://ppt.cc/UJEn
: data set:
: test.csv
: http://ppt.cc/dfrK
: 麻烦各位了!