出差回來了，回個文給你。
## 建立一組虛假資料並進行 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
: 麻煩各位了!