出差回来了，回个文给你。 ## 建立一组虚假资料并进行 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/cn.aspx?n=bbs/R_Language/M.1422369395.A.D52.html