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

※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 125.230.76.152 ※ 文章網址: https://webptt.com/m.aspx?n=bbs/R_Language/M.1422369395.A.D52.html