Nothing
## ------------------------------------------------------------------------
library(isni)
data(sos)
sos[sample(nrow(sos),10),]
## ------------------------------------------------------------------------
ymodel= sexact ~ gender*faculty
summary(glm(ymodel,family=binomial, data=sos))
## ------------------------------------------------------------------------
sos.isni<-isniglm(ymodel, family=binomial, data=sos)
sos.isni
## ------------------------------------------------------------------------
summary(sos.isni)
## ------------------------------------------------------------------------
ygmodel <- sexact | is.na(sexact) ~ gender*faculty | gender *faculty
summary(isniglm(ygmodel, family=binomial, data=sos))
## ------------------------------------------------------------------------
gender <- c(0,0,1,1,0,0,1,1)
faculty <- c(0,0,0,0,1,1,1,1)
gender <- factor(gender, levels = c(0, 1), labels =c("male", "female"))
faculty <- factor(faculty, levels = c(0, 1), labels =c("other", "mdv"))
SAcount <- c(NA, 1277, NA, 1247, NA, 126, NA, 152)
total <- c(1189,1710,978,1657,68,215,73,246)
sosgrp <- data.frame(gender=gender, faculty=faculty, SAcount=SAcount, total=total)
ymodel <- SAcount/total ~gender*faculty
sosgrp.isni<-isniglm(ymodel, family=binomial, data=sosgrp, weight=total)
summary(sosgrp.isni)
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.