小编Rfa*_*fan的帖子

有没有办法从"glmer"对象中获得"边际效应"

我正在使用估计随机效应logit模型glmer,我想报告自变量的边际效应.对于glm模型,包mfx有助于计算边际效应.glmer对象有任何包或函数吗?

谢谢你的帮助.

下面给出了可重现的例子

mydata <- read.csv("http://www.ats.ucla.edu/stat/data/binary.csv")
mydata$rank <- factor(mydata$rank) #creating ranks
id <- rep(1:ceiling(nrow(mydata)/2), times=c(2)) #creating ID variable
mydata <- cbind(mydata,data.frame(id,stringsAsFactors=FALSE)) 
set.seed(12345)
mydata$ran <- runif(nrow(mydata),0,1) #creating a random variable

library(lme4)
cfelr <- glmer(admit ~ (1 | id) + rank + gpa + ran + gre, data=mydata ,family = binomial)
summary(cfelr)
Run Code Online (Sandbox Code Playgroud)

r lme4 marginal-effects

21
推荐指数
2
解决办法
3210
查看次数

在每组数据后插入一个空行

这个问题已被要求excel.如何在一组数据后自动插入空行

我想知道是否有一个函数在R中做同样的事情.

例:

group <- c("a","b","b","c","c","c","d","d","d","d")
xvalue <- c(16:25)
yvalue <- c(1:10)
df <- data.frame(cbind(group,xvalue,yvalue))
   group xvalue yvalue
1      a     16      1
2      b     17      2
3      b     18      3
4      c     19      4
5      c     20      5
6      c     21      6
7      d     22      7
8      d     23      8
9      d     24      9
10     d     25     10
Run Code Online (Sandbox Code Playgroud)

我想在每个组之后有一个空行,以便我可以手动检查条目

   group xvalue yvalue
1      a     16      1
2
3      b     17      2
4      b     18      3
5
6      c     19 …
Run Code Online (Sandbox Code Playgroud)

r

8
推荐指数
2
解决办法
6582
查看次数

标签 统计

r ×2

lme4 ×1

marginal-effects ×1