Nothing
## ----comment=NA---------------------------------------------------------------
require(moonBook) # for use of data radial
fit=lm(NTAV~age,data=radial)
summary(fit)
## ----comment=NA,message=FALSE,fig.width=5,fig.height=5------------------------
require(ggplot2)
ggplot(radial,aes(y=NTAV,x=age))+geom_point()+geom_smooth(method="lm")
## ----comment=NA,message=FALSE-------------------------------------------------
require(ggiraph)
require(ggiraphExtra)
require(plyr)
ggPredict(fit,se=TRUE,interactive=TRUE)
## ----comment=NA---------------------------------------------------------------
fit1=lm(NTAV~age+sex,data=radial)
summary(fit1)
## ----comment=NA,message=FALSE,fig.width=5,fig.height=5------------------------
equation1=function(x){coef(fit1)[2]*x+coef(fit1)[1]}
equation2=function(x){coef(fit1)[2]*x+coef(fit1)[1]+coef(fit1)[3]}
ggplot(radial,aes(y=NTAV,x=age,color=sex))+geom_point()+
stat_function(fun=equation1,geom="line",color=scales::hue_pal()(2)[1])+
stat_function(fun=equation2,geom="line",color=scales::hue_pal()(2)[2])
## -----------------------------------------------------------------------------
ggPredict(fit1,se=TRUE,interactive=TRUE)
## ----comment=NA---------------------------------------------------------------
fit2=lm(NTAV~age*DM,data=radial)
summary(fit2)
## ----comment=NA,message=FALSE,fig.width=5,fig.height=5------------------------
ggplot(radial,aes(y=NTAV,x=age,color=factor(DM)))+geom_point()+stat_smooth(method="lm",se=FALSE)
## -----------------------------------------------------------------------------
ggPredict(fit2,colorAsFactor = TRUE,interactive=TRUE)
## ----comment=NA---------------------------------------------------------------
fit3=lm(NTAV~age*weight,data=radial)
summary(fit3)
## ----comment=NA,message=FALSE,fig.width=5,fig.height=5------------------------
ggplot(radial,aes(y=NTAV,x=age,color=weight))+geom_point()+stat_smooth(method="lm",se=FALSE)
## -----------------------------------------------------------------------------
ggPredict(fit3,interactive=TRUE)
## ----comment=NA,message=FALSE,fig.width=5,fig.height=5------------------------
fit4=lm(NTAV~age*weight*HBP,data=radial)
summary(fit4)
## -----------------------------------------------------------------------------
ggPredict(fit4,interactive = TRUE)
## ----comment=NA,message=FALSE-------------------------------------------------
require(TH.data) # for use data GBSG2
fit5=glm(cens~pnodes*horTh,data=GBSG2,family=binomial)
summary(fit5)
## -----------------------------------------------------------------------------
ggPredict(fit5,se=TRUE,interactive=TRUE,digits=3)
## ----comment=NA,message=FALSE-------------------------------------------------
fit6=glm(cens~pnodes+horTh,data=GBSG2,family=binomial)
summary(fit6)
## -----------------------------------------------------------------------------
ggPredict(fit6,se=TRUE,interactive=TRUE,digits=3)
## ----comment=NA,message=FALSE-------------------------------------------------
fit7=glm(cens~pnodes*age,data=GBSG2,family=binomial)
summary(fit7)
## ----warning=FALSE------------------------------------------------------------
ggPredict(fit7,interactive=TRUE)
## ----warning=FALSE------------------------------------------------------------
ggPredict(fit7,interactive=TRUE,colorn=100,jitter=FALSE)
## ----comment=NA,message=FALSE-------------------------------------------------
fit8=glm(cens~pnodes*age*horTh,data=GBSG2,family=binomial)
summary(fit8)
## ----warning=FALSE------------------------------------------------------------
ggPredict(fit8,interactive=TRUE,colorn=100,jitter=FALSE)
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.